WRF | 为什么wrfout中经向风和纬向风的格点数不一样

文摘   2024-10-13 16:05   北京  

WRF | 为什么wrfout中经向风和纬向风的格点数不一样

个人信息

公众号:气python风雨

Image Name

关注我获取更多学习资料,第一时间收到我的Python学习资料,也可获取我的联系方式沟通合作

温馨提示

由于可视化代码过长隐藏,可点击运行Fork查看
若没有成功加载可视化图,点击运行可以查看
ps:隐藏代码在【代码已被隐藏】所在行,点击所在行,可以看到该行的最右角,会出现个三角形,点击查看即可

项目目标

本项目旨在解答为什么wrfout中经向风和纬向风的格点数不一样

楔子

今日偷得浮生半日闲,突然一位读者私信问问题

他说在处理wrfout文件时碰到奇怪的现象

Image Name
Image Name

为什么U和V的格点会不一致,而且维度变量名字有区别

答案

真相只有一个,那就是荒川c网格

Image Name

WRF的网格分布如上,U和V并非是位于格点正中,他们会错开半个网格

south_north_stag或者west_east_stag则是交错网格的格点数

将网格平铺是这样的

Image Name

质量网格

  • • 压力、温度、湿度等与质量相关的量在网格单元中心计算,标记为“x”点。

  • • 这些“x”点统称为“质量网格”。

U网格

  • • 水平风的u分量在网格单元的左右侧中心计算,标记为绿色虚线。

  • • 这些点统称为“U网格”。

  • • U网格在y方向上的点数与质量网格相同,在x方向上多一个点。

V网格

  • • 水平风的v分量在网格单元的上下侧中心计算,标记为红色条。

  • • 这些点统称为“V网格”。

  • • V网格在x方向上的点数与质量网格相同,在y方向上多一个点。

交错网格

  • • 定义质量网格单元角的点由示意图中的黑点表示。

  • • 这些点统称为“交错网格”。

  • • WRF交错网格在质量网格的每个方向上都有一个额外的点。

如何处理

在python方面,我们经常用wrf-python的getvar函数直接获取插值后的ua和va

ua = getvar(ncfile, "ua")  
va = getvar(ncfile, "va")  

它们格点是一致的,直接使用即可

其他语言例如matlab则需要自行插值(博主对matlab了解不多,大家如有其他方法也可评论区说出)

参考

  1. 1. https://www.researchgate.net/figure/The-Arakawa-C-staggered-grid-used-in-WRF_fig2_258420398

  2. 2. https://amps-backup.ucar.edu/information/configuration/wrf_grid_structure.html

  3. 3. https://bbs.06climate.com/forum.php?mod=viewthread&tid=98890


气python风雨
主要发一些涉及大气科学的Python文章与个人学习备忘录
 最新文章