anchor.diff:          real, allocatable, dimension(:,:,:) :: du_dt,dv_dt,dp_dt,dth_dt,dqv_dt,dqc_dt
           ! sh, lh, and pblh fields are only 2d. 
           ! These are only used with LSM option 1 and are derived from forcing file
  -        real, allocatable, dimension(:,:) :: dsh_dt,dlh_dt,dpblh_dt
  +        real, allocatable, dimension(:,:)   :: dsh_dt,dlh_dt,dpblh_dt
           ! change in shortwave and longwave at surface if read from forcing
  -        real, allocatable, dimension(:,:) :: dsw_dt, dlw_dt
  +        real, allocatable, dimension(:,:)   :: dsw_dt, dlw_dt
           ! change in sst if read from forcing file
  -        real, allocatable, dimension(:,:) :: dsst_dt
  -        ! store the low resolution versionf of terrain and atmospheric elevations
  -        real,allocatable,dimension(:,:)::lowres_terrain
  -        real,allocatable,dimension(:,:,:)::lowres_z
  +        real, allocatable, dimension(:,:)   :: dsst_dt
  +        
  +        ! store the low resolution versiof of terrain and atmospheric elevations
  +        real,allocatable,dimension(:,:)     :: lowres_terrain
  +        real,allocatable,dimension(:,:,:)   :: lowres_z
           ! store the full high-res 3D grid for the next time step to compute dXdt fields
           ! includes high res versions of low res terrain and z
           type(domain_type)::next_domain
  
positive.code: !! pbl_height    = Height of the planetary boundary layer   [m]
  !!
  !!      ---- Radiation variables ----
  !! cloudfrac        = Cloud fraction                            [0-1]
  !! swdown       = Shortwave down at land surface            [W/m^2]
  !! lwdown       = Longwave down at land surface             [W/m^2]
  !! lwup         = Lonwave up from the land surface          [W/m^2]
  !!
  !! ---- Land Surface variables ---- 
  !!   3D fields ---- NX x NZ x NY
  !! soil_t       = 3D Soil temperature                       [K]
  !! soil_vwc     = 3D Soil volumetric water content          [m^3/m^3]
  !! 
  
negative.code:     ! store Advection options
      ! ------------------------------------------------
      type adv_options_type
          logical :: boundary_buffer          ! buffer to smooth a bit to cut down on oscillations at the border if FCT is not used
          logical :: flux_corrected_transport ! use Flux Corrected Transport (FCT) to maintain stability and prevent any wild oscllations
          integer :: mpdata_order             ! accuracy order for MP_DATA advection scheme. 
      end type adv_options_type
  
  
  