DustEmission

INTERFACE:

   subroutine DustEmission (lbp, ubp, num_nolakep, filter_nolakep)

DESCRIPTION:

Dust mobilization. This code simulates dust mobilization due to wind from the surface into the lowest atmospheric layer On output flx_mss_vrt_dst(ndst) is the surface dust emission (kg/m4#4/s) [+ = to atm]. Source: C. Zender's dust model

USES:

    use shr_const_mod, only : SHR_CONST_RHOFW

ARGUMENTS:

     implicit none
     integer, intent(in) :: lbp, ubp                    ! pft bounds
     integer, intent(in) :: num_nolakep                 ! number of column non
     integer, intent(in) :: filter_nolakep(num_nolakep) ! pft filter for non-l

LOCAL VARIABLES:

   local pointers to implicit in arguments
     integer , pointer :: pcolumn(:)         ! pft's column index
     integer , pointer :: plandunit(:)       ! pft's landunit index
     integer , pointer :: pgridcell(:)       ! pft's gridcell index
     integer , pointer :: ityplun(:)         ! landunit type
     real(r8), pointer :: tlai(:)            ! one-sided leaf area index, no b
     real(r8), pointer :: tsai(:)            ! one-sided stem area index, no b
     real(r8), pointer :: frac_sno(:)        ! fraction of ground covered by s
     real(r8), pointer :: gwc_thr(:)         ! threshold gravimetric soil mois
     real(r8), pointer :: forc_rho(:)        ! density (kg/m**3)
     real(r8), pointer :: fv(:)              ! friction velocity (m/s) (for du
     real(r8), pointer :: u10(:)             ! 10-m wind (m/s) (created for du
     real(r8), pointer :: mbl_bsn_fct(:)     ! basin factor
     real(r8), pointer :: mss_frc_cly_vld(:) ! [frc] Mass fraction clay limite
     real(r8), pointer :: h2osoi_vol(:,:)    ! volumetric soil water (0<=h2oso
     real(r8), pointer :: h2osoi_liq(:,:)    ! liquid soil water (kg/m2)
     real(r8), pointer :: h2osoi_ice(:,:)    ! frozen soil water (kg/m2)
     real(r8), pointer :: watsat(:,:)        ! saturated volumetric soil water
 
   local pointers to implicit out arguments
     real(r8), pointer :: flx_mss_vrt_dst(:,:)   ! surface dust emission (kg/m
     real(r8), pointer :: flx_mss_vrt_dst_tot(:) ! total dust flux into atmosp

REVISION HISTORY:

   Created by Sam Levis
   Migrated to new data structures by Peter Thornton and Mariana Vertenstein
   !Created by Peter Thornton and Mariana Vertenstein