Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

Plot moist potential temperature at 500 hPa

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

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
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.

View attachment