 Timestamp:
 Nov 25, 2019 12:04:48 PM (4 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/chemistry_model_mod.f90
r4304 r4306 27 27 !  28 28 ! $Id$ 29 ! Corretion for r4304 commit 30 ! 31 ! 4304 20191125 10:43:03Z banzhafs 29 32 ! Precision cleanup in drydepo_aero_zhang_vd subroutine 30 33 ! … … 5683 5686 5684 5687 ! 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 ) 5686 5689 5687 5690 END FUNCTION sedimentation_velocity … … 5716 5719 REAL(wp), PARAMETER :: grav = 9.80665_wp !< acceleration of gravity (m/s2) 5717 5720 5718 REAL(wp), PARAMETER :: beta = 2.0_wp5719 5721 REAL(wp), PARAMETER :: epsilon0 = 3.0_wp 5720 5722 REAL(wp), PARAMETER :: kb = 1.38066E23_wp … … 5757 5759 ! and sticking efficiency R (1 = no rebound) 5758 5760 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 ) 5760 5762 Einterc = 0.0_wp 5761 5763 Reffic = 1.0_wp 5762 5764 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 ) 5764 5766 Einterc = 0.0_wp 5765 5767 Reffic = exp( Stokes**0.5_wp ) 5766 5768 ELSE 5767 5769 stokes = vs1 * ustar / ( grav * A_lu(luc) * 1.0E3_wp ) 5768 Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1.0E3_wp ) )**2 .0_wp5770 Einterc = 0.5_wp * ( partsize / (A_lu(luc) * 1.0E3_wp ) )**2 5769 5771 Reffic = exp( Stokes**0.5_wp ) 5770 5772 END IF … … 5774 5776 ! 5775 5777 ! determine impaction efficiency: 5776 Eimpac = ( stokes / ( alfa_lu(luc) + stokes ) )** beta5778 Eimpac = ( stokes / ( alfa_lu(luc) + stokes ) )**2 5777 5779 ! 5778 5780 ! sedimentation resistance:
Note: See TracChangeset
for help on using the changeset viewer.