Project

General

Profile

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

    (1-2/2)