Page 1 of 1

Plot moist potential temperature at 500 hPa

Posted: Wed Oct 30, 2019 1:56 am
by bruceyoung01
Do we have NCL or python program to interpolate MPAS history data into specific geopotential height, e.g., 500 hPa? Or other relative easy way to interpolate it. I know we have some variables interpolated at specific geopotential height from diag files, e.g., temperature_500hPa. Thank you

Re: Plot moist potential temperature at 500 hPa

Posted: Thu Nov 07, 2019 5:23 pm
by mcurry
Here's an example Python script that interpolates qv to 500 hPA for all nCells. Reading in a history file, it uses Scipy's interp1d to create an interpolation function of Pressure and QV for all the grid levels of a single cell. From that function we can request our desired pressure level, and receive the corresponding QV value.

Since Theta, QV and Pressure are all dimensioned by Time, nCells, nVertLevels, you shouldn't have much difficulty converting this script for Theta (Or any other similarly dimensioned variable). However, feel free to let us know if you have a problem or a question.

interp_qv.py
Interpolate QV to 500 hPa
(1.28 KiB) Downloaded 38 times