Changeset 4292 for palm/trunk

Timestamp:

Nov 11, 2019 1:04:50 PM (5 years ago)

Author:

banzhafs

Message:

Bugfix for r4290

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Bugfix for r4290
+!
+! 4290 2019-11-11 12:06:14Z banzhafs $
 ! Bugfix in sedimentation resistance calculation in drydepo_aero_zhang_vd subroutine
 !
-!4273 2019-10-24 13:40:54Z monakurppa
+! 4273 2019-10-24 13:40:54Z monakurppa
 ! Add logical switches nesting_chem and nesting_offline_chem (both .TRUE.
 ! by default)
@@ -5756 +5759 @@
        Reffic = exp( -Stokes**0.5_wp )
     ELSE
-       stokes = vs1 * ustar / (grav * A_lu(luc) * 1.e-3)
-       Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1e-3 ) )**2
+       stokes = vs1 * ustar / (grav * A_lu(luc) * 1.E-3)
+       Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1E-3 ) )**2
        Reffic = exp( -Stokes**0.5_wp )
     END IF
@@ -5768 +5771 @@
 !
 !-- sedimentation resistance:
-    rs = 1.0_wp / ( epsilon0 * MAX( 1.0e-5, ustar ) * ( Ebrown + Eimpac + Einterc ) * Reffic )
+    rs = 1.0_wp / ( epsilon0 * MAX( 1.0E-5_wp, ustar ) * ( Ebrown + Eimpac + Einterc ) * Reffic )
 
 !-- deposition velocity according to Seinfeld and Pandis (2006; eq 19.7):