anchor.diff:          ! these variables are required for any boundary/forcing file type
           call io_read2d(options%boundary_files(1),options%latvar,boundary%lat)
           call io_read2d(options%boundary_files(1),options%lonvar,boundary%lon)
  +        call convert_longitudes(boundary%lon)
           call io_read2d(options%boundary_files(1),options%ulat,boundary%u_geo%lat)
           call io_read2d(options%boundary_files(1),options%ulon,boundary%u_geo%lon)
  +        call convert_longitudes(boundary%u_geo%lon)
           call io_read2d(options%boundary_files(1),options%vlat,boundary%v_geo%lat)
           call io_read2d(options%boundary_files(1),options%vlon,boundary%v_geo%lon)
  +        call convert_longitudes(boundary%v_geo%lon)
           call io_read2d(options%boundary_files(1),options%hgtvar,boundary%terrain)
           
           ! read in the vertical coordinate
  
positive.code: !
  ! ------------------------------------------------------------------------------------------
  
  !   Allow running over a sub-domain, by removing the outer N grid cells from all sides of the domain (lat,lon,terrain)
      subroutine remove_edges(domain,edgesize)
          implicit none
  
negative.code:         call io_read2d(options%ext_wind_files(1),options%hgt_hi,bc%ext_winds%terrain,1)
          call io_read2d(options%ext_wind_files(1),options%latvar,bc%ext_winds%lat)
          call io_read2d(options%ext_wind_files(1),options%lonvar,bc%ext_winds%lon)
          call convert_longitudes(bc%ext_winds%lon)
          call io_read2d(options%ext_wind_files(1),options%ulat_hi,bc%ext_winds%u_geo%lat)
          call io_read2d(options%ext_wind_files(1),options%ulon_hi,bc%ext_winds%u_geo%lon)
          call convert_longitudes(bc%ext_winds%u_geo%lon)
          call io_read2d(options%ext_wind_files(1),options%vlat_hi,bc%ext_winds%v_geo%lat)
          call io_read2d(options%ext_wind_files(1),options%vlon_hi,bc%ext_winds%v_geo%lon)
          call convert_longitudes(bc%ext_winds%v_geo%lon)
          
          write(*,*) "Setting up ext wind geoLUTs"
          call geo_LUT(bc%next_domain%u_geo, bc%ext_winds%u_geo)
          call geo_LUT(bc%next_domain%v_geo, bc%ext_winds%v_geo)
  