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

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

_______________________________________________
pyngl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/pyngl-talk
Received on Tue May 26 2009 - 23:50:35 MDT