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

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



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