anchor.diff:          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)
  
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%init_conditions_file,options%hgt_hi,domain%terrain,1)
          call io_read2d(options%init_conditions_file,options%lat_hi,domain%lat,1)
          call io_read2d(options%init_conditions_file,options%lon_hi,domain%lon,1)
          call convert_longitudes(domain%lon)
          call io_read2d(options%init_conditions_file,options%ulat_hi,domain%u_geo%lat,1)
          call io_read2d(options%init_conditions_file,options%ulon_hi,domain%u_geo%lon,1)
          call convert_longitudes(domain%u_geo%lon)
          call io_read2d(options%init_conditions_file,options%vlat_hi,domain%v_geo%lat,1)
          call io_read2d(options%init_conditions_file,options%vlon_hi,domain%v_geo%lon,1)
          call convert_longitudes(domain%v_geo%lon)
          
          if (options%landvar/="") then
              call io_read2d(options%init_conditions_file,options%landvar,domain%landmask,1)
              where(domain%landmask==0) domain%landmask=kLC_WATER
  