In [1]:

    

%matplotlib inline

from pyradarmet.BeamBlock import BeamBlock as BeamBlock


    

In [2]:

    

# The first example is the placement for the SPOL radar during the NAME campaign (2006?)
# 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_spol = BeamBlock(radar_lon=-107., radar_lat=24., radar_antenna_height=10., 
               elev_angle=0.8, beamwidth=.91, 
               range_dist=150., azimuthal_angle=0.)


    

    




/Users/nickguy/software/python/anaconda/lib/python2.7/site-packages/pyradarmet-0.1.1-py2.7.egg/pyradarmet/geometry.py:318: RuntimeWarning: invalid value encountered in sqrt
  Numer = (y * np.sqrt(a**2 - y**2)) + (a**2 * np.arcsin(y/a)) + (np.pi * a**2 /2.)
/Users/nickguy/software/python/anaconda/lib/python2.7/site-packages/pyradarmet-0.1.1-py2.7.egg/pyradarmet/geometry.py:318: RuntimeWarning: divide by zero encountered in divide
  Numer = (y * np.sqrt(a**2 - y**2)) + (a**2 * np.arcsin(y/a)) + (np.pi * a**2 /2.)
/Users/nickguy/software/python/anaconda/lib/python2.7/site-packages/pyradarmet-0.1.1-py2.7.egg/pyradarmet/geometry.py:318: RuntimeWarning: invalid value encountered in arcsin
  Numer = (y * np.sqrt(a**2 - y**2)) + (a**2 * np.arcsin(y/a)) + (np.pi * a**2 /2.)

In [3]:

    

# Now Plot the 3 panel display
# All the of the keywords are shown for interest, we'll use default colormaps
bb_spol.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)


    

    




/Users/nickguy/software/python/anaconda/lib/python2.7/site-packages/matplotlib/colors.py:583: RuntimeWarning: invalid value encountered in less
  cbook._putmask(xa, xa < 0.0, -1)






    

In [4]:

    

# Now the next example is the NPOL radar deployment near Alamosa in 2014
bb_npol = BeamBlock(radar_lon=-105.86, radar_lat=37.44, radar_antenna_height=4., 
               elev_angle=0.9, beamwidth=.88, 
               range_dist=150., azimuthal_angle=0.)


    

In [5]:

    

# Now Plot the 3 panel display
# All the of the keywords are shown for interest, we'll use default colormaps
bb_npol.plot_3panel(beam_blockage='complete',
                    range_rings=(30., 60., 90.),
                    saveFig=False)


    

In [ ]: