Ignore:
Timestamp:
Apr 17, 2018 10:27:57 AM (4 years ago)
Author:
kanani
Message:

Fixes for radiative transfer model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_3d_model.f90

    r2938 r2977  
    2525! -----------------
    2626! $Id$
     27! Implement changes from branch radiation (r2948-2971) with minor modifications
     28! (moh.hefny):
     29! - set radiation_interactions according to the existence of urban/land vertical
     30!   surfaces and trees to activiate RTM
     31! - set average_radiation to TRUE if RTM is activiated
     32!
     33! 2938 2018-03-27 15:52:42Z suehring
    2734! - Revise Inifor initialization for geostrophic wind components
    2835! - Initialize synthetic turbulence generator in case of Inifor initialization 
     
    506513
    507514    USE radiation_model_mod,                                                   &
    508         ONLY:  radiation_init, radiation, radiation_control, radiation_scheme, &
     515        ONLY:  average_radiation,                                              &
     516               radiation_init, radiation, radiation_control, radiation_scheme, &
    509517               radiation_calc_svf, radiation_write_svf,                        &
    510518               radiation_interaction, radiation_interactions,                  &
    511519               radiation_interaction_init, radiation_read_svf,                 &
    512                radiation_presimulate_solar_pos
     520               radiation_presimulate_solar_pos, radiation_interactions_on
    513521   
    514522    USE random_function_mod
     
    534542    USE surface_mod,                                                           &
    535543        ONLY :  init_surface_arrays, init_surfaces, surf_def_h, surf_lsm_h,    &
    536                 surf_usm_h, get_topography_top_index_ji
     544                surf_usm_h, get_topography_top_index_ji, vertical_surfaces_exist
    537545   
    538546    USE transpose_indices
     
    23572365!-- If required, set chemical emissions
    23582366!-- (todo(FK): This should later on be CALLed time-dependently in init_3d_model)
    2359     IF ( air_chemistry ) THEN
     2367    IF ( air_chemistry )  THEN
    23602368       CALL chem_emissions
    23612369    ENDIF
     
    23642372!-- Initialize radiation processes
    23652373    IF ( radiation )  THEN
     2374!
     2375!--    Activate radiation_interactions according to the existence of vertical surfaces and/or trees.
     2376!--    The namelist parameter radiation_interactions_on can override this behavior.
     2377!--    (This check cannot be performed in check_parameters, because vertical_surfaces_exist is first set in
     2378!--    init_surface_arrays.)
     2379       IF ( radiation_interactions_on )  THEN
     2380          IF ( vertical_surfaces_exist  .OR.  plant_canopy )  THEN
     2381             radiation_interactions    = .TRUE.
     2382             average_radiation         = .TRUE.
     2383          ELSE
     2384             radiation_interactions_on = .FALSE.   !< reset namelist parameter: no interactions
     2385                                                   !< calculations necessary in case of flat surface
     2386          ENDIF
     2387       ELSEIF ( vertical_surfaces_exist  .OR.  plant_canopy )  THEN
     2388          message_string = 'radiation_interactions_on is set to .FALSE. although '     // &
     2389                           'vertical surfaces and/or trees exist. The model will run ' // &
     2390                           'without RTM (no shadows, no radiation reflections)'
     2391          CALL message( 'init_3d_model', 'PA0348', 0, 1, 0, 6, 0 )
     2392       ENDIF
    23662393!
    23672394!--    If required, initialize radiation interactions between surfaces
Note: See TracChangeset for help on using the changeset viewer.