CDO

Hi,
I am trying to regrid CMIP5 thetao data to a regular 2.5x2.5 grid. However, there is something wrong with the Russia model(INMCM4.0).

cdo intlevel,10,20,30,40,50,60,70,80,90,100,110,120,130,140,150,160,170,180,190,200 -selname,thetao thetao_Omon_inmcm4_historical_r1i1p1_197001-197412.nc interim.nc
cdo remapbil,mygrid interim.nc thetao_197001-197412.nc

cdo intlevel: Started child process "selname,thetao thetao_Omon_inmcm4_historical_r1i1p1_197001-197412.nc (pipe1.1)".
Warning (cdfScanVarAttributes) : NetCDF: Variable not found - areacello
cdo(2) selname: Processed 293760000 values from 3 variables over 60 timesteps ( 11.80s )
cdo intlevel: Processed 293760000 values from 1 variable over 60 timesteps ( 11.80s )
cdo remapbil: SCRIP bilinear remapping from curvilinear (360x340) to lonlat (144x49) grid, with source mask (0)
cdo remapbil: Processed 146880000 values from 1 variable over 60 timesteps ( 2.54s )

There seem no errors, but I find that the datas in the output(thetao_197001-197412.nc) are all missing values. I also checked the original input, and it seems normal (both normal values and missing values).
I don't know what is wrong.
Could anyone help me with this.

A small size data "thetao_Omon_inmcm4_historical_r1i1p1_197001-197002.lite.nc" has been uploaded.

And a header of thetao_Omon_inmcm4_historical_r1i1p1_197001-197412.nc:
netcdf thetao_Omon_inmcm4_historical_r1i1p1_197001-197412 {
dimensions:
time = UNLIMITED ; // (60 currently)
lev = 40 ;
rlat = 340 ;
rlon = 360 ;
bnds = 2 ;
vertices = 4 ;
variables:
double time(time) ;
time:bounds = "time_bnds" ;
time:units = "days since 1850-1-1" ;
time:calendar = "365_day" ;
time:axis = "T" ;
time:long_name = "time" ;
time:standard_name = "time" ;
double time_bnds(time, bnds) ;
double lev(lev) ;
lev:bounds = "lev_bnds" ;
lev:axis = "Z" ;
lev:positive = "up" ;
lev:long_name = "ocean sigma coordinate" ;
lev:standard_name = "ocean_sigma_coordinate" ;
lev:formula = "z(n,k,j,i) = eta(n,j,i) + sigma(k)*(depth(j,i)+eta(n,j,i))" ;
lev:formula_terms = "sigma: lev eta: eta depth: depth" ;
double lev_bnds(lev, bnds) ;
lev_bnds:formula = "z(n,k,j,i) = eta(n,j,i) + sigma(k)*(depth(j,i)+eta(n,j,i))" ;
lev_bnds:standard_name = "ocean_sigma_coordinate" ;
lev_bnds:units = "" ;
lev_bnds:formula_terms = "sigma: lev_bnds eta: eta depth: depth" ;
float eta(time, rlat, rlon) ;
eta:long_name = "Sea Surface Height" ;
eta:units = "m" ;
eta:cell_methods = "time: mean" ;
float depth(rlat, rlon) ;
depth:long_name = "Sea Floor Depth" ;
depth:comment = "Ocean bathymetry." ;
depth:units = "m" ;
double rlat(rlat) ;
rlat:units = "degrees" ;
rlat:axis = "Y" ;
rlat:long_name = "latitude in rotated pole grid" ;
rlat:standard_name = "grid_latitude" ;
double rlon(rlon) ;
rlon:units = "degrees" ;
rlon:axis = "X" ;
rlon:long_name = "longitude in rotated pole grid" ;
rlon:standard_name = "grid_longitude" ;
int rotated_latitude_longitude ;
rotated_latitude_longitude:grid_mapping_name = "rotated_latitude_longitude" ;
rotated_latitude_longitude:grid_north_pole_latitude = 70. ;
rotated_latitude_longitude:grid_north_pole_longitude = 100. ;
rotated_latitude_longitude:north_pole_grid_longitude = -80. ;
float lat(rlat, rlon) ;
lat:standard_name = "latitude" ;
lat:long_name = "latitude coordinate" ;
lat:units = "degrees_north" ;
lat:bounds = "lat_vertices" ;
float lon(rlat, rlon) ;
lon:standard_name = "longitude" ;
lon:long_name = "longitude coordinate" ;
lon:units = "degrees_east" ;
lon:bounds = "lon_vertices" ;
float lat_vertices(rlat, rlon, vertices) ;
lat_vertices:units = "degrees_north" ;
float lon_vertices(rlat, rlon, vertices) ;
lon_vertices:units = "degrees_east" ;
float thetao(time, lev, rlat, rlon) ;
thetao:standard_name = "sea_water_potential_temperature" ;
thetao:long_name = "Sea Water Potential Temperature" ;
thetao:units = "K" ;
thetao:original_name = "thetao" ;
thetao:cell_methods = "time: mean (interval: 1 month)" ;
thetao:cell_measures = "area: areacello volume: volcello" ;
thetao:history = "2010-09-01T15:38:09Z altered by CMOR: Reordered dimensions, original order: time rlat rlon lev." ;
thetao:missing_value = 1.e+20f ;
thetao:_FillValue = 1.e+20f ;
thetao:associated_files = "baseURL: http://cmip-pcmdi.llnl.gov/CMIP5/dataLocation gridspecFile: gridspec_fx_inmcm4_historical_r0i0p0.nc areacello: areacello_fx_inmcm4_historical_r0i0p0.nc volcello: volcello_fx_inmcm4_historical_r0i0p0.nc" ;
thetao:grid_mapping = "rotated_latitude_longitude" ;
thetao:coordinates = "lat lon" ;

// global attributes:
:institution = "INM (Institute for Numerical Mathematics, Moscow, Russia)" ;
:institute_id = "INM" ;
:experiment_id = "historical" ;
:source = "inmcm4 (2009)" ;
:model_id = "inmcm4" ;
:forcing = "N/A" ;
:parent_experiment_id = "piControl" ;
:branch_time = 0. ;
:contact = "Evgeny Volodin, volodin@inm.ras.ru,INM RAS, Gubkina 8, Moscow, 119333 Russia,+7-495-9383904" ;
:history = "Output from /data5/volodin/20C 2010-09-01T15:38:10Z CMOR rewrote data to comply with CF standards and CMIP5 requirements." ;
:comment = "no comments" ;
:references = "Volodin, Diansky, Gusev 2010. Climate model INMCM4.0. Izvestia RAS. Atmospheric and oceanic physics, V.46, N4, in print." ;
:initialization_method = 1 ;
:physics_version = 1 ;
:tracking_id = "a630b82d-8d22-485d-b549-dfa503078b72" ;
:product = "output" ;
:experiment = "Historical" ;
:frequency = "mon" ;
:creation_date = "2010-09-01T15:38:10Z" ;
:Conventions = "CF-1.4" ;
:project_id = "CMIP5" ;
:table_id = "Table Omon (12 May 2010) f2afe576fb73a3a11aaa3cc8f2e62cf3" ;
:title = "inmcm4 model output prepared for CMIP5 Historical" ;
:parent_experiment = "pre-industrial control" ;
:modeling_realm = "ocean" ;
:realization = 1 ;
:cmor_version = "2.0.0" ;