Changeset 2977 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Apr 17, 2018 10:27:57 AM (6 years ago)

Author:

kanani

Message:

Fixes for radiative transfer model

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (5 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Implement changes from branch radiation (r2948-2971) with minor modifications
+! (moh.hefny):
+! - set radiation_interactions according to the existence of urban/land vertical
+!   surfaces and trees to activiate RTM
+! - set average_radiation to TRUE if RTM is activiated
+! 
+! 2938 2018-03-27 15:52:42Z suehring
 ! - Revise Inifor initialization for geostrophic wind components
 ! - Initialize synthetic turbulence generator in case of Inifor initialization  
@@ -506 +513 @@
 
     USE radiation_model_mod,                                                   &
-        ONLY:  radiation_init, radiation, radiation_control, radiation_scheme, &
+        ONLY:  average_radiation,                                              &
+               radiation_init, radiation, radiation_control, radiation_scheme, &
                radiation_calc_svf, radiation_write_svf,                        &
                radiation_interaction, radiation_interactions,                  &
                radiation_interaction_init, radiation_read_svf,                 &
-               radiation_presimulate_solar_pos
+               radiation_presimulate_solar_pos, radiation_interactions_on
     
     USE random_function_mod 
@@ -534 +542 @@
     USE surface_mod,                                                           &
         ONLY :  init_surface_arrays, init_surfaces, surf_def_h, surf_lsm_h,    &
-                surf_usm_h, get_topography_top_index_ji
+                surf_usm_h, get_topography_top_index_ji, vertical_surfaces_exist
    
     USE transpose_indices
@@ -2357 +2365 @@
 !-- If required, set chemical emissions
 !-- (todo(FK): This should later on be CALLed time-dependently in init_3d_model)
-    IF ( air_chemistry ) THEN
+    IF ( air_chemistry )  THEN
        CALL chem_emissions
     ENDIF
@@ -2364 +2372 @@
 !-- Initialize radiation processes
     IF ( radiation )  THEN
+!
+!--    Activate radiation_interactions according to the existence of vertical surfaces and/or trees.
+!--    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
+             radiation_interactions    = .TRUE.
+             average_radiation         = .TRUE.
+          ELSE
+             radiation_interactions_on = .FALSE.   !< reset namelist parameter: no interactions 
+                                                   !< calculations necessary in case of flat surface
+          ENDIF
+       ELSEIF ( vertical_surfaces_exist  .OR.  plant_canopy )  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)'
+          CALL message( 'init_3d_model', 'PA0348', 0, 1, 0, 6, 0 )
+       ENDIF
 !
 !--    If required, initialize radiation interactions between surfaces