Re: New Ticket - [netCDFDecoders !FEG-733021]: Problem reading ECMWF Grib file

Yes, we are not supporting quasi grids at this point.

Robb Kambic wrote:
> On Tue, 19 Dec 2006, address@hidden wrote:
>
>
> > New Ticket: Problem reading ECMWF Grib file
>
>
> this is a quasi-regular Grib1 grid, we decided not to do the work since
> it's Grib1.  are we going to bring this up?  last quasi work was 3-5 days
> effort.  i commented out the print in the code so it won't come up in IDV
> and did another grib jar ci.
>
> robb...
>
>
>
>
>
>
> > Hi Robb-
> >
> > I have a GRIB file from ECMWF that does not give the correct
> > navigation:
> >
> > /upc/share/testdata/grid/grib/grib1/data/ana_Z_20031006_00_ml.grib
> >
> > When I run Grib1Dump on this, it prints out:
> >
> > $ java  ucar.grib.grib1.Grib1Dump ana_Z_20031006_00_ml.grib
> > GDS  numPV[ 0 ] = 0
> > GDS  numPV[ 1 ] = 0
> > GDS  numPV[ 2 ] = 16916
> > GDS  numPV[ 3 ] = 0
> > GDS  numPV[ 4 ] = 16934
> > GDS  numPV[ 5 ] = 27875
> > GDS  numPV[ 6 ] = 16959
> > GDS  numPV[ 7 ] = -23082
> > GDS  numPV[ 8 ] = 16991
> > GDS  numPV[ 9 ] = -23792
> > GDS  numPV[ 10 ] = 17030
> > GDS  numPV[ 11 ] = 31674
> > GDS  numPV[ 12 ] = 17076
> > GDS  numPV[ 13 ] = -27240
> > GDS  numPV[ 14 ] = 17130
> > GDS  numPV[ 15 ] = -14480
> > GDS  numPV[ 16 ] = 17170
> > GDS  numPV[ 17 ] = -22545
> > GDS  numPV[ 18 ] = 17175
> > GDS  numPV[ 19 ] = 24461
> > GDS  numPV[ 20 ] = 17181
> > GDS  numPV[ 21 ] = 2532
> > GDS  numPV[ 22 ] = 17187
> > GDS  numPV[ 23 ] = -1429
> > GDS  numPV[ 24 ] = 17196
> > GDS  numPV[ 25 ] = -27779
> > GDS  numPV[ 26 ] = 17207
> > GDS  numPV[ 27 ] = 14993
> > GDS  numPV[ 28 ] = 17220
> > GDS  numPV[ 29 ] = 27995
> > GDS  numPV[ 30 ] = 17236
> > GDS  numPV[ 31 ] = -14440
> > GDS  numPV[ 32 ] = 17257
> > GDS  numPV[ 33 ] = 2623
> > GDS  numPV[ 34 ] = 17282
> > GDS  numPV[ 35 ] = 9314
> > GDS  numPV[ 36 ] = 17313
> > GDS  numPV[ 37 ] = 15928
> > GDS  numPV[ 38 ] = 17351
> > GDS  numPV[ 39 ] = -14657
> > GDS  numPV[ 40 ] = 17399
> > GDS  numPV[ 41 ] = -31574
> > GDS  numPV[ 42 ] = 17427
> > GDS  numPV[ 43 ] = 10933
> > GDS  numPV[ 44 ] = 17431
> > GDS  numPV[ 45 ] = -32251
> > GDS  numPV[ 46 ] = 17436
> > GDS  numPV[ 47 ] = -30046
> > GDS  numPV[ 48 ] = 17442
> > GDS  numPV[ 49 ] = 15630
> > GDS  numPV[ 50 ] = 17448
> > GDS  numPV[ 51 ] = -30688
> > GDS  numPV[ 52 ] = 17455
> > GDS  numPV[ 53 ] = 11634
> > GDS  numPV[ 54 ] = 17461
> > GDS  numPV[ 55 ] = -12461
> > GDS  numPV[ 56 ] = 17468
> > GDS  numPV[ 57 ] = 5070
> > GDS  numPV[ 58 ] = 17473
> > GDS  numPV[ 59 ] = -19591
> > GDS  numPV[ 60 ] = 17478
> > GDS  numPV[ 61 ] = 32047
> > GDS  numPV[ 62 ] = 17482
> > GDS  numPV[ 63 ] = 21426
> > GDS  numPV[ 64 ] = 17485
> > GDS  numPV[ 65 ] = 11036
> > GDS  numPV[ 66 ] = 17486
> > GDS  numPV[ 67 ] = -460
> > GDS  numPV[ 68 ] = 17487
> > GDS  numPV[ 69 ] = -12835
> > GDS  numPV[ 70 ] = 17487
> > GDS  numPV[ 71 ] = -24453
> > GDS  numPV[ 72 ] = 17486
> > GDS  numPV[ 73 ] = -32409
> > GDS  numPV[ 74 ] = 17484
> > GDS  numPV[ 75 ] = -32684
> > GDS  numPV[ 76 ] = 17481
> > GDS  numPV[ 77 ] = -20288
> > GDS  numPV[ 78 ] = 17478
> > GDS  numPV[ 79 ] = 10588
> > GDS  numPV[ 80 ] = 17474
> > GDS  numPV[ 81 ] = 888
> > GDS  numPV[ 82 ] = 17469
> > GDS  numPV[ 83 ] = 23155
> > GDS  numPV[ 84 ] = 17464
> > GDS  numPV[ 85 ] = 19232
> > GDS  numPV[ 86 ] = 17458
> > GDS  numPV[ 87 ] = -3272
> > GDS  numPV[ 88 ] = 17453
> > GDS  numPV[ 89 ] = 28866
> > GDS  numPV[ 90 ] = 17447
> > GDS  numPV[ 91 ] = -7808
> > GDS  numPV[ 92 ] = 17442
> > GDS  numPV[ 93 ] = 25179
> > GDS  numPV[ 94 ] = 17437
> > GDS  numPV[ 95 ] = 3790
> > GDS  numPV[ 96 ] = 17432
> > GDS  numPV[ 97 ] = 81
> > GDS  numPV[ 98 ] = 17427
> > GDS  numPV[ 99 ] = 19911
> > GDS  numPV[ 100 ] = 17392
> > GDS  numPV[ 101 ] = -20835
> > GDS  numPV[ 102 ] = 17332
> > GDS  numPV[ 103 ] = 31525
> > GDS  numPV[ 104 ] = 17280
> > GDS  numPV[ 105 ] = -902
> > GDS  numPV[ 106 ] = 17238
> > GDS  numPV[ 107 ] = -24934
> > GDS  numPV[ 108 ] = 17205
> > GDS  numPV[ 109 ] = 30154
> > GDS  numPV[ 110 ] = 17181
> > GDS  numPV[ 111 ] = 13654
> > GDS  numPV[ 112 ] = 17106
> > GDS  numPV[ 113 ] = 25814
> > GDS  numPV[ 114 ] = 16961
> > GDS  numPV[ 115 ] = -7259
> > GDS  numPV[ 116 ] = 16757
> > GDS  numPV[ 117 ] = -7610
> > GDS  numPV[ 118 ] = 0
> > GDS  numPV[ 119 ] = 0
> > GDS  numPV[ 120 ] = 0
> > GDS  numPV[ 121 ] = 0
> > GDS  numPV[ 122 ] = 0
> > GDS  numPV[ 123 ] = 0
> > GDS  numPV[ 124 ] = 0
> > GDS  numPV[ 125 ] = 0
> > GDS  numPV[ 126 ] = 0
> > GDS  numPV[ 127 ] = 0
> > GDS  numPV[ 128 ] = 0
> > GDS  numPV[ 129 ] = 0
> > GDS  numPV[ 130 ] = 0
> > GDS  numPV[ 131 ] = 0
> > GDS  numPV[ 132 ] = 0
> > GDS  numPV[ 133 ] = 0
> > GDS  numPV[ 134 ] = 0
> > GDS  numPV[ 135 ] = 0
> > GDS  numPV[ 136 ] = 0
> > GDS  numPV[ 137 ] = 0
> > GDS  numPV[ 138 ] = 0
> > GDS  numPV[ 139 ] = 0
> > GDS  numPV[ 140 ] = 0
> > GDS  numPV[ 141 ] = 0
> > GDS  numPV[ 142 ] = 0
> > GDS  numPV[ 143 ] = 0
> > GDS  numPV[ 144 ] = 0
> > GDS  numPV[ 145 ] = 0
> > GDS  numPV[ 146 ] = 0
> > GDS  numPV[ 147 ] = 0
> > GDS  numPV[ 148 ] = 0
> > GDS  numPV[ 149 ] = 0
> > GDS  numPV[ 150 ] = 0
> > GDS  numPV[ 151 ] = 0
> > GDS  numPV[ 152 ] = 0
> > GDS  numPV[ 153 ] = 0
> > GDS  numPV[ 154 ] = 0
> > GDS  numPV[ 155 ] = 0
> > GDS  numPV[ 156 ] = 0
> > GDS  numPV[ 157 ] = 0
> > GDS  numPV[ 158 ] = 0
> > GDS  numPV[ 159 ] = 0
> > GDS  numPV[ 160 ] = 0
> > GDS  numPV[ 161 ] = 0
> > GDS  numPV[ 162 ] = 0
> > GDS  numPV[ 163 ] = 0
> > GDS  numPV[ 164 ] = 0
> > GDS  numPV[ 165 ] = 0
> > GDS  numPV[ 166 ] = 0
> > GDS  numPV[ 167 ] = 0
> > GDS  numPV[ 168 ] = 0
> > GDS  numPV[ 169 ] = 0
> > GDS  numPV[ 170 ] = 15695
> > GDS  numPV[ 171 ] = -32329
> > GDS  numPV[ 172 ] = 15902
> > GDS  numPV[ 173 ] = 15596
> > GDS  numPV[ 174 ] = 15990
> > GDS  numPV[ 175 ] = -2748
> > GDS  numPV[ 176 ] = 16148
> > GDS  numPV[ 177 ] = -12304
> > GDS  numPV[ 178 ] = 16173
> > GDS  numPV[ 179 ] = -23504
> > GDS  numPV[ 180 ] = 16212
> >
> > For a particular grid:
> >
> > --------------------------------------------------------------------
> >                         Header : GRIB1
> >                     Discipline : 0 Meteorological Products
> >                   GRIB Edition : 1
> >                    GRIB length : 130916
> >             Originating Center : 98 European Center for Medium-Range
> > Weather For
> > ecasts (RSMC)
> >         Originating Sub-Center : 0 WMO Secretariat
> >             Product Definition : 0 Forecast/Uninitialized
> > Analysis/Image Product
> >
> >             Parameter Category : -1 Meteorological Parameters
> >                 Parameter Name : 131 U U velocity
> >                Parameter Units : m s^-1
> >                 Reference Time : 2003-10-06T00:00:00Z
> >                     Time Units : hour
> >           Time Range Indicator : product valid at RT + P1
> >                    Time 1 (P1) : 0
> >                    Time 2 (P2) : 0
> >        Generating Process Type : 203 Unknown
> >                     Level Type : 109 hybrid level
> >                  Level Value 1 : 60.0
> >                  Level Value 2 : 0.0
> >                     GDS Exists : true
> >                     BMS Exists : false
> >          Number of data points : 65160
> >                      Grid Name : Latitude/Longitude Grid
> >                      Grid Shape: 0 spherical
> >          Spherical earth radius: 6367.47
> >                             Nx : 360
> >                             Ny : 181
> >                            La1 : 90.0
> >                            Lo1 : 0.0
> >   Resolution & Component flags : 128
> >                            La2 : -90.0
> >                            Lo2 : 359.0
> >                             Dx : -9999.0
> >                             Dy : 1.0
> >                  Scanning mode : 0
> >
> > notice that the Dx is -9999.0.  This ends up creating non-sensical
> > longitudes.
> >
> > So, one issue is the debug messages printing out and the other is
> > the bad longitude values.  This shows up in the tools UI as well as
> > the IDV.  In the Tools UI, a dump of the longitude variable gives:
> >
> >
> >    double lon(lon=360);
> >      :units = "degrees_east";
> >      :long_name = "longitude coordinate";
> >      :standard_name = "longitude";
> >      :grid_spacing = "-9999.0 degrees_east";
> >      :_CoordinateAxisType = "Lon";
> >
> >  data:
> >
> >   {0.0, -9999.0, -19998.0, -29997.0, -39996.0, -49995.0, -59994.0,
> > -69993.0, -79992.0, -89991.0, -99990.0, -109989.0, -119988.0, -129987.0,
> > -139986.0, -149985.0, -159984.0, -169983.0, -179982.0, -189981.0,
> > -199980.0, -209979.0, -219978.0, -229977.0, -239976.0, -249975.0,
> > -259974.0, -269973.0, -279972.0, -289971.0, -299970.0, -309969.0,
> > -319968.0, -329967.0, -339966.0, -349965.0, -359964.0, -369963.0,
> > -379962.0, -389961.0, -399960.0, -409959.0, -419958.0, -429957.0,
> > -439956.0, -449955.0, -459954.0, -469953.0, -479952.0, -489951.0,
> > -499950.0, -509949.0, -519948.0, -529947.0, -539946.0, -549945.0,
> > -559944.0, -569943.0, -579942.0, -589941.0, -599940.0, -609939.0,
> > -619938.0, -629937.0, -639936.0, -649935.0, -659934.0, -669933.0,
> > -679932.0, -689931.0, -699930.0, -709929.0, -719928.0, -729927.0,
> > -739926.0, -749925.0, -759924.0, -769923.0, -779922.0, -789921.0,
> > -799920.0, -809919.0, -819918.0, -829917.0, -839916.0, -849915.0,
> > -859914.0, -869913.0, -879912.0, -889911.0, -899910.0, -909909.0,
> > -919908.0, -929907.0, -939906.0, -949905.0, -959904.0, -969903.0,
> > -979902.0, -989901.0, -999900.0, -1009899.0, -1019898.0, -1029897.0,
> > -1039896.0, ....................
> >
> >
> > There is another file in the directory (ana_Z_20031006_00.grib) from
> > ECMWF which is on pressure coordinates. It does not have this problem.
> >
> > Don
> > *************************************************************
> > Don Murray                               UCAR Unidata Program
> > address@hidden                        P.O. Box 3000
> > (303) 497-8628                              Boulder, CO 80307
> > http://www.unidata.ucar.edu/staff/donm
> > *************************************************************
> >
> >
> >
> >
> > Ticket Details
> > ===================
> > Ticket ID: FEG-733021
> > Department: Support netCDF Decoders
> > Priority: Normal
> > Status: Open
> > Link: 
> > http://www.unidata.ucar.edu/esupport/staff/index.php?_m=tickets&_a=viewticket&ticketid=2889
>
>
> ===============================================================================
> Robb Kambic                                Unidata Program Center
> Software Engineer III                      Univ. Corp for Atmospheric Research
> address@hidden             WWW: http://www.unidata.ucar.edu/
> ===============================================================================