CDO

I am using RCM netCDF data on a rotated pole grid. When ever I apply any CDO function to the data two grid related variables disappear, making plotting (in IDL or NCL) not possible.

Has anyone else had similar problems or have any solution ideas?

see below the file header before and after applying a CDO tool:

before . . . .

netcdf cahpaa.pm.6190.05216.rr8 {
dimensions:
time = 360 ;
bounds2 = 2 ;
grid_latitude = 192 ;
grid_longitude = 192 ;
variables:
double time(time) ;
time:units = "days since 1860-12-1" ;
time:calendar = "360_day" ;
time:standard_name = "time" ;
time:axis = "T" ;
time:bounds = "time_bounds" ;
double time_bounds(time, bounds2) ;
float grid_latitude(grid_latitude) ;
grid_latitude:units = "degrees" ;
grid_latitude:standard_name = "grid_latitude" ;
grid_latitude:axis = "Y" ;
grid_latitude:bounds = "grid_latitude_bounds" ;
float grid_latitude_bounds(grid_latitude, bounds2) ;
float grid_longitude(grid_longitude) ;
grid_longitude:units = "degrees" ;
grid_longitude:standard_name = "grid_longitude" ;
grid_longitude:axis = "X" ;
grid_longitude:bounds = "grid_longitude_bounds" ;
float grid_longitude_bounds(grid_longitude, bounds2) ;
float latitude(grid_latitude, grid_longitude) ;
latitude:units = "degrees_north" ;
latitude:standard_name = "latitude" ;
float longitude(grid_latitude, grid_longitude) ;
longitude:units = "degrees_east" ;
longitude:standard_name = "longitude" ;
char rotated_latitude_longitude ;
rotated_latitude_longitude:grid_north_pole_latitude = 75.f ;
rotated_latitude_longitude:grid_north_pole_longitude = 289.f ;
rotated_latitude_longitude:grid_mapping_name = "rotated_latitude_longitude" ;
float precipitation_flux(time, grid_latitude, grid_longitude) ;
precipitation_flux:_FillValue = -1.073742e+09f ;
precipitation_flux:long_name = "TOTAL PRECIPITATION RATE KG/M2/S" ;
precipitation_flux:standard_name = "precipitation_flux" ;
precipitation_flux:units = "kg m-2 s-1" ;
precipitation_flux:coordinates = "latitude longitude" ;
precipitation_flux:grid_mapping = "rotated_latitude_longitude" ;
precipitation_flux:cell_methods = "time: mean" ;

// global attributes:
:Conventions = "CF-1.5" ;
:runid = "cahpa" ;
:stash_code = 5216 ;
:submodel = 1 ;
:lbproc = 128 ;
:history = "2012-09-27 11:08:44 UTC Aggregated by cf-python v0.9.3.3" ;

after . . .

netcdf cahpaa.pm.6190.05216.rr8.mmday {
dimensions:
x = 192 ;
y = 192 ;
time = UNLIMITED ; // (360 currently)
nb2 = 2 ;
variables:
float longitude(y, x) ;
longitude:standard_name = "longitude" ;
longitude:long_name = "longitude" ;
longitude:units = "degrees_east" ;
longitude:_CoordinateAxisType = "Lon" ;
float latitude(y, x) ;
latitude:standard_name = "latitude" ;
latitude:long_name = "latitude" ;
latitude:units = "degrees_north" ;
latitude:_CoordinateAxisType = "Lat" ;
double time(time) ;
time:standard_name = "time" ;
time:bounds = "time_bnds" ;
time:units = "days since 1860-12-01 00:00:00" ;
time:calendar = "360_day" ;
double time_bnds(time, nb2) ;
time_bnds:units = "days since 1860-12-01 00:00:00" ;
time_bnds:calendar = "360_day" ;
float precipitation_flux(time, y, x) ;
precipitation_flux:standard_name = "precipitation_flux" ;
precipitation_flux:long_name = "TOTAL PRECIPITATION RATE KG/M2/S" ;
precipitation_flux:units = "kg m-2 s-1" ;
precipitation_flux:coordinates = "longitude latitude" ;
precipitation_flux:_FillValue = -1.073742e+09f ;
precipitation_flux:cell_methods = "time: mean" ;

// global attributes:
:CDI = "Climate Data Interface version 1.5.4 (http://code.zmaw.de/projects/cdi)" ;
:Conventions = "CF-1.5" ;
:history = "Mon Oct 01 15:03:37 2012: cdo mulc,86400 cahpaa.pm.6190.05216.rr8.nc cahpaa.pm.6190.05216.rr8.mmday.nc\n",
"2012-09-27 11:08:44 UTC Aggregated by cf-python v0.9.3.3" ;
:runid = "cahpa" ;
:stash_code = 5216 ;
:submodel = 1 ;
:lbproc = 128 ;
:CDO = "Climate Data Operators version 1.5.4 (http://code.zmaw.de/projects/cdo)" ;
}

Thanks very much!

Grace