[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

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/
===============================================================================