anchor.diff:              if (urban_category==-1) urban_category = 13
               if (ice_category==-1)   ice_category = 15
               if (water_category==-1) water_category = 17
  +
           elseif (trim(LU_Categories)=="USGS") then
               if (urban_category==-1) urban_category = 1
               if (ice_category==-1)   ice_category = -1
               if (water_category==-1) water_category = 16
  +
           elseif (trim(LU_Categories)=="USGS-RUC") then
               if (urban_category==-1) urban_category = 1
               if (ice_category==-1)   ice_category = 24
               if (water_category==-1) water_category = 16
               ! also note, lakes_category = 28
               write(*,*) "WARNING: not handling lake category (28)"
  +
           elseif (trim(LU_Categories)=="MODI-RUC") then
               if (urban_category==-1) urban_category = 13
               if (ice_category==-1)   ice_category = 15
               if (water_category==-1) water_category = 17
               ! also note, lakes_category = 21
               write(*,*) "WARNING: not handling lake category (21)"
  +
           elseif (trim(LU_Categories)=="NLCD40") then
               if (urban_category==-1) urban_category = 13
               if (ice_category==-1)   ice_category = 15 ! and 22?
  
positive.code: 
  
  contains
  
      !> -------------------------------
      !! Initialize an options object
  
negative.code:                 this%lsm_options    = this%lsm_options[1]
                  this%bias_options   = this%bias_options[1]
              endif
              sync all
          end do
  
          call collect_physics_requests(this)
  
      end subroutine init
  
      !> -------------------------------
      !! Call all physics driver var_request routines
      !!
      !! var_request routines allow physics modules to requested
      !! which variables they need to have allocated, advected, and written in restart files
      !!
      !! -------------------------------
      subroutine collect_physics_requests(options)
          type(options_t) :: options
  
  