import pyart
from matplotlib import pyplot as plt
%matplotlib inline
filename = 'data/nsaxsaprrhiC1.a1.20131203.141936.nc'

radar = pyart.io.read(filename)

dir(radar)

radar.azimuth.keys()

radar.azimuth['standard_name']

radar.azimuth['data']

f = plt.figure(figsize=[15,8])
plt.plot(radar.time['data'], radar.azimuth['data'] )
plt.xlabel(radar.time['standard_name'] + ' (' + radar.time['units'] + ')')
plt.ylabel(radar.azimuth['standard_name'] + ' (' + radar.azimuth['units'] + ')')

print radar.range['data'].min(), radar.range['data'].max(), radar.range['units']
f = plt.figure(figsize=[15,8])
plt.plot(radar.time['data'], radar.elevation['data'] )
plt.xlabel(radar.time['standard_name'] + ' (' + radar.time['units'] + ')')
plt.ylabel(radar.elevation['standard_name'] + ' (' + radar.elevation['units'] + ')')

for mykey in radar.metadata.keys():
    print mykey, ': ', radar.metadata[mykey]

radar.sweep_end_ray_index['data']
f = plt.figure(figsize=[15,8])
for i in range(len(radar.sweep_end_ray_index['data'])):
    start_index = radar.sweep_start_ray_index['data'][i]
    end_index = radar.sweep_end_ray_index['data'][i]
    plt.plot(radar.time['data'][start_index:end_index], 
             radar.elevation['data'][start_index:end_index], 
             label = 'Sweep number '+ str(radar.sweep_number['data'][i]))
plt.legend()
plt.xlabel(radar.time['standard_name'] + ' (' + radar.time['units'] + ')')
plt.ylabel(radar.elevation['standard_name'] + ' (' + radar.elevation['units'] + ')')    

print radar.fields.keys()
print ""
for mykey in radar.fields.keys():
    print mykey,':', radar.fields[mykey]['standard_name'] + ' (' + radar.fields[mykey]['units'] + ')'

f = plt.figure(figsize=[15,8])
my_pc = plt.pcolormesh(radar.range['data'], radar.time['data'],
                       radar.fields['reflectivity_horizontal']['data'])
plt.xlabel(radar.range['standard_name'] + ' (' + radar.range['units'] + ')')
plt.ylabel(radar.time['standard_name'] + ' (' + radar.time['units'] + ')')
cb = plt.colorbar(mappable = my_pc)
cb.set_label(radar.fields['reflectivity_horizontal']['standard_name'] +\
             ' (' + radar.fields['reflectivity_horizontal']['units'] + ')')
    


f = plt.figure(figsize=[15,8])
start_index = radar.sweep_start_ray_index['data'][0]
end_index = radar.sweep_end_ray_index['data'][0]
my_pc = plt.pcolormesh(radar.range['data'], radar.time['data'][start_index:end_index],
                       radar.fields['reflectivity_horizontal']['data'][start_index:end_index, :])
plt.xlabel(radar.range['standard_name'] + ' (' + radar.range['units'] + ')')
plt.ylabel(radar.time['standard_name'] + ' (' + radar.time['units'] + ')')
cb = plt.colorbar(mappable = my_pc)
cb.set_label(radar.fields['reflectivity_horizontal']['standard_name'] +\
             ' (' + radar.fields['reflectivity_horizontal']['units'] + ')')
    


radar.info('compact')   # see what happens with 'standard' or 'full'

!radar_info --compact data/nsaxsaprrhiC1.a1.20131203.141936.nc