Remapbil -> generic grid type (2 Grid Coordinates available: generic and lonlat)
Added by Anja Thielke almost 8 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 almost 8 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 almost 8 years ago
Hi Ralf,
thank you for your reply. This is exactly what I needed.
Greetings
Anja