anchor.diff:          else
               call interpolate_in_x(domain%lon,domain%u_geo%lon)
           endif
  -        call convert_longitudes(domain%u_geo%lon)
  +
  +        call standardize_coordinates(domain%u_geo)
   
           if (options%vlat_hi/="") then
               call io_read(options%init_conditions_file,options%vlat_hi,domain%v_geo%lat,1)
  
positive.code:         else
              call interpolate_in_y(domain%lon,domain%v_geo%lon)
          endif
          call convert_longitudes(domain%v_geo%lon)
  
          if (options%landvar/="") then
              call io_read(options%init_conditions_file,options%landvar,domain%landmask,1)
  
negative.code:         call io_read(options%init_conditions_file,options%hgt_hi,domain%terrain,1)
          call io_read(options%init_conditions_file,options%lat_hi,domain%lat,1)
          call io_read(options%init_conditions_file,options%lon_hi,domain%lon,1)
          call standardize_coordinates(domain)
  
          !  because u/vlat/lon_hi are optional, we calculated them by interplating from the mass lat/lon if necessary
          if (options%ulat_hi/="") then
  