Skip to content

Commit

Permalink
DOC: Add xradar gridding example
Browse files Browse the repository at this point in the history
  • Loading branch information
mgrover1 committed Oct 5, 2023
1 parent 7bc1922 commit d69b02c
Showing 1 changed file with 40 additions and 0 deletions.
40 changes: 40 additions & 0 deletions examples/xradar/plot_grid_xradar.py
Original file line number Diff line number Diff line change
@@ -0,0 +1,40 @@
"""
==================================
Plot a PPI Using Xradar and Py-ART
==================================
An example which uses xradar and Py-ART to grid a PPI file.
"""

# Author: Max Grover ([email protected])
# License: BSD 3 clause


import xradar as xd

import pyart
from pyart.testing import get_test_data

# Locate the test data and read in using xradar
filename = get_test_data("swx_20120520_0641.nc")
tree = xd.io.open_cfradial1_datatree(filename)

# Give the tree Py-ART radar methods
radar = pyart.xradar.Xradar(tree)

# Grid using 11 vertical levels, and 101 horizontal grid cells at a resolution on 1 km
grid = pyart.map.grid_from_radars(
(radar,),
grid_shape=(11, 101, 101),
grid_limits=(
(0.0, 10_000),
(-50_000.0, 50_000.0),
(-50_000, 50_000.0),
),
)

display = pyart.graph.GridMapDisplay(grid)
display.plot_grid(
"reflectivity_horizontal", level=0, vmin=-20, vmax=60, cmap="ChaseSpectral"
)

0 comments on commit d69b02c

Please sign in to comment.