Plot moist potential temperature at 500 hPa

Questions, issues, or suggestions related to post-processing of model output.
Post Reply
bruceyoung01
Posts: 56
Joined: Mon Sep 10, 2018 6:06 pm

Plot moist potential temperature at 500 hPa

Post by bruceyoung01 » 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

mcurry
Posts: 47
Joined: Mon Oct 29, 2018 5:33 pm
Location: Boulder, Co

Re: Plot moist potential temperature at 500 hPa

Post by mcurry » 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.

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

Post Reply

Return to “Post-processing”