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

Last change on this file since 4784 was 4768, checked in by suehring, 3 years ago

Enable 3D data output also with 64-bit precision

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