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

Timestamp:

Jun 1, 2017 2:12:31 PM (8 years ago)

Author:

maronga

Message:

revised number of soil layers. added support for RRTMG runs with dry atmosphere/soil

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
-!
+! Allow for RRTMG runs without humidity/cloud physics
+! 
+! 2233 2017-05-30 18:08:54Z suehring
+! 
 ! 2232 2017-05-30 17:47:52Z suehring
 ! Adjustments to new topography concept
@@ -1415 +1418 @@
                                                   / 1000.0_wp )**0.286_wp
 
-             DO k = nzb+1, nzt+1
-                rrtm_tlay(0,k) = pt(k,j,i) * ( (hyp(k) ) / 100000.0_wp         &
-                                 )**0.286_wp + l_d_cp * ql(k,j,i)
-                rrtm_h2ovmr(0,k) = mol_mass_air_d_wv * (q(k,j,i) - ql(k,j,i))
-
-             ENDDO
+
+             IF ( cloud_physics )  THEN
+                DO k = nzb+1, nzt+1
+                   rrtm_tlay(0,k) = pt(k,j,i) * ( (hyp(k) ) / 100000.0_wp      &
+                                    )**0.286_wp + l_d_cp * ql(k,j,i)
+                   rrtm_h2ovmr(0,k) = mol_mass_air_d_wv * (q(k,j,i) - ql(k,j,i))
+                ENDDO
+             ELSE
+                DO k = nzb+1, nzt+1
+                   rrtm_tlay(0,k) = pt(k,j,i) * ( (hyp(k) ) / 100000.0_wp      &
+                                    )**0.286_wp
+                   rrtm_h2ovmr(0,k) = 0.0_wp
+                ENDDO
+             ENDIF
 
 !
@@ -1455 +1466 @@
              rrtm_icld   = 0
 
-             DO k = nzb+1, nzt+1
-                rrtm_cliqwp(0,k) =  ql(k,j,i) * 1000.0_wp *                    &
-                                    (rrtm_plev(0,k) - rrtm_plev(0,k+1))        &
-                                    * 100.0_wp / g 
-
-                IF ( rrtm_cliqwp(0,k) > 0.0_wp )  THEN
-                   rrtm_cldfr(0,k) = 1.0_wp
-                   IF ( rrtm_icld == 0 )  rrtm_icld = 1
-
-!
-!--                Calculate cloud droplet effective radius
-                   IF ( cloud_physics )  THEN
-                      rrtm_reliq(0,k) = 1.0E6_wp * ( 3.0_wp * ql(k,j,i)        &
-                                        * rho_surface                          &
-                                        / ( 4.0_wp * pi * nc_const * rho_l )   &
-                                        )**0.33333333333333_wp                 &
-                                        * EXP( LOG( sigma_gc )**2 )
-
-                   ELSEIF ( cloud_droplets )  THEN
-                      number_of_particles = prt_count(k,j,i)
-
-                      IF (number_of_particles <= 0)  CYCLE
-                      particles => grid_particles(k,j,i)%particles(1:number_of_particles)
-                      s_r2 = 0.0_wp
-                      s_r3 = 0.0_wp
-
-                      DO  n = 1, number_of_particles
-                         IF ( particles(n)%particle_mask )  THEN
-                            s_r2 = s_r2 + particles(n)%radius**2 * &
-                                   particles(n)%weight_factor
-                            s_r3 = s_r3 + particles(n)%radius**3 * &
-                                   particles(n)%weight_factor
-                         ENDIF
-                      ENDDO
-
-                      IF ( s_r2 > 0.0_wp )  rrtm_reliq(0,k) = s_r3 / s_r2
-
+             IF ( cloud_physics )  THEN
+                DO k = nzb+1, nzt+1
+                   rrtm_cliqwp(0,k) =  ql(k,j,i) * 1000.0_wp *                 &
+                                       (rrtm_plev(0,k) - rrtm_plev(0,k+1))     &
+                                       * 100.0_wp / g 
+
+                   IF ( rrtm_cliqwp(0,k) > 0.0_wp )  THEN
+                      rrtm_cldfr(0,k) = 1.0_wp
+                      IF ( rrtm_icld == 0 )  rrtm_icld = 1
+
+!
+!--                   Calculate cloud droplet effective radius
+                      IF ( cloud_physics )  THEN
+                         rrtm_reliq(0,k) = 1.0E6_wp * ( 3.0_wp * ql(k,j,i)        &
+                                           * rho_surface                          &
+                                           / ( 4.0_wp * pi * nc_const * rho_l )   &
+                                           )**0.33333333333333_wp                 &
+                                           * EXP( LOG( sigma_gc )**2 )
+
+                      ELSEIF ( cloud_droplets )  THEN
+                         number_of_particles = prt_count(k,j,i)
+
+                         IF (number_of_particles <= 0)  CYCLE
+                         particles => grid_particles(k,j,i)%particles(1:number_of_particles)
+                         s_r2 = 0.0_wp
+                         s_r3 = 0.0_wp
+
+                         DO  n = 1, number_of_particles
+                            IF ( particles(n)%particle_mask )  THEN
+                               s_r2 = s_r2 + particles(n)%radius**2 * &
+                                      particles(n)%weight_factor
+                               s_r3 = s_r3 + particles(n)%radius**3 * &
+                                      particles(n)%weight_factor
+                            ENDIF
+                         ENDDO
+
+                         IF ( s_r2 > 0.0_wp )  rrtm_reliq(0,k) = s_r3 / s_r2
+
+                      ENDIF
+
+!
+!--                   Limit effective radius
+                      IF ( rrtm_reliq(0,k) > 0.0_wp )  THEN
+                         rrtm_reliq(0,k) = MAX(rrtm_reliq(0,k),2.5_wp)
+                         rrtm_reliq(0,k) = MIN(rrtm_reliq(0,k),60.0_wp)
+                     ENDIF
                    ENDIF
-
-!
-!--                Limit effective radius
-                   IF ( rrtm_reliq(0,k) > 0.0_wp )  THEN
-                      rrtm_reliq(0,k) = MAX(rrtm_reliq(0,k),2.5_wp)
-                      rrtm_reliq(0,k) = MIN(rrtm_reliq(0,k),60.0_wp)
-                  ENDIF
-                ENDIF
-             ENDDO
+                ENDDO
+             ENDIF
 
 !