In [1]:
%matplotlib inline
from pyradarmet.BeamBlock import BeamBlock as BeamBlock
In [6]:
# This test is to see about the implementation of combining multiple DEM files
# when the point occurs close to an edge.
# Line 1: Supply the radar latitude, longitude, (No altitude means it is calculated)
# Line 2: Elevation angle of interest and beam width of radar
# Line 3: Distance from radar and azimuthal angle for profile of interest
bb_test_edge = BeamBlock(radar_lon=-100.5, radar_lat=39.9, radar_antenna_height=10.,
elev_angle=0.8, beamwidth=.91,
range_dist=150., azimuthal_angle=0.)
In [7]:
# Now Plot the 3 panel display
# All the of the keywords are shown for interest, we'll use default colormaps
bb_test_edge.plot_3panel(beam_blockage='complete',
lon_plot_range=2., lat_plot_range=2.,
terrain_cmap=None, beam_blockage_cmap=None,
range_rings=(50., 100.),
saveFig=False)
In [12]:
# Now the next example is the NPOL radar deployment near Alamosa in 2014
bb_test2 = BeamBlock(radar_lon=-60., radar_lat=9.5, radar_antenna_height=4.,
elev_angle=0.9, beamwidth=.88,
range_dist=150., azimuthal_angle=0.)
In [13]:
# Now Plot the 3 panel display
# All the of the keywords are shown for interest, we'll use default colormaps
bb_test2.plot_3panel(beam_blockage='complete',
range_rings=(30., 60., 90.),
saveFig=False)
In [ ]: