Changeset 2604 for palm

Timestamp:

Nov 6, 2017 1:29:00 PM (6 years ago)

Author:

schwenkel

Message:

last commit documented

File:

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

palm/trunk/SOURCE/radiation_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! bugfix for calculation of effective radius using morrison microphysics
+! 
+! 2601 2017-11-02 16:22:46Z scharf
 ! added emissivity to namelist
 ! 
@@ -207 +210 @@
 
     USE microphysics_mod,                                                      &
-        ONLY:  nc_const, sigma_gc
+        ONLY:  na_init, nc_const, sigma_gc
 
 #if defined ( __netcdf )
@@ -1556 +1559 @@
 !--                   Calculate cloud droplet effective radius
                       IF ( cloud_physics )  THEN
-                         nc_rad = MERGE( nc(k,j,i), nc_const, microphysics_morrison )
+!
+!--                      Calculete effective droplet radius. In case of using 
+!--                      cloud_scheme = 'morrison' and a non reasonable number 
+!--                      of cloud droplets the inital aerosol number  
+!--                      concentration is considered.
+                         IF ( microphysics_morrison ) THEN
+                            IF ( nc(k,j,i) > 1.0E-20_wp ) THEN
+                               nc_rad = nc(k,j,i)
+                            ELSE
+                               nc_rad = na_init
+                            ENDIF 
+                         ELSE
+                            nc_rad = nc_const
+                         ENDIF  
+
                          rrtm_reliq(0,k) = 1.0E6_wp * ( 3.0_wp * ql(k,j,i)     &
                                            * rho_surface                       &