Changeset 4292 for palm


Ignore:
Timestamp:
Nov 11, 2019 1:04:50 PM (5 years ago)
Author:
banzhafs
Message:

Bugfix for r4290

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/chemistry_model_mod.f90

    r4290 r4292  
    2727! -----------------
    2828! $Id$
     29! Bugfix for r4290
     30!
     31! 4290 2019-11-11 12:06:14Z banzhafs $
    2932! Bugfix in sedimentation resistance calculation in drydepo_aero_zhang_vd subroutine
    3033!
    31 !4273 2019-10-24 13:40:54Z monakurppa
     34! 4273 2019-10-24 13:40:54Z monakurppa
    3235! Add logical switches nesting_chem and nesting_offline_chem (both .TRUE.
    3336! by default)
     
    57565759       Reffic = exp( -Stokes**0.5_wp )
    57575760    ELSE
    5758        stokes = vs1 * ustar / (grav * A_lu(luc) * 1.e-3)
    5759        Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1e-3 ) )**2
     5761       stokes = vs1 * ustar / (grav * A_lu(luc) * 1.E-3)
     5762       Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1E-3 ) )**2
    57605763       Reffic = exp( -Stokes**0.5_wp )
    57615764    END IF
     
    57685771!
    57695772!-- sedimentation resistance:
    5770     rs = 1.0_wp / ( epsilon0 * MAX( 1.0e-5, ustar ) * ( Ebrown + Eimpac + Einterc ) * Reffic )
     5773    rs = 1.0_wp / ( epsilon0 * MAX( 1.0E-5_wp, ustar ) * ( Ebrown + Eimpac + Einterc ) * Reffic )
    57715774
    57725775!-- deposition velocity according to Seinfeld and Pandis (2006; eq 19.7):
Note: See TracChangeset for help on using the changeset viewer.