Example pages containing:
tips |
resources |
functions/procedures
NCL Graphics: WRF Interpolations and Cross Sections
wrf_user_interp_line is used to the interpolation across a horizontal line, using an angle of 90 and a pivot point in the lat/lon plane. The result is an array that has been interpolated at a particular latitude value, for a series of longitude values.
Three plots are created and overlaid on a single map plot:
- contour line plot of height
- filled contour plot of wind speed
- wind bars (vector) plot of U and V
By setting the special linecoords attribute to True, the lat/lon values representing the cross section are returned, and used to label the X axis.
A Python version of this projection is available here.
A different type of plot is created for each level, using wrf_contour, wrf_vector, and wrf_map_overlays. This results in some nice titles above each plot and an informational label below each plot.
This script was written to handle WRF files that may have multiple timesteps. This particular file only had one time step.
This time, however, start and end lat/lon points are used to define the line for the interpolation in the lat/lon plane, resulting in an interpolation across a series of lat, lon pairs.
The XY plot displays the lat, lon pairs on the X axis. As an added bonus, a map plot is added as an annotation of the XY plot, showing the lat/lon cross-section line.