source: palm/trunk/SOURCE/user_module.f90 @ 4504

Last change on this file since 4504 was 4504, checked in by raasch, 4 years ago

file re-formatted to follow the PALM coding standard, hint for setting rmask arrays added

  • Property svn:keywords set to Id
File size: 40.9 KB
RevLine 
[1873]1!> @file user_module.f90
[4497]2!--------------------------------------------------------------------------------------------------!
[2696]3! This file is part of the PALM model system.
[1036]4!
[4497]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.
[1036]8!
[4497]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.
[1036]12!
[4497]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/>.
[1036]15!
[4360]16! Copyright 1997-2020 Leibniz Universitaet Hannover
[4497]17!--------------------------------------------------------------------------------------------------!
[1036]18!
[4497]19!
[484]20! Current revisions:
[211]21! -----------------
[1683]22!
[4504]23!
[1321]24! Former revisions:
25! -----------------
26! $Id: user_module.f90 4504 2020-04-20 12:11:24Z raasch $
[4504]27! hint for setting rmask arrays added
28!
29! 4497 2020-04-15 10:20:51Z raasch
[4497]30! file re-formatted to follow the PALM coding standard
31!
32! 4495 2020-04-13 20:11:20Z raasch
[4495]33! restart data handling with MPI-IO added
[4497]34!
[4495]35! 4360 2020-01-07 11:25:50Z suehring
[4497]36! Introduction of wall_flags_total_0, which currently sets bits based on static topography
37! information used in wall_flags_static_0
38!
[4346]39! 4329 2019-12-10 15:46:36Z motisi
[4329]40! Renamed wall_flags_0 to wall_flags_static_0
[4497]41!
[4329]42! 4287 2019-11-01 14:50:20Z raasch
[4287]43! reading of namelist file and actions in case of namelist errors revised so that statement labels
[4497]44! and goto statements are not required any more; this revision also removes a previous bug which
45! appeared when the namelist has been commented out in the namelist file
46!
[4287]47! 4182 2019-08-22 15:20:23Z scharf
[4182]48! Corrected "Former revisions" section
[4497]49!
[4182]50! 3986 2019-05-20 14:08:14Z Giersch
[3986]51! Redundant integration of control parameters in user_rrd_global removed
[4497]52!
[3986]53! 3911 2019-04-17 12:26:19Z knoop
[3911]54! Bugfix: added before_prognostic_equations case in user_actions
[4497]55!
[3911]56! 3768 2019-02-27 14:35:58Z raasch
[3768]57! variables commented + statements added to avoid compiler warnings about unused variables
58!
59! 3767 2019-02-27 08:18:02Z raasch
[3767]60! unused variable for file index removed from rrd-subroutines parameter list
[4497]61!
[3767]62! 3747 2019-02-16 15:15:23Z gronemeier
[3747]63! Add routine user_init_arrays
[4497]64!
[3747]65! 3703 2019-01-29 16:43:53Z knoop
[2894]66! An example for a user defined global variable has been added (Giersch)
[1321]67!
[4182]68! Revision 1.1  1998/03/24 15:29:04  raasch
69! Initial revision
70!
71!
[211]72! Description:
73! ------------
[4497]74!> Declaration of user-defined variables. This module may only be used in the user-defined routines
75!> (contained in user_interface.f90).
76!--------------------------------------------------------------------------------------------------!
[1682]77 MODULE user
[211]78
[3687]79    USE arrays_3d
80
81    USE control_parameters
82
83    USE cpulog
84
85    USE indices
86
[1320]87    USE kinds
[211]88
[3687]89    USE pegrid
90
91    USE statistics
92
93    USE surface_mod
94
[1320]95    IMPLICIT NONE
96
[4497]97    INTEGER(iwp) ::  dots_num_palm      !<
98    INTEGER(iwp) ::  dots_num_user = 0  !<
99    INTEGER(iwp) ::  user_idummy        !<
[1320]100
[4497]101    LOGICAL ::  user_module_enabled = .FALSE.  !<
102
103    REAL(wp) ::  user_rdummy  !<
104
[211]105!
106!-- Sample for user-defined output
[4497]107!    REAL(wp) :: global_parameter  !< user defined global parameter
[2894]108!
[4497]109!    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  u2      !< user defined array
110!    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  u2_av   !< user defined array
111!    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  ustvst  !< user defined array
[211]112
113    SAVE
114
[3701]115    PRIVATE
[3687]116
117!
118!- Public functions
[4497]119    PUBLIC                                                                                         &
120       user_actions,                                                                               &
121       user_check_data_output,                                                                     &
122       user_check_data_output_pr,                                                                  &
123       user_check_data_output_ts,                                                                  &
124       user_check_parameters,                                                                      &
125       user_data_output_2d,                                                                        &
126       user_data_output_3d,                                                                        &
127       user_define_netcdf_grid,                                                                    &
128       user_header,                                                                                &
129       user_init,                                                                                  &
130       user_init_arrays,                                                                           &
131       user_last_actions,                                                                          &
132       user_parin,                                                                                 &
133       user_rrd_global,                                                                            &
134       user_rrd_local,                                                                             &
135       user_statistics,                                                                            &
136       user_3d_data_averaging,                                                                     &
137       user_wrd_global,                                                                            &
138       user_wrd_local
[3687]139
[4497]140
[3687]141!
142!- Public parameters, constants and initial values
[4497]143   PUBLIC                                                                                          &
[3701]144      user_module_enabled
[3687]145
146    INTERFACE user_parin
147       MODULE PROCEDURE user_parin
148    END INTERFACE user_parin
149
150    INTERFACE user_check_parameters
151       MODULE PROCEDURE user_check_parameters
152    END INTERFACE user_check_parameters
153
[3700]154    INTERFACE user_check_data_output_ts
155       MODULE PROCEDURE user_check_data_output_ts
156    END INTERFACE user_check_data_output_ts
157
[3687]158    INTERFACE user_check_data_output_pr
159       MODULE PROCEDURE user_check_data_output_pr
160    END INTERFACE user_check_data_output_pr
161
162    INTERFACE user_check_data_output
163       MODULE PROCEDURE user_check_data_output
164    END INTERFACE user_check_data_output
165
[3701]166    INTERFACE user_define_netcdf_grid
167       MODULE PROCEDURE user_define_netcdf_grid
168    END INTERFACE user_define_netcdf_grid
169
[3687]170    INTERFACE user_init
171       MODULE PROCEDURE user_init
172    END INTERFACE user_init
173
[3747]174    INTERFACE user_init_arrays
175       MODULE PROCEDURE user_init_arrays
176    END INTERFACE user_init_arrays
177
[3687]178    INTERFACE user_header
179       MODULE PROCEDURE user_header
180    END INTERFACE user_header
181
182    INTERFACE user_actions
183       MODULE PROCEDURE user_actions
184       MODULE PROCEDURE user_actions_ij
185    END INTERFACE user_actions
186
187    INTERFACE user_3d_data_averaging
188       MODULE PROCEDURE user_3d_data_averaging
189    END INTERFACE user_3d_data_averaging
190
191    INTERFACE user_data_output_2d
192       MODULE PROCEDURE user_data_output_2d
193    END INTERFACE user_data_output_2d
194
195    INTERFACE user_data_output_3d
196       MODULE PROCEDURE user_data_output_3d
197    END INTERFACE user_data_output_3d
198
199    INTERFACE user_statistics
200       MODULE PROCEDURE user_statistics
201    END INTERFACE user_statistics
202
203    INTERFACE user_rrd_global
[4495]204       MODULE PROCEDURE user_rrd_global_ftn
205       MODULE PROCEDURE user_rrd_global_mpi
[3687]206    END INTERFACE user_rrd_global
207
208    INTERFACE user_rrd_local
209       MODULE PROCEDURE user_rrd_local
210    END INTERFACE user_rrd_local
211
212    INTERFACE user_wrd_global
213       MODULE PROCEDURE user_wrd_global
214    END INTERFACE user_wrd_global
215
216    INTERFACE user_wrd_local
217       MODULE PROCEDURE user_wrd_local
218    END INTERFACE user_wrd_local
219
[3701]220    INTERFACE user_last_actions
221       MODULE PROCEDURE user_last_actions
222    END INTERFACE user_last_actions
[3687]223
[3701]224
[3687]225 CONTAINS
226
227
[4497]228!--------------------------------------------------------------------------------------------------!
[3687]229! Description:
230! ------------
231!> Parin for &user_parameters for user module
[4497]232!--------------------------------------------------------------------------------------------------!
[3687]233 SUBROUTINE user_parin
234
[4497]235    CHARACTER (LEN=80) ::  line  !< string containing the last line read from namelist file
[3687]236
[4497]237    INTEGER(iwp) ::  i          !<
238    INTEGER(iwp) ::  io_status  !< status after reading the namelist file
239    INTEGER(iwp) ::  j          !<
[3687]240
241
[4287]242    NAMELIST /user_parameters/                                                                     &
243       data_output_masks_user,                                                                     &
244       data_output_pr_user,                                                                        &
245       data_output_user,                                                                           &
246       region
[3687]247
248!
[3768]249!-- Next statement is to avoid compiler warnings about unused variables. Please remove in case
250!-- that you are using them.
251    IF ( dots_num_palm == 0  .OR.  dots_num_user == 0  .OR.  user_idummy == 0  .OR.                &
252         user_rdummy == 0.0_wp )  CONTINUE
253
254!
[4497]255!-- Set revision number of this default interface version. It will be checked within the main
256!-- program (palm). Please change the revision number in case that the current revision does not
257!-- match with previous revisions (e.g. if routines have been added/deleted or if parameter lists
258!-- in subroutines have been changed).
[4495]259    user_interface_current_revision = 'r4495'
[3687]260
261!
[4497]262!-- Position the namelist-file at the beginning (it has already been opened in parin), and try to
263!-- read (find) a namelist named "user_parameters".
[3687]264    REWIND ( 11 )
[4287]265    READ( 11, user_parameters, IOSTAT=io_status )
[3687]266
[4287]267!
268!-- Actions depending on the READ status
269    IF ( io_status == 0 )  THEN
270!
271!--    User namelist found and correctly read. Set default module switch to true. This activates
272!--    calls of the user-interface subroutines.
273       user_module_enabled = .TRUE.
[3687]274
[4287]275    ELSEIF ( io_status > 0 )  THEN
276!
277!--    User namelist was found, but contained errors. Print an error message containing the line
278!--    that caused the problem
279       BACKSPACE( 11 )
280       READ( 11 , '(A)') line
281       CALL parin_fail_message( 'user_parameters', line )
[3701]282
[4287]283    ENDIF
[3687]284
285!
[4497]286!-- Determine the number of user-defined profiles and append them to the standard data output
287!-- (data_output_pr)
[3701]288    IF ( user_module_enabled )  THEN
[3687]289       IF ( data_output_pr_user(1) /= ' ' )  THEN
290          i = 1
291          DO WHILE ( data_output_pr(i) /= ' '  .AND.  i <= 100 )
292             i = i + 1
293          ENDDO
294          j = 1
295          DO WHILE ( data_output_pr_user(j) /= ' '  .AND.  j <= 100 )
296             data_output_pr(i) = data_output_pr_user(j)
297             max_pr_user_tmp   = max_pr_user_tmp + 1
298             i = i + 1
299             j = j + 1
300          ENDDO
301       ENDIF
302    ENDIF
303
[3701]304
[3687]305 END SUBROUTINE user_parin
306
307
[4497]308!--------------------------------------------------------------------------------------------------!
[3687]309! Description:
310! ------------
311!> Check &userpar control parameters and deduce further quantities.
[4497]312!--------------------------------------------------------------------------------------------------!
[3687]313 SUBROUTINE user_check_parameters
314
[4497]315!
316!-- Here the user may add code to check the validity of further &userpar control parameters or
317!-- deduce further quantities.
[3687]318
319
320 END SUBROUTINE user_check_parameters
321
322
[4497]323!--------------------------------------------------------------------------------------------------!
[3687]324! Description:
325! ------------
[3700]326!> Set module-specific timeseries units and labels
[4497]327!--------------------------------------------------------------------------------------------------!
[3700]328 SUBROUTINE user_check_data_output_ts( dots_max, dots_num, dots_label, dots_unit )
329
[4497]330    INTEGER(iwp),      INTENT(IN)     ::  dots_max  !<
331    INTEGER(iwp),      INTENT(INOUT)  ::  dots_num  !<
[3700]332
[4497]333    CHARACTER(LEN=*), DIMENSION(dots_max), INTENT(INOUT)  ::  dots_label  !<
334    CHARACTER(LEN=*), DIMENSION(dots_max), INTENT(INOUT)  ::  dots_unit   !<
[3700]335
[3768]336!
337!-- Next line is to avoid compiler warning about unused variables. Please remove.
338    IF ( dots_num == 0  .OR.  dots_label(1)(1:1) == ' '  .OR.  dots_unit(1)(1:1) == ' ' )  CONTINUE
[3700]339
[3768]340!
[4497]341!-- Sample for user-defined time series:
342!-- For each time series quantity you have to give a label and a unit, which will be used for the
343!-- NetCDF file. They must not contain more than seven characters. The value of dots_num has to be
344!-- increased by the number of new time series quantities. Its old value has to be stored in
345!-- dots_num_palm. See routine user_statistics on how to calculate and output these quantities.
[3700]346
347!    dots_num_palm = dots_num
348
349!    dots_num = dots_num + 1
350!    dots_num_user = dots_num_user + 1
351!    dots_label(dots_num) = 'abs_umx'
352!    dots_unit(dots_num)  = 'm/s'
353
354!    dots_num = dots_num + 1
355!    dots_num_user = dots_num_user + 1
356!    dots_label(dots_num) = 'abs_vmx'
357!    dots_unit(dots_num)  = 'm/s'
358
359
360 END SUBROUTINE user_check_data_output_ts
361
362
[4497]363!--------------------------------------------------------------------------------------------------!
[3700]364! Description:
365! ------------
[4497]366!> Set the unit of user defined profile output quantities. For those variables not recognized by the
367!> user, the parameter unit is set to "illegal", which tells the calling routine that the
368!> output variable is not defined and leads to a program abort.
369!--------------------------------------------------------------------------------------------------!
[3700]370 SUBROUTINE user_check_data_output_pr( variable, var_count, unit, dopr_unit )
[3687]371
372
373    USE profil_parameter
374
375
[4497]376    CHARACTER (LEN=*) ::  unit      !<
377    CHARACTER (LEN=*) ::  variable  !<
[3700]378    CHARACTER (LEN=*) ::  dopr_unit !< local value of dopr_unit
[3687]379
[4497]380!    INTEGER(iwp) ::  user_pr_index  !<
381    INTEGER(iwp) ::  var_count      !<
[3687]382
[3768]383!
384!-- Next line is to avoid compiler warning about unused variables. Please remove.
385    IF ( unit(1:1) == ' '  .OR.  dopr_unit(1:1) == ' '  .OR.  var_count == 0 )  CONTINUE
386
[3687]387    SELECT CASE ( TRIM( variable ) )
388
389!
390!--    Uncomment and extend the following lines, if necessary.
[4497]391!--    Add additional CASE statements depending on the number of quantities for which profiles are
392!--    to be calculated. The respective calculations to be performed have to be added in routine
393!--    user_statistics. The quantities are (internally) identified by a user-profile-number
394!--    (see variable "user_pr_index" below). The first user-profile must be assigned the number
395!--    "pr_palm+1", the second one "pr_palm+2", etc. The respective user-profile-numbers have also
396!--    to be used in routine user_statistics!
397!       CASE ( 'u*v*' )                      ! quantity string as given in data_output_pr_user
[3687]398!          user_pr_index = pr_palm + 1
399!          dopr_index(var_count)  = user_pr_index    ! quantities' user-profile-number
[3700]400!          dopr_unit = 'm2/s2'  ! quantity unit
401!          unit = dopr_unit
[4497]402!          hom(:,2,user_pr_index,:) = SPREAD( zu, 2, statistic_regions+1 )
403!                                            ! grid on which the quantity is defined (use zu or zw)
404!
[3687]405
406       CASE DEFAULT
407          unit = 'illegal'
408
409    END SELECT
410
411
412 END SUBROUTINE user_check_data_output_pr
413
414
[4497]415!--------------------------------------------------------------------------------------------------!
[3687]416! Description:
417! ------------
[4497]418!> Set the unit of user defined output quantities. For those variables not recognized by the user,
419!> the parameter unit is set to "illegal", which tells the calling routine that the output variable
420!> is not defined and leads to a program abort.
421!--------------------------------------------------------------------------------------------------!
[3687]422 SUBROUTINE user_check_data_output( variable, unit )
423
424
[4497]425    CHARACTER (LEN=*) ::  unit      !<
426    CHARACTER (LEN=*) ::  variable  !<
[3687]427
428
429    SELECT CASE ( TRIM( variable ) )
430
431!
432!--    Uncomment and extend the following lines, if necessary
433!       CASE ( 'u2' )
434!          unit = 'm2/s2'
435!
436!       CASE ( 'u*v*' )
437!          unit = 'm2/s2'
438!
439       CASE DEFAULT
440          unit = 'illegal'
441
442    END SELECT
443
444
445 END SUBROUTINE user_check_data_output
446
447
[4497]448!--------------------------------------------------------------------------------------------------!
[3687]449! Description:
450! ------------
[3747]451!> Initialize user-defined arrays
[4497]452!--------------------------------------------------------------------------------------------------!
[3747]453 SUBROUTINE user_init_arrays
454
455
[3768]456!    INTEGER(iwp) :: i       !< loop index
457!    INTEGER(iwp) :: j       !< loop index
458!    INTEGER(iwp) :: region  !< index for loop over statistic regions
[3747]459
460!
461!-- Allocate user-defined arrays and set flags for statistic regions.
462!-- Sample for user-defined output
463!    ALLOCATE( u2(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
464!    ALLOCATE( ustvst(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
465
466!
467!-- Example for defining a statistic region:
[4504]468!-- ATTENTION: rmask = 0 is required at the ghost boundaries to guarantee correct statistic
469!--            evaluations (otherwise ghost points would be counted twice). This setting has
470!--            already been cared for in routine init_3d_model. Please don't set the ghost points
471!--            /= 0. i.e. run the following loop only over nxl,nxr and nys,nyn.
[3747]472!     IF ( statistic_regions >= 1 )  THEN
473!        region = 1
[4497]474!
[3747]475!        rmask(:,:,region) = 0.0_wp
476!        DO  i = nxl, nxr
477!           IF ( i >= INT( 0.25 * nx ) .AND. i <= INT( 0.75 * nx ) )  THEN
478!              DO  j = nys, nyn
479!                 IF ( i >= INT( 0.25 * ny ) .AND. i <= INT( 0.75 * ny ) )  THEN
480!                    rmask(j,i,region) = 1.0_wp
481!                 ENDIF
482!              ENDDO
483!           ENDIF
484!        ENDDO
[4497]485!
[3747]486!     ENDIF
487
488 END SUBROUTINE user_init_arrays
489
490
[4497]491!--------------------------------------------------------------------------------------------------!
[3747]492! Description:
493! ------------
[3687]494!> Execution of user-defined initializing actions
[4497]495!--------------------------------------------------------------------------------------------------!
[3687]496 SUBROUTINE user_init
497
498
[4497]499!    CHARACTER(LEN=20) :: field_char  !<
[3687]500!
501!-- Here the user-defined initializing actions follow:
502!-- Sample for user-defined output
[3747]503!    ustvst = 0.0_wp
[3687]504
505
506 END SUBROUTINE user_init
507
508
[4497]509!--------------------------------------------------------------------------------------------------!
[3687]510! Description:
511! ------------
[4497]512!> Set the grids on which user-defined output quantities are defined. Allowed values for grid_x are
513!> "x" and "xu", for grid_y "y" and "yv", and for grid_z "zu" and "zw".
514!--------------------------------------------------------------------------------------------------!
[3700]515 SUBROUTINE user_define_netcdf_grid( variable, found, grid_x, grid_y, grid_z )
516
517
518    CHARACTER (LEN=*) ::  grid_x     !<
519    CHARACTER (LEN=*) ::  grid_y     !<
520    CHARACTER (LEN=*) ::  grid_z     !<
521    CHARACTER (LEN=*) ::  variable   !<
522
523    LOGICAL ::  found   !<
524
525
526    SELECT CASE ( TRIM( variable ) )
527
528!
529!--    Uncomment and extend the following lines, if necessary
530!       CASE ( 'u2', 'u2_xy', 'u2_xz', 'u2_yz' )
531!          found  = .TRUE.
532!          grid_x = 'xu'
533!          grid_y = 'y'
534!          grid_z = 'zu'
535
536!       CASE ( 'u*v*', 'u*v*_xy', 'u*v*_xz', 'u*v*_yz' )
537!          found  = .TRUE.
538!          grid_x = 'x'
539!          grid_y = 'y'
540!          grid_z = 'zu'
541
542       CASE DEFAULT
543          found  = .FALSE.
544          grid_x = 'none'
545          grid_y = 'none'
546          grid_z = 'none'
547
548    END SELECT
549
550
551 END SUBROUTINE user_define_netcdf_grid
552
553
554
555
[4497]556!--------------------------------------------------------------------------------------------------!
[3700]557! Description:
558! ------------
[3687]559!> Print a header with user-defined information.
[4497]560!--------------------------------------------------------------------------------------------------!
[3687]561 SUBROUTINE user_header( io )
562
563
[4497]564    INTEGER(iwp) ::  i   !<
565    INTEGER(iwp) ::  io  !<
[3687]566
567!
[4497]568!-- If no user-defined variables are read from the namelist-file, no information will be printed.
[3701]569    IF ( .NOT. user_module_enabled )  THEN
[3687]570       WRITE ( io, 100 )
571       RETURN
572    ENDIF
573
574!
575!-- Printing the information.
576    WRITE ( io, 110 )
577
578    IF ( statistic_regions /= 0 )  THEN
579       WRITE ( io, 200 )
580       DO  i = 0, statistic_regions
581          WRITE ( io, 201 )  i, region(i)
582       ENDDO
583    ENDIF
584
585!
586!-- Format-descriptors
587100 FORMAT (//' *** no user-defined variables found'/)
[4497]588110 FORMAT (//1X,78('#') // ' User-defined variables and actions:' /                               &
589            ' -----------------------------------'//)
[3687]590200 FORMAT (' Output of profiles and time series for following regions:' /)
591201 FORMAT (4X,'Region ',I1,':   ',A)
592
593
594 END SUBROUTINE user_header
595
596
[4497]597!--------------------------------------------------------------------------------------------------!
[3687]598! Description:
599! ------------
600!> Call for all grid points
[4497]601!--------------------------------------------------------------------------------------------------!
[3687]602 SUBROUTINE user_actions( location )
603
604
[4497]605    CHARACTER(LEN=*) ::  location  !<
[3687]606
[4497]607!    INTEGER(iwp) ::  i  !<
608!    INTEGER(iwp) ::  j  !<
609!    INTEGER(iwp) ::  k  !<
[3687]610
611    CALL cpu_log( log_point(24), 'user_actions', 'start' )
612
613!
[4497]614!-- Here the user-defined actions follow. No calls for single grid points are allowed at locations
615!-- before and after the timestep, since these calls are not within an i,j-loop
[3687]616    SELECT CASE ( location )
617
618       CASE ( 'before_timestep' )
619!
620!--       Enter actions to be done before every timestep here
621
[3911]622       CASE ( 'before_prognostic_equations' )
623!
624!--       Enter actions to be done before all prognostic equations here
[3687]625
626       CASE ( 'after_integration' )
627!
[4497]628!--       Enter actions to be done after every time integration (before data output)
[3687]629!--       Sample for user-defined output:
630!          DO  i = nxlg, nxrg
631!             DO  j = nysg, nyng
632!                DO  k = nzb, nzt
633!                   u2(k,j,i) = u(k,j,i)**2
634!                ENDDO
635!             ENDDO
636!          ENDDO
637!          DO  i = nxlg, nxr
638!             DO  j = nysg, nyn
639!                DO  k = nzb, nzt+1
640!                   ustvst(k,j,i) =  &
[4497]641!                      ( 0.5_wp * ( u(k,j,i) + u(k,j,i+1) ) - hom(k,1,1,0) ) *                      &
[3687]642!                      ( 0.5_wp * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,0) )
643!                ENDDO
644!             ENDDO
645!          ENDDO
646
647
648       CASE ( 'after_timestep' )
649!
650!--       Enter actions to be done after every timestep here
651
652
653       CASE ( 'u-tendency' )
654!
655!--       Enter actions to be done in the u-tendency term here
656
657
658       CASE ( 'v-tendency' )
659
660
661       CASE ( 'w-tendency' )
662
663
664       CASE ( 'pt-tendency' )
665
666
667       CASE ( 'sa-tendency' )
668
669
670       CASE ( 'e-tendency' )
671
672
673       CASE ( 'q-tendency' )
674
675
676       CASE ( 's-tendency' )
677
678
679       CASE DEFAULT
[3911]680          CONTINUE
[3687]681
682    END SELECT
683
684    CALL cpu_log( log_point(24), 'user_actions', 'stop' )
685
686 END SUBROUTINE user_actions
687
688
[4497]689!--------------------------------------------------------------------------------------------------!
[3687]690! Description:
691! ------------
692!> Call for grid point i,j
[4497]693!--------------------------------------------------------------------------------------------------!
[3687]694 SUBROUTINE user_actions_ij( i, j, location )
695
696
[4497]697    CHARACTER(LEN=*) ::  location  !<
[3687]698
[4497]699    INTEGER(iwp) ::  i  !<
700    INTEGER(iwp) ::  j  !<
[3687]701
702!
703!-- Here the user-defined actions follow
704    SELECT CASE ( location )
705
706       CASE ( 'u-tendency' )
[3768]707
[3687]708!
[3768]709!--       Next line is to avoid compiler warning about unused variables. Please remove.
710          IF ( i == 0  .OR.  j == 0 )  CONTINUE
711
712!
[3687]713!--       Enter actions to be done in the u-tendency term here
714
715
716       CASE ( 'v-tendency' )
717
718
719       CASE ( 'w-tendency' )
720
721
722       CASE ( 'pt-tendency' )
723
724
725       CASE ( 'sa-tendency' )
726
727
728       CASE ( 'e-tendency' )
729
730
731       CASE ( 'q-tendency' )
732
733
734       CASE ( 's-tendency' )
735
736
737       CASE DEFAULT
[3911]738          CONTINUE
[3687]739
740    END SELECT
741
742 END SUBROUTINE user_actions_ij
743
744
[4497]745!--------------------------------------------------------------------------------------------------!
[3687]746! Description:
747! ------------
[4497]748!> Sum up and time-average user-defined output quantities as well as allocate the array necessary
749!> for storing the average.
750!--------------------------------------------------------------------------------------------------!
[3687]751 SUBROUTINE user_3d_data_averaging( mode, variable )
752
753
[4497]754    CHARACTER(LEN=*) ::  mode      !<
755    CHARACTER(LEN=*) ::  variable  !<
[3687]756
[4497]757!    INTEGER(iwp) ::  i  !<
758!    INTEGER(iwp) ::  j  !<
759!    INTEGER(iwp) ::  k  !<
[3687]760
761    IF ( mode == 'allocate' )  THEN
762
763       SELECT CASE ( TRIM( variable ) )
764
765!
766!--       Uncomment and extend the following lines, if necessary.
[4497]767!--       The arrays for storing the user defined quantities (here u2_av) have to be declared and
768!--       defined by the user!
[3687]769!--       Sample for user-defined output:
770!          CASE ( 'u2' )
771!             IF ( .NOT. ALLOCATED( u2_av ) )  THEN
772!                ALLOCATE( u2_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
773!             ENDIF
774!             u2_av = 0.0_wp
775
776          CASE DEFAULT
777             CONTINUE
778
779       END SELECT
780
781    ELSEIF ( mode == 'sum' )  THEN
782
783       SELECT CASE ( TRIM( variable ) )
784
785!
786!--       Uncomment and extend the following lines, if necessary.
[4497]787!--       The arrays for storing the user defined quantities (here u2 and u2_av) have to be declared
788!--       and defined by the user!
[3687]789!--       Sample for user-defined output:
790!          CASE ( 'u2' )
[4497]791!             IF ( ALLOCATED( u2_av ) )  THEN
[3687]792!                DO  i = nxlg, nxrg
793!                   DO  j = nysg, nyng
794!                      DO  k = nzb, nzt+1
795!                         u2_av(k,j,i) = u2_av(k,j,i) + u2(k,j,i)
796!                      ENDDO
797!                   ENDDO
798!                ENDDO
799!             ENDIF
800
801          CASE DEFAULT
802             CONTINUE
803
804       END SELECT
805
806    ELSEIF ( mode == 'average' )  THEN
807
808       SELECT CASE ( TRIM( variable ) )
809
810!
811!--       Uncomment and extend the following lines, if necessary.
[4497]812!--       The arrays for storing the user defined quantities (here u2_av) have to be declared and
813!--       defined by the user!
[3687]814!--       Sample for user-defined output:
815!          CASE ( 'u2' )
[4497]816!             IF ( ALLOCATED( u2_av ) )  THEN
[3687]817!                DO  i = nxlg, nxrg
818!                   DO  j = nysg, nyng
819!                      DO  k = nzb, nzt+1
820!                         u2_av(k,j,i) = u2_av(k,j,i) / REAL( average_count_3d, KIND=wp )
821!                      ENDDO
822!                   ENDDO
823!                ENDDO
824!             ENDIF
825
826       END SELECT
827
828    ENDIF
829
830
831 END SUBROUTINE user_3d_data_averaging
832
833
[4497]834!--------------------------------------------------------------------------------------------------!
[3687]835! Description:
836! ------------
[4497]837!> Resorts the user-defined output quantity with indices (k,j,i) to a temporary array with indices
838!> (i,j,k) and sets the grid on which it is defined. Allowed values for grid are "zu" and "zw".
839!--------------------------------------------------------------------------------------------------!
[3687]840 SUBROUTINE user_data_output_2d( av, variable, found, grid, local_pf, two_d, nzb_do, nzt_do )
841
842
[4497]843    CHARACTER(LEN=*) ::  grid      !<
844    CHARACTER(LEN=*) ::  variable  !<
[3687]845
[4497]846    INTEGER(iwp) ::  av      !< flag to control data output of instantaneous or time-averaged data
847!    INTEGER(iwp) ::  i       !< grid index along x-direction
848!    INTEGER(iwp) ::  j       !< grid index along y-direction
849!    INTEGER(iwp) ::  k       !< grid index along z-direction
850!    INTEGER(iwp) ::  m       !< running index surface elements
851    INTEGER(iwp) ::  nzb_do  !< lower limit of the domain (usually nzb)
852    INTEGER(iwp) ::  nzt_do  !< upper limit of the domain (usually nzt+1)
[3687]853
[4497]854    LOGICAL      ::  found  !<
855    LOGICAL      ::  two_d  !< flag parameter that indicates 2D variables (horizontal cross sections)
[3687]856
[3768]857!    REAL(wp) ::  fill_value = -999.0_wp    !< value for the _FillValue attribute
[3687]858
[4497]859    REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf  !<
[3687]860
[3768]861!
862!-- Next line is to avoid compiler warning about unused variables. Please remove.
863    IF ( av == 0  .OR.  local_pf(nxl,nys,nzb_do) == 0.0_wp  .OR.  two_d )  CONTINUE
[3687]864
[3768]865
[3687]866    found = .TRUE.
867
868    SELECT CASE ( TRIM( variable ) )
869
870!
871!--    Uncomment and extend the following lines, if necessary.
[4497]872!--    The arrays for storing the user defined quantities (here u2 and u2_av) have to be declared
873!--    and defined by the user!
[3687]874!--    Sample for user-defined output:
875!       CASE ( 'u2_xy', 'u2_xz', 'u2_yz' )
876!          IF ( av == 0 )  THEN
877!             DO  i = nxl, nxr
878!                DO  j = nys, nyn
879!                   DO  k = nzb_do, nzt_do
880!                      local_pf(i,j,k) = u2(k,j,i)
881!                   ENDDO
882!                ENDDO
883!             ENDDO
884!          ELSE
[4497]885!             IF ( .NOT. ALLOCATED( u2_av ) )  THEN
[3687]886!                ALLOCATE( u2_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
887!                u2_av = REAL( fill_value, KIND = wp )
888!             ENDIF
889!             DO  i = nxl, nxr
890!                DO  j = nys, nyn
891!                   DO  k = nzb_do, nzt_do
892!                      local_pf(i,j,k) = u2_av(k,j,i)
893!                   ENDDO
894!                ENDDO
895!             ENDDO
896!          ENDIF
897!
898!          grid = 'zu'
899!
[4497]900!--    In case two-dimensional surface variables are output, the user has to access related
901!--    surface-type. Uncomment and extend following lines appropriately (example output of vertical
902!--    surface momentum flux of u-component). Please note, surface elements can be distributed over
903!--    several data types, depending on their respective surface properties.
[3687]904!       CASE ( 'usws_xy' )
905!          IF ( av == 0 )  THEN
906!
907!--           Horizontal default-type surfaces
908!             DO  m = 1, surf_def_h(0)%ns
909!                i = surf_def_h(0)%i(m)
910!                j = surf_def_h(0)%j(m)
911!                local_pf(i,j,1) = surf_def_h(0)%usws(m)
912!             ENDDO
913!
914!--           Horizontal natural-type surfaces
915!             DO  m = 1, surf_lsm_h%ns
916!                i = surf_lsm_h%i(m)
917!                j = surf_lsm_h%j(m)
918!                local_pf(i,j,1) = surf_lsm_h%usws(m)
919!             ENDDO
920!
921!--           Horizontal urban-type surfaces
922!             DO  m = 1, surf_usm_h%ns
923!                i = surf_usm_h%i(m)
924!                j = surf_usm_h%j(m)
925!                local_pf(i,j,1) = surf_usm_h%usws(m)
926!             ENDDO
927!          ENDIF
928!
929!          grid = 'zu'
[4497]930!--
[3687]931
932
933       CASE DEFAULT
934          found = .FALSE.
935          grid  = 'none'
936
937    END SELECT
938
939
940 END SUBROUTINE user_data_output_2d
941
942
[4497]943!--------------------------------------------------------------------------------------------------!
[3687]944! Description:
945! ------------
[4497]946!> Resorts the user-defined output quantity with indices (k,j,i) to a temporary array with indices
947!> (i,j,k).
948!--------------------------------------------------------------------------------------------------!
[3687]949 SUBROUTINE user_data_output_3d( av, variable, found, local_pf, nzb_do, nzt_do )
950
951
[4497]952    CHARACTER(LEN=*) ::  variable  !<
[3687]953
[4497]954    INTEGER(iwp) ::  av     !<
955!    INTEGER(iwp) ::  i      !<
956!    INTEGER(iwp) ::  j      !<
957!    INTEGER(iwp) ::  k      !<
[3687]958    INTEGER(iwp) ::  nzb_do !< lower limit of the data output (usually 0)
959    INTEGER(iwp) ::  nzt_do !< vertical upper limit of the data output (usually nz_do3d)
960
[4497]961    LOGICAL      ::  found  !<
[3687]962
[4497]963!    REAL(wp) ::  fill_value = -999.0_wp  !< value for the _FillValue attribute
[3687]964
[4497]965    REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) ::  local_pf  !<
[3687]966
[3768]967!
968!-- Next line is to avoid compiler warning about unused variables. Please remove.
969    IF ( av == 0  .OR.  local_pf(nxl,nys,nzb_do) == 0.0_wp )  CONTINUE
[3687]970
[3768]971
[3687]972    found = .TRUE.
973
974    SELECT CASE ( TRIM( variable ) )
975
976!
977!--    Uncomment and extend the following lines, if necessary.
[4497]978!--    The arrays for storing the user defined quantities (here u2 and u2_av) have to be declared
979!--    and defined by the user!
[3687]980!--    Sample for user-defined output:
981!       CASE ( 'u2' )
982!          IF ( av == 0 )  THEN
983!             DO  i = nxl, nxr
984!                DO  j = nys, nyn
985!                   DO  k = nzb_do, nzt_do
986!                      local_pf(i,j,k) = u2(k,j,i)
987!                   ENDDO
988!                ENDDO
989!             ENDDO
990!          ELSE
[4497]991!             IF ( .NOT. ALLOCATED( u2_av ) )  THEN
[3687]992!                ALLOCATE( u2_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
993!                u2_av = REAL( fill_value, KIND = wp )
994!             ENDIF
995!             DO  i = nxl, nxr
996!                DO  j = nys, nyn
997!                   DO  k = nzb_do, nzt_do
998!                      local_pf(i,j,k) = u2_av(k,j,i)
999!                   ENDDO
1000!                ENDDO
1001!             ENDDO
1002!          ENDIF
1003!
1004
1005       CASE DEFAULT
1006          found = .FALSE.
1007
1008    END SELECT
1009
1010
1011 END SUBROUTINE user_data_output_3d
1012
1013
[4497]1014!--------------------------------------------------------------------------------------------------!
[3687]1015! Description:
1016! ------------
[4497]1017!> Calculation of user-defined statistics, i.e. horizontally averaged profiles and time series.
1018!> This routine is called for every statistic region sr defined by the user, but at least for the
1019!> region "total domain" (sr=0). See section 3.5.4 on how to define, calculate, and output user
1020!> defined quantities.
1021!--------------------------------------------------------------------------------------------------!
[3687]1022 SUBROUTINE user_statistics( mode, sr, tn )
1023
1024
[4497]1025    CHARACTER(LEN=*) ::  mode  !<
1026!    INTEGER(iwp) ::  i   !<
1027!    INTEGER(iwp) ::  j   !<
1028!    INTEGER(iwp) ::  k   !<
1029    INTEGER(iwp) ::  sr  !<
1030    INTEGER(iwp) ::  tn  !<
[3687]1031
[4497]1032!    REAL(wp), DIMENSION(:), ALLOCATABLE ::  ts_value_l  !<
[3687]1033
[3768]1034!
1035!-- Next line is to avoid compiler warning about unused variables. Please remove.
1036    IF ( sr == 0  .OR.  tn == 0 )  CONTINUE
1037
[3687]1038    IF ( mode == 'profiles' )  THEN
1039
1040!
[4497]1041!--    Sample on how to calculate horizontally averaged profiles of user-defined quantities. Each
1042!--    quantity is identified by the index "pr_palm+#" where "#" is an integer starting from 1.
1043!--    These user-profile-numbers must also be assigned to the respective strings given by
1044!--    data_output_pr_user in routine user_check_data_output_pr.
[3687]1045!       !$OMP DO
1046!       DO  i = nxl, nxr
1047!          DO  j = nys, nyn
1048!             DO  k = nzb+1, nzt
1049!!
[4497]1050!!--             Sample on how to calculate the profile of the resolved-scale horizontal momentum
1051!!--             flux u*v*
1052!                sums_l(k,pr_palm+1,tn) = sums_l(k,pr_palm+1,tn) +                                  &
1053!                                         ( 0.5_wp * ( u(k,j,i) + u(k,j,i+1) ) - hom(k,1,1,sr) ) *  &
1054!                                         ( 0.5_wp * ( v(k,j,i) + v(k,j+1,i) ) - hom(k,1,2,sr) ) *  &
1055!                                         rmask(j,i,sr) * MERGE( 1.0_wp, 0.0_wp,                    &
1056!                                         BTEST( wall_flags_total_0(k,j,i), 0 ) )
[3687]1057!!
[4497]1058!!--             Further profiles can be defined and calculated by increasing the second index of
1059!!--             array sums_l (replace ... appropriately)
1060!                sums_l(k,pr_palm+2,tn) = sums_l(k,pr_palm+2,tn) + ...   * rmask(j,i,sr)
[3687]1061!             ENDDO
1062!          ENDDO
1063!       ENDDO
1064
1065    ELSEIF ( mode == 'time_series' )  THEN
1066
[3700]1067
1068!       ALLOCATE ( ts_value_l(dots_num_user) )
[3687]1069!
1070!--    Sample on how to add values for the user-defined time series quantities.
[4497]1071!--    These have to be defined before in routine user_init. This sample creates two time series for
1072!--    the absolut values of the horizontal velocities u and v.
[3687]1073!       ts_value_l = 0.0_wp
[3700]1074!       ts_value_l(1) = ABS( u_max )
1075!       ts_value_l(2) = ABS( v_max )
[3687]1076!
1077!--     Collect / send values to PE0, because only PE0 outputs the time series.
[4497]1078!--     CAUTION: Collection is done by taking the sum over all processors. You may have to normalize
1079!--              this sum, depending on the quantity that you like to calculate. For serial runs,
1080!--              nothing has to be done.
1081!--     HINT: If the time series value that you are calculating has the same value on all PEs, you
1082!--           can omit the MPI_ALLREDUCE call and assign ts_value(dots_num_palm+1:,sr) = ts_value_l directly.
[3687]1083!#if defined( __parallel )
1084!       IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
[4497]1085!       CALL MPI_ALLREDUCE( ts_value_l(1), ts_value(dots_num_palm+1,sr), dots_num_user, MPI_REAL,   &
1086!                           MPI_MAX, comm2d, ierr )
[3687]1087!#else
[3700]1088!       ts_value(dots_num_palm+1:dots_num_palm+dots_num_user,sr) = ts_value_l
[3687]1089!#endif
1090
1091    ENDIF
1092
1093 END SUBROUTINE user_statistics
1094
1095
[4497]1096!--------------------------------------------------------------------------------------------------!
[3687]1097! Description:
1098! ------------
[4495]1099!> Read module-specific global restart data (Fortran binary format).
[4497]1100!--------------------------------------------------------------------------------------------------!
[4495]1101 SUBROUTINE user_rrd_global_ftn( found )
[3687]1102
1103
[4497]1104    LOGICAL, INTENT(OUT)  ::  found  !<
[3687]1105
1106
[3768]1107    found = .TRUE.
[3687]1108
1109
[3768]1110    SELECT CASE ( restart_string(1:length) )
[3687]1111
[3768]1112       CASE ( 'global_paramter' )
1113!          READ ( 13 )  global_parameter
[3687]1114
[3768]1115       CASE DEFAULT
[4497]1116
[3768]1117          found = .FALSE.
[3687]1118
[3768]1119    END SELECT
[3687]1120
1121
[4495]1122 END SUBROUTINE user_rrd_global_ftn
[3687]1123
1124
[4497]1125!--------------------------------------------------------------------------------------------------!
[3687]1126! Description:
1127! ------------
[4495]1128!> Read module-specific global restart data (MPI-IO).
[4497]1129!--------------------------------------------------------------------------------------------------!
[4495]1130 SUBROUTINE user_rrd_global_mpi
1131
1132!    CALL rrd_mpi_io( 'global_parameter', global_parameter )
1133    CONTINUE
1134
1135 END SUBROUTINE user_rrd_global_mpi
1136
1137
[4497]1138!--------------------------------------------------------------------------------------------------!
[4495]1139! Description:
1140! ------------
[4497]1141!> Reading processor specific restart data from file(s) that has been defined by the user. Subdomain
1142!> index limits on file are given by nxl_on_file, etc. Indices nxlc, etc. indicate the range of
1143!> gridpoints to be mapped from the subdomain on file (f) to the subdomain of the current PE (c).
1144!> They have been calculated in routine rrd_local.
1145!--------------------------------------------------------------------------------------------------!
1146 SUBROUTINE user_rrd_local( k, nxlf, nxlc, nxl_on_file, nxrf, nxrc, nxr_on_file, nynf, nync,       &
1147                            nyn_on_file, nysf, nysc, nys_on_file, tmp_3d, found )
[3687]1148
1149
[3768]1150    INTEGER(iwp) ::  idum            !<
1151    INTEGER(iwp) ::  k               !<
1152    INTEGER(iwp) ::  nxlc            !<
1153    INTEGER(iwp) ::  nxlf            !<
1154    INTEGER(iwp) ::  nxl_on_file     !<
1155    INTEGER(iwp) ::  nxrc            !<
1156    INTEGER(iwp) ::  nxrf            !<
1157    INTEGER(iwp) ::  nxr_on_file     !<
1158    INTEGER(iwp) ::  nync            !<
1159    INTEGER(iwp) ::  nynf            !<
1160    INTEGER(iwp) ::  nyn_on_file     !<
1161    INTEGER(iwp) ::  nysc            !<
1162    INTEGER(iwp) ::  nysf            !<
1163    INTEGER(iwp) ::  nys_on_file     !<
[3687]1164
[4497]1165    LOGICAL, INTENT(OUT)  ::  found  !<
[3687]1166
[4497]1167    REAL(wp), DIMENSION(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,nxl_on_file-nbgp:nxr_on_file+nbgp) :: tmp_3d  !<
[3687]1168
1169!
[3768]1170!-- Next line is to avoid compiler warning about unused variables. Please remove.
1171    idum = k + nxlc + nxlf + nxrc + nxrf + nync + nynf + nysc + nysf +                             &
1172           INT( tmp_3d(nzb,nys_on_file,nxl_on_file) )
1173
1174!
[3687]1175!-- Here the reading of user-defined restart data follows:
1176!-- Sample for user-defined output
1177
[3768]1178    found = .TRUE.
[3687]1179
[3768]1180    SELECT CASE ( restart_string(1:length) )
[3687]1181
[3768]1182       CASE ( 'u2_av' )
[4497]1183!          IF ( .NOT. ALLOCATED( u2_av ) )  THEN
[3768]1184!               ALLOCATE( u2_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
1185!          ENDIF
1186!          IF ( k == 1 )  READ ( 13 )  tmp_3d
[4497]1187!             u2_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) =                                    &
1188!             tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
[3687]1189!
[3768]1190       CASE DEFAULT
[3687]1191
[3768]1192          found = .FALSE.
[3687]1193
[3768]1194    END SELECT
[3687]1195
[3768]1196 END SUBROUTINE user_rrd_local
[3687]1197
1198
[4497]1199!--------------------------------------------------------------------------------------------------!
[3687]1200! Description:
1201! ------------
[4497]1202!> Writes global and user-defined restart data into binary file(s) for restart runs.
1203!--------------------------------------------------------------------------------------------------!
[3768]1204 SUBROUTINE user_wrd_global
[3687]1205
[4495]1206    IF ( TRIM( restart_data_format_output ) == 'fortran_binary' )  THEN
[3687]1207
[4495]1208!       CALL wrd_write_string( 'global_parameter' )
1209!       WRITE ( 14 )  global_parameter
1210
1211    ELSEIF ( TRIM( restart_data_format_output ) == 'mpi' )  THEN
1212
1213!    CALL rrd_mpi_io( 'global_parameter', global_parameter )
1214
1215    ENDIF
1216
[3768]1217 END SUBROUTINE user_wrd_global
[3687]1218
1219
[4497]1220!--------------------------------------------------------------------------------------------------!
[3687]1221! Description:
1222! ------------
[4497]1223!> Writes processor specific and user-defined restart data into binary file(s) for restart runs.
1224!--------------------------------------------------------------------------------------------------!
[3768]1225 SUBROUTINE user_wrd_local
[3687]1226
1227!
1228!-- Here the user-defined actions at the end of a job follow.
1229!-- Sample for user-defined output:
[4495]1230    IF ( TRIM( restart_data_format_output ) == 'fortran_binary' )  THEN
[3687]1231
[4495]1232!       IF ( ALLOCATED( u2_av ) )  THEN
1233!          CALL wrd_write_string( 'u2_av' )
1234!          WRITE ( 14 )  u2_av
1235!       ENDIF
1236
1237    ELSEIF ( TRIM( restart_data_format_output ) == 'mpi' )  THEN
1238
1239!       IF ( ALLOCATED( u2_av ) )  CALL wrd_mpi_io( 'u2_av', u2_av )
1240
1241    ENDIF
1242
[3768]1243 END SUBROUTINE user_wrd_local
[3687]1244
1245
[4497]1246!--------------------------------------------------------------------------------------------------!
[3687]1247! Description:
1248! ------------
1249!> Execution of user-defined actions at the end of a job.
[4497]1250!--------------------------------------------------------------------------------------------------!
[3687]1251 SUBROUTINE user_last_actions
1252
1253!
1254!-- Here the user-defined actions at the end of a job might follow.
1255
1256
1257 END SUBROUTINE user_last_actions
1258
[211]1259 END MODULE user
Note: See TracBrowser for help on using the repository browser.