Changeset 4418 for palm/trunk

Timestamp:

Feb 21, 2020 9:41:13 AM (5 years ago)

Author:

raasch

Message:

bugfix for raindrop number adjustment

File:

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

palm/trunk/SOURCE/bulk_cloud_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! bugfix for raindrop number adjustment
+! 
+! 4370 2020-01-10 14:00:44Z raasch
 ! vector directives added to force vectorization on Intel19 compiler
 ! 
@@ -228 +231 @@
     REAL(wp) ::  w_precipitation = 9.65_wp  !< maximum terminal velocity (m s-1)
     REAL(wp) ::  x0 = 2.6E-10_wp           !< separating drop mass (kg)
-    REAL(wp) ::  xamin = 5.24E-19_wp       !< average aerosol mass (kg) (~ 0.05µm)
+!    REAL(wp) ::  xamin = 5.24E-19_wp       !< average aerosol mass (kg) (~ 0.05µm)
     REAL(wp) ::  xcmin = 4.18E-15_wp       !< minimum cloud drop size (kg) (~ 1µm)
     REAL(wp) ::  xrmin = 2.6E-10_wp        !< minimum rain drop size (kg)
@@ -3140 +3143 @@
              ELSE
                 IF ( nc(k,j,i) * xcmin > qc(k,j,i) * hyrho(k) )  THEN
-                   nc(k,j,i) = qc(k,j,i) * hyrho(k) / xamin * flag
+                   nc(k,j,i) = qc(k,j,i) * hyrho(k) / xcmin * flag
                 ENDIF
              ENDIF