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

Timestamp:

Jan 21, 2019 1:02:11 AM (5 years ago)

Author:

knoop

Message:

Some interface calls moved to module_interface + cleanup

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Some interface calls moved to module_interface + cleanup
+! 
+! 3667 2019-01-10 14:26:24Z schwenkel
 ! Modified check for rrtmg input files
 ! 
@@ -567 +570 @@
                message_string, plant_canopy, pt_surface,                       &
                rho_surface, simulated_time, spinup_time, surface_pressure,     &
+               read_svf, write_svf,                                            &
                time_since_reference_point, urban_surface, varnamelength
 
@@ -630 +634 @@
         ONLY:  get_topography_top_index, get_topography_top_index_ji,          &
                ind_pav_green, ind_veg_wall, ind_wat_win,                       &
-               surf_lsm_h, surf_lsm_v, surf_type, surf_usm_h, surf_usm_v
+               surf_lsm_h, surf_lsm_v, surf_type, surf_usm_h, surf_usm_v,      &
+               vertical_surfaces_exist
 
     IMPLICIT NONE
@@ -1739 +1744 @@
        INTEGER(iwp) ::  ind_type  !< running index for subgrid-surface tiles
 #endif
+
+!
+!--    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 
+!--    via sky-view factors. This must be done before radiation is initialized.
+       IF ( radiation_interactions )  CALL radiation_interaction_init
+
+!
+!--    Initialize radiation model
+       CALL location_message( 'initializing radiation model', .FALSE. )
 
 !
@@ -2671 +2704 @@
        CALL init_date_and_time
 
+       CALL location_message( 'finished', .TRUE. )
+
+!
+!--    Find all discretized apparent solar positions for radiation interaction.
+       IF ( radiation_interactions )  CALL radiation_presimulate_solar_pos
+
+!
+!--    If required, read or calculate and write out the SVF 
+       IF ( radiation_interactions .AND. read_svf)  THEN
+!
+!--       Read sky-view factors and further required data from file
+          CALL location_message( '    Start reading SVF from file', .FALSE. )
+          CALL radiation_read_svf()
+          CALL location_message( '    Reading SVF from file has finished', .TRUE. )
+
+       ELSEIF ( radiation_interactions .AND. .NOT. read_svf)  THEN
+!
+!--       calculate SFV and CSF
+          CALL location_message( '    Start calculation of SVF', .FALSE. )
+          CALL radiation_calc_svf()
+          CALL location_message( '    Calculation of SVF has finished', .TRUE. )
+       ENDIF
+
+       IF ( radiation_interactions .AND. write_svf)  THEN
+!
+!--       Write svf, csf svfsurf and csfsurf data to file
+          CALL location_message( '    Start writing SVF in file', .FALSE. )
+          CALL radiation_write_svf()
+          CALL location_message( '    Writing SVF in file has finished', .TRUE. )
+       ENDIF
+
+!
+!--    Adjust radiative fluxes. In case of urban and land surfaces, also 
+!--    call an initial interaction.
+       IF ( radiation_interactions )  THEN
+          CALL radiation_interaction
+       ENDIF
+
        RETURN
 
@@ -5838 +5909 @@
            !WRITE(message_string, '(a,f6.3)') 'Precomputing effective box optical ' &
            !    // 'depth using prototype leaf area density = ', prototype_lad
-           !CALL message('usm_init_urban_surface', 'PA0520', 0, 0, -1, 6, 0)
+           !CALL message('radiation_interaction_init', 'PA0520', 0, 0, -1, 6, 0)
        ENDIF