initGridcells

INTERFACE:

   subroutine initGridcells (vegxy, wtxy)

DESCRIPTION:

   Initialize sub-grid mapping and allocates space for derived type
   hierarchy.  For each land gridcell determine landunit, column and
   pft properties.  Note that ngcells, nlunits, ncols and npfts are
   per-processor totals here and are currently not used for anything other
   than placeholders.  Determine if there are any vegetated landunits and
   if so---the weight of the vegetated landunit relative to the gridcell
   The first landunit contains all the vegetated patches (if any) For now,
   the vegetated patches will all be gathered on a single landunit, with
   each vegetated type having its own column on that landunit.  The special
   patches (urban, lake, wetland, glacier) each get their own landunit
   having a single column and one non-vegetated pfts

USES:

     use decompMod    , only : get_proc_bounds, get_proc_global, &
                               get_gcell_info, get_gcell_xyind, &
                               get_sn_index
     use shr_const_mod, only : SHR_CONST_PI

ARGUMENTS:

     implicit none
     integer , intent(in) :: vegxy(lsmlon,lsmlat,maxpatch) ! PFT type
     real(r8), intent(in) :: wtxy(lsmlon,lsmlat,maxpatch)  ! subgrid patch
                                                           ! weights

REVISION HISTORY:

   Created by Peter Thornton and Mariana Vertenstein