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: 
  
  
      !> -------------------------------
      !! Reads the name of the options file from the commandline
      !!
      !! If no options file was specified, defaults to icar_options.nml
      !!
      !! -------------------------------
      function get_options_file() result(options_file)
          implicit none
          character(len=MAXFILELENGTH) ::options_file
          integer :: error
          logical :: file_exists
  
          if (command_argument_count()>0) then
              call get_command_argument(1,options_file, status=error)
              if (error>0) then
                  options_file="icar_options.nml"
              elseif (error==-1) then
                  write(*,*) "Options filename = ", trim(options_file), " ...<cutoff>"
                  write(*,*) "Maximum filename length = ", MAXFILELENGTH
                  stop "ERROR: options filename too long"
              endif
          else
              options_file="icar_options.nml"
          endif
          INQUIRE(file=trim(options_file), exist=file_exists)
          if (.not.file_exists) then
              write(*,*) "Using options file = ", trim(options_file)
              stop "Options file does not exist. "
  
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
  
  