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

Timestamp:

Jul 25, 2018 3:19:37 PM (7 years ago)

Author:

suehring

Message:

Bugfix in radiation forcing in case of RRTMG; further bugfix in output of surface variables in case of overhanging structures

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Bugfix, map signle-column radiation forcing profiles on top of any topography
+! 
+! 3156 2018-07-19 16:30:54Z knoop
 ! Bugfix: replaced usage of the pt array with the surf%pt_surface array
 !
@@ -2997 +3000 @@
 
 !
-!--          Save heating rates (convert from K/d to K/h)
-             DO k = nzb+1, nzt+1
-                rad_lw_hr(k,:,:)     = rrtm_lwhr(0,k)  * d_hours_day
-                rad_lw_cs_hr(k,:,:)  = rrtm_lwhrc(0,k) * d_hours_day
+!--          Save heating rates (convert from K/d to K/h). 
+!--          Further, even though an aggregated radiation is computed, map
+!--          signle-column profiles on top of any topography, in order to 
+!--          obtain correct near surface radiation heating/cooling rates.
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   k_topo = get_topography_top_index_ji( j, i, 's' )
+                   DO k = k_topo+1, nzt+1
+                      rad_lw_hr(k,j,i)     = rrtm_lwhr(0,k)  * d_hours_day
+                      rad_lw_cs_hr(k,j,i)  = rrtm_lwhrc(0,k) * d_hours_day
+                   ENDDO
+                ENDDO
              ENDDO