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

Timestamp:

Oct 30, 2018 6:08:55 PM (5 years ago)

Author:

kanani

Message:

from branch resler@3462: add MRT shaping function (radiation_model_mod), use basic constants (biometeorology_mod), adjust precision to wp (biometeorology_pt_mod)

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! From branch resler@3462, pavelkrc:
+! add MRT shaping function for human
+! 
+! 3449 2018-10-29 19:36:56Z suehring
 ! New RTM version 3.0: (Pavel Krc, Jaroslav Resler, ICS, Prague)
 !   - Interaction of plant canopy with LW radiation
@@ -902 +906 @@
     LOGICAL                                        ::  mrt_skip_roof = .TRUE.             !< do not calculate MRT above roof surfaces
     LOGICAL                                        ::  mrt_include_sw = .TRUE.            !< should MRT calculation include SW radiation as well?
+    LOGICAL                                        ::  mrt_geom_human = .TRUE.            !< MRT direction weights simulate human body instead of a sphere
     INTEGER(iwp)                                   ::  nrefsteps = 3                      !< number of reflection steps to perform
     REAL(wp), PARAMETER                            ::  ext_coef = 0.6_wp                  !< extinction coefficient (a.k.a. alpha)
@@ -1172 +1177 @@
 !-- Public variables and constants / NEEDS SORTING
     PUBLIC albedo, albedo_type, decl_1, decl_2, decl_3, dots_rad, dt_radiation,&
-           emissivity, force_radiation_call, lat, lon,                         &
+           emissivity, force_radiation_call, lat, lon, mrt_geom_human,         &
            mrt_include_sw, mrt_nlevels, mrtbl, mrtinsw, mrtinlw, nmrtbl,       &
            rad_net_av, radiation, radiation_scheme, rad_lw_in,                 &
@@ -2942 +2947 @@
                                   constant_albedo, dt_radiation, emissivity,    &
                                   lw_radiation, max_raytracing_dist,            &
-                                  min_irrf_value, mrt_include_sw, mrt_nlevels,  &
+                                  min_irrf_value, mrt_geom_human,               &
+                                  mrt_include_sw, mrt_nlevels,                  &
                                   mrt_skip_roof, net_radiation, nrefsteps,      &
                                   plant_lw_interact, rad_angular_discretization,&
@@ -2957 +2963 @@
                                   constant_albedo, dt_radiation, emissivity,    &
                                   lw_radiation, max_raytracing_dist,            &
-                                  min_irrf_value, mrt_include_sw, mrt_nlevels,  &
+                                  min_irrf_value, mrt_geom_human,               &
+                                  mrt_include_sw, mrt_nlevels,                  &
                                   mrt_skip_roof, net_radiation, nrefsteps,      &
                                   plant_lw_interact, rad_angular_discretization,&
@@ -6553 +6560 @@
            zbdry(:) = (/( zn0+REAL(izn,wp)*zns, izn=0, nzn )/)
            vffrac0(:) = (COS(zbdry(0:nzn-1)) - COS(zbdry(1:nzn))) / 2._wp / REAL(naz, wp)
+           !
+           !--Modify direction weights to simulate human body (lower weight for top-down)
+           IF ( mrt_geom_human )  THEN
+              vffrac0(:) = vffrac0(:) * MAX(0._wp, SIN(zcent(:))*0.88_wp + COS(zcent(:))*0.12_wp)
+              vffrac0(:) = vffrac0(:) / (SUM(vffrac0) * REAL(naz, wp))
+           ENDIF
 
            DO  imrt = 1, nmrtbl