Changeset 2544 for palm/trunk/SOURCE/radiation_model_mod.f90

Timestamp:

Oct 13, 2017 6:09:32 PM (7 years ago)

Author:

maronga

Message:

introduced new module date_and_time_mod

File:

• palm/trunk/SOURCE/radiation_model_mod.f90 (modified) (11 diffs)

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Moved date and time quantitis to separate module date_and_time_mod 
+! 
+! 2512 2017-10-04 08:26:59Z raasch
 ! upper bounds of cross section and 3d output changed from nx+1,ny+1 to nx,ny
 ! no output of ghost layer data
@@ -183 +186 @@
                surface_pressure, time_since_reference_point
 
+    USE date_and_time_mod,                                                     &
+        ONLY:  calc_date_and_time, d_hours_day, d_seconds_hour, day_of_year,   &
+               time_utc
+
     USE indices,                                                               &
         ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
@@ -262 +269 @@
 
     INTEGER(iwp) :: albedo_type  = 9999999, & !< Albedo surface type
-                    day,                    & !< current day of the year
-                    day_init     = 172,     & !< day of the year at model start (21/06)
                     dots_rad     = 0          !< starting index for timeseries output 
 
@@ -276 +281 @@
 
 
-    REAL(wp), PARAMETER :: d_seconds_hour  = 0.000277777777778_wp,  & !< inverse of seconds per hour (1/3600)
-                           d_hours_day    = 0.0416666666667_wp,     & !< inverse of hours per day (1/24)
-                           sigma_sb       = 5.67037321E-8_wp,       & !< Stefan-Boltzmann constant
+    REAL(wp), PARAMETER :: sigma_sb       = 5.67037321E-8_wp,       & !< Stefan-Boltzmann constant
                            solar_constant = 1368.0_wp                 !< solar constant at top of atmosphere
 
@@ -297 +300 @@
                 skip_time_do_radiation = 0.0_wp, & !< Radiation model is not called before this time
                 sky_trans,                       & !< sky transmissivity
-                time_radiation = 0.0_wp,         & !< time since last call of radiation code
-                time_utc,                        & !< current time in UTC
-                time_utc_init = 43200.0_wp         !< UTC time at model start (noon)
+                time_radiation = 0.0_wp            !< time since last call of radiation code
+
 
     REAL(wp), DIMENSION(0:0) ::  zenith,         & !< cosine of solar zenith angle
@@ -562 +564 @@
            rad_sw_in_av, rad_sw_out, rad_sw_out_av, rad_sw_cs_hr,              &
            rad_sw_cs_hr_av, rad_sw_hr, rad_sw_hr_av, sigma_sb,                 &
-           skip_time_do_radiation, time_radiation, unscheduled_radiation_calls,&
-           zenith, calc_zenith, sun_direction, sun_dir_lat, sun_dir_lon,       &
-           day_init, time_utc_init
+           skip_time_do_radiation, solar_constant, time_radiation,             &
+           unscheduled_radiation_calls, zenith, calc_zenith, sun_direction,    &
+           sun_dir_lat, sun_dir_lon
 
 
@@ -1360 +1362 @@
        IF ( radiation .AND. radiation_scheme /= 'constant' )  THEN
           WRITE ( io, 1 )  lambda
-          WRITE ( io, 2 )  day_init, time_utc_init
        ENDIF
 
@@ -1367 +1368 @@
 
  1 FORMAT ('    Geograph. longitude            :   lambda = ',F4.1,' degr')
- 2 FORMAT ('    Day of the year at model start :   day_init = ',I3             &
-            /'    UTC time at model start        :   time_utc_init = ',F7.1' s')
  3 FORMAT (//' Radiation model information:'/                                  &
               ' ----------------------------'/)
@@ -1401 +1400 @@
        NAMELIST /radiation_par/   albedo, albedo_type, albedo_lw_dir,          &
                                   albedo_lw_dif, albedo_sw_dir, albedo_sw_dif, &
-                                  constant_albedo, day_init, dt_radiation,     &
+                                  constant_albedo, dt_radiation,               &
                                   lambda, lw_radiation, net_radiation,         &
                                   radiation_scheme, skip_time_do_radiation,    &
-                                  sw_radiation, time_utc_init,                 &
-                                  unscheduled_radiation_calls
+                                  sw_radiation, unscheduled_radiation_calls
        
        line = ' '
@@ -1630 +1628 @@
                rrtm_ch4vmr     , rrtm_n2ovmr  , rrtm_o2vmr , rrtm_asdir(:,j,i),&
                rrtm_asdif(:,j,i), rrtm_aldir(:,j,i), rrtm_aldif(:,j,i), zenith,&
-               0.0_wp          , day          , solar_constant,   rrtm_inflgsw,&
+               0.0_wp          , day_of_year  , solar_constant,   rrtm_inflgsw,&
                rrtm_iceflgsw   , rrtm_liqflgsw, rrtm_cldfr , rrtm_sw_taucld   ,&
                rrtm_sw_ssacld  , rrtm_sw_asmcld, rrtm_sw_fsfcld, rrtm_cicewp  ,&
@@ -1693 +1691 @@
 !--    Calculate current day and time based on the initial values and simulation
 !--    time
-       day = day_init + INT(FLOOR( (time_utc_init + time_since_reference_point)&
-                               / 86400.0_wp ), KIND=iwp)
-       time_utc = MOD((time_utc_init + time_since_reference_point), 86400.0_wp)
-
+       CALL calc_date_and_time
 
 !
 !--    Calculate solar declination and hour angle   
-       declination = ASIN( decl_1 * SIN(decl_2 * REAL(day, KIND=wp) - decl_3) )
+       declination = ASIN( decl_1 * SIN(decl_2 * REAL(day_of_year, KIND=wp) - decl_3) )
        hour_angle  = 2.0_wp * pi * (time_utc / 86400.0_wp) + lon - pi