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

[IDV #OUE-658238]: problem selecting spatial subset of ECMWF grib file

Hi HP-

> In fact, after double-checking there is no difference btw XP and Linux.
> My particular grib files fails spatial subsetting in both cases. The
> failing file happens to be quite large (60MB). Might be size is linked
> to the problem? I will send you the file via YouSendIt.

Thanks for the file.  I'm not understanding the problem here.  The
grid goes from 45N 0E (UL) to 5N 65E (LR) with a spacing of .25 degrees
(261x161 pts):

--------------------------------------------------------------------
                        Header : GRIB1
                    Discipline : 0 Meteorological Products
                  GRIB Edition : 1
                   GRIB length : 84150
            Originating Center : 98 European Center for Medium-Range
ecasts (RSMC)
        Originating Sub-Center : 0 WMO Secretariat
            Product Definition : 0 Forecast/Uninitialized Analysis/I
            Parameter Category : -1 Meteorological Parameters
                Parameter Name : 129 Z Geopotential
               Parameter Units : m2 s-2
                Reference Time : 2009-02-01T00:00:00Z
                    Time Units : hour
          Time Range Indicator : product valid at RT + P1
                   Time 1 (P1) : 0
                   Time 2 (P2) : 0
       Generating Process Type : 133 Unknown
                    Level Type : 100 isobaric
                 Level Value 1 : 200.0
                 Level Value 2 : 0.0
                    GDS Exists : true
                    BMS Exists : false
         Number of data points : 42021
                     Grid Name : Latitude/Longitude Grid
                     Grid Shape: 0 spherical
         Spherical earth radius: 6367.47
                            Nx : 261
                            Ny : 161
                           La1 : 45.0
                           Lo1 : 0.0
  Resolution & Component flags : 128
                         Winds : True
                           La2 : 5.0
                           Lo2 : 65.0
                            Dx : 0.25
                            Dy : 0.25
                 Scanning mode : 0

That would put the grid bounding box (blue rectangle) where it shows on 
the screen shot.  You can only subset within that region.
When I try to subset within that box, it seems to work for me with
2.6u2 and 2.7a1.

Can you send me a bundle pointing to that file on your local system
with a grid subset that doesn't work the way you think it should?

> From cutting trees, HP

Another ridge heading our way next week.  Back to spring weather.
We need snow!

Don
 
> Unidata IDV Support wrote:
> > Hi HP-
> >
> >
> >> Full Name: HansPeter Roesli
> >> Email Address: address@hidden
> >> Organization: EUMETSAT
> >> Package Version: 2.6u2 build date:@DATE@
> >> Operating System: Linux
> >> Hardware: Java: home: /home/user/IDV_2.6u2/jre version: 1.6.0_07 j3d:1.5.2 
> >> fcs (build4)
> >> Description of problem: Hi -
> >>
> >> It is quite some time I did not try to spatially subset a ECMWF grib
> >> file. I discover now to my surprise under Linux (not under XP) that it
> >> does not work properly anymore. It is only possible to subset whithin a
> >> blue rectanlge, and even here the subended crop area behaves in a
> >> strange way.
> >>
> >> Look a the attached screen capture. The problem is also with 2.7a1 of
> >> course.
> >>
> >
> > There should be no difference between Linux and XP.  I can believe there is
> > a difference between 2.6uX and 2.7a1 - we made some changes to the map
> > selection panel.
> >
> > Can you send us a sample file that we can test with?
> >
> > Thanks.
> >
> > Don
> >
> >
> > Ticket Details
> > ===================
> > Ticket ID: OUE-658238
> > Department: Support IDV
> > Priority: Normal
> > Status: Open
> >
> >
> >
> 
> 


Ticket Details
===================
Ticket ID: OUE-658238
Department: Support IDV
Priority: Normal
Status: Open