INTERFACE:
subroutine DustDryDep (lbp, ubp)DESCRIPTION:
Determine Turbulent dry deposition for dust. Calculate the turbulent component of dust dry deposition, (the turbulent deposition velocity through the lowest atmospheric layer. CAM will calculate the settling velocity through the whole atmospheric column. The two calculations will determine the dust dry deposition flux to the surface. Note: Same process should occur over oceans. For the coupled CESM, we may find it more efficient to let CAM calculate the turbulent dep velocity over all surfaces. This would require passing the aerodynamic resistance, ram(1), and the friction velocity, fv, from the land to the atmosphere component. In that case, dustini need not calculate particle diamter (dmt_vwr) and particle density (dns_aer). Source: C. Zender's dry deposition code
!USES
use shr_const_mod, only : SHR_CONST_PI, SHR_CONST_RDAIR, SHR_CONST_BOLTZ use clm_atmlnd , only : clm_a2l
ARGUMENTS:
implicit none
integer, intent(in) :: lbp, ubp ! pft bounds
!LOCAL VARIABLES
local pointers to implicit in arguments
integer , pointer :: pgridcell(:) ! pft's gridcell index
real(r8), pointer :: pwtgcell(:) ! weight of pft relative to corresponding gridcell
real(r8), pointer :: forc_t(:) ! atm temperature (K)
real(r8), pointer :: forc_pbot(:) ! atm pressure (Pa)
real(r8), pointer :: forc_rho(:) ! atm density (kg/m**3)
real(r8), pointer :: fv(:) ! friction velocity (m/s)
real(r8), pointer :: ram1(:) ! aerodynamical resistance (s/m)
real(r8), pointer :: vlc_trb(:,:) ! Turbulent deposn velocity (m/s)
real(r8), pointer :: vlc_trb_1(:) ! Turbulent deposition velocity 1
real(r8), pointer :: vlc_trb_2(:) ! Turbulent deposition velocity 2
real(r8), pointer :: vlc_trb_3(:) ! Turbulent deposition velocity 3
real(r8), pointer :: vlc_trb_4(:) ! Turbulent deposition velocity 4
!REVISION HISTORY
Created by Sam Levis
!LOCAL VARIABLES