Changeset 4227 for palm/trunk/SOURCE/urban_surface_mod.f90
- Timestamp:
- Sep 10, 2019 6:04:34 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/urban_surface_mod.f90
r4214 r4227 28 28 ! ----------------- 29 29 ! $Id$ 30 ! implement new palm_date_time_mod 31 ! 32 ! 4214 2019-09-02 15:57:02Z suehring 30 33 ! Bugfix, missing initialization and clearing of soil-moisture tendency 31 34 ! (J.Resler) … … 222 225 ONLY: cpu_log, log_point, log_point_s 223 226 224 USE date_and_time_mod, &225 ONLY: time_utc_init226 227 227 USE grid_variables, & 228 228 ONLY: dx, dy, ddx, ddy, ddx2, ddy2 … … 236 236 USE kinds 237 237 238 USE palm_date_time_mod, & 239 ONLY: get_date_time, seconds_per_hour 240 238 241 USE pegrid 239 242 … … 8305 8308 !-- calculation of actual profile coefficient 8306 8309 !-- ??? check time_since_reference_point ??? 8307 dtime = mod(simulated_time + time_utc_init, 24.0_wp*3600.0_wp) 8308 dhour = INT(dtime/3600.0_wp) 8310 CALL get_date_time( time_since_reference_point, hour=dhour, second_of_day=dtime ) 8309 8311 8310 8312 !-- TO_DO: activate, if testcase is available … … 8319 8321 !-- given to anthropogenic heat aheat*acoef (W*m-2) 8320 8322 !-- linear interpolation of coeficient 8321 acoef = (REAL(dhour+1,wp)-dtime/ 3600.0_wp)*aheatprof(k,dhour) + &8322 (dtime/ 3600.0_wp-REAL(dhour,wp))*aheatprof(k,dhour+1)8323 acoef = (REAL(dhour+1,wp)-dtime/seconds_per_hour)*aheatprof(k,dhour) + & 8324 (dtime/seconds_per_hour-REAL(dhour,wp))*aheatprof(k,dhour+1) 8323 8325 IF ( aheat(k,j,i) > 0.0_wp ) THEN 8324 8326 !
Note: See TracChangeset
for help on using the changeset viewer.