anchor.diff:          call io_read(options%ext_wind_files(1),options%hgt_hi,bc%ext_winds%terrain,1)
           call io_read(options%ext_wind_files(1),options%latvar,bc%ext_winds%lat)
           call io_read(options%ext_wind_files(1),options%lonvar,bc%ext_winds%lon)
  -        call convert_longitudes(bc%ext_winds%lon)
  +        call standardize_coordinates(bc%ext_winds)
   
           ! ulat_hi and vlat_hi are optional, if they are not supplied interpolate the mass lat/lon grid
           if (options%ulat_hi/="") then
  
positive.code:         else
              call interpolate_in_x(bc%ext_winds%lon,bc%ext_winds%u_geo%lon)
          endif
          call convert_longitudes(bc%ext_winds%u_geo%lon)
  
          if (options%vlat_hi/="") then
              call io_read(options%ext_wind_files(1),options%vlat_hi,bc%ext_winds%v_geo%lat,1)
  
negative.code:         ! these variables are required for any boundary/forcing file type
          call io_read(options%boundary_files(1),options%latvar,boundary%lat)
          call io_read(options%boundary_files(1),options%lonvar,boundary%lon)
          call standardize_coordinates(boundary)
          call io_read(options%boundary_files(1),options%ulat,boundary%u_geo%lat)
          call io_read(options%boundary_files(1),options%ulon,boundary%u_geo%lon)
          call standardize_coordinates(boundary%u_geo)
          call io_read(options%boundary_files(1),options%vlat,boundary%v_geo%lat)
          call io_read(options%boundary_files(1),options%vlon,boundary%v_geo%lon)
          call standardize_coordinates(boundary%v_geo)
          call io_read(options%boundary_files(1),options%hgtvar,boundary%terrain)
  
          ! read in the vertical coordinate
  