CDO

Original from pmsoares@fc.ul.pt:

I’ve been trying to use the cdo to remapbil some RCM’s results, but in two model cases the cdo gives an error. For example, for the PROMES model output, I do the following:

cat > mygrid << EOF

gridtype = lonlat

xsize = 87         # nx

ysize = 46         # ny

xfirst = -11.6     # initial lon value

xinc = 0.2         # dx

yfirst = 35.2      # initial lat value

yinc = 0.2         # dy

EOF

cdo remapbil,mygrid UCLM-PROMES_CTL_ERA40_DM_25km_1961-1970_pr.nc UCLM-PROMES_CTL_ERA40_DM_25km_1961-1970_pr_i.nc

and i get the following error:

Error (cdoDefineGrid) : Read grid from mygrid failed!

When I look the ncdump –v of the files I see:

[11:13:39 wrf@tornado Dados] $ ncdump -h UCLM-PROMES_CTL_ERA40_DM_25km_1961-1970_pr.nc

netcdf UCLM-PROMES_CTL_ERA40_DM_25km_1961-1970_pr {

dimensions:

        time = UNLIMITED ; // (3652 currently)

        y = 189 ;

        x = 201 ;

variables:

        float pr(time, y, x) ;

                pr:missing_value = 1.e+30f ;

                pr:long_name = "Precipitation" ;

                pr:standard_name = "precipitation_flux" ;

                pr:cell_methods = "time:mean" ;

                pr:_FillValue = 1.e+30f ;

                pr:units = "kg m-2 s-1" ;

        float time(time) ;

                time:units = "days since 1950-01-01 00:00:00" ;

                time:calendar = "gregorian" ;

                time:title = "Time" ;

                time:long_name = "time" ;

                time:time_origin = "01-JAN-1950" ;

                time:bounds = "time_bnds" ;

                time:actual_range = "1950-01-01 00:00:00 2001-01-01 00:00:00" ;

                time:delta_t = "0000-00-01 00:00:00" ;

        double x(x) ;

                x:axis = "X" ;

                x:standard_name = "projection_x_coordinate" ;

                x:long_name = "x coordinate of projection" ;

                x:units = "km" ;

        double y(y) ;

                y:axis = "Y" ;

                y:standard_name = "projection_y_coordinate" ;

                y:long_name = "y coordinate of projection" ;

                y:units = "km" ;

        float lat(y, x) ;

                lat:units = "degrees_north" ;

                lat:long_name = "Latitude" ;

                lat:nav_model = "Default grid" ;

                lat:standard_name = "Latitude" ;

        float lon(y, x) ;

                lon:units = "degrees_east" ;

                lon:long_name = "Longitude" ;

                lon:nav_model = "Default grid" ;

                lon:standard_name = "Longitude" ;

        float Lambert_Conformal ;

                Lambert_Conformal:grid_mapping_name = "lambert_conformal_conic" ;

                Lambert_Conformal:longitude_of_central_meridian = 30. ;

                Lambert_Conformal:latitude_of_projection_origin = 50. ;

                Lambert_Conformal:cone_type = "secant" ;

                Lambert_Conformal:northern_parallel = "60.0N" ;

                Lambert_Conformal:southern_parallel = "30.0N" ;

// global attributes:

                :Conventions = "CF-1.0" ;

                :institution = "UCLM" ;

                :source = "PROMES model 2005" ;

                :experiment = "Control simulation" ;

}

Could you help me to understand what is the problem?

Gratefull.

Pedro