Home

Context Navigation

source: palm/trunk/SOURCE/gust_mod.f90 @ 4748

Last change on this file since 4748 was 4646, checked in by raasch, 4 years ago
files re-formatted to follow the PALM coding standard
Property svn:keywords set to `Id`
File size: 25.8 KB

Line
1	!> @file gust_mod.f90
2	!--------------------------------------------------------------------------------------------------!
3	! This file is part of the PALM model system.
4	!
5	! PALM is free software: you can redistribute it and/or modify it under the terms of the GNU General
6	! Public License as published by the Free Software Foundation, either version 3 of the License, or
7	! (at your option) any later version.
8	!
9	! PALM is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the
10	! implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General
11	! Public License for more details.
12	!
13	! You should have received a copy of the GNU General Public License along with PALM. If not, see
14	! <http://www.gnu.org/licenses/>.
15	!
16	! Copyright 1997-2020 Leibniz Universitaet Hannover
17	!--------------------------------------------------------------------------------------------------!
18	!
19	! Current revisions:
20	! -----------------
21	!
22	!
23	! Former revisions:
24	! -----------------
25	! $Id: gust_mod.f90 4646 2020-08-24 16:02:40Z pavelkrc $
26	! file re-formatted to follow the PALM coding standard
27	!
28	! 4535 2020-05-15 12:07:23Z raasch
29	! bugfix for restart data format query
30	!
31	! 4517 2020-05-03 14:29:30Z raasch
32	! added restart with MPI-IO for reading local arrays
33	!
34	! 4495 2020-04-13 20:11:20Z raasch
35	! restart data handling with MPI-IO added
36	!
37	! 4360 2020-01-07 11:25:50Z suehring
38	! CASE statement for dummy variable u2_av in gust_rrd_local changed to avoid unintended
39	! interdependencies with user-defined variables
40	!
41	! 3837 2019-03-28 16:55:58Z knoop
42	! unused variable for file index removed from rrd-subroutines parameter list
43	!
44	! 3725 2019-02-07 10:11:02Z raasch
45	! dummy statement modified to avoid compiler warnings about unused variables
46	!
47	! 3685 2019-01-21 01:02:11Z knoop
48	! Some interface calls moved to module_interface + cleanup
49	!
50	! 3665 2019-01-10 08:28:24Z raasch
51	! dummy statements added to avoid compiler warnings about unused variables
52	!
53	! 3655 2019-01-07 16:51:22Z knoop
54	! Bugfix: domain bounds of local_pf corrected
55	!
56	!
57	! Interfaces concerning data output updated
58	!
59	!
60	! renamed gust_par to gust_parameters
61	!
62	!
63	! Initial interface definition
64	!
65	!
66	! Description:
67	! ------------
68	!> Gust model.
69	!>
70	!> @todo This is just a dummy module. The actual module ist not released yet.
71	!--------------------------------------------------------------------------------------------------!
72	MODULE gust_mod
73
74	USE control_parameters, &
75	ONLY: restart_data_format_output
76
77	USE indices, &
78	ONLY: nxl, nxlg, nxr, nxrg, nys, nysg, nyn, nyng, nzb, nzt
79
80	USE kinds
81
82	! USE restart_data_mpi_io_mod, &
83	! ONLY:
84
85	IMPLICIT NONE
86
87	INTEGER(iwp) :: idum !< dummy variable used to avoid compiler warnings about unused variables
88
89	LOGICAL :: dummy_logical = .FALSE. !< switch to avoid compiler warnings about unused variables
90	LOGICAL :: gust_module_enabled = .FALSE. !< switch, if the entire module is used at all
91
92	SAVE
93
94	PRIVATE
95
96	!
97	!-- Public functions
98	PUBLIC &
99	gust_parin, &
100	gust_check_parameters, &
101	gust_check_data_output_pr, &
102	gust_check_data_output, &
103	gust_init_arrays, &
104	gust_init, &
105	gust_define_netcdf_grid, &
106	gust_header, &
107	gust_actions, &
108	gust_prognostic_equations, &
109	gust_swap_timelevel, &
110	gust_3d_data_averaging, &
111	gust_data_output_2d, &
112	gust_data_output_3d, &
113	gust_statistics, &
114	gust_rrd_global, &
115	gust_wrd_global, &
116	gust_rrd_local, &
117	gust_wrd_local
118	!
119	!-- Public parameters, constants and initial values
120	PUBLIC &
121	gust_module_enabled
122
123
124	INTERFACE gust_parin
125	MODULE PROCEDURE gust_parin
126	END INTERFACE gust_parin
127
128	INTERFACE gust_check_parameters
129	MODULE PROCEDURE gust_check_parameters
130	END INTERFACE gust_check_parameters
131
132	INTERFACE gust_check_data_output_pr
133	MODULE PROCEDURE gust_check_data_output_pr
134	END INTERFACE gust_check_data_output_pr
135
136	INTERFACE gust_check_data_output
137	MODULE PROCEDURE gust_check_data_output
138	END INTERFACE gust_check_data_output
139
140	INTERFACE gust_init_arrays
141	MODULE PROCEDURE gust_init_arrays
142	END INTERFACE gust_init_arrays
143
144	INTERFACE gust_init
145	MODULE PROCEDURE gust_init
146	END INTERFACE gust_init
147
148	INTERFACE gust_define_netcdf_grid
149	MODULE PROCEDURE gust_define_netcdf_grid
150	END INTERFACE gust_define_netcdf_grid
151
152	INTERFACE gust_header
153	MODULE PROCEDURE gust_header
154	END INTERFACE gust_header
155
156	INTERFACE gust_actions
157	MODULE PROCEDURE gust_actions
158	MODULE PROCEDURE gust_actions_ij
159	END INTERFACE gust_actions
160
161	INTERFACE gust_prognostic_equations
162	MODULE PROCEDURE gust_prognostic_equations
163	MODULE PROCEDURE gust_prognostic_equations_ij
164	END INTERFACE gust_prognostic_equations
165
166	INTERFACE gust_swap_timelevel
167	MODULE PROCEDURE gust_swap_timelevel
168	END INTERFACE gust_swap_timelevel
169
170	INTERFACE gust_3d_data_averaging
171	MODULE PROCEDURE gust_3d_data_averaging
172	END INTERFACE gust_3d_data_averaging
173
174	INTERFACE gust_data_output_2d
175	MODULE PROCEDURE gust_data_output_2d
176	END INTERFACE gust_data_output_2d
177
178	INTERFACE gust_data_output_3d
179	MODULE PROCEDURE gust_data_output_3d
180	END INTERFACE gust_data_output_3d
181
182	INTERFACE gust_statistics
183	MODULE PROCEDURE gust_statistics
184	END INTERFACE gust_statistics
185
186	INTERFACE gust_rrd_global
187	MODULE PROCEDURE gust_rrd_global_ftn
188	MODULE PROCEDURE gust_rrd_global_mpi
189	END INTERFACE gust_rrd_global
190
191	INTERFACE gust_wrd_global
192	MODULE PROCEDURE gust_wrd_global
193	END INTERFACE gust_wrd_global
194
195	INTERFACE gust_rrd_local
196	MODULE PROCEDURE gust_rrd_local_ftn
197	MODULE PROCEDURE gust_rrd_local_mpi
198	END INTERFACE gust_rrd_local
199
200	INTERFACE gust_wrd_local
201	MODULE PROCEDURE gust_wrd_local
202	END INTERFACE gust_wrd_local
203
204	CONTAINS
205
206
207	!--------------------------------------------------------------------------------------------------!
208	! Description:
209	! ------------
210	!> Parin for &gust_parameters for gust module
211	!--------------------------------------------------------------------------------------------------!
212	SUBROUTINE gust_parin
213
214
215	IMPLICIT NONE
216
217	CHARACTER (LEN=80) :: line !< dummy string that contains the current line of the parameter file
218
219	NAMELIST /gust_parameters/ &
220	gust_module_enabled
221
222	line = ' '
223	!
224	!-- Try to find gust module package
225	REWIND ( 11 )
226	line = ' '
227	DO WHILE ( INDEX( line, '&gust_parameters' ) == 0 )
228	READ ( 11, '(A)', END=10 ) line
229	ENDDO
230	BACKSPACE ( 11 )
231	!
232	!-- Read user-defined namelist
233	READ ( 11, gust_parameters )
234	!
235	!-- Set flag that indicates that the gust module is switched on
236	gust_module_enabled = .TRUE.
237
238	10 CONTINUE
239
240
241	END SUBROUTINE gust_parin
242
243
244	!--------------------------------------------------------------------------------------------------!
245	! Description:
246	! ------------
247	!> Check parameters routine for gust module
248	!--------------------------------------------------------------------------------------------------!
249	SUBROUTINE gust_check_parameters
250
251
252	IMPLICIT NONE
253
254
255	END SUBROUTINE gust_check_parameters
256
257
258	!--------------------------------------------------------------------------------------------------!
259	! Description:
260	! ------------
261	!> Check data output of profiles for gust module
262	!--------------------------------------------------------------------------------------------------!
263	SUBROUTINE gust_check_data_output_pr( variable, var_count, unit, dopr_unit )
264
265
266	IMPLICIT NONE
267
268	CHARACTER (LEN=*) :: dopr_unit !< local value of dopr_unit
269	CHARACTER (LEN=*) :: unit !<
270	CHARACTER (LEN=*) :: variable !<
271
272	INTEGER(iwp) :: var_count !<
273
274	!
275	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
276	IF ( dummy_logical ) idum = LEN( unit ) + LEN( variable ) + LEN( dopr_unit ) + var_count
277
278	END SUBROUTINE gust_check_data_output_pr
279
280	!--------------------------------------------------------------------------------------------------!
281	! Description:
282	! ------------
283	!> Check data output for gust module
284	!--------------------------------------------------------------------------------------------------!
285	SUBROUTINE gust_check_data_output( var, unit )
286
287
288	IMPLICIT NONE
289
290	CHARACTER (LEN=*) :: unit !<
291	CHARACTER (LEN=*) :: var !<
292
293	!
294	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
295	IF ( dummy_logical ) idum = LEN( var ) + LEN( unit )
296
297	END SUBROUTINE gust_check_data_output
298
299
300	!--------------------------------------------------------------------------------------------------!
301	! Description:
302	! ------------
303	!> Allocate gust module arrays and define pointers
304	!--------------------------------------------------------------------------------------------------!
305	SUBROUTINE gust_init_arrays
306
307
308	IMPLICIT NONE
309
310
311	END SUBROUTINE gust_init_arrays
312
313
314	!--------------------------------------------------------------------------------------------------!
315	! Description:
316	! ------------
317	!> Initialization of the gust module
318	!--------------------------------------------------------------------------------------------------!
319	SUBROUTINE gust_init
320
321
322	IMPLICIT NONE
323
324
325	END SUBROUTINE gust_init
326
327
328	!--------------------------------------------------------------------------------------------------!
329	!
330	! Description:
331	! ------------
332	!> Subroutine defining appropriate grid for netcdf variables.
333	!> It is called out from subroutine netcdf.
334	!--------------------------------------------------------------------------------------------------!
335	SUBROUTINE gust_define_netcdf_grid( var, found, grid_x, grid_y, grid_z )
336
337
338	IMPLICIT NONE
339
340	CHARACTER (LEN=*), INTENT(IN) :: grid_x !<
341	CHARACTER (LEN=*), INTENT(IN) :: grid_y !<
342	CHARACTER (LEN=*), INTENT(IN) :: grid_z !<
343	CHARACTER (LEN=*), INTENT(IN) :: var !<
344
345	LOGICAL, INTENT(IN) :: found !<
346
347	!
348	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
349	IF ( found ) idum = LEN( var ) + LEN( grid_x ) + LEN( grid_y ) + LEN( grid_z )
350
351	END SUBROUTINE gust_define_netcdf_grid
352
353
354	!--------------------------------------------------------------------------------------------------!
355	! Description:
356	! ------------
357	!> Header output for gust module
358	!--------------------------------------------------------------------------------------------------!
359	SUBROUTINE gust_header ( io )
360
361
362	IMPLICIT NONE
363
364	INTEGER(iwp), INTENT(IN) :: io !< Unit of the output file
365
366	!
367	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
368	IF ( dummy_logical ) idum = io
369
370	END SUBROUTINE gust_header
371
372
373	!--------------------------------------------------------------------------------------------------!
374	! Description:
375	! ------------
376	!> Call for all grid points
377	!--------------------------------------------------------------------------------------------------!
378	SUBROUTINE gust_actions( location )
379
380
381	IMPLICIT NONE
382
383	CHARACTER (LEN=*) :: location !<
384
385	!
386	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
387	IF ( dummy_logical ) idum = LEN( location )
388
389	END SUBROUTINE gust_actions
390
391
392	!--------------------------------------------------------------------------------------------------!
393	! Description:
394	! ------------
395	!> Call for grid point i,j
396	!--------------------------------------------------------------------------------------------------!
397	SUBROUTINE gust_actions_ij( i, j, location )
398
399
400	IMPLICIT NONE
401
402	CHARACTER (LEN=*) :: location
403
404	INTEGER(iwp) :: i
405	INTEGER(iwp) :: j
406
407	!
408	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
409	IF ( dummy_logical ) idum = i + j + LEN( location )
410
411	END SUBROUTINE gust_actions_ij
412
413
414	!--------------------------------------------------------------------------------------------------!
415	! Description:
416	! ------------
417	!> Call for all grid points
418	!--------------------------------------------------------------------------------------------------!
419	SUBROUTINE gust_prognostic_equations()
420
421	!
422	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
423	IF ( dummy_logical ) idum = 1
424
425	END SUBROUTINE gust_prognostic_equations
426
427
428	!--------------------------------------------------------------------------------------------------!
429	! Description:
430	! ------------
431	!> Call for grid point i,j
432	!--------------------------------------------------------------------------------------------------!
433	SUBROUTINE gust_prognostic_equations_ij( i, j, i_omp_start, tn )
434
435
436	INTEGER(iwp), INTENT(IN) :: i !< grid index in x-direction
437	INTEGER(iwp), INTENT(IN) :: i_omp_start !< first loop index of i-loop in prognostic_equations
438	INTEGER(iwp), INTENT(IN) :: j !< grid index in y-direction
439	INTEGER(iwp), INTENT(IN) :: tn !< task number of openmp task
440
441	!
442	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
443	IF ( dummy_logical ) idum = i + j + i_omp_start + tn
444
445	END SUBROUTINE gust_prognostic_equations_ij
446
447
448	!--------------------------------------------------------------------------------------------------!
449	! Description:
450	! ------------
451	!> Swapping of timelevels
452	!--------------------------------------------------------------------------------------------------!
453	SUBROUTINE gust_swap_timelevel ( mod_count )
454
455
456	IMPLICIT NONE
457
458	INTEGER, INTENT(IN) :: mod_count
459
460	!
461	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
462	IF ( dummy_logical ) idum = mod_count
463
464	END SUBROUTINE gust_swap_timelevel
465
466
467	!--------------------------------------------------------------------------------------------------!
468	!
469	! Description:
470	! ------------
471	!> Subroutine for averaging 3D data
472	!--------------------------------------------------------------------------------------------------!
473	SUBROUTINE gust_3d_data_averaging( mode, variable )
474
475
476	IMPLICIT NONE
477
478	CHARACTER (LEN=*) :: mode !<
479	CHARACTER (LEN=*) :: variable !<
480
481	!
482	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
483	IF ( dummy_logical ) idum = LEN( mode ) + LEN( variable )
484
485	END SUBROUTINE gust_3d_data_averaging
486
487	!--------------------------------------------------------------------------------------------------!
488	!
489	! Description:
490	! ------------
491	!> Subroutine defining 2D output variables
492	!--------------------------------------------------------------------------------------------------!
493	SUBROUTINE gust_data_output_2d( av, variable, found, grid, mode, local_pf, two_d, nzb_do, &
494	nzt_do, fill_value )
495
496
497	IMPLICIT NONE
498
499	CHARACTER (LEN=*), INTENT(INOUT) :: grid !< name of vertical grid
500	CHARACTER (LEN=*), INTENT(IN) :: mode !< either 'xy', 'xz' or 'yz'
501	CHARACTER (LEN=*), INTENT(IN) :: variable !< name of variable
502
503	INTEGER(iwp), INTENT(IN) :: av !< flag for (non-)average output
504	INTEGER(iwp), INTENT(IN) :: nzb_do !< vertical output index (bottom)
505	INTEGER(iwp), INTENT(IN) :: nzt_do !< vertical output index (top)
506
507	LOGICAL, INTENT(INOUT) :: found !< flag if output variable is found
508	LOGICAL, INTENT(INOUT) :: two_d !< flag parameter that indicates 2D variables (horizontal cross sections)
509
510	REAL(wp), INTENT(IN) :: fill_value !< value for the _FillValue attribute
511
512	REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) :: local_pf !< local
513	!< array to which output data is resorted to
514
515	!
516	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
517	IF ( found .AND. two_d ) THEN
518	idum = av + LEN( variable ) + LEN( grid // mode ) + local_pf(nxl,nys,nzb_do) + fill_value
519	ENDIF
520
521	END SUBROUTINE gust_data_output_2d
522
523
524	!--------------------------------------------------------------------------------------------------!
525	!
526	! Description:
527	! ------------
528	!> Subroutine defining 3D output variables
529	!--------------------------------------------------------------------------------------------------!
530	SUBROUTINE gust_data_output_3d( av, variable, found, local_pf, fill_value, nzb_do, nzt_do )
531
532
533	IMPLICIT NONE
534
535	CHARACTER (LEN=*), INTENT(IN) :: variable !< name of variable
536
537	INTEGER(iwp), INTENT(IN) :: av !< flag for (non-)average output
538	INTEGER(iwp), INTENT(IN) :: nzb_do !< lower limit of the data output (usually 0)
539	INTEGER(iwp), INTENT(IN) :: nzt_do !< vertical upper limit of the data output (usually nz_do3d)
540
541	LOGICAL, INTENT(INOUT) :: found !< flag if output variable is found
542
543	REAL(wp), INTENT(IN) :: fill_value !< value for the _FillValue attribute
544
545	REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do), INTENT(INOUT) :: local_pf !< local
546	!< array to which output data is resorted to
547
548	!
549	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
550	IF ( found ) idum = av + LEN( variable ) + fill_value + local_pf(nxl,nys,nzb_do)
551
552	END SUBROUTINE gust_data_output_3d
553
554
555	!--------------------------------------------------------------------------------------------------!
556	! Description:
557	! ------------
558	!> This routine computes profile and timeseries data for the gust module.
559	!--------------------------------------------------------------------------------------------------!
560	SUBROUTINE gust_statistics( mode, sr, tn, dots_max )
561
562
563	IMPLICIT NONE
564
565	CHARACTER (LEN=*) :: mode !<
566
567	INTEGER(iwp) :: dots_max !<
568	INTEGER(iwp) :: sr !<
569	INTEGER(iwp) :: tn !<
570
571	!
572	!-- Next line is just to avoid compiler warnings about unused variables. You may remove it.
573	IF ( dummy_logical ) idum = dots_max + sr + tn + LEN( mode )
574
575	END SUBROUTINE gust_statistics
576
577
578	!--------------------------------------------------------------------------------------------------!
579	! Description:
580	! ------------
581	!> Read module-specific global restart data (Fortran binary format).
582	!--------------------------------------------------------------------------------------------------!
583	SUBROUTINE gust_rrd_global_ftn( found )
584
585
586	USE control_parameters, &
587	ONLY: length, restart_string
588
589
590	IMPLICIT NONE
591
592	LOGICAL, INTENT(OUT) :: found
593
594
595	found = .TRUE.
596
597
598	SELECT CASE ( restart_string(1:length) )
599
600	CASE ( 'global_paramter' )
601	! READ ( 13 ) global_parameter
602
603	CASE DEFAULT
604
605	found = .FALSE.
606
607	END SELECT
608
609
610	END SUBROUTINE gust_rrd_global_ftn
611
612
613	!--------------------------------------------------------------------------------------------------!
614	! Description:
615	! ------------
616	!> Read module-specific global restart data (MPI-IO).
617	!--------------------------------------------------------------------------------------------------!
618	SUBROUTINE gust_rrd_global_mpi
619
620
621	! CALL rrd_mpi_io( 'global_parameter', global_parameter )
622	! READ ( 13 ) global_parameter
623
624	END SUBROUTINE gust_rrd_global_mpi
625
626
627	!--------------------------------------------------------------------------------------------------!
628	! Description:
629	! ------------
630	!> Read module-specific local restart data arrays (Fortran binary format).
631	!--------------------------------------------------------------------------------------------------!
632	SUBROUTINE gust_rrd_local_ftn( k, nxlf, nxlc, nxl_on_file, nxrf, nxrc, nxr_on_file, nynf, nync,&
633	nyn_on_file, nysf, nysc, nys_on_file, tmp_2d, tmp_3d, found )
634
635
636	USE control_parameters
637
638	USE indices
639
640	USE kinds
641
642	USE pegrid
643
644
645	IMPLICIT NONE
646
647	INTEGER(iwp) :: k !<
648	INTEGER(iwp) :: nxlc !<
649	INTEGER(iwp) :: nxlf !<
650	INTEGER(iwp) :: nxl_on_file !<
651	INTEGER(iwp) :: nxrc !<
652	INTEGER(iwp) :: nxrf !<
653	INTEGER(iwp) :: nxr_on_file !<
654	INTEGER(iwp) :: nync !<
655	INTEGER(iwp) :: nynf !<
656	INTEGER(iwp) :: nyn_on_file !<
657	INTEGER(iwp) :: nysc !<
658	INTEGER(iwp) :: nysf !<
659	INTEGER(iwp) :: nys_on_file !<
660
661	LOGICAL, INTENT(OUT) :: found
662
663	REAL(wp), DIMENSION(nys_on_file-nbgp:nyn_on_file+nbgp,nxl_on_file-nbgp:nxr_on_file+nbgp) :: tmp_2d !<
664	REAL(wp), DIMENSION(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,nxl_on_file-nbgp:nxr_on_file+nbgp) :: tmp_3d !<
665
666
667	!
668	!-- Next lins are just to avoid compiler warnings about unused variables in case of empty user interface routine.
669	!-- You may remove them.
670	IF ( dummy_logical ) THEN
671	idum = k + nxlc + nxlf + nxrc + nxrf + nync + nynf + nysc + nysf + &
672	tmp_2d(nys_on_file,nxl_on_file) + tmp_3d(nzb,nys_on_file,nxl_on_file)
673	ENDIF
674
675	!
676	!-- Here the reading of user-defined restart data follows:
677	!-- Sample for user-defined output
678	found = .TRUE.
679
680	SELECT CASE ( restart_string(1:length) )
681
682	CASE ( '.......' )
683	! IF ( .NOT. ALLOCATED( u2_av ) ) THEN
684	! ALLOCATE( u2_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
685	! ENDIF
686	! IF ( k == 1 ) READ ( 13 ) tmp_3d
687	! u2_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
688	! tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
689	!
690	CASE DEFAULT
691
692	found = .FALSE.
693
694	END SELECT
695
696
697	END SUBROUTINE gust_rrd_local_ftn
698
699
700	!--------------------------------------------------------------------------------------------------!
701	! Description:
702	! ------------
703	!> Read module-specific local restart data arrays (MPI-IO).
704	!--------------------------------------------------------------------------------------------------!
705	SUBROUTINE gust_rrd_local_mpi
706
707
708	! CALL rrd_mpi_io( 'local_array', local_array )
709
710	END SUBROUTINE gust_rrd_local_mpi
711
712
713	!--------------------------------------------------------------------------------------------------!
714	! Description:
715	! ------------
716	!> This routine writes the respective restart data for the gust module.
717	!--------------------------------------------------------------------------------------------------!
718	SUBROUTINE gust_wrd_global
719
720
721	IMPLICIT NONE
722
723
724	IF ( TRIM( restart_data_format_output ) == 'fortran_binary' ) THEN
725
726	! CALL wrd_write_string( 'global_parameter' )
727	! WRITE ( 14 ) global_parameter
728
729	! IF ( ALLOCATED( inflow_damping_factor ) ) THEN
730	! CALL wrd_write_string( 'inflow_damping_factor' )
731	! WRITE ( 14 ) inflow_damping_factor
732	! ENDIF
733
734	ELSEIF ( restart_data_format_output(1:3) == 'mpi' ) THEN
735
736	! CALL wrd_mpi_io( 'global_parameter', global_parameter )
737	! IF ( ALLOCATED( inflow_damping_factor ) ) THEN
738	! CALL wrd_mpi_io_global_array( 'inflow_damping_factor', inflow_damping_factor )
739	! ENDIF
740
741	ENDIF
742
743	END SUBROUTINE gust_wrd_global
744
745
746	!--------------------------------------------------------------------------------------------------!
747	! Description:
748	! ------------
749	!> This routine writes the respective restart data for the gust module.
750	!--------------------------------------------------------------------------------------------------!
751	SUBROUTINE gust_wrd_local
752
753
754	IMPLICIT NONE
755
756
757	IF ( TRIM( restart_data_format_output ) == 'fortran_binary' ) THEN
758
759	! IF ( ALLOCATED( u2_av ) ) THEN
760	! CALL wrd_write_string( 'u2_av' )
761	! WRITE ( 14 ) u2_av
762	! ENDIF
763
764	ELSEIF ( restart_data_format_output(1:3) == 'mpi' ) THEN
765
766	! IF ( ALLOCATED( u2_av ) ) CALL wrd_mpi_io( 'u2_av', u2_av )
767
768	ENDIF
769
770	END SUBROUTINE gust_wrd_local
771
772
773
774	END MODULE gust_mod

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

| Impressum | ©Leibniz Universität Hannover |