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

Timestamp:

Oct 16, 2017 12:41:56 PM (7 years ago)

Author:

schwenkel

Message:

extended by cloud_droplets option

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! extended by cloud_droplets option, minor bugfix and correct calculation of
+! cloud droplet number concentration
+! 
+! 2544 2017-10-13 18:09:32Z maronga
 ! Moved date and time quantitis to separate module date_and_time_mod 
 ! 
@@ -173 +177 @@
  
     USE arrays_3d,                                                             &
-        ONLY:  dzw, hyp, pt, q, ql, zu, zw
+        ONLY:  dzw, hyp, nc, pt, q, ql, zu, zw
 
     USE cloud_parameters,                                                      &
@@ -183 +187 @@
     USE control_parameters,                                                    &
         ONLY:  cloud_droplets, cloud_physics, g, initializing_actions,         &
-               large_scale_forcing, lsf_surf, phi, pt_surface, rho_surface,    &
-               surface_pressure, time_since_reference_point
+               large_scale_forcing, lsf_surf, microphysics_morrison, phi,      & 
+               pt_surface, rho_surface, surface_pressure,                      &
+               time_since_reference_point
 
     USE date_and_time_mod,                                                     &
@@ -1252 +1257 @@
              rad_lw_out(0,j,i) = emis(j,i) * sigma_sb * (pt(k,j,i) * exn)**4
 
-             IF ( cloud_physics )  THEN
+             IF ( cloud_physics  .OR.  cloud_droplets )  THEN
                 pt1 = pt(k+1,j,i) + l_d_cp / exn1 * ql(k+1,j,i)
                 rad_lw_in(0,j,i)  = 0.8_wp * sigma_sb * (pt1 * exn1)**4
@@ -1303 +1308 @@
              exn1 = (hyp(k+1) / 100000.0_wp )**0.286_wp
 
-             IF ( cloud_physics )  THEN
+             IF ( cloud_physics  .OR.  cloud_droplets )  THEN
                 pt1 = pt(k+1,j,i) + l_d_cp / exn1 * ql(k+1,j,i)
                 rad_lw_in(0,j,i)  = 0.8_wp * sigma_sb * (pt1 * exn1)**4
@@ -1450 +1455 @@
        INTEGER(iwp) :: i, j, k, n !< loop indices
 
-       REAL(wp)     ::  s_r2, &   !< weighted sum over all droplets with r^2
-                        s_r3      !< weighted sum over all droplets with r^3
+       REAL(wp)     ::  nc_rad, &    !< number concentration of cloud droplets
+                        s_r2,   &    !< weighted sum over all droplets with r^2
+                        s_r3         !< weighted sum over all droplets with r^3
 
 !
@@ -1490 +1496 @@
                    rrtm_h2ovmr(0,k) = mol_mass_air_d_wv * (q(k,j,i) - ql(k,j,i))
                 ENDDO
+             ELSEIF ( cloud_droplets ) 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) 
+                ENDDO     
              ELSE
                 DO k = nzb+1, nzt+1
@@ -1531 +1543 @@
              rrtm_icld   = 0
 
-             IF ( cloud_physics )  THEN
+             IF ( cloud_physics  .OR.  cloud_droplets )  THEN
                 DO k = nzb+1, nzt+1
                    rrtm_cliqwp(0,k) =  ql(k,j,i) * 1000.0_wp *                 &
@@ -1544 +1556 @@
 !--                   Calculate cloud droplet effective radius
                       IF ( cloud_physics )  THEN
+                         nc_rad = MERGE( nc(k,j,i), nc_const, microphysics_morrison )
                          rrtm_reliq(0,k) = 1.0E6_wp * ( 3.0_wp * ql(k,j,i)     &
                                            * rho_surface                       &
-                                           / ( 4.0_wp * pi * nc_const * rho_l )&
+                                           / ( 4.0_wp * pi * nc_rad * rho_l )&
                                            )**0.33333333333333_wp              &
                                            * EXP( LOG( sigma_gc )**2 )