Changeset 4291 for palm/trunk/SOURCE

Timestamp:

Nov 11, 2019 12:36:54 PM (5 years ago)

Author:

moh.hefny

Message:

RTM is ON in case of biometeorology even if there is no vertical surfaces and/or trees to enable biommet-output in flat surfaces domain

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -27 +27 @@
 ! Former revisions:
 ! -----------------
+!
 ! $Id$
+! Enabled RTM in case of biometeorology even if there is no vertical 
+! surfaces or 3D vegetation in the domain
+! 
+! 4286 2019-10-30 16:01:14Z resler
 ! - Fix wrong treating of time_rad during interpolation in radiation_model_mod
 ! - Fix wrong checks of time_rad from dynamic driver in radiation_model_mod
 ! - Add new directional model of human body for MRT: ellipsoid
-! 
+!
 ! 4271 2019-10-23 10:46:41Z maronga
 ! Bugfix: missing parentheses in calculation of snow albedo
@@ -230 +235 @@
 
     USE control_parameters,                                                    &
-        ONLY:  cloud_droplets, coupling_char,                                  &
+        ONLY:  biometeorology, cloud_droplets, coupling_char,                  &
                debug_output, debug_output_timestep, debug_string,              &
                dt_3d,                                                          &
@@ -1539 +1544 @@
        IF ( debug_output )  CALL debug_message( 'radiation_init', 'start' )
 !
-!--    Activate radiation_interactions according to the existence of vertical surfaces and/or trees.
+!--    Activate radiation_interactions according to the existence of vertical surfaces and/or trees
+!      or if biometeorology output is required for flat surfaces.
 !--    The namelist parameter radiation_interactions_on can override this behavior.
 !--    (This check cannot be performed in check_parameters, because vertical_surfaces_exist is first set in 
 !--    init_surface_arrays.)
        IF ( radiation_interactions_on )  THEN
-          IF ( vertical_surfaces_exist  .OR.  plant_canopy )  THEN
+          IF ( vertical_surfaces_exist  .OR.  plant_canopy  .OR.  biometeorology )  THEN
              radiation_interactions    = .TRUE.
              average_radiation         = .TRUE.
@@ -1551 +1557 @@
                                                    !< calculations necessary in case of flat surface
           ENDIF
-       ELSEIF ( vertical_surfaces_exist  .OR.  plant_canopy )  THEN
+       ELSEIF ( vertical_surfaces_exist  .OR.  plant_canopy  .OR.  biometeorology )  THEN
           message_string = 'radiation_interactions_on is set to .FALSE. although '     // &
-                           'vertical surfaces and/or trees exist. The model will run ' // &
-                           'without RTM (no shadows, no radiation reflections)'
+                           'vertical surfaces and/or trees or biometeorology exist '   // & 
+                           'is ON. The model will run without RTM (no shadows, no '    // &
+                           'radiation reflections)'
           CALL message( 'init_3d_model', 'PA0348', 0, 1, 0, 6, 0 )
        ENDIF