Re: Grabbing ECMWF ensemble values for a station (37.6200, -122.4000)

*our dev config:*

Linux chef 2.6.22.14-72.fc6 #1 SMP Wed Nov 21 14:10:25 EST 2007 x86_64
x86_64 x86_64 GNU/Linux
gcc (GCC) 4.1.2 20070626 (Red Hat 4.1.2-13)
GNU Fortran (GCC) 4.1.2 20070626 (Red Hat 4.1.2-13)
Python 2.5.2

*our production config:*

Linux prod 2.6.21.7-2.fc8xen #1 SMP Fri Feb 15 12:34:28 EST 2008 x86_64
GNU/Linux
gcc (GCC) 4.2.4 (Ubuntu 4.2.4-1ubuntu3)
GNU Fortran (GCC) 4.2.4 (Ubuntu 4.2.4-1ubuntu3)
Python 2.5.2

Any sample code to extract the all the 51 ensembles (min temp, max temp and
precip) for station (lat, long) is appreciated.

Thanks
Satish

On Tue, May 26, 2009 at 10:50 PM, David Brown <dbrown_at_ucar.edu> wrote:

> Hi Satish,
> This file does indeed contain ensemble data. It uses ECMWF local extension
> number 30. Support for this extension was added to the NIO library fairly
> recently. It is not part of the current release of PyNIO. We can send you a
> current version if you contact Mary Haley with the results on your system of
> executing:
> uname -a
> gcc --version
> gfortran --version
>
> -dave
>
> On May 26, 2009, at 11:24 PM, satish viswanatham wrote:
>
> Thank you Dave for helping me out. I uploaded the following file:
>
>
> W1E_20090526_0000_20090526_0000_1
>
>
> On Tue, May 26, 2009 at 9:16 PM, David Brown <*dbrown_at_ucar.edu*<dbrown_at_ucar.edu>>
> wrote:
>
> Hi Satish,
>
> The print of the file contents does not show an ensemble dimension. If you
> are sure that the file contents ensemble data, then your version of PyNIO
> does not recognize the ECMWF local extension that defines this file as an
> ensemble file. Can you send instructions detailing how to get this file, or
> else transfer it using ftp:
>
>
> ftp *ftp.cgd.ucar.edu* <http://ftp.cgd.ucar.edu/>
>
> login as anonymous with your email address as password
>
> cd incoming
>
> put <the-file-name>
>
>
> Then send me the name of the file you transferred.
>
> Thanks,
>
> -dave
>
>
>
>
>
>
> On May 26, 2009, at 6:53 PM, satish viswanatham wrote:
>
>
>
> Hi,
>
>
> I am trying to use PyNIO library to extract precip (forecast) information
> for few stations. Looks like there are two variables namely: TP_GDS0_SFC and
> TP_GDS0_SFC_1.
>
>
> f =
> Nio.open_file('/export/disk0/grb-files/W1E_20090526_0000_20090530_1200_1.grb2')
>
> print f
>
> print f.variables['TP_GDS0_SFC'].shape #(49, 53)
>
> print f.variables['TP_GDS0_SFC_1'].shape # (71, 141)
>
>
> How do I get hold of 51 ensembles for lat and long (37.6200,-122.4000)? Any
> sample code to extract that?
>
>
> Thanks
>
>
> Nio file:
> /export/disk0/grb-files/W1E_20090526_0000_20090530_1200_1.grb2
>
> global attributes:
>
> dimensions:
>
> g0_lat_0 = 49
>
> g0_lon_1 = 53
>
> g0_lat_2 = 71
>
> g0_lon_3 = 141
>
> variables:
>
> float VAR_121_GDS0_SFC [ g0_lat_0, g0_lon_1 ]
>
> center : European Center for Medium-Range Weather Forecasts -
> Reading
>
> long_name : Unknown Variable Name
>
> units : unknown
>
> _FillValue : -999
>
> level_indicator : 1
>
> gds_grid_type : 0
>
> parameter_table_version : 128
>
> parameter_number : 121
>
> forecast_time : 108
>
> forecast_time_units : hours
>
> initial_time : 05/26/2009 (00:00)
>
> float VAR_121_GDS0_SFC_1 [ g0_lat_2, g0_lon_3 ]
>
> center : European Center for Medium-Range Weather Forecasts -
> Reading
>
> long_name : Unknown Variable Name
>
> units : unknown
>
> _FillValue : -999
>
> level_indicator : 1
>
> gds_grid_type : 0
>
> parameter_table_version : 128
>
> parameter_number : 121
>
> forecast_time : 108
>
> forecast_time_units : hours
>
> initial_time : 05/26/2009 (00:00)
>
> float VAR_122_GDS0_SFC [ g0_lat_0, g0_lon_1 ]
>
> center : European Center for Medium-Range Weather Forecasts -
> Reading
>
> long_name : Unknown Variable Name
>
> units : unknown
>
> _FillValue : -999
>
> level_indicator : 1
>
> gds_grid_type : 0
>
> parameter_table_version : 128
>
> parameter_number : 122
>
> forecast_time : 108
>
> forecast_time_units : hours
>
> initial_time : 05/26/2009 (00:00)
>
> float VAR_122_GDS0_SFC_1 [ g0_lat_2, g0_lon_3 ]
>
> center : European Center for Medium-Range Weather Forecasts -
> Reading
>
> long_name : Unknown Variable Name
>
> units : unknown
>
> _FillValue : -999
>
> level_indicator : 1
>
> gds_grid_type : 0
>
> parameter_table_version : 128
>
> parameter_number : 122
>
> forecast_time : 108
>
> forecast_time_units : hours
>
> initial_time : 05/26/2009 (00:00)
>
> float TP_GDS0_SFC [ g0_lat_0, g0_lon_1 ]
>
> center : European Center for Medium-Range Weather Forecasts -
> Reading
>
> long_name : Total precipitation
>
> units : m
>
> _FillValue : -999
>
> level_indicator : 1
>
> gds_grid_type : 0
>
> parameter_table_version : 128
>
> parameter_number : 228
>
> forecast_time : 108
>
> forecast_time_units : hours
>
> initial_time : 05/26/2009 (00:00)
>
> float TP_GDS0_SFC_1 [ g0_lat_2, g0_lon_3 ]
>
> center : European Center for Medium-Range Weather Forecasts -
> Reading
>
> long_name : Total precipitation
>
> units : m
>
> _FillValue : -999
>
> level_indicator : 1
>
> gds_grid_type : 0
>
> parameter_table_version : 128
>
> parameter_number : 228
>
> forecast_time : 108
>
> forecast_time_units : hours
>
> initial_time : 05/26/2009 (00:00)
>
> float g0_lat_2 [ g0_lat_2 ]
>
> long_name : latitude
>
> GridType : Cylindrical Equidistant Projection Grid
>
> units : degrees_north
>
> Dj : 0.5
>
> Di : 0.5
>
> Lo2 : -65
>
> La2 : 25
>
> Lo1 : -135
>
> La1 : 60
>
> float g0_lon_3 [ g0_lon_3 ]
>
> long_name : longitude
>
> GridType : Cylindrical Equidistant Projection Grid
>
> units : degrees_east
>
> Dj : 0.5
>
> Di : 0.5
>
> Lo2 : -65
>
> La2 : 25
>
> Lo1 : -135
>
> La1 : 60
>
> float g0_lat_0 [ g0_lat_0 ]
>
> long_name : latitude
>
> GridType : Cylindrical Equidistant Projection Grid
>
> units : degrees_north
>
> Dj : 0.5
>
> Di : 0.5
>
> Lo2 : 15
>
> La2 : 36
>
> Lo1 : -11
>
> La1 : 60
>
> float g0_lon_1 [ g0_lon_1 ]
>
> long_name : longitude
>
> GridType : Cylindrical Equidistant Projection Grid
>
> units : degrees_east
>
> Dj : 0.5
>
> Di : 0.5
>
> Lo2 : 15
>
> La2 : 36
>
> Lo1 : -11
>
> La1 : 60
>
>
> _______________________________________________
>
> pyngl-talk mailing list
>
> List instructions, subscriber options, unsubscribe:
>
> *http://mailman.ucar.edu/mailman/listinfo/pyngl-talk***<http://mailman.ucar.edu/mailman/listinfo/pyngl-talk>
>
>
>
>

_______________________________________________
pyngl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/pyngl-talk
Received on Wed May 27 2009 - 08:32:14 MDT