【python进行wrfout后处理基础绘图】安装wrf-python、读取wrfout、读取高度z变量场并插值到500hPa_z、输出nc文件、绘制等值线标值、副高区域填色

---

---

• 资料：台风“菲特“fitow模拟结果文件，https://blog.youkuaiyun.com/nice_clever/article/details/127340492#comments_24201637
• 必要python包：netCDF4、wrf-puthon

【anaconda安装wrf-python 】

conda install -c conda-forge wrf-python

本文主要介绍python对wrfout结果文件的初步后处理操作，以及基础绘图。
wrfout后处理包括：【读取wrfout文件、读取wrfout文件中变量metadata及数据、对高度场进行500hPa插值、输出nc文件】
基础绘图操作包括：【设置投影和范围、绘制等值线contour和等值线标值、副高区域填色contourf】
仅展示初步评估模拟的效果，若精美绘制需要进一步的设置、细化。

---

效果图

以下代码运行后输出plot()：

一、 python读取wrfout

import numpy as np
import xarray as xr
from netCDF4 import Dataset
from wrf import getvar, ALL_TIMES, interplevel
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as cfeature# ### 使用 cat 方法合并多个文件

wrfin = Dataset('./wrfout_d01_2013-10-05_00_00_00')
# print(wrfin)
# 提取位势高度和压力场
z = getvar(wrfin, 'z')  # 提取WRF netCDF 变量 # model height
p = getvar