【2025开箱即用】Cygwin WRF4.6.1 For Windows 绿色免安装版

kermit

WRF数值模式大多需要在高性能Linux集群上运行，但有时我们也想要在自己的个人电脑上学习WRF或者做一些简单的天气模拟。以往我们在Windows环境下运行WRF大多使用容器或者虚拟机或者子系统。但它们往往不够便携，还需要花时间配置各种环境，对于初学者的门槛也很高。

为了让WRF在Windows下运行更简便，Kermit Chen开发了 Piggy_Packages ，它基于Cygwin 深度定制，解压即用。它不仅可以运行WRF，还可以运行NCL、CDO、NCO、Ncview、Panoply、MET、UPP、ESMF、RTTOV等各种Linux下好用的气象软件。

——————系统需求——————

64位版本的Windows 10 (22H2) 、Windows 11

ARM版本的Windows 11；

请注意，随着Windows 10支持周期的结束，对于 Windows 10 的支持将很快被移除。

如果您目前使用的是Windows 10，安装Windows Terminal将会获得更佳的使用体验，请在此处下载：https://aka.ms/terminal

Windows 11 已经内置了Windows Terminal，您无需额外安装。

～～～～～～～～～～～～～～～～～～～～～～～
下载地址【没有气象家园账号的同学，可以加入Piggy气象科学交流QQ群852983715获取下载地址喔】

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

～～～～～～～～～～～～～～～～～～～～～～～
简要教程～～
下面以运行WPS为例带大家体验一下Cygwin WRF4.6.1
首先下载并解压缩 Piggy_Packages 然后右键管理员身份运行Piggy-Bash.sh 或者 Piggy-Zsh.sh
Piggy-Bash.sh 提供一种简约经典的交互方式；而Piggy-Zsh.sh 提供一种现代炫酷的交互方式。
（Win10 用户使用Piggy-Zsh.sh需要安装Windows Terminal，下载地址见上文，否则会很丑）
首次启动时会进行环境部署，需要几分钟时间。

Piggy_Packages Version 2025.1 的WRF数值模式安装在/pp_model/WRF-4.6.1

现在我们切换到切换到该目录：

1. cd /pp_model/WRF-4.6.1

复制代码

目录结构如下：

1. /pp_model/WRF-4.6.1
   
2. ├── download_geog.sh # 高精度地形数据下载脚本
   
3. ├── RTTOV # 快速辐射传输模式（v12.1），用于WRFDA卫星资料同化
   
4. ├── WPS # WRF前处理系统
   
5. ├── WPS_GEOG # 静态地形数据（低精度）
   
6. ├── WRF -> WRF-REAL # 指向 WRF-REAL 的软链接
   
7. ├── WRF-Chem # WRF化学模式
   
8. ├── WRFDA # WRF同化系统
   
9. ├── WRF-FIRE # WRF火灾模式
   
10. ├── WRF-Hydro # WRF水文模式
    
11. ├── WRFPLUS # WRF切线性和伴随模式，用于敏感性试验和WRFDA四维变分
    
12. ├── WRF-REAL # WRF 中尺度大气模式
    
13. └── WRF-VF # WRF 中尺度大气模式（启用vortex following）

复制代码

请注意，WRF目录下的RTTOV为WRFDA兼容版本，我们在/pp_module/module/gcc-15.1.0/RTTOV-14.0目录下提供了一个更新版本的RTTOV。
1）地形数据下载：
Piggy_Packages 已经内置了低精度地形数据，对于学习使用已经足够，但是如果您想要更好的模拟结果，您可以通过download_geog.sh脚本下载高精度（解压后约29GB）和可选的完整地形数据（解压后约200GB）。

1. ./download_geog.sh

复制代码

2）WPS和WRF运行
让我们返回HOME目录

1. cd ~

复制代码

然后将WPS和WRF拷贝过来：

1. cp -r /pp_model/WRF-4.6.1/WPS ~/
   
2. cp -r /pp_model/WRF-4.6.1/WRF-REAL ~/WRF

复制代码

进入WPS目录：

1. cd WPS

复制代码

用记事本或VS Code打开namelist.wps: (请注意，虽然它的后缀是wps，但是请勿使用金山WPS Office 或者Microsoft Word打开)

1. open -e namelist.wps

复制代码

提示：如果您安装了VS Code，并在安装时选择了将VS Code添加到PATH，上述命令将默认使用VS Code 打开；否则将使用记事本打开。

清空文件中原有内容，将以下内容粘贴进去：

1. &share
   
2. wrf_core = 'ARW',
   
3. max_dom = 3,
   
4. start_date = '2025-05-01_00:00:00','2025-05-01_00:00:00','2025-05-01_00:00:00'
   
5. end_date = '2025-05-01_18:00:00','2025-05-01_18:00:00','2025-05-01_18:00:00'
   
6. interval_seconds = 10800,
   
7. io_form_geogrid = 2,
   
8. /
   
9. 
   
10. &geogrid
    
11. parent_id = 1, 1, 2,
    
12. parent_grid_ratio = 1, 3, 3,
    
13. i_parent_start = 1, 40, 17,
    
14. j_parent_start = 1, 26, 21,
    
15. e_we = 81, 49, 55,
    
16. e_sn = 57, 49, 55,
    
17. geog_data_res = 'lowres','lowres','lowres',
    
18. dx = 81000,
    
19. dy = 81000,
    
20. map_proj = 'lambert',
    
21. ref_lat = 35,
    
22. ref_lon = 110,
    
23. truelat1 = 30,
    
24. truelat2 = 60,
    
25. stand_lon = 116.397,
    
26. geog_data_path = '/pp_model/WRF-4.6.1/WPS_GEOG/'
    
27. /
    
28. 
    
29. &ungrib
    
30. out_format = 'WPS',
    
31. prefix = 'FILE',
    
32. /
    
33. 
    
34. &metgrid
    
35. fg_name = 'FILE'
    
36. io_form_metgrid = 2,
    
37. /
    

复制代码

【可选】如果您已经下载了高精度地形数据，您可以将：

1. geog_data_res = 'lowres','lowres','lowres',

复制代码

更改为：

1. geog_data_res = 'highres','highres','highres',

复制代码

更改完之后，您可以输入以下命令可视化WRF模拟区域：

1. ncl util/plotgrids_new.ncl

复制代码

(点击图像以关闭图窗)

接下来我们运行geogrid.exe将地形数据插值到模式网格上：

1. mpirun -np $NUMBER_OF_PROCESSORS ./geogrid.exe

复制代码

$NUMBER_OF_PROCESSORS 为您的CPU核心数量，也可以任意指定不高于您CPU核心数量的值（下同），例如：

1. mpirun -np 4 ./geogrid.exe

复制代码

下面我们用ungrib.exe解码和提取背景场资料。首先从 https://pan.baidu.com/s/1uiSDnEcsSwERz2l6sEGKOQ?pwd=1234 下载一些FNL背景场文件；
然后输入一下命令在资源管理器中打开HOME目录：

1. open ~

复制代码

打开之后，将下载的FNL背景场文件拷贝进去。

使用链接脚本将背景场文件链接到WPS所在目录：

1. ./link_grib.csh ~/fnl_*.grib2

复制代码

由于我们使用的FNL背景场，它和GFS同源，因此我们需要适用于GFS的Vtable，将它链接到当前目录下。Vtable的作用在于告诉ungrib.exe需要提取哪些变量，以及这些变量的GRIB编码信息。

1. ln -s ungrib/Variable_Tables/Vtable.GFS Vtable

复制代码

然后我们运行ungrib.exe进行背景场数据的解码和提取。

1. ./ungrib.exe

复制代码

请注意，ungrib.exe无法使用mpirun并行运行。
运行完成后将会生成FILE:开头的文件。

metgrid.exe的作用是将ungrib.exe提取出的背景场信息水平插值到模式网格上，并且与geogrid.exe插值的地形数据结合，生成模式的输入数据。

1. mpirun -np $NUMBER_OF_PROCESSORS ./metgrid.exe

复制代码

运行成功后，将会生成met_em.开头的文件。

让我们切换到WRF运行目录：

1. cd ~/WRF/run

复制代码

将metgrid.exe的输出数据链接到此处：

1. ln -s ~/WPS/met_em.* ./

复制代码

然后我们编辑WRF的配置文件namelist.input：

1. open -e namelist.input

复制代码

清除文件中的所有内容，然后将下方内容粘贴进去：

1. &time_control
   
2. run_days                            = 0,
   
3. run_hours                           = 18,
   
4. run_minutes                         = 0,
   
5. run_seconds                         = 0,
   
6. start_year                          = 2025, 2025, 2025,
   
7. start_month                         = 05,   05,   05,
   
8. start_day                           = 01,   01,   01,
   
9. start_hour                          = 00,   00,   00,
   
10. end_year                            = 2025, 2025, 2025,
    
11. end_month                           = 05,   05,   05,
    
12. end_day                             = 01,   01,   01,
    
13. end_hour                            = 18,   18,   18,
    
14. interval_seconds                    = 10800
    
15. input_from_file                     = .true.,.true.,.true.,
    
16. history_interval                    = 60,  60,  60,
    
17. frames_per_outfile                  = 1, 1, 1,
    
18. restart                             = .false.,
    
19. restart_interval                    = 7200,
    
20. io_form_history                     = 2
    
21. io_form_restart                     = 2
    
22. io_form_input                       = 2
    
23. io_form_boundary                    = 2
    
24. /
    
25. 
    
26. &domains
    
27. time_step                           = 90,
    
28. time_step_fract_num                 = 0,
    
29. time_step_fract_den                 = 1,
    
30. max_dom                             = 3,
    
31. e_we                                = 81, 49, 55,
    
32. e_sn                                = 57, 49, 55,
    
33. e_vert                              = 45, 45, 45,
    
34. dzstretch_s                         = 1.1
    
35. p_top_requested                     = 5000,
    
36. num_metgrid_levels                  = 34,
    
37. num_metgrid_soil_levels             = 4,
    
38. dx                                  = 81000,
    
39. dy                                  = 81000,
    
40. grid_id                             = 1,     2,     3,
    
41. parent_id                           = 1, 1, 2,
    
42. i_parent_start                      = 1, 40, 17,
    
43. j_parent_start                      = 1, 26, 21,
    
44. parent_grid_ratio                   = 1, 3, 3,
    
45. parent_time_step_ratio              = 1, 1, 1,
    
46. feedback                            = 1,
    
47. smooth_option                       = 0
    
48. /
    
49. 
    
50. &physics
    
51. physics_suite                       = 'CONUS'
    
52. mp_physics                          = -1,    -1,
    
53. cu_physics                          = -1,    -1,
    
54. ra_lw_physics                       = -1,    -1,
    
55. ra_sw_physics                       = -1,    -1,
    
56. bl_pbl_physics                      = -1,    -1,
    
57. sf_sfclay_physics                   = -1,    -1,
    
58. sf_surface_physics                  = -1,    -1,
    
59. radt                                = 15,    15,
    
60. bldt                                = 0,     0,
    
61. cudt                                = 0,     0,
    
62. icloud                              = 1,
    
63. num_land_cat                        = 21,
    
64. sf_urban_physics                    = 0,     0,
    
65. fractional_seaice                   = 1,
    
66. /
    
67. 
    
68. &fdda
    
69. /
    
70. 
    
71. &dynamics
    
72. hybrid_opt                          = 2,
    
73. w_damping                           = 0,
    
74. diff_opt                            = 2,      2,
    
75. km_opt                              = 4,      4,
    
76. diff_6th_opt                        = 0,      0,
    
77. diff_6th_factor                     = 0.12,   0.12,
    
78. base_temp                           = 290.
    
79. damp_opt                            = 3,
    
80. zdamp                               = 5000.,  5000.,
    
81. dampcoef                            = 0.2,    0.2,
    
82. khdif                               = 0,      0,
    
83. kvdif                               = 0,      0,
    
84. non_hydrostatic                     = .true., .true.,
    
85. moist_adv_opt                       = 1,      1,
    
86. scalar_adv_opt                      = 1,      1,
    
87. gwd_opt                             = 1,      0,
    
88. /
    
89. 
    
90. &bdy_control
    
91. spec_bdy_width                      = 5,
    
92. specified                           = .true.
    
93. /
    
94. 
    
95. &grib2
    
96. /
    
97. 
    
98. &namelist_quilt
    
99. nio_tasks_per_group = 0,
    
100. nio_groups = 1,
     
101. /
     

复制代码

real.exe用于将metgrid.exe的输出在垂直方向上插值到模式网格，生成真实的三维大气结构场。

1. mpirun -np $NUMBER_OF_PROCESSORS ./real.exe

复制代码

提示：由于我们的d03区域很小，使用的核心数建议不超过16，否则会出现区域不够分割的错误。
real.exe运行后将会生成WRF的初始场wrfinput_*和边界场wrfbdy_*文件。

下面正式开始运行WRF数值模式：

1. mpirun -np $NUMBER_OF_PROCESSORS ./wrf.exe

复制代码

WRF数值模式运行是一个非常耗时的工作，运行上述命令之后，您可以开启一个新的Piggy_Packages 终端，显示WRF运行进度。
开启新终端之后，输入：

1. tail -f ~/WRF/run/rsl.out.0000

复制代码

WRF运行结束后，将会生成wrfout开头的结果文件。

3）WRF运行结果的可视化
wrfout开头的结果文件为NetCDF格式，因此，您可以使用Ncview或者Panoply来可视化这些文件。

1. ncview wrfout_d01_2025-05-01_18:00:00

复制代码

1. Panoply wrfout_d01_2025-05-01_18:00:00

复制代码

帖子推荐：【原创资源】全国超高精度海拔高程数据
http://bbs.06climate.com/forum.php?mod=viewthread&tid=104420&fromuid=113111
如想了解 Piggy_Packages 的更多功能，请查看这篇帖子：
Piggy_Packages 气象科研常用软件集合
http://bbs.06climate.com/forum.php?mod=viewthread&tid=104814&fromuid=113111

——————技术支持——————
~欢迎加入Piggy气象科学交流QQ群852983715学习交流~

                               
登录/注册后可看大图

kermit

本教程旨在帮助初学者运行并理解WRF数据前处理、数值预报、资料同化的完整流程，所描述的运行流程不是最简单的，但能够尽量避免遇到难以处理的问题。建议初学者按步骤操作，运行成功之后再尝试自己或者别的教程的方法。

1.    Piggy Packages气象工具包

欢迎使用 Piggy_Packages！

Piggy_Packages 是一个专为气象科研人员和爱好者打造的、高度集成的软件包合集。它基于 Cygwin 进行了深度定制和优化，旨在 Windows 操作系统环境下提供一个“解压即用”的便捷科研平台。本软件包致力于消除在 Windows 系统中安装和配置复杂气象软件的障碍，让用户能够快速投入到实际的科研工作中。

Piggy_Packages 目前支持 Windows 10 （22H2）和 Windows 11。请注意，随着Windows 10支持周期的结束，对于 Windows 10 的支持将很快被移除。

如果您目前使用的是Windows 10，安装Windows Terminal将会获得更佳的使用体验，请在此处下载：https://aka.ms/terminal
Windows 11 已经内置了Windows Terminal，您无需额外安装。

本文所用到的全部数据下载地址：http://nui.st/4fqm

2.    WPS数据前处理系统的运行

假设现在是2022年4月1日06时（世界时），我们想要预报如图所示的区域在接下来24小时内的逐小时天气要素变化。

由于现在刚刚到2022年4月1日06时，虽然我们可以获取到这个时刻全球不同站点的气象观测，但是各种地面、高空、卫星、雷达以及微波辐射计、飞机报等观测资料时空维度不一，有的区域观测资料密集而有的区域观测资料稀疏，质量也参差不齐，无法直接处理成均匀的格点数据。我们的预报需要这个时刻的气象格点资料，最简便的方法是从全球模式的预报场获得。

为了完成24小时的逐小时天气要素预报，我们已获取了全球模式于2022年4月1日00时起报的6、9、12、15、18、21、34、27、30小时预报场。其中6小时预报场，即6小时前对当前时刻的全球大气预报，作为本次预报的初始场，9-30小时预报场作为预报过程中的边界条件（因为区域模式不关注区域外发生了什么，但是区域外的大尺度天气系统又足以对区域内的天气造成影响，因此区域模式预报时需要全球模式的预报场为其提供边界条件）。

首先将 WPS程序包拷贝到当前目录：

1. cp -r /pp_model/WRF-4.6.1/WPS ./

复制代码

进入WPS目录：

1. cd WPS

复制代码

打开当前目录：

1. open

复制代码

将刚刚下载的GFS全球预报场资料拷贝至此。

在运行WPS之前，需要编辑namelist.wps 文件，该文件记录了WPS运行时的各种参数。

打开namelist.wps ：

1. open -e namelist.wps

复制代码

清空所有内容，并且将以下内容粘贴进去：

1. &share
   
2. wrf_core = 'ARW',
   
3. max_dom = 3,
   
4. start_date = '2022-04-01_06:00:00','2022-04-01_06:00:00','2022-04-01_06:00:00'
   
5. end_date = '2022-04-02_06:00:00','2022-04-02_06:00:00','2022-04-02_06:00:00'
   
6. interval_seconds = 10800,
   
7. io_form_geogrid = 2,
   
8. /
   
9. 
   
10. &geogrid
    
11. parent_id = 1, 1, 2,
    
12. parent_grid_ratio = 1, 3, 3,
    
13. i_parent_start = 1, 40, 17,
    
14. j_parent_start = 1, 26, 21,
    
15. e_we = 81, 49, 55,
    
16. e_sn = 57, 49, 55,
    
17. geog_data_res = 'lowres','lowres','lowres',
    
18. dx = 81000,
    
19. dy = 81000,
    
20. map_proj = 'lambert',
    
21. ref_lat = 35,
    
22. ref_lon = 110,
    
23. truelat1 = 30,
    
24. truelat2 = 60,
    
25. stand_lon = 116.397,
    
26. geog_data_path = '/pp_model/WRF-4.6.1/WPS_GEOG/'
    
27. /
    
28. 
    
29. &ungrib
    
30. out_format = 'WPS',
    
31. prefix = 'FILE',
    
32. /
    
33. 
    
34. &metgrid
    
35. fg_name = 'FILE'
    
36. io_form_metgrid = 2,
    
37. /
    

复制代码

然后生成模拟区域的地形数据，运行geogrid.exe

1. ./geogrid.exe

复制代码

成功之后将生成geo_em.d01.nc、geo_em.d02.nc、geo_em.d03.nc文件。

接下来从GRIB格式的GFS全球预报场数据中提取模式需要的气象要素场：

1. ./link_grib.csh gfs.0p25.2022*
   
2. ln -sf ungrib/Variable_Tables/Vtable.GFS Vtable
   
3. ./ungrib.exe

复制代码

成功之后将生成FILE:2022-04-01_06、FILE:2022-04-01_09等文件

接下来将气象要素场插值到地形格点上：

1. ./metgrid.exe

复制代码

成功之后将生成met_em.d01.2022-04-01_06:00:00.nc等文件

至此整个数据前处理流程运行完毕。

3.    REAL生成模式的初始场和边界

回到Home目录：

1. cd

复制代码

将 WRF程序包拷贝到当前目录：

1. cp -r /pp_model/WRF-4.6.1/WRF-REAL ~/WRF

复制代码

切换到WRF运行目录：

1. cd WRF/run

复制代码

real.exe的作用是将metgrid.exe生成的气象要素场处理为模式的初始场和边界条件，其运行参数储存在namelist.input 文件中。初学者可以直接使用我们提供的namelist.input文件。

打开当前目录：

1. open

复制代码

将我们提供的namelist.input文件拷贝至此，覆盖原文件。

将metgrid.exe生成的met_em文件拷贝（对于初学者来说比较保险）到当前目录：

1. cp ../../WPS/met_em.* .

复制代码

运行real.exe

1. ./real.exe

复制代码

运行完成后生成三层嵌套网格的初始场wrfinput_d01、wrfinput_d02、wrfinput_d03和最外层网格边界条件wrfbdy_d01 （内层网格的边界条件由外层网格提供，因此只有最外层网格需要全球模式为其提供边界条件）

4.    数值模式WRF的运行

现在到了激动人心的时刻了，让我们来运行数值模式WRF、占卜未来天气！

1. ./wrf.exe

复制代码

或者，如果你的电脑CPU和内存充足，你可以多核运行WRF：

1. mpirun -np $NUMBER_OF_PROCESSORS ./wrf.exe #[使用全部CPU核心]
   
2. mpirun -np 2 ./wrf.exe #[使用2个CPU核心]

复制代码

等待其运行完成。

如果在等待过程中想要查看进度，可以另外开启一个终端窗口，输入：

1. tail -f ~/WRF/run/rsl.out.0000

复制代码

运行完成后，将能看到每个小时的wrfout文件。

5.    绘图查看预报结果

您可以使用NCVIEW或者PANOPLY打开wrfout文件，直接查看各个变量绘图，或者使用NCL官网提供的脚本绘图。

1. ncview wrfout_d01_2022-04-01_07:00:00

复制代码

6.    WRFDA资料同化

我们知道，模式预报效果很大程度取决于初始场的准确度。但是还记得我们给模式提供的初始场是从哪里来的吗？

是的，它来自于6小时之前全球模式的预报结果。既然是预报，就一定会存在误差！那么我们要怎么修正预报的误差呢？这时，我们就要用到当前时刻的观测资料。

我们知道，全球模式预报场尽管存在误差，但毕竟是整齐的网格化数据，能够基本反映全球大气真实状况；而观测资料数量众多，但空间分布疏密不均，质量参差不齐，由于设备故障或算法问题导致的数据错误或不准确时常发生。我们通过资料同化的方法充分利用优质的观测资料、按权重采信准确度一般但又有一定利用价值的观测资料，以修正模式初始场；而明显不可信的观测资料一般直接剔除。相关理论知识在此不予展开，请感兴趣的同学自行学习。

回到Home目录：

1. cd

复制代码

新建WRFDA运行目录：

1. mkdir -p WRFDA/d01 WRFDA/d02 WRFDA/d03 WRFDA/obs

复制代码

然后将所需的namelist.input拷贝到对应目录，d01还需要拷贝更新边界条件的配置文件parame.in，观测资料拷贝到obs目录。

依次对三层嵌套的初始场进行同化：

1. #d01:
   
2. cd d01
   
3. cp ~/WRF/run/wrfinput_d01 fg #背景场
   
4. cp ../obs/prepbufr.gdas.20220401.t06z.nr.48h ob.bufr #常规观测
   
5. cp ../obs/gdas.1bamua.t06z.20220401.bufr amsua.bufr #卫星AMSUA观测
   
6. cp /pp_model/WRF-4.6.1/WRFDA/var/run/be.dat.cv3 be.dat #背景误差协方差矩阵
   
7. cp /pp_model/WRF-4.6.1/WRFDA/run/LANDUSE.TBL .
   
8. # cp /pp_model/WRF-4.6.1/WRFDA/var/run/VarBC.in .  ##如果导致WRFDA卡死，请删除VarBC.in文件
   
9. cp /pp_model/WRF-4.6.1/WRFDA/var/run/radiance_info .
   
10. ln -sf /pp_model/WRF-4.6.1/RTTOV/rtcoef_rttov12/rttov7pred54L rttov_coeffs
    
11. /pp_model/WRF-4.6.1/WRFDA/var/build/da_wrfvar.exe
    
12. #d01同化之后需要更新边界条件
    
13. cp ~/WRF/run/wrfbdy_d01 .
    
14. /pp_model/WRF-4.6.1/WRFDA/var/build/da_update_bc.exe
    
15. cd ..
    
16. #d02:
    
17. cd d02
    
18. cp ~/WRF/run/wrfinput_d02 fg #背景场
    
19. cp ../obs/prepbufr.gdas.20220401.t06z.nr.48h ob.bufr #常规观测
    
20. cp ../obs/gdas.1bamua.t06z.20220401.bufr amsua.bufr #卫星AMSUA观测
    
21. cp /pp_model/WRF-4.6.1/WRFDA/var/run/be.dat.cv3 be.dat #背景误差协方差矩阵
    
22. cp /pp_model/WRF-4.6.1/WRFDA/run/LANDUSE.TBL .
    
23. cp /pp_model/WRF-4.6.1/WRFDA/var/run/VarBC.in .
    
24. cp /pp_model/WRF-4.6.1/WRFDA/var/run/radiance_info .
    
25. ln -sf /pp_model/WRF-4.6.1/RTTOV/rtcoef_rttov12/rttov7pred54L rttov_coeffs
    
26. /pp_model/WRF-4.6.1/WRFDA/var/build/da_wrfvar.exe
    
27. cd ..
    
28. #d03:
    
29. cd d03
    
30. cp ~/WRF/run/wrfinput_d03 fg #背景场
    
31. cp ../obs/prepbufr.gdas.20220401.t06z.nr.48h ob.bufr #常规观测
    
32. cp ../obs/gdas.1bamua.t06z.20220401.bufr amsua.bufr #卫星AMSUA观测
    
33. cp /pp_model/WRF-4.6.1/WRFDA/var/run/be.dat.cv3 be.dat #背景误差协方差矩阵
    
34. cp /pp_model/WRF-4.6.1/WRFDA/run/LANDUSE.TBL .
    
35. cp /pp_model/WRF-4.6.1/WRFDA/var/run/VarBC.in .
    
36. cp /pp_model/WRF-4.6.1/WRFDA/var/run/radiance_info .
    
37. ln -sf /pp_model/WRF-4.6.1/RTTOV/rtcoef_rttov12/rttov7pred54L rttov_coeffs
    
38. /pp_model/WRF-4.6.1/WRFDA/var/build/da_wrfvar.exe
    
39. cd ..

复制代码

至此WRFDA资料同化的全部流程已经完成。在科研工作中，我们还需要对同化结果进行细致的诊断分析，以评估资料和同化算法的质量。感兴趣的同学可以自行学习。

7.    使用同化后的初始场进行数值预报

我们通过WRFDA同化系统，得到了更加准确的模式初始场，我们将这些初始场拷贝到WRF运行路径下，即可进行数值预报试验。

1. cd ~/WRF/run
   
2. #删除之前的初始场和预报结果
   
3. rm -rf wrfbdy_d01 wrfinput_d0* wrfout_*
   
4. #拷贝同化后的初始场和边界条件
   
5. cp ~/WRFDA/d01/wrfbdy_d01 .
   
6. cp ~/WRFDA/d01/wrfvar_output wrfinput_d01
   
7. cp ~/WRFDA/d02/wrfvar_output wrfinput_d02
   
8. cp ~/WRFDA/d03/wrfvar_output wrfinput_d03
   
9. #运行WRF
   
10. mpirun -np $NUMBER_OF_PROCESSORS ./wrf.exe

复制代码

kermit

由于版本更新，本帖已被移除！

792161805

------ ERROR while reading namelist dynamics ------
Maybe here?:      scalar_adv_opt                      = 1,      1,
Maybe here?:      gwd_opt                             = 1,      0,
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:   10851
ERRORS while reading one or more namelists from namelist.input.
大佬，报错显示这个是什么原因啊

清宵逝水

666新手拿来学习太方便了

kermit

sunxupeng 发表于 2022-5-5 17:25
您好，我已经查询到了此网址，但是我发现有一个问题不知道是不是我找寻数据错误的问题（我找的网址https: ...

rda.ucar.edu一般用来下载历史数据 最新数据还是GFS官网比较好

kermit

sunxupeng 发表于 2022-5-5 09:01
您好，是否是这个网站呢？
https://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/

可以的
不过我们一般用这个网站 https://rda.ucar.edu

kermit

sunxupeng 发表于 2022-5-4 17:41
请问预报场数据是从哪里获取的？

GFS官网或者RDA

Ronin.

赞赞大佬，非常厉害

Running

谢谢楼主，非常游泳，期待已久。

挚zhi

666 来学习一下！！

单纯xiao可爱

小炜真的好厉害！！

灭火器

可怕，能不能稍微透露下怎么开发的

AllenLu

不知道用起来怎么样，期待教程，感谢楼主分享

最好的明天

好厉害啊，感觉会很有用

2223xia

我又来了，学习下楼主大大分享的好东西

si18

期待之后教程 感激不尽

topmad

我虚拟机的序列号要保不住了

Ifare

强，点赞