Changeset 4306 for palm


Ignore:
Timestamp:
Nov 25, 2019 12:04:48 PM (5 years ago)
Author:
banzhafs
Message:

Correction to r4304 commit in chemistry_model_mod

File:
1 edited

Legend:

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

    r4304 r4306  
    2727! -----------------
    2828! $Id$
     29! Corretion for r4304 commit
     30!
     31! 4304 2019-11-25 10:43:03Z banzhafs
    2932! Precision clean-up in drydepo_aero_zhang_vd subroutine
    3033!
     
    56835686
    56845687!-- sedimentation velocity
    5685     vs = rhopart * ( partsize**2.0_wp ) * grav * slipcor / ( 18.0_wp * visc )
     5688    vs = rhopart * ( partsize**2 ) * grav * slipcor / ( 18.0_wp * visc )
    56865689
    56875690 END FUNCTION sedimentation_velocity
     
    57165719    REAL(wp), PARAMETER ::  grav     = 9.80665_wp             !< acceleration of gravity (m/s2)
    57175720
    5718     REAL(wp), PARAMETER ::  beta     = 2.0_wp
    57195721    REAL(wp), PARAMETER ::  epsilon0 = 3.0_wp
    57205722    REAL(wp), PARAMETER ::  kb       = 1.38066E-23_wp
     
    57575759!-- and sticking efficiency R (1 = no rebound)
    57585760    IF ( luc == ilu_ice .OR. nwet==9 .OR. luc == ilu_water_sea .OR. luc == ilu_water_inland )  THEN
    5759        stokes = vs1 * ustar**2.0_wp / ( grav * kinvisc )
     5761       stokes = vs1 * ustar**2 / ( grav * kinvisc )
    57605762       Einterc = 0.0_wp
    57615763       Reffic = 1.0_wp
    57625764    ELSE IF ( luc == ilu_other .OR. luc == ilu_desert )  THEN     !<tundra of desert
    5763        stokes = vs1 * ustar**2.0_wp / ( grav * kinvisc )
     5765       stokes = vs1 * ustar**2 / ( grav * kinvisc )
    57645766       Einterc = 0.0_wp
    57655767       Reffic = exp( -Stokes**0.5_wp )
    57665768    ELSE
    57675769       stokes = vs1 * ustar / ( grav * A_lu(luc) * 1.0E-3_wp )
    5768        Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1.0E-3_wp ) )**2.0_wp
     5770       Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1.0E-3_wp ) )**2
    57695771       Reffic = exp( -Stokes**0.5_wp )
    57705772    END IF
     
    57745776!
    57755777!-- determine impaction efficiency:
    5776     Eimpac = ( stokes / ( alfa_lu(luc) + stokes ) )**beta
     5778    Eimpac = ( stokes / ( alfa_lu(luc) + stokes ) )**2
    57775779!
    57785780!-- sedimentation resistance:
Note: See TracChangeset for help on using the changeset viewer.