Vertical cross-section plot of radial wind

[MR. VIJAY VISHWAKARMA]

Hello everyone,

I want to plot vertical cross-section of radial wind ( pressure or height on y-axis and radius on x-axis). I have a WRF simulated output dataset in netcdf format which include pressue, u-wind, v-wind, time, lat, lon, lev(pressure from 1000 to 50 hpa), and other variables.

As I am currently beginning to start working on wrf-python, please suggest how I should proceed to find a solution for this plot. Please share your view or any working script. 

Thank you.

Sincerely,

Vijay

[Marco Miani]

Hi Vijay,

There is more than one way to skin a cat:

1. https://geocat-examples.readthedocs.io/en/latest/gallery/WRF/NCL_wrf_interp_3.html#sphx-glr-gallery-wrf-ncl-wrf-interp-3-py

This example illustrates the use of (python based) geocat tool, specifically for slicing your 3D data over the vertical coordinate.

2. If you are familiar with NCL (NCAR) you can have a look at this:
https://www.ncl.ucar.edu/Applications/wrfinterp.shtml

3. https://wrf-python.readthedocs.io/en/latest/plot.html
This is also a very good start, rich of examples to get inspired from.

In general, and for all 3 of the above mentioned, you can extract vertical profiles (planes) by either providing start and end coordinate point, OR via pivot point+angle.

The example 3, shows how to label your xaxis (distance along section) with both lat, lon.

Cheers

Marco