Anonymous
Использование getvar для получения MSLP из файла wrfout дает ошибку
Сообщение
Anonymous » 15 ноя 2024, 23:43
Просто небольшой вопрос. Я пытаюсь извлечь ветры U и V мощностью 850 гПа вместе с MSLP из моего файла wrfout. Вот мой текущий код:
Код: Выделить всё
import numpy as np
import pandas as pd
import xarray as xr
from netCDF4 import Dataset
from wrf import getvar, interplevel, latlon_coords, ALL_TIMES
fileName = "wrfout1_d01_2024-11-14_12_00_00_5-day_fcst_rain"
wrf_file = Dataset(fileName)
# Extract the p, u, v, and mslp
p = getvar(wrf_file, "pressure", ALL_TIMES)
ua = getvar(wrf_file, "ua", ALL_TIMES, units="kt")
va = getvar(wrf_file, "va", ALL_TIMES, units="kt")
slp = getvar(wrf_file, "slp", ALL_TIMES, units="hPa")
# Interpolate u and v winds to 850 hPa
u_850 = interplevel(ua, p, 850)
v_850 = interplevel(va, p, 850)
# Get the lat/lon coordinates
lats, lons = latlon_coords(u_850)
Строка slp выдает ошибку KeyError: «QVAPOR».
Есть идеи?
Подробнее здесь:
https://stackoverflow.com/questions/791 ... s-an-error
1731703398
Anonymous
Просто небольшой вопрос. Я пытаюсь извлечь ветры U и V мощностью 850 гПа вместе с MSLP из моего файла wrfout. Вот мой текущий код: [code]import numpy as np import pandas as pd import xarray as xr from netCDF4 import Dataset from wrf import getvar, interplevel, latlon_coords, ALL_TIMES fileName = "wrfout1_d01_2024-11-14_12_00_00_5-day_fcst_rain" wrf_file = Dataset(fileName) # Extract the p, u, v, and mslp p = getvar(wrf_file, "pressure", ALL_TIMES) ua = getvar(wrf_file, "ua", ALL_TIMES, units="kt") va = getvar(wrf_file, "va", ALL_TIMES, units="kt") slp = getvar(wrf_file, "slp", ALL_TIMES, units="hPa") # Interpolate u and v winds to 850 hPa u_850 = interplevel(ua, p, 850) v_850 = interplevel(va, p, 850) # Get the lat/lon coordinates lats, lons = latlon_coords(u_850) [/code] Строка slp выдает ошибку KeyError: «QVAPOR». Есть идеи? Подробнее здесь: [url]https://stackoverflow.com/questions/79192052/using-getvar-to-obtain-mslp-from-a-wrfout-file-gives-an-error[/url]