In [8]:

    

import pyart
from matplotlib import pyplot as plt
import numpy as np
from matplotlib import rc
%matplotlib inline


    

In [30]:

    

#radar = pyart.io.read('/data/radar/nsaxsaprppiC1.a1.20140209.232402.nc')
radar = pyart.io.read('/data/radar/nsaxsaprppiC1.a1.20140201.184802.nc')


    

In [31]:

    

print radar.fields.keys()
print radar.fixed_angle['data']
print radar.range['data'].max()
print radar.longitude['data']
print radar.latitude['data']


    

    




[u'reflectivity_horizontal', u'cross_correlation_ratio', u'normalized_coherent_power', u'total_power', u'mean_doppler_velocity', u'doppler_spectrum_width', u'differential_reflectivity', u'specific_differential_phase', u'differential_phase']
[  0.           0.49987793   0.99975586   1.6973877    2.49938965
   3.49914551   4.49890137   5.49865723   6.49841309   7.49816895
   8.99780273  10.99731445  13.50219727  16.50146484  20.00061035  45.        ]
49950.0
[-156.66809997]
[ 71.32494003]

In [45]:

    

font = {'size' : 16}
rc('font', **font)
f = plt.figure(figsize=[15,8])
myd = pyart.graph.RadarMapDisplay(radar)
myd.plot_ppi_map('reflectivity_horizontal', 1, vmin=-20, vmax=20, 
                 min_lon=-157.1, max_lon=-156, min_lat=71.2, max_lat=71.6,
                 lon_lines = np.arange(-158, -154, .2), 
                 lat_lines = np.arange(69, 72, .1), resolution = 'h',
                 auto_range=False)
ax1 = plt.gca()
x0,x1 = ax1.get_xlim()
y0,y1 = ax1.get_ylim()
ax1.set_aspect((x1-x0)/(y1-y0) - .15)


    

In [56]:

    

grid = pyart.map.grid_from_radars(
    (radar,),
    grid_shape=(21, 201, 201),
    grid_limits=((0, 5000.), (-25000., 25000.), (-25000., 25000.)),
    fields=radar.fields.keys(), refl_field='reflectivity_horizontal')


    

In [62]:

    

# panel sizes
map_panel_axes = [0.05, 0.05, .4, .80]
x_cut_panel_axes = [0.55, 0.10, .4, .30]
y_cut_panel_axes = [0.55, 0.50, .4, .30]
colorbar_panel_axes = [0.05, 0.90, .4, .03]

#min_lon=-157.1, max_lon=-156, min_lat=71.2, max_lat=71.6

# parameters
level = 5
vmin = -20
vmax = 20
lat = radar.latitude['data']
lon = radar.longitude['data']

fig = plt.figure(figsize=[15, 8])

grid_display = pyart.graph.GridMapDisplay(grid)


# panel 1, basemap, radar reflectivity and NARR overlay
ax1 = fig.add_axes(map_panel_axes)
grid_display.plot_basemap(resolution = 'h',
                       lon_lines = np.arange(-158, -154, .2), 
                       lat_lines = np.arange(69, 72, .1))

grid_display.plot_grid('reflectivity_horizontal', 
                       level=level, vmin=vmin, vmax=vmax)

grid_display.plot_crosshairs(lon=lon, lat=lat)

#grid_display.plot_crosshairs()

# colorbar
cbax = fig.add_axes(colorbar_panel_axes)
grid_display.plot_colorbar(cax=cbax)

# panel 2, longitude slice.
ax2 = fig.add_axes(x_cut_panel_axes)
grid_display.plot_longitude_slice('reflectivity_horizontal', 
                                  vmin=vmin, vmax=vmax, lon=lon, lat=lat)
ax2.set_xlabel('Distance from radar (km)')

# panel 3, latitude slice
ax3 = fig.add_axes(y_cut_panel_axes)
grid_display.plot_latitude_slice('reflectivity_horizontal',
                                 vmin=vmin, vmax=vmax, lon=lon, lat=lat)

x0,x1 = ax1.get_xlim()
y0,y1 = ax1.get_ylim()
ax1.set_aspect((x1-x0)/(y1-y0))


    

In [63]:

    

pyart.io.write_grid('/data/radar/nsa_grid.nc', grid)


    

In [ ]: