Project

General

Profile

Grid definition

Added by Serena ILLIG almost 14 years ago

Hello,

I would like to interpolate ECMWF Grib files onto my model Grid.
I am trying to define a simple grid that is regular in longitude and slihtly irregular in latitudes.
I use the PINGO convention (see below) to define it. I also tried the CDO convention but I do not know if is is correct to choose 'lonlat' as gridtype, and the results is the same.
My problem is that I get a regularly spaced latitudes.
Here are the CDO commands I execute :

cdo -s -R remapbil,bigben_grid4cdo /data/Observations/ERA_INTERIM/Grib/EI_ecmwf_200001.grib temp.gribR
cdo -s -f nc copy -setreftime,19000101,000000,days -setpartab,ecmwf temp.gribR JUNK/temp_EI_ecmwf_Y2000M1.nc

Thanks for your help,

Serena


Grid Description File
(Comments start at non digit characters and end at end of line)
First part: The dimensions.
114 226 = Number of longitudes and latitudes
Second part: The listed longitudes.
2 means equidistant longitudes
-1.500000 -1.333333 = Most western and second most western longitude
Third part: The listed latitudes
226 means all latitudes are given in the following list:
-29.000000 -28.854230 -28.708255 -28.562076 -28.415692 -28.269106 -28.122317 -27.975326 -27.828135 -27.680743 -27.533151 -27.385361 -27.237372 -27.089186 -26.940802 -26.792223 -26.643449 -26.494480 -26.345317 -26.195961 -26.046413 -25.896673 -25.746742 -25.596622 -25.446312 -25.295814 -25.145128 -24.994256 -24.843198 -24.691954 -24.540526 -24.388915 -24.237121 -24.085145 -23.932989 -23.780652 -23.628136 -23.475442 -23.322570 -23.169521 -23.016297 -22.862898 -22.709325 -22.555580 -22.401662 -22.247572 -22.093313 -21.938884 -21.784287 -21.629523 -21.474591 -21.319495 -21.164233 -21.008809 -20.853221 -20.697472 -20.541562 -20.385492 -20.229264 -20.072878 -19.916335 -19.759636 -19.602783 -19.445776 -19.288617 -19.131306 -18.973844 -18.816233 -18.658473 -18.500566 -18.342513 -18.184314 -18.025971 -17.867485 -17.708857 -17.550088 -17.391179 -17.232131 -17.072946 -16.913624 -16.754166 -16.594575 -16.434850 -16.274993 -16.115005 -15.954887 -15.794641 -15.634267 -15.473767 -15.313141 -15.152392 -14.991519 -14.830525 -14.669411 -14.508177 -14.346825 -14.185356 -14.023771 -13.862072 -13.700260 -13.538335 -13.376299 -13.214154 -13.051900 -12.889539 -12.727072 -12.564501 -12.401825 -12.239048 -12.076169 -11.913191 -11.750114 -11.586940 -11.423669 -11.260304 -11.096846 -10.933296 -10.769654 -10.605923 -10.442104 -10.278197 -10.114205 -9.950128 -9.785969 -9.621727 -9.457405 -9.293004 -9.128524 -8.963969 -8.799338 -8.634633 -8.469855 -8.305006 -8.140087 -7.975100 -7.810045 -7.644924 -7.479739 -7.314490 -7.149180 -6.983809 -6.818379 -6.652891 -6.487347 -6.321747 -6.156094 -5.990389 -5.824632 -5.658826 -5.492971 -5.327070 -5.161123 -4.995132 -4.829099 -4.663024 -4.496908 -4.330755 -4.164564 -3.998338 -3.832076 -3.665782 -3.499457 -3.333101 -3.166716 -3.000304 -2.833866 -2.667403 -2.500917 -2.334409 -2.167881 -2.001333 -1.834768 -1.668187 -1.501591 -1.334982 -1.168360 -1.001728 -0.835087 -0.668438 -0.501783 -0.335122 -0.168459 -0.001793 0.164874 0.331540 0.498204 0.664864 0.831520 0.998169 1.164810 1.331442 1.498064 1.664674 1.831270 1.997852 2.164417 2.330965 2.497494 2.664002 2.830489 2.996952 3.163391 3.329803 3.496189 3.662545 3.828871 3.995166 4.161428 4.327655 4.493846 4.660001 4.826116 4.992192 5.158227 5.324218 5.490166 5.656068 5.821923 5.987730 6.153488 6.319194 6.484848 6.650448 6.815993 6.981482 7.146913


Replies (2)

RE: Grid definition - Added by Uwe Schulzweida almost 14 years ago

Hi Serena,
It is not possible to store irregular spaced lonlat grids in GRIB format.
So, you have either to store the result of the remapping in netCDF format

cdo -s -R -f nc remapbil,bigben_grid4cdo /data/Observations/ERA_INTERIM/Grib/EI_ecmwf_200001.grib temp.nc
cdo -s -f nc setreftime,19000101,000000,days -setpartab,ecmwf temp.nc JUNK/temp_EI_ecmwf_Y2000M1.nc
or to overwrite the "wrong" grid description of the resulting GRIB file
cdo -s -R remapbil,bigben_grid4cdo /data/Observations/ERA_INTERIM/Grib/EI_ecmwf_200001.grib temp.gribR
cdo -s -f nc setreftime,19000101,000000,days -setpartab,ecmwf -setgrid,bigben_grid4cdo temp.gribR JUNK/temp_EI_ecmwf_Y2000M1.nc
Regards,
Uwe

RE: Grid definition - Added by Serena ILLIG almost 14 years ago

Thanks a lot.
It's good to know that Grib do not like irregular spaced lonlat grids...
I also manage to do everything in one command and the output grid remain irregular.

Thanks for your help,

Serena

    (1-2/2)