WRFOUT 绘制台站探空图与简单分析
前言
实际应用中探空图可以分析所在区域的动热力特征，是预报员的好朋友
而在WRF应用中可以将其作为模式是否准确的检验工具
下面进行WRFOUT数据的探空图绘制
导入库
#
#库
from wrf import uvmet, to_np, getvar, interplevel, smooth2d, get_cartopy, cartopy_xlim, cartopy_ylim, latlon_coords,ll_to_xy
import numpy as np
from netCDF4 import Dataset
import matplotlib.pyplot as plt
import matplotlib as m
import metpy.calc as mpcalc
from metpy.plots import Hodograph, SkewT
from metpy.units import units
import metpy.calc as mpcalc
from metpy.cbook import get_test_data
from metpy.plots import add_metpy_logo, SkewT
from metpy.units import units
import matplotlib.pyplot as plt
from mpl_toolkits.axes_grid1.inset_locator import inset_axes
读数据
#读取WRF输出文件
wrfin = Dataset('/home/mw/input/wrfout3385/wrfout_d02_2022-07-14_0700.nc')
#指定要提取的经纬度坐标点
lat_lon =
#将经纬度坐标转换为模型坐标系（x, y）
x_y = ll_to_xy(wrfin, lat_lon, lat_lon)
#提取所需变量数据
p = getvar(wrfin, "pressure", timeidx=0)[:, x_y, x_y] * units.hPa
T = getvar(wrfin, "tc", timeidx=0)[:, x_y, x_y] * units.degC
Td = getvar(wrfin, "td", timeidx=0)[:, x_y, x_y] * units.degC
u = getvar(wrfin, "ua", timeidx=0)[:, x_y, x_y] * units('m/s')
v = getvar(wrfin, "va", timeidx=0)[:, x_y, x_y] * units('m/s')
h = getvar(wrfin, "height", timeidx=0)[:, x_y, x_y]
# 计算风向和风速
wind_dir = mpcalc.wind_direction(u, v)
wind_speed = mpcalc.wind_speed(u, v)
简单绘图
#创建9*9英寸，100dpi的画布
fig = plt.figure(figsize=(9, 9), dpi=100)
#在画布上添加SkewT对象并设置旋转角度为45度
skew = SkewT(fig, rotation=45)
#绘制温度和露点温度线
skew.plot(p, T, 'r')
skew.plot(p, Td, 'g')
#绘制风羽标志
skew.plot_barbs(p, u, v)
#设置y轴范围为1050hPa到100hPa，x轴范围为0到40度
skew.ax.set_ylim(1050, 100)
skew.ax.set_xlim(-30, 40)
#添加特殊线，如干绝热线、湿绝热线和混合比线
skew.plot_dry_adiabats()
skew.plot_moist_adiabats()
skew.plot_mixing_lines()
# 创建高空风速分析图
ax_hod = inset_axes(skew.ax, '30%', '30%', loc=1)
h = Hodograph(ax_hod, component_range=70.)
h.add_grid(increment=20)
h.plot_colormapped(u, v, wind_speed)# Plot a line colored by wind speed
#设置x轴和y轴标签
skew.ax.set_xlabel('温度 (°C)')
skew.ax.set_ylabel('压强 (hPa)')
Text(0, 0.5, '压强 (hPa)')
缝合绘图
层结曲线（红色），露点曲线（绿色），LCL(黑点)
CIN 蓝色阴影 CAPE 红色阴影
干绝热线（红色虚线）湿绝热线（蓝色虚线）等饱和混合比线（绿色虚线）
小结
1、要注意单位
2、简单分析：CAPE >CIN-- 不稳定 风向逆转 --冷平流 露点曲线：整层大气较干 层结曲线：温度垂直递减率大
番外 计算CAPE CIN
cape, cin = mpcalc.cape_cin(p, T, Td, prof, which_lfc='bottom', which_el='top')
print(cape, cin)
1360.9249456932523 joule / kilogram -104.38382141033901 joule /
完整文件与代码在此
文章来源于微信公众号：气ython风雨
页:
[1]