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

Timestamp:

Aug 8, 2019 11:26:00 AM (5 years ago)

Author:

suehring

Message:

indoor_model_mod: Revision of indoor model; urban_surface_mod: parameters for waste heat from cooling and heating are introduced to building data base; initialization of building data base moved to an earlier point of time before indoor model will be initialized; radiation_model_mod: minor improvement in some comments; synthetic_turbulence_generator: unused variable removed

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Comments added
+! 
+! 4134 2019-08-02 18:39:57Z suehring
 ! Bugfix in formatted write statement
 ! 
@@ -3569 +3572 @@
 
        IF ( average_radiation ) THEN
-
+!
+!--       Determine minimum topography top index. 
           k_topo_l = MINVAL( get_topography_top_index( 's' ) )
 #if defined( __parallel )
@@ -3705 +3709 @@
                              rrtm_lw_tauaer(0:0,k_topo+1:nzt_rad+1,1:nbndlw+1)
           
-!              CALL rrtmg_lw( 1, nzt_rad      , rrtm_icld    , rrtm_idrv      ,&
-!              rrtm_play       , rrtm_plev    , rrtm_tlay    , rrtm_tlev      ,&
-!              rrtm_tsfc       , rrtm_h2ovmr  , rrtm_o3vmr   , rrtm_co2vmr    ,&
-!              rrtm_ch4vmr     , rrtm_n2ovmr  , rrtm_o2vmr   , rrtm_cfc11vmr  ,&
-!              rrtm_cfc12vmr   , rrtm_cfc22vmr, rrtm_ccl4vmr , rrtm_emis      ,&
-!              rrtm_inflglw    , rrtm_iceflglw, rrtm_liqflglw, rrtm_cldfr     ,&
-!              rrtm_lw_taucld  , rrtm_cicewp  , rrtm_cliqwp  , rrtm_reice     ,& 
-!              rrtm_reliq      , rrtm_lw_tauaer,                               &
-!              rrtm_lwuflx     , rrtm_lwdflx  , rrtm_lwhr  ,                   &
-!              rrtm_lwuflxc    , rrtm_lwdflxc , rrtm_lwhrc ,                   &
-!              rrtm_lwuflx_dt  ,  rrtm_lwuflxc_dt )
-
              CALL rrtmg_lw( 1,                                                 &                                        
                             nzt_rad-k_topo,                                    &
@@ -3809 +3801 @@
              rrtm_sw_asmaer_dum = rrtm_sw_asmaer(0:0,k_topo+1:nzt_rad+1,1:nbndsw+1)
              rrtm_sw_ecaer_dum  = rrtm_sw_ecaer(0:0,k_topo+1:nzt_rad+1,1:naerec+1)
-!              CALL rrtmg_sw( 1, nzt_rad      , rrtm_icld     , rrtm_iaer      ,&
-!              rrtm_play      , rrtm_plev     , rrtm_tlay     , rrtm_tlev      ,&
-!              rrtm_tsfc      , rrtm_h2ovmr   , rrtm_o3vmr    , rrtm_co2vmr    ,&
-!              rrtm_ch4vmr    , rrtm_n2ovmr   , rrtm_o2vmr    , rrtm_asdir     ,&
-!              rrtm_asdif     , rrtm_aldir    , rrtm_aldif    , zenith         ,&
-!              0.0_wp         , day_of_year   , solar_constant, rrtm_inflgsw   ,&
-!              rrtm_iceflgsw  , rrtm_liqflgsw , rrtm_cldfr    , rrtm_sw_taucld ,&
-!              rrtm_sw_ssacld , rrtm_sw_asmcld, rrtm_sw_fsfcld, rrtm_cicewp    ,&
-!              rrtm_cliqwp    , rrtm_reice    , rrtm_reliq    , rrtm_sw_tauaer ,&
-!              rrtm_sw_ssaaer , rrtm_sw_asmaer, rrtm_sw_ecaer , rrtm_swuflx    ,&
-!              rrtm_swdflx    , rrtm_swhr     , rrtm_swuflxc  , rrtm_swdflxc   ,&
-!              rrtm_swhrc     , rrtm_dirdflux , rrtm_difdflux )
+
              CALL rrtmg_sw( 1,                                                 &
                             nzt_rad-k_topo,                                    &
@@ -3879 +3860 @@
   
 !
-!--          Save fluxes:
-!--          - whole domain
+!--          Save radiation fluxes for the entire depth of the model domain
              DO k = nzb, nzt+1
                 rad_sw_in(k,:,:)  = rrtm_swdflx(0,k)
                 rad_sw_out(k,:,:) = rrtm_swuflx(0,k)
              ENDDO
-!--          - direct and diffuse SW at urban-surface-layer (required by RTM)
+!--          Save direct and diffuse SW radiation at the surface (required by RTM)
              rad_sw_in_dir(:,:) = rrtm_dirdflux(0,k_topo)
              rad_sw_in_diff(:,:) = rrtm_difdflux(0,k_topo)