DustDryDep

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 CCSM, 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

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 :: 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