收藏本站 劰载中...网站公告 | 吾爱海洋论坛交流QQ群:835383472

[其他] WRF系列教程5:WPS前处理的一些利器

[复制链接]

WRF模式模拟中,前处理WPS有许多有用的小工具,本篇文章列出几个常用的小工具,对一些特殊用途或者诊断错误或有帮助。

1 查看GRIB文件

ungrib.exe解压的是GRIB格式的数据,实际使用中,可能会遇到下载变量不全的情况,如果想要检验,可以使用以下工具打印grib文件的气象信息、变量、层次和日期,进行检验。

位于WPS/util文件夹:

g1print.exe  # grib1格式

g2print.exe   # grib2格式

' R6 ~0 T% `/ }  q; U
  • 9 Y0 N& s7 k( w5 a

    8 C9 J' W# l: n* ]& W

[C] 纯文本查看 复制代码
./g2print.exe gfs.t12z.pgrb2.0p50.f000


! |7 v1 ]7 ~: l& v: p& Q' l; s2 X4 K

8dab38651c84eb9d6cba99d54a130fc8.png

2 将气象数据转化为中间格式

metgrid.exe步骤,读取中间格式(Intermediate Format)的气象数据,一般由ungrib.exe步骤得到,而ungrib解压的数据为grib格式,并且需要提供Vtable。一些其他全球模式,比如耦合模式为NC格式,如果需要用这些模式数据驱动WRF,那么就需要将其转为中间格式,可以使用官方提供的Fortran代码(文末的参考链接)。需要根据具体应用做些修改并编译才能使用。

7ddfbccc0fe3d9bfc6d4478e8b00c187.jpeg

3 查看中间格式

ungrib.exe解压得到的中间格式或许不全,如果想查看,可以通过rd_intermediate.exe程序进行查看,会打印中间格式的文件信息。

% a/ P, V9 c! F5 g- C6 c, W

  • 1 X; F. u: ~8 o8 T5 ?- O: s

    8 [% p2 Q; j3 Y# I* b$ E

[C] 纯文本查看 复制代码
./rd_intermediate.exe FILE\:2021-12-28_12


$ s& _( _: m2 l$ K( Y) @/ y

e58cce76a886ea26c302ace394d450dd.png

注: 如果只看到最后SUCCESSFUL这句话,说明中间格式的文件是空的,查看其大小会显示为0,说明ungrib.exe步骤实际没成功(虽然可能也没报错)。

4 中间格式转为nc格式

int2nc.exe 可以将中间格式转为NC格式, ungrib.exe + int2nc.exe结合,可以将grib格式转为NC格式,如果有此需求,也算曲线救国路线。

+ L6 s7 G) @8 Z

  • ! _1 q; z) ~8 N) c' z2 |! m( _3 R

    ! u1 v7 ^' X1 W# s0 o

[C] 纯文本查看 复制代码
./int2nc.exe FILE:2021-12-28_12


  N8 {# I5 V% Y5 q, n# X; Q

得到FILE:2021-12-28_12.nc

5 绘制模拟区域

plotgrids_new.ncl脚本可以读取namelist.wps的设置,绘制设置的模拟区域,这一步骤可以检查设置模拟区域是否正确。

" v! o9 ~7 X) Y% O& ]5 C
  • ' E2 c1 y! O3 q: z8 c5 T. z
    , l1 A  b1 r. [( b$ Q

[C] 纯文本查看 复制代码
ncl  util/plotgrids_new.ncl


" {& u0 U0 |7 S$ Z( }

da9805877c08175d30060854692e0ce8.png

6 其他工具

还有一些其他小工具,虽然没用过,这里列举一下,读者有需要的可以自行探究。

avg_tsfc.exe

计算中间格式数据的日平均陆面温度,起止日期为namelist.wps中share部分设置的日期。

mod_levs.exe

去除intermediate format files格式数据中的某些层,使得数据层与namelist.wps中mod_levs层一致,mod_levs部分的press_pa 表示需要保留的气压层。

calc_ecmwf_p.exe

real.exe步骤需要气象场处于相同垂直层的3D气压场和位势高度场。calc_ecmwf_p.exe可进行垂直插值,将ECMWF的sigma层数据构建成需要的气象场。


0 q5 L) |3 w2 w1 Z% p0 t. v

参考链接:

游客,如果您要查看本帖隐藏内容请回复

2 B, v' k1 N6 V  A1 n
回复

举报 使用道具

相关帖子

全部回帖
暂无回帖,快来参与回复吧
懒得打字?点击右侧快捷回复 【吾爱海洋论坛发文有奖】
您需要登录后才可以回帖 登录 | 立即注册
细雨静聆
活跃在2024-12-26
快速回复 返回顶部 返回列表