Changeset 4227 for palm/trunk/SOURCE/time_integration_spinup.f90

Timestamp:

Sep 10, 2019 6:04:34 PM (5 years 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:

• palm/trunk/SOURCE/time_integration_spinup.f90 (modified) (8 diffs)

palm/trunk/SOURCE/time_integration_spinup.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! implement new palm_date_time_mod
+! 
+! 4223 2019-09-10 09:20:47Z gronemeier
 ! Corrected "Former revisions" section
 ! 
@@ -79 +82 @@
         ONLY:  cpu_log, log_point_s
 
-    USE date_and_time_mod,                                                     &
-        ONLY: day_of_year_init, time_utc_init
-
     USE indices,                                                               &
         ONLY:  nbgp, nzb, nzt, nysg, nyng, nxlg, nxrg
 
-
     USE land_surface_model_mod,                                                &
         ONLY:  lsm_energy_balance, lsm_soil_model, lsm_swap_timelevel
+
     USE pegrid
 
@@ -94 +94 @@
 
     USE kinds
+
+    USE palm_date_time_mod,                                                    &
+        ONLY:  get_date_time, seconds_per_hour
 
     USE radiation_model_mod,                                                   &
@@ -132 +135 @@
 
     INTEGER(iwp) :: current_timestep_number_spinup = 0  !< number if timestep during spinup
+    INTEGER(iwp) :: day_of_year                         !< day of the year
   
     LOGICAL :: run_control_header_spinup = .FALSE.  !< flag parameter for steering whether the header information must be output
 
-    REAL(wp) ::  pt_spinup   !< temporary storage of temperature
-    REAL(wp) ::  dt_save     !< temporary storage for time step
+    REAL(wp) ::  pt_spinup      !< temporary storage of temperature
+    REAL(wp) ::  dt_save        !< temporary storage for time step
+    REAL(wp) ::  second_of_day  !< second of the day
                   
     REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::  pt_save  !< temporary storage of temperature
@@ -234 +239 @@
 !--       shifted by one hour to simulate a lag between air temperature and
 !--       incoming radiation
+          CALL get_date_time( simulated_time - spinup_time - seconds_per_hour, &
+                              day_of_year=day_of_year,                         &
+                              second_of_day=second_of_day                      )
+
           pt_spinup = spinup_pt_mean + spinup_pt_amplitude                     &
-             * solar_angle (time_utc_init + time_since_reference_point - 3600.0)
+                                     * solar_angle(day_of_year, second_of_day)
 
 !
@@ -546 +555 @@
 !-- Returns the cosine of the solar zenith angle at a given time. This routine
 !-- is similar to that for calculation zenith (see radiation_model_mod.f90)
-    FUNCTION solar_angle( local_time ) 
+    !> @todo Load function calc_zenith of radiation model instead of
+    !>       rewrite the function here.
+    FUNCTION solar_angle( day_of_year, second_of_day ) 
 
        USE basic_constants_and_equations_mod,                                  &
@@ -559 +570 @@
 
 
-       REAL(wp) ::  solar_angle     !< cosine of the solar zenith angle
-
-       REAL(wp) ::  day             !< day of the year
-       REAL(wp) ::  declination     !< solar declination angle
-       REAL(wp) ::  hour_angle      !< solar hour angle
-       REAL(wp) ::  time_utc        !< current time in UTC
-       REAL(wp), INTENT(IN) :: local_time
-!
-!--    Calculate current day and time based on the initial values and simulation
-!--    time
-
-       day = day_of_year_init + INT(FLOOR( local_time / 86400.0_wp ), KIND=iwp)
-       time_utc = MOD(local_time, 86400.0_wp)
-
-
+       INTEGER(iwp), INTENT(IN) ::  day_of_year  !< day of the year
+
+       REAL(wp)             ::  declination      !< solar declination angle
+       REAL(wp)             ::  hour_angle       !< solar hour angle
+       REAL(wp), INTENT(IN) ::  second_of_day    !< current time of the day in UTC
+       REAL(wp)             ::  solar_angle      !< cosine of the solar zenith angle
 !
 !--    Calculate solar declination and hour angle   
-       declination = ASIN( decl_1 * SIN(decl_2 * REAL(day, KIND=wp) - decl_3) )
-       hour_angle  = 2.0_wp * pi * (time_utc / 86400.0_wp) + lon - pi
+       declination = ASIN( decl_1 * SIN(decl_2 * REAL(day_of_year, KIND=wp) - decl_3) )
+       hour_angle  = 2.0_wp * pi * (second_of_day / 86400.0_wp) + lon - pi
 
 !
@@ -584 +586 @@
                      * COS(hour_angle)
 
-
     END FUNCTION solar_angle