Changeset 4227 for palm/trunk/SOURCE/biometeorology_mod.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/biometeorology_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/biometeorology_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! implement new palm_date_time_mod
+! 
+! 4223 2019-09-10 09:20:47Z gronemeier
 ! Corrected "Former revisions" section
 ! 
@@ -140 +143 @@
                surface_pressure
 
-    USE date_and_time_mod,                                                     &
-        ONLY:  calc_date_and_time, day_of_year, time_utc
-
     USE grid_variables,                                                        &
         ONLY:  ddx, dx, ddy, dy
@@ -155 +155 @@
         ONLY:  netcdf_data_input_uvem, uvem_projarea_f, uvem_radiance_f,       &
                uvem_irradiance_f, uvem_integration_f, building_obstruction_f
+
+    USE palm_date_time_mod,                                                    &
+        ONLY:  get_date_time
 !
 !-- Import radiation model to obtain input for mean radiant temperature
@@ -4226 +4229 @@
 !-- UVEM specific subroutines
 
-!---------------------------------------------------------------------------------------------------------------------!
+!--------------------------------------------------------------------------------------------------!
 ! Description:
 ! ------------
 !> Module-specific routine for new module
-!---------------------------------------------------------------------------------------------------------------------!
+!--------------------------------------------------------------------------------------------------!
  SUBROUTINE uvem_solar_position
-    
-    USE date_and_time_mod,                                                                                            &
-       ONLY:  calc_date_and_time, day_of_year, time_utc 
-    
-    USE control_parameters,                                                                                           &
-       ONLY:  latitude, longitude   
+
+    USE control_parameters,                                                                        &
+       ONLY:  latitude, longitude, time_since_reference_point
 
     IMPLICIT NONE
-    
-   
-    REAL(wp) ::  alpha       = 0.0_wp   !< solar azimuth angle in radiant    
-    REAL(wp) ::  doy_r       = 0.0_wp   !< real format of day_of_year           
-    REAL(wp) ::  declination = 0.0_wp   !< declination
-    REAL(wp) ::  dtor        = 0.0_wp   !< factor to convert degree to radiant
-    REAL(wp) ::  js          = 0.0_wp   !< parameter for solar position calculation
-    REAL(wp) ::  lat         = 52.39_wp !< latitude
-    REAL(wp) ::  lon         = 9.7_wp   !< longitude        
-    REAL(wp) ::  thetar      = 0.0_wp   !< angle for solar zenith angle calculation
-    REAL(wp) ::  thetasr     = 0.0_wp   !< angle for solar azimuth angle calculation    
-    REAL(wp) ::  zgl         = 0.0_wp   !< calculated exposure by direct beam   
-    REAL(wp) ::  woz         = 0.0_wp   !< calculated exposure by diffuse radiation
-    REAL(wp) ::  wsp         = 0.0_wp   !< calculated exposure by direct beam   
-   
-
-    CALL calc_date_and_time
-    doy_r = real(day_of_year)   
+
+    INTEGER(iwp) ::  day_of_year = 0       !< day of year
+
+    REAL(wp) ::  alpha         = 0.0_wp    !< solar azimuth angle in radiant    
+    REAL(wp) ::  declination   = 0.0_wp    !< declination
+    REAL(wp) ::  dtor          = 0.0_wp    !< factor to convert degree to radiant
+    REAL(wp) ::  js            = 0.0_wp    !< parameter for solar position calculation
+    REAL(wp) ::  lat           = 52.39_wp  !< latitude
+    REAL(wp) ::  lon           = 9.7_wp    !< longitude
+    REAL(wp) ::  second_of_day = 0.0_wp    !< current second of the day
+    REAL(wp) ::  thetar        = 0.0_wp    !< angle for solar zenith angle calculation
+    REAL(wp) ::  thetasr       = 0.0_wp    !< angle for solar azimuth angle calculation    
+    REAL(wp) ::  zgl           = 0.0_wp    !< calculated exposure by direct beam   
+    REAL(wp) ::  woz           = 0.0_wp    !< calculated exposure by diffuse radiation
+    REAL(wp) ::  wsp           = 0.0_wp    !< calculated exposure by direct beam   
+
+
+    CALL get_date_time( time_since_reference_point, &
+                        day_of_year = day_of_year, second_of_day = second_of_day )
     dtor = pi / 180.0_wp
     lat = latitude
@@ -4263 +4264 @@
 !
 !-- calculation of js, necessary for calculation of equation of time (zgl) :
-    js=  72.0_wp * ( doy_r + ( time_utc / 86400.0_wp ) ) / 73.0_wp 
+    js=  72.0_wp * ( REAL( day_of_year, KIND=wp ) + ( second_of_day / 86400.0_wp ) ) / 73.0_wp 
 !
 !-- calculation of equation of time (zgl):
@@ -4270 +4271 @@
 !
 !-- calculation of apparent solar time woz:
-    woz = ( ( time_utc / 3600.0_wp ) - ( 4.0_wp * ( 15.0_wp - lon ) ) / 60.0_wp ) + ( zgl / 60.0_wp )
+    woz = ( ( second_of_day / 3600.0_wp ) - ( 4.0_wp * ( 15.0_wp - lon ) ) / 60.0_wp ) + ( zgl / 60.0_wp )
 !
 !-- calculation of hour angle (wsp):