Changeset 4358 for palm/trunk/SOURCE/nesting_offl_mod.f90
 Timestamp:
 Dec 30, 2019 1:03:09 PM (3 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/nesting_offl_mod.f90
r4346 r4358 25 25 !  26 26 ! $Id$ 27 ! Bugfix, time coordinate is relative to origin_time rather than to 00:00:00 28 ! UTC. 29 ! 30 ! 4346 20191218 11:55:56Z motisi 27 31 ! Introduction of wall_flags_total_0, which currently sets bits based on static 28 32 ! topography information used in wall_flags_static_0 … … 213 217 pids_id 214 218 215 USE palm_date_time_mod, &216 ONLY: get_date_time217 218 219 USE pegrid 219 220 … … 304 305 305 306 REAL(wp) :: fac_dt !< interpolation factor 306 REAL(wp) :: time_utc_init !< time in secondsofday of origin_date_time307 307 REAL(wp) :: zi_ribulk = 0.0_wp !< boundarylayer depth according to bulk Richardson criterion, i.e. the height where Ri_bulk exceeds the critical 308 308 !< bulk Richardson number of 0.2 … … 433 433 ! Check if dynamic driver data input is required. 434 434 IF ( nest_offl%time(nest_offl%tind_p) <= & 435 MAX( time_since_reference_point, 0.0_wp) + time_utc_init .OR.&435 MAX( time_since_reference_point, 0.0_wp) .OR. & 436 436 .NOT. nest_offl%init ) THEN 437 437 CONTINUE … … 447 447 ! 448 448 ! Obtain time index for current point in time. Note, the time coordinate 449 ! in the input file is relative to time_utc_init. Since time_since_...450 ! is negativ when spinup is used, use MAX function to obtain correct451 ! time at the beginning.452 nest_offl%tind = MINLOC( ABS( nest_offl%time  ( & 453 time_utc_init +&454 MAX( time_since_reference_point, 0.0_wp) )&449 ! in the input file is always relative to the initial time in UTC, i.e. 450 ! the time coordinate always starts at 0.0 even if the initial UTC is e.g. 451 ! 7200.0. Further, since time_since_reference_point is negativ here when 452 ! spinup is applied, use MAX function to obtain correct time index. 453 nest_offl%tind = MINLOC( ABS( nest_offl%time  & 454 MAX( time_since_reference_point, 0.0_wp) & 455 455 ), DIM = 1 )  1 456 456 nest_offl%tind_p = nest_offl%tind + 1 … … 1582 1582 ! t+dt can slightly exceed time(tind_p) before boundary data is updated 1583 1583 ! again. 1584 fac_dt = ( time_ utc_init + time_since_reference_point&1584 fac_dt = ( time_since_reference_point & 1585 1585  nest_offl%time(nest_offl%tind) + dt_3d ) / & 1586 1586 ( nest_offl%time(nest_offl%tind_p)  nest_offl%time(nest_offl%tind) ) … … 1886 1886 1887 1887 INTEGER(iwp) :: n !< running index for chemical species 1888 1889 !1890 ! Get time_utc_init from origin_date_time1891 CALL get_date_time( 0.0_wp, second_of_day = time_utc_init )1892 1888 1893 1889 ! Allocate arrays for geostrophic wind components. Arrays will
Note: See TracChangeset
for help on using the changeset viewer.