CDO

Dear All, I have a netcdf file which has originally stretched rotated grid. I converted it into rotated grid with eliminate stretching coefficient by rotated grid transform function of MATLAB. The original file structure is:

@netcdf ARPEGE4.5_CTL_CNRM-CM3_DM_50km_1950-2000_pr {
dimensions:
rgrid = 17859 ;
time = UNLIMITED ; // (18628 currently)
variables:
float longitude(rgrid) ;
longitude:standard_name = "longitude" ;
longitude:long_name = "longitude" ;
longitude:units = "degrees_east" ;
float latitude(rgrid) ;
latitude:standard_name = "latitude" ;
latitude:long_name = "latitude" ;
latitude:units = "degrees_north" ;
int rgrid(rgrid) ;
rgrid:compress = "latitude longitude" ;
double time(time) ;
time:units = "days since 1950-01-01 00:00" ;
time:standard_name = "time" ;
time:long_name = "time" ;
time:calendar = "standard" ;
float pr(time, rgrid) ;
pr:standard_name = "precipitation_flux" ;
pr:long_name = "Precipitation" ;
pr:units = "kg m-2 s-1" ;
pr:coordinates = "lon lat" ;
pr:grid_mapping = "rotated_stretched_latitude_longitude" ;
pr:rotated_north_pole_latitude = "40." ;
pr:rotated_north_pole_longitude = "12." ;
pr:stretching_coefficient = "2.5" ;
pr:cell_method = "time: mean" ;

// global attributes:
:institution = "Meteo-France/CNRM" ;
:Conventions = "CF-1.0" ;
:conventionsURL = "http://www.cgd.ucar.edu/cms/eaton/cf-metadata/index.html" ;
:source = "ARPEGE-Climat v4.6" ;
:project_id = "CECILIA" ;
:experiment_id = "EB9" ;
:realization = "0a" ;
:experiment_ref = "CNRM RCM experiments EB9 EE1" ;
:creation_date = "2008-07-08 17:52:45" ;
:filename = "ARPEGE4.5_CTL_CNRM-CM3_DM_50km_195001_pr.nc" ;
:history = "Fri Jul 11 10:12:51 2008: ncatted -O -a history,global,a,c,(ncrcat 195001 ... 200012) new/ARPEGE4.5_CTL_CNRM-CM3_DM_50km_1950-2000_pr.nc\n",
"(ncrcat 195001 ... 200012)" ;
}
@
and the new one is:

netcdf ARPEGE_1961-1990_pr {
dimensions:
rgrid = 17859 ;
time = 10957 ;
lon = 17859 ;
lat = 17859 ;
variables:
float lon(rgrid) ;
lon:standard_name = "longitude" ;
lon:long_name = "longitude" ;
lon:units = "degrees_east" ;
float lat(rgrid) ;
lat:standard_name = "latitude" ;
lat:long_name = "latitude" ;
lat:units = "degrees_north" ;
int rgrid(rgrid) ;
rgrid:compress = "latitude longitude" ;
double time(time) ;
time:units = "days since 1950-01-01 00:00" ;
time:standard_name = "time" ;
time:long_name = "time" ;
time:calendar = "standard" ;
float pr(time, rgrid) ;
pr:standard_name = "precipitation_flux" ;
pr:long_name = "Precipitation" ;
pr:units = "kg m-2 s-1" ;
pr:coordinates = "lon lat" ;
pr:grid_mapping = "rotated_stretched_latitude_longitude" ;
pr:rotated_north_pole_latitude = "40." ;
pr:rotated_north_pole_longitude = "12." ;
pr:cell_method = "time: mean" ;
}


I would like to remap this new netcdf file by remapcon operator, but it does not have grid cell corner. Is there a way to determine grid cell corner values?

Thank you very much for your help in advance!