使用wrfpython计算WRF输出的月度统计数据

Question

我使用wrf python来读取WRFOUT文件

wrflist = [Dataset("wrfout_d01_2012-12-01_00:00:00"),
           Dataset("wrfout_d01_2012-12-02_00:00:00"),
           Dataset("wrfout_d01_2012-12-03_00:00:00"),
           Dataset("wrfout_d01_2012-12-04_00:00:00"),
           Dataset("wrfout_d01_2012-12-05_00:00:00")]

p_join = getvar(wrflist, "P", timeidx=ALL_TIMES, method="join")

print(p_join)

我有包含一个月每小时输出的每日文件。如何计算特定变量的月平均值。与cdo monmean相当的东西

Answer 1

如果您不介意基于CDO的响应，我只需要计算每个文件的平均时间

for file in wrfout_d01_2012-12-??_00:00:00 ; do 
    cdo timmean $file ${file}_monmean.nc 
done

然后你可以通过使用

cdo ensmean ${file}_monmean.nc out.nc 

然而！这仅对WRF文件中的瞬时字段有效。对于通量，您需要使用总和而不是平均值，和/或进行一些单位转换。我倾向于使用cdo selvar提取所需的通量/非通量变量列表，然后只处理这些变量（还因为标准WRF输出包含通常不需要的ton）