Remapbil -> generic grid type (2 Grid Coordinates available: generic and lonlat)
Added by Anja Thielke over 6 years ago
Hi,
I've got an issue to remap my data.
When I use
cdo remapbil,../../../mygrid.txt tos.nc sst.nc
I get the error
cdo remapbil (Abort): Unsupported grid type: generic
But when checking the data (sinfov), I find that there are lon and lat values just as expected, but as a 2nd grid coordinate?
File format : netCDF4 classic ZIP -1 : Institut Source Ttype Levels Num Points Num Dtype : Parameter ID 1 : unknown ICOADS instant 1 1 2 1 F64 : -1 2 : unknown ICOADS instant 1 1 16020 2 F32z : -2 Grid coordinates : 1 : generic : points=2 2 : lonlat : points=16020 (180x89) lon : 0 to 358 by 2 degrees_east circular lat : 88 to -88 by -2 degrees_north Vertical coordinates : 1 : surface : levels=1 Time coordinate : 1942 steps RefTime = 1800-01-01 00:00:00 Units = days Calendar = standard YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss YYYY-MM-DD hh:mm:ss 1854-01-01 00:00:00 1854-02-01 00:00:00 1854-03-01 00:00:00 1854-04-01 00:00:00 1854-05-01 00:00:00 1854-06-01 00:00:00 1854-07-01 00:00:00 1854-08-01 00:00:00
Is it possible to chose the lonlat grid instead of the generic grid?
the lon lat variables are also available in the data
variables: float lat(lat) ; lat:units = "degrees_north" ; lat:long_name = "Latitude" ; lat:actual_range = 88.f, -88.f ; lat:standard_name = "latitude" ; lat:axis = "Y" ; lat:coordinate_defines = "center" ; float lon(lon) ; lon:units = "degrees_east" ; lon:long_name = "Longitude" ; lon:actual_range = 0.f, 358.f ; lon:standard_name = "longitude" ; lon:axis = "X" ; lon:coordinate_defines = "center" ; double time_bnds(time, nbnds) ; time_bnds:long_name = "Time Boundaries" ; double time(time) ; time:units = "days since 1800-1-1 00:00:00" ; time:long_name = "Time" ; time:delta_t = "0000-01-00 00:00:00" ; time:avg_period = "0000-01-00 00:00:00" ; time:prev_avg_period = "0000-00-07 00:00:00" ; time:standard_name = "time" ; time:axis = "T" ; time:actual_range = 19723., 78800. ; float sst(time, lat, lon) ; sst:long_name = "Monthly Means of Sea Surface Temperature" ; sst:units = "degC" ; sst:var_desc = "Sea Surface Temperature" ; sst:level_desc = "Surface" ; sst:statistic = "Mean" ; sst:missing_value = -9.96921e+36f ; sst:actual_range = -1.8f, 33.95f ; sst:valid_range = -5.f, 40.f ; sst:dataset = "NOAA Extended Reconstructed SST V4" ; sst:parent_stat = "Individual Values" ; sst:coordinates = "lon lat" ;
Cheers,
Anja
Replies (2)
RE: Remapbil -> generic grid type (2 Grid Coordinates available: generic and lonlat) - Added by Ralf Mueller over 6 years ago
I recommend
cdo remapbil,../../../mygrid.txt -selname,sst tos.nc sst.nc
backgroud: CDO recognizes the variable time_bnds
as a data variable living on 2-point-grid, which cannot be interpolated. Thats why you need to select the variabe(s) to be interpolated instead of using the whole file.
hth
ralf
RE: Remapbil -> generic grid type (2 Grid Coordinates available: generic and lonlat) - Added by Anja Thielke over 6 years ago
Hi Ralf,
thank you for your reply. This is exactly what I needed.
Greetings
Anja