Ignore:
Timestamp:
Sep 10, 2019 6:04:34 PM (22 months ago)
Author:
gronemeier
Message:

implement new palm_date_time_mod; replaced namelist parameters time_utc_init and day_of_year_init by origin_date_time

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r4214 r4227  
    2828! -----------------
    2929! $Id$
     30! implement new palm_date_time_mod
     31!
     32! 4214 2019-09-02 15:57:02Z suehring
    3033! Bugfix, missing initialization and clearing of soil-moisture tendency
    3134! (J.Resler)
     
    222225        ONLY:  cpu_log, log_point, log_point_s
    223226
    224     USE date_and_time_mod,                                                     &
    225         ONLY:  time_utc_init
    226 
    227227    USE grid_variables,                                                        &
    228228        ONLY:  dx, dy, ddx, ddy, ddx2, ddy2
     
    236236    USE kinds
    237237             
     238    USE palm_date_time_mod,                                                    &
     239        ONLY:  get_date_time, seconds_per_hour
     240
    238241    USE pegrid
    239242       
     
    83058308!--        calculation of actual profile coefficient
    83068309!--        ??? 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 )
    83098311
    83108312!--         TO_DO: activate, if testcase is available
     
    83198321!--                    given to anthropogenic heat aheat*acoef (W*m-2)
    83208322!--                    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)
    83238325                        IF ( aheat(k,j,i) > 0.0_wp )  THEN
    83248326!
Note: See TracChangeset for help on using the changeset viewer.