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

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