使用MetPy计算WRF地转风与非地转风遇问题求助

一、警告原因说明 #

• 多时间坐标警告：WRF输出的元数据包含冗余时间字段，MetPy无法确定唯一时间坐标，属于WRF格式非CF标准导致的兼容性问题。
• 水平维度识别警告：WRF用south_north/west_east命名维度，和MetPy默认识别的y/x不匹配，因此MetPy只能默认假设维度顺序为[Y,X]。
• 弧度输入警告：WRF纬度存储为十进制度数，但MetPy计算科里奥利参数时，若坐标单位未正确标记，会误将度数当作弧度处理（度数范围0-90远大于π/2≈1.57弧度），触发警告。

二、非地转风计算错误原因 #

geo_wind_u/v出现无量纲单位，是因为MetPy未从WRF的z变量中正确提取带单位的纬度数据，导致科里奥利参数计算异常，最终地转风单位失效，无法和带m/s单位的ua/va做减法运算。

三、修正后的代码 #

import numpy as np
from netCDF4 import Dataset
import metpy.calc as mpcalc
from metpy.units import units
from metpy.xarray import assign_latitude_longitude
from wrf import getvar

# 打开WRF文件
filename = "wrfout_d01_2016-10-04_12:00:00"
ncfile = Dataset(filename)

# 提取变量并添加标准CF坐标，让MetPy识别维度
z = getvar(ncfile, "z", units="m").pipe(assign_latitude_longitude)
ua = getvar(ncfile, "ua", units="m s-1").pipe(assign_latitude_longitude)
va = getvar(ncfile, "va", units="m s-1").pipe(assign_latitude_longitude)

# 平滑位势高度
z = mpcalc.smooth_gaussian(z, 3)

# 显式指定水平维度，避免MetPy猜测；同时确保纬度单位正确（度转弧度由MetPy自动处理）
geo_wind_u, geo_wind_v = mpcalc.geostrophic_wind(z, latitude=z.lat)

# 计算非地转风
ageo_wind_u = ua - geo_wind_u
ageo_wind_v = va - geo_wind_v

关键改动说明 #

1. assign_latitude_longitude：给WRF变量添加标准CF风格的lat/lon坐标，让MetPy能正确识别水平维度和坐标信息，消除维度相关警告。
2. 显式传入latitude参数：强制geostrophic_wind使用正确的纬度数据（带单位），确保科里奥利参数计算正确，避免单位异常和弧度警告。
3. 保留单位传递：所有变量通过getvar指定单位后，MetPy的计算会自动维持单位一致性，解决无量纲和m/s的维度冲突。

内容的提问来源于stack exchange，提问作者n0rthern.dancer