Changeset 4795 for palm/trunk/SOURCE/nesting_offl_mod.f90

Timestamp:

Nov 25, 2020 3:55:14 PM (4 years ago)

Author:

suehring

Message:

mesoscale nesting: bugfix in obtaining the correct timestamp in case of restart runs; virtual measurements: add missing control flags

File:

• palm/trunk/SOURCE/nesting_offl_mod.f90 (modified) (4 diffs)

palm/trunk/SOURCE/nesting_offl_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in obtaining the correct timestamp in case of restart runs
+!
+! 4724 2020-10-06 17:20:39Z suehring
 ! - Enable LOD=1 input of boundary conditions
 ! - Minor bugfix - add missing initialization of the top boundary
@@ -300 +303 @@
        INTEGER(iwp) ::  nzu               !< number of vertical levels on scalar grid in dynamic input file
        INTEGER(iwp) ::  nzw               !< number of vertical levels on w grid in dynamic input file
-       INTEGER(iwp) ::  tind              !< time index for reference time in mesoscale-offline nesting
-       INTEGER(iwp) ::  tind_p            !< time index for following time in mesoscale-offline nesting
+       INTEGER(iwp) ::  tind = 0          !< time index for reference time in mesoscale-offline nesting
+       INTEGER(iwp) ::  tind_p = 0        !< time index for following time in mesoscale-offline nesting
 
        LOGICAL      ::  init = .FALSE.    !< flag indicating that offline nesting is already initialized
@@ -502 +505 @@
        ENDIF
 !
-!--    CPU measurement
+!--    Start of CPU measurement
        CALL cpu_log( log_point_s(86), 'NetCDF input forcing', 'start' )
-
 !
 !--    Obtain time index for current point in time. Note, the time coordinate
@@ -511 +513 @@
 !--    7200.0. Further, since time_since_reference_point is negativ here when 
 !--    spinup is applied, use MAX function to obtain correct time index. 
-       nest_offl%tind = MINLOC( ABS( nest_offl%time -                                              &
-                                     MAX( time_since_reference_point, 0.0_wp)                      &
-                                   ), DIM = 1 ) - 1
+       nest_offl%tind = MINLOC( ABS( nest_offl%time - MAX( time_since_reference_point, 0.0_wp) ),  &
+                                DIM = 1 ) - 1
+!
+!--    Note, in case of restart runs, the time index for the boundary data may indicate a time
+!--    in the future. This needs to be checked and corrected.
+       IF ( TRIM( initializing_actions ) == 'read_restart_data'  .AND.                             &
+            nest_offl%time(nest_offl%tind) > time_since_reference_point )  THEN
+          nest_offl%tind = nest_offl%tind - 1
+       ENDIF
        nest_offl%tind_p = nest_offl%tind + 1
 !