- Timestamp:
- Nov 11, 2019 1:04:50 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/chemistry_model_mod.f90
r4290 r4292 27 27 ! ----------------- 28 28 ! $Id$ 29 ! Bugfix for r4290 30 ! 31 ! 4290 2019-11-11 12:06:14Z banzhafs $ 29 32 ! Bugfix in sedimentation resistance calculation in drydepo_aero_zhang_vd subroutine 30 33 ! 31 ! 4273 2019-10-24 13:40:54Z monakurppa34 ! 4273 2019-10-24 13:40:54Z monakurppa 32 35 ! Add logical switches nesting_chem and nesting_offline_chem (both .TRUE. 33 36 ! by default) … … 5756 5759 Reffic = exp( -Stokes**0.5_wp ) 5757 5760 ELSE 5758 stokes = vs1 * ustar / (grav * A_lu(luc) * 1. e-3)5759 Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1 e-3 ) )**25761 stokes = vs1 * ustar / (grav * A_lu(luc) * 1.E-3) 5762 Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1E-3 ) )**2 5760 5763 Reffic = exp( -Stokes**0.5_wp ) 5761 5764 END IF … … 5768 5771 ! 5769 5772 !-- sedimentation resistance: 5770 rs = 1.0_wp / ( epsilon0 * MAX( 1.0 e-5, ustar ) * ( Ebrown + Eimpac + Einterc ) * Reffic )5773 rs = 1.0_wp / ( epsilon0 * MAX( 1.0E-5_wp, ustar ) * ( Ebrown + Eimpac + Einterc ) * Reffic ) 5771 5774 5772 5775 !-- deposition velocity according to Seinfeld and Pandis (2006; eq 19.7):
Note: See TracChangeset
for help on using the changeset viewer.