Page 1 of 1
Plot moist potential temperature at 500 hPa
Posted: Wed Oct 30, 2019 1:56 am
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
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.