Establishment

INTERFACE:

   subroutine Establishment(lbp, ubp, lbg, ubg)

DESCRIPTION:

   Calculates establishment of new pfts
   Called once per year

USES:

     use clmtype
     use clm_varpar   , only : numpft
     use clm_varcon   , only : istsoil
     use pftvarcon    , only : noveg
     use shr_const_mod, only : SHR_CONST_CDAY, SHR_CONST_PI, SHR_CONST_TKFRZ

ARGUMENTS:

     implicit none
     integer , intent(in) :: lbp, ubp         ! pft bounds
     integer , intent(in) :: lbg, ubg         ! gridcell bounds

CALLED FROM:

   subroutine lpj in module DGVMMod

REVISION HISTORY:

   Author: Sam Levis (adapted from Stephen Sitch's LPJ subr. establishment)
   3/4/02, Peter Thornton: Migrated to new data structures.

LOCAL VARIABLES:

   local pointers to implicit in arguments
     integer , pointer :: plandunit(:)     ! landunit of corresponding pft
     integer , pointer :: pgridcell(:)     ! gridcell of corresponding pft
     integer , pointer :: ltype(:)         ! landunit type for corresponding p
     real(r8), pointer :: wtgcell(:)       ! pft weight relative to grid cell
     real(r8), pointer :: tmomin20(:)      ! 20-yr running mean of tmomin
     real(r8), pointer :: agdd20(:)        ! 20-yr running mean of agdd
     real(r8), pointer :: agddtw(:)        ! accumulated growing degree days a
     real(r8), pointer :: prec365(:)       ! 365-day running mean of tot. prec
     real(r8), pointer :: sla(:)           ! ecophys const - sp. leaf area [m2
     logical , pointer :: tree(:)          ! ecophys const - true=> tree is pr
     real(r8), pointer :: crownarea_max(:) ! ecophys const - tree maximum crow
     real(r8), pointer :: lm_sapl(:)       ! ecophys const - leaf mass of sapl
     real(r8), pointer :: sm_sapl(:)       ! ecophys const - stem mass of sapl
     real(r8), pointer :: hm_sapl(:)       ! ecophys const - heartwood mass of
     real(r8), pointer :: rm_sapl(:)       ! ecophys const - root mass of sapi
     real(r8), pointer :: reinickerp(:)    ! ecophys const - parameter in allo
     real(r8), pointer :: wooddens(:)      ! ecophys const - wood density (gC/
     real(r8), pointer :: latosa(:)        ! ecophys const - ratio of leaf are
     real(r8), pointer :: allom1(:)        ! ecophys const - parameter in allo
     real(r8), pointer :: allom2(:)        ! ecophys const - parameter in allo
     real(r8), pointer :: allom3(:)        ! ecophys const - parameter in allo
     real(r8), pointer :: tcmin(:)         ! ecophys const - minimum coldest m
     real(r8), pointer :: tcmax(:)         ! ecophys const - maximum coldest m
     real(r8), pointer :: gddmin(:)        ! ecophys const - minimum growing d
   local pointers to implicit in/out arguments
     integer , pointer :: ivt(:)           ! vegetation type for this pft
     logical , pointer :: present(:)       ! true=> PFT present in patch
     real(r8), pointer :: nind(:)          ! number of individuals (#/m**2)
     real(r8), pointer :: lm_ind(:)        ! individual leaf mass
     real(r8), pointer :: sm_ind(:)        ! individual sapwood mass
     real(r8), pointer :: hm_ind(:)        ! individual heartwood mass
     real(r8), pointer :: rm_ind(:)        ! individual root mass
     real(r8), pointer :: litterag(:)      ! above ground litter
     real(r8), pointer :: litterbg(:)      ! below ground litter
   local pointers to implicit out arguments
     real(r8), pointer :: fpcgrid(:)       ! foliar projective cover on gridce
     real(r8), pointer :: htop(:)          ! canopy top (m)
     real(r8), pointer :: lai_ind(:)       ! LAI per individual
     real(r8), pointer :: crownarea(:)     ! area that each individual tree ta