Project

General

Profile

netcdf CF-1.0 coordinate system for projected CRS not supported

Added by Cary Lynch about 9 years ago

I am having a problem using CDO commands on netcdf files that have coordinate information in the CRS CF-1.0 format.
The file attributes are as follows:
netcdf pr_bcc-csm1-1_r1i1p1_historical_1950_1969_daily {
dimensions:
lat = 444 ;
lon = 922 ;
time = 7305 ;
crs = 1 ;
variables:
double lat(lat) ;
lat:long_name = "latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
lat:axis = "Y" ;
lat:description = "Latitude of the center of the grid cell" ;
double lon(lon) ;
lon:long_name = "longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
lon:axis = "X" ;
lon:description = "Longitude of the center of the grid cell" ;
float time(time) ;
time:units = "days since 1900-01-01 00:00:00" ;
time:calendar = "gregorian" ;
time:description = "days since 1900-01-01" ;
float precipitation(time, lat, lon) ;
precipitation:_FillValue = -9999.f ;
precipitation:long_name = "Precipitation" ;
precipitation:units = "mm" ;
precipitation:grid_mapping = "crs" ;
precipitation:standard_name = "precipitation" ;
precipitation:cell_methods = "time: sum(interval: 24 hours)" ;
precipitation:coordinates = "time lon lat" ;
int crs(crs) ;
crs:grid_mapping_name = "latitude_longitude" ;
crs:longitude_of_prime_meridian = 0. ;
crs:semi_major_axis = 6378137. ;
crs:inverse_flattening = 298.257223563 ;

I am unsure of the proper sequence commands to use to re-project this data. I want to keep the nlat and nlon the same, can anyone help with the remapbil command sequence to make this data compatible with CDOs?