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
Received on Tue May 26 2009 - 18:53:27 MDT
This archive was generated by hypermail 2.2.0 : Mon Jun 01 2009 - 07:48:18 MDT