from pylab import *
import netCDF4
import pandas as pd
nc = netCDF4.Dataset('gdas.44512.2010040100.nc')
print netCDF4.chartostring(nc.variables['hdr_sid'][-1,:])
PRDA2
nobs,dum = shape(nc.variables['obs_arr'])
obs_arr = nc.variables['obs_arr'][:,:]
hdr_typ = netCDF4.chartostring(nc.variables['hdr_typ'][:,:])
hdr_sid = netCDF4.chartostring(nc.variables['hdr_sid'][:,:])
hdr_vld = netCDF4.chartostring(nc.variables['hdr_vld'][:,:])
obs_qty = netCDF4.chartostring(nc.variables['obs_qty'][:,:])
hdr_arr = nc.variables['hdr_arr'][:,:]
id=int8(obs_arr[:,0])
d={}
d['station_id'] = hdr_sid[id]
d['lon'] = hdr_arr[id,0]
d['lat'] = hdr_arr[id,1]
d['elev'] = hdr_arr[id,2]
d['time'] = hdr_vld[id]
d['grib_code'] = obs_arr[:,1]
d['press'] = obs_arr[:,2]
d['height'] = obs_arr[:,3]
d['obs_val'] = obs_arr[:,4]
d['quality'] = obs_qty[:]
df = pd.DataFrame(d)
df['year'] = df['time'].str[0:4]
df['mon'] = df['time'].str[4:6]
df['day'] = df['time'].str[6:8]
df['hour'] = df['time'].str[9:11]
del df['time']
print df.ix[0]
elev 15 grib_code 11 height 14.99732 lat -69.96 lon 41.66 obs_val 284.75 press 1008 quality 2 station_id 74494 mon 04 year 2010 day 01 hour 00 Name: 0, dtype: object