Changeset 4306

Timestamp:

Nov 25, 2019 12:04:48 PM (5 years ago)

Author:

banzhafs

Message:

Correction to r4304 commit in chemistry_model_mod

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Corretion for r4304 commit
+!
+! 4304 2019-11-25 10:43:03Z banzhafs
 ! Precision clean-up in drydepo_aero_zhang_vd subroutine
 !
@@ -5683 +5686 @@
 
 !-- sedimentation velocity
-    vs = rhopart * ( partsize**2.0_wp ) * grav * slipcor / ( 18.0_wp * visc )
+    vs = rhopart * ( partsize**2 ) * grav * slipcor / ( 18.0_wp * visc )
 
  END FUNCTION sedimentation_velocity
@@ -5716 +5719 @@
     REAL(wp), PARAMETER ::  grav     = 9.80665_wp             !< acceleration of gravity (m/s2)
 
-    REAL(wp), PARAMETER ::  beta     = 2.0_wp
     REAL(wp), PARAMETER ::  epsilon0 = 3.0_wp
     REAL(wp), PARAMETER ::  kb       = 1.38066E-23_wp
@@ -5757 +5759 @@
 !-- and sticking efficiency R (1 = no rebound)
     IF ( luc == ilu_ice .OR. nwet==9 .OR. luc == ilu_water_sea .OR. luc == ilu_water_inland )  THEN
-       stokes = vs1 * ustar**2.0_wp / ( grav * kinvisc )
+       stokes = vs1 * ustar**2 / ( grav * kinvisc )
        Einterc = 0.0_wp
        Reffic = 1.0_wp
     ELSE IF ( luc == ilu_other .OR. luc == ilu_desert )  THEN     !<tundra of desert
-       stokes = vs1 * ustar**2.0_wp / ( grav * kinvisc )
+       stokes = vs1 * ustar**2 / ( grav * kinvisc )
        Einterc = 0.0_wp
        Reffic = exp( -Stokes**0.5_wp )
     ELSE
        stokes = vs1 * ustar / ( grav * A_lu(luc) * 1.0E-3_wp )
-       Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1.0E-3_wp ) )**2.0_wp
+       Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1.0E-3_wp ) )**2
        Reffic = exp( -Stokes**0.5_wp )
     END IF
@@ -5774 +5776 @@
 !
 !-- determine impaction efficiency:
-    Eimpac = ( stokes / ( alfa_lu(luc) + stokes ) )**beta
+    Eimpac = ( stokes / ( alfa_lu(luc) + stokes ) )**2
 !
 !-- sedimentation resistance: