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

Timestamp:

Dec 30, 2019 1:03:09 PM (4 years ago)

Author:

suehring

Message:

Bugfix, time coordinate is relative to origin_time rather than to 00:00:00 UTC; Refine post-initialization check for realistically inital values of mixing ratio. Give an error message for faulty initial values, but only a warning in a restart run.

File:

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

palm/trunk/SOURCE/nesting_offl_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix, time coordinate is relative to origin_time rather than to 00:00:00 
+! UTC.
+! 
+! 4346 2019-12-18 11:55:56Z motisi
 ! Introduction of wall_flags_total_0, which currently sets bits based on static
 ! topography information used in wall_flags_static_0
@@ -213 +217 @@
                pids_id
 
-    USE palm_date_time_mod,                                                    &
-        ONLY:  get_date_time
-
     USE pegrid
 
@@ -304 +305 @@
 
     REAL(wp) ::  fac_dt              !< interpolation factor
-    REAL(wp) ::  time_utc_init       !< time in seconds-of-day of origin_date_time
     REAL(wp) ::  zi_ribulk = 0.0_wp  !< boundary-layer depth according to bulk Richardson criterion, i.e. the height where Ri_bulk exceeds the critical
                                      !< bulk Richardson number of 0.2
@@ -433 +433 @@
 !--    Check if dynamic driver data input is required.
        IF ( nest_offl%time(nest_offl%tind_p) <=                                &
-            MAX( time_since_reference_point, 0.0_wp) + time_utc_init  .OR.     &
+            MAX( time_since_reference_point, 0.0_wp)  .OR.                     &
             .NOT.  nest_offl%init )  THEN
           CONTINUE
@@ -447 +447 @@
 !
 !--    Obtain time index for current point in time. Note, the time coordinate
-!--    in the input file is relative to time_utc_init. Since time_since_...
-!--    is negativ when spinup is used, use MAX function to obtain correct 
-!--    time at the beginning. 
-       nest_offl%tind = MINLOC( ABS( nest_offl%time - (                        &
-                                     time_utc_init +                           &
-                                     MAX( time_since_reference_point, 0.0_wp) )&
+!--    in the input file is always relative to the initial time in UTC, i.e.
+!--    the time coordinate always starts at 0.0 even if the initial UTC is e.g.
+!--    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_p = nest_offl%tind + 1
@@ -1582 +1582 @@
 !--    t+dt can slightly exceed time(tind_p) before boundary data is updated 
 !--    again. 
-       fac_dt = ( time_utc_init + time_since_reference_point                   &
+       fac_dt = ( time_since_reference_point                                   &
                 - nest_offl%time(nest_offl%tind) + dt_3d ) /                   &
            ( nest_offl%time(nest_offl%tind_p) - nest_offl%time(nest_offl%tind) )
@@ -1886 +1886 @@
            
        INTEGER(iwp) ::  n !< running index for chemical species
-
-!
-!--    Get time_utc_init from origin_date_time
-       CALL get_date_time( 0.0_wp, second_of_day = time_utc_init )
 
 !--    Allocate arrays for geostrophic wind components. Arrays will