Changeset 3859 for palm/trunk/SOURCE

Timestamp:

Apr 3, 2019 8:30:31 PM (5 years ago)

Author:

maronga

Message:

comments in radiation model updated, minor bugfix in palm_csd

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -28 +28 @@
 ! -----------------
 ! $Id$
+! Added some descriptions
+! 
+! 3847 2019-04-01 14:51:44Z suehring
 ! Implement check for dt_radiation (must be > 0)
 ! 
@@ -831 +834 @@
 
 !
-!-- Flag parameters for RRTMGS (should not be changed)
+!-- Flag parameters to be passed to RRTMG (should not be changed until ice phase in clouds is allowed)
     INTEGER(iwp), PARAMETER :: rrtm_idrv     = 1, & !< flag for longwave upward flux calculation option (0,1)
                                rrtm_inflglw  = 2, & !< flag for lw cloud optical properties (0,1,2)
@@ -3492 +3495 @@
 
 !
-!--       Avoid temperature/humidity jumps at the top of the LES domain by 
-!--       linear interpolation from nzt+2 to nzt+7
+!--       Avoid temperature/humidity jumps at the top of the PALM domain by 
+!--       linear interpolation from nzt+2 to nzt+7. Jumps are induced by
+!--       discrepancies between the values in the  domain and those above that
+!--       are prescribed in RRTMG
           DO k = nzt+2, nzt+7
              rrtm_tlay(0,k) = rrtm_tlay(0,nzt+1)                            &
@@ -3507 +3512 @@
           ENDDO
 
-!--       Linear interpolate to zw grid
+!--       Linear interpolate to zw grid. Loop reaches one level further up
+!--       due to the staggered grid in RRTMG
           DO k = nzb+2, nzt+8
              rrtm_tlev(0,k)   = rrtm_tlay(0,k-1) + (rrtm_tlay(0,k) -        &
@@ -4614 +4620 @@
 ! Description:
 ! ------------
-!> Read trace gas data from file
+!> Read trace gas data from file and convert into trace gas paths / volume
+!> mixing ratios. If a user-defined input file is provided it needs to follow
+!> the convections used in RRTMG (see respective netCDF files shipped with 
+!> RRTMG)
 !------------------------------------------------------------------------------!
     SUBROUTINE read_trace_gas_data
@@ -4638 +4647 @@
                        id_var    !< NetCDf id ot the absorber
 
-       REAL(wp) :: p_mls_l, p_mls_u, p_wgt_l, p_wgt_u, p_mls_m
+       REAL(wp) :: p_mls_l, &    !< pressure lower limit for interpolation
+                   p_mls_u, &    !< pressure upper limit for interpolation
+                   p_wgt_l, &    !< pressure weight lower limit for interpolation
+                   p_wgt_u, &    !< pressure weight upper limit for interpolation
+                   p_mls_m       !< mean pressure between upper and lower limits