In [1]:
#first we do some imports and check the version of Py-ART for consistency
import pyart
from matplotlib import pyplot as plt
import netCDF4
import numpy as np
from copy import deepcopy
import copy
import scipy
from matplotlib import rc
from netCDF4 import num2date
import h5py
%matplotlib inline
print pyart.__version__
In [2]:
filename = '/Volumes/LaCie/data/radar/kurnel/SCollis_KurnellHDF5_20140902_1800-1900/5402_20140902_182402.vol.h5'
radar = pyart.aux_io.read_odim_h5(filename, file_field_names=True)
In [14]:
print radar.fields.keys()
print radar.fixed_angle['data']
In [19]:
myd = pyart.graph.RadarMapDisplay(radar)
fig = plt.figure(figsize = [18,10])
myd.plot_ppi_map( 'DBZH', vmin = -8, vmax = 64,
resolution = 'h', max_lat = -33,
min_lat = -35, min_lon = 150, max_lon = 153)
m = myd.basemap
m.drawparallels(np.linspace(-32, -35, 11),labels=[1,0,0,0])
m.drawmeridians(np.linspace(150, 153, 11),labels=[0,0,0,1])
m.drawrivers()
m.drawmapscale(150.8, -33.2, 151, -33.2, 100, barstyle='fancy', fontcolor='k', fillcolor1='b', fillcolor2='k')
Out[19]:
In [16]:
myd.plot_ppi_map?
In [ ]: