Changeset 4290 for palm/trunk

Timestamp:

Nov 11, 2019 12:06:14 PM (5 years ago)

Author:

banzhafs

Message:

Bugfix in sedimentation resistance calculation in drydepo_aero_zhang_vd subroutine

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Bugfix in sedimentation resistance calculation in drydepo_aero_zhang_vd subroutine
+!
+!4273 2019-10-24 13:40:54Z monakurppa
 ! Add logical switches nesting_chem and nesting_offline_chem (both .TRUE.
 ! by default)
@@ -5765 +5768 @@
 !
 !-- sedimentation resistance:
-    rs = 1.0_wp / ( epsilon0 * ustar * ( Ebrown + Eimpac + Einterc ) * Reffic )
+    rs = 1.0_wp / ( epsilon0 * MAX( 1.0e-5, ustar ) * ( Ebrown + Eimpac + Einterc ) * Reffic )
 
 !-- deposition velocity according to Seinfeld and Pandis (2006; eq 19.7):