WRF案例运行-Nested Model Runs

已于 2022-05-20 08:58:07 修改
阅读量1.8k 收藏 23
点赞数 4
分类专栏： WRF 文章标签： linux ubuntu
于 2022-04-27 17:36:07 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_40505953/article/details/124440905
版权

1、说明

   本次案例是利用WRF进行Nested Model Runs的嵌套模拟，下面的实例来自官网Nested Model Runs。如果还没安装和配置好WRF,可以先看我这篇博客WRF模型配置。如果之前没有做过单个域WRF模拟的同学，建议可以先试试单个域的WRF模拟WRF案例运行-Single Domain Case。

2、Nested Model Runs的基本流程

   WRF使用基本上会涉及到嵌套，因为嵌套可以在局地使用分辨率更高的数据来节约计算资源，同时达到提高局地模拟精度的目的。嵌套主要有两种方式：单向嵌套和双向嵌套。
   2.1 双向嵌套
   双向嵌套是一种在不同网格分辨率下同时运行多个域并相互通信的方法。粗域为嵌套域提供边界值，然后嵌套域将其计算反馈给粗域。(feedback=1)
   2.1.1 Two-way nested run, with one input file（双向嵌套，一个输入文件）
   本例的预处理步骤类似于单个域设置。当执行多个嵌套域的时候，唯一的区别发生在wrf.exe执行阶段。嵌套的子区域不需要另外准备初始场数据（比如Static Geographic Data和 Grided Meteorological Data），嵌套子区域的静态数据和气象数据都是由粗域插值而成。流程如下所示。
   优势：嵌套区域的模拟时间可以晚于粗域。
   劣势：嵌套区域的模拟精度可能降低，影响模拟质量。

#在namelist.input中
input_from_file=.true.,.false.,

在这里插入图片描述
   2.1.2 Two-way nested run, with two input files（双向嵌套，两个输入文件）
   在wrf.exe运行时，有两个输入文件wrfinput_do1和wrfinput_do2，这种嵌套方式也有两种方法，主要的区别在于：
   （1）嵌套区域都有静态地理数据和气象数据，都来自wrfinput_do2。推荐使用此方法
   （2）嵌套区域的静态地理数据来自wrfinput_do2，但是气象数据由粗域插值而成。当嵌套域的模拟时间晚于粗域的时候，可以使用该方法。
在这里插入图片描述
上面流程图中，A：WPS可以设置为嵌套域生成多个met_em.d02.*文件，但只需要初始时间。B：到目前为止，上述情况（1）和（2）的步骤是一样的。C：嵌套域将始终从粗略域获取其边界条件，因此不会创建wrfbdy_d02。需要生成嵌套域和粗域的geo_em.dox,粗区域的met_em.do1(所有时刻)，嵌套区域met_em.do2(初始时刻)，母区域和子区域的wrfinput_dox。
   2.2 单向嵌套
   单向嵌套是一种不同网格分辨率的多个域同时运行的方法，但唯一的通信来自粗域，它为嵌套区域提供边界值。但是嵌套域不会将计算结果返回给粗域。 此方法的运行方式与 双向嵌套相同，但需在 namelist.input 文件中设置 feedback=0。
   2.2.1 使用ndown.exe嵌套
   仅在已经有一个域已经运行了很长时间（几个月或几年）的时候，然后
决定使用更高分辨率来对嵌套区域进行研究，当觉得如果重新在来成本太高，就可以使用ndown.exe。

3、Nested Domain (2-way, 2 Input Files)

   本次实例使用双向嵌套，两个输入文件
   3.1 Geogrid
   （1）确保在WPS文件夹下

test@test:~/Build_WRF/WPS-4.1$

   （2）编辑namelist.wps

test@test:~/Build_WRF/WPS-4.1$ nano namelist.wps

对namelist.wps做以下改变：

&share
 wrf_core = 'ARW',
 max_dom = 2,
 start_date = '2016-10-06_00:00:00','2016-10-06_00:00:00',
 end_date   = '2016-10-08_00:00:00','2016-10-06_00:00:00',
 interval_seconds = 21600
 io_form_geogrid = 2,
/

&geogrid
 parent_id         =   1,1,
 parent_grid_ratio =   1,3,
 i_parent_start    =   1,23,
 j_parent_start    =   1,46,
 e_we              =  91,88,
 e_sn              =  100,79,
geog_data_res = 'default','default',
 dx = 27000,
 dy = 27000,
 map_proj = 'mercator',
 ref_lat   = 28,
 ref_lon   = -79,
 truelat1  = 30.0,
 truelat2  = 60.0,
 stand_lon = -79,
 geog_data_path = '/home/test/Build_WRF/WPS_GEOG/'
/
&ungrib
 out_format = 'WPS',
 prefix = 'FILE',
/

&metgrid
 fg_name = 'FILE'
 io_form_metgrid = 2,
/

如果你还不知道静态地理如何如下载，且不知道放在哪个文件夹，请看先看完这篇博客WRF案例运行-Single Domain Case，再进行下一步的操作。
   （3）在运行geogrib.exe之前，可以先确保这个域在正确的位置

test@test:~/Build_WRF/WPS-4.1$ ncl util/plotgrids_new.ncl

在这里插入图片描述
   （4）运行geogrib.exe

test@test:~/Build_WRF/WPS-4.1$ ./geogrid.exe

运行完后，看到Successful completion of geogrid，说明运行成功了，并且在WPS中生成以下两个文件。
在这里插入图片描述
   3.2 Geogrid
   （1）在WPS和WRF同一个文件夹下建立DATA文件夹

test@test:~/Build_WRF/WPS-4.1$ cd ..
test@test:~/Build_WRF$ mkdir DATA

在这里插入图片描述
   （2）下载Matthew case study data数据，放入刚刚创立的DATA文件夹中，并进行解压。

test@test:~/Build_WRF$ cd DATA/
test@test:~/Build_WRF/DATA$ tar -xf matthew_1deg.tar.gz

在这里插入图片描述
   （3）连接 GFS Vtable：

test@test:~/Build_WRF/DATA$ cd ..
test@test:~/Build_WRF$ cd WPS-4.1/
test@test:~/Build_WRF/WPS-4.1$ ln -sf ungrib/Variable_Tables/Vtable.GFS Vtable

   （4）通过link_grib.csh脚本来连接GRIB data

test@test:~/Build_WRF/WPS-4.1$ ./link_grib.csh ../DATA/matthew/fnl

在WPS文件夹中可以看到以下文件：
在这里插入图片描述
GRIBFILE.AAA -> …/DATA/matthew/fnl_20161006_00_00.grib2
GRIBFILE.AAB -> …/DATA/matthew/fnl_20161006_06_00.grib2
GRIBFILE.AAC -> …/DATA/matthew/fnl_20161006_12_00.grib2
GRIBFILE.AAD -> …/DATA/matthew/fnl_20161006_18_00.grib2
GRIBFILE.AAE -> …/DATA/matthew/fnl_20161007_00_00.grib2
GRIBFILE.AAF -> …/DATA/matthew/fnl_20161007_06_00.grib2
GRIBFILE.AAG -> …/DATA/matthew/fnl_20161007_12_00.grib2
GRIBFILE.AAH -> …/DATA/matthew/fnl_20161007_18_00.grib2
GRIBFILE.AAI -> …/DATA/matthew/fnl_20161008_00_00.grib2
   （5）运行ungrib.exe

test@test:~/Build_WRF/WPS-4.1$ ./ungrib.exe

在WPS文件加下产生以下文件：
在这里插入图片描述
   3.3 Metgrid
   将输入数据插值到模型域中

test@test:~/Build_WRF/WPS-4.1$ ./metgrid.exe

wps文件夹中将会产生以下文件：
met_em.d01.2016-10-06_00:00:00.nc
met_em.d01.2016-10-06_06:00:00.nc
met_em.d01.2016-10-06_12:00:00.nc
met_em.d01.2016-10-06_18:00:00.nc
met_em.d01.2016-10-07_00:00:00.nc
met_em.d01.2016-10-07_06:00:00.nc
met_em.d01.2016-10-07_12:00:00.nc
met_em.d01.2016-10-07_18:00:00.nc
met_em.d01.2016-10-08_00:00:00.nc
met_em.d02.2016-10-06_00:00:00.nc
在这里插入图片描述
   3.4 Real and WRF
   使用嵌套区域都有静态地理数据和气象数据都有数据的方法。
   （1）确保运行在WRF的test的em_real文件夹下

test@test:~/Build_WRF/WPS-4.1$ cd ..
test@test:~/Build_WRF$ cd WRF-4.1.2/
test@test:~/Build_WRF/WRF-4.1.2$ cd test
test@test:~/Build_WRF/WRF-4.1.2/test$ cd em_real
test@test:~/Build_WRF/WRF-4.1.2/test/em_real$

   （2）编辑namelist.input
   本次实例namelist.input中进行如下设置：

&time_control
 run_days                            = 0,
 run_hours                           = 48,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2016,2016,
 start_month                         = 10,10,
 start_day                           = 06,06,
 start_hour                          = 00,00,
 end_year                            = 2016,2016,
 end_month                           = 10,10,
 end_day                             = 08,08,
 end_hour                            = 00,00,
 interval_seconds                    = 21600
 input_from_file                     = .true.,.true.
 history_interval                    = 180,60
 frames_per_outfile                  = 1000,1000,
 restart                             = .false.,
 restart_interval                    = 5000,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 /

 &domains
 time_step                           = 150,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 2,
 e_we                                = 91,88,
 e_sn                                = 100,79
 e_vert                              = 35,35
 p_top_requested                     = 5000,
 num_metgrid_levels                  = 32,
 num_metgrid_soil_levels             = 4,
 dx                                  = 27000,9000
 dy                                  = 27000,9000
 grid_id                             = 1,2,
 parent_id                           = 1,1,
 i_parent_start                      = 1,23,
 j_parent_start                      = 1,46,
 parent_grid_ratio                   = 1,3,
 parent_time_step_ratio              = 1,3,
 feedback                            = 1,
 smooth_option                       = 0,
 sfcp_to_sfcp = .true.
 /

 &physics
 physics_suite                       = 'CONUS'
 mp_physics                          = 8,8,
 cu_physics                          = 3,3,
 ra_lw_physics                       = 4,4,
 ra_sw_physics                       = 4,4,
 bl_pbl_physics                      = 2,2,
 sf_sfclay_physics                   = 2,2,
 sf_surface_physics                  = 2,2,
 radt                                = 10,10,
 bldt                                = 0,0,
 cudt                                = 5,
 icloud                              = 1,
 num_land_cat                        = 21,
 sf_urban_physics                    = 0,0,
 /

 &fdda
 /

 &dynamics 
 w_damping                           = 0,
 diff_opt                            = 2,2,
 km_opt                              = 4,4,
 diff_6th_opt                        = 0,
 diff_6th_factor                     = 0.12,
 base_temp                           = 290.
 damp_opt                            = 3,
 zdamp                               = 5000.,5000.,
 dampcoef                            = 0.,0.,
 khdif                               = 0,0,
 kvdif                               = 0,0,
 non_hydrostatic                     = .true.,.true.,
 moist_adv_opt                       = 1,1,    
 scalar_adv_opt                      = 1,1,    
 /

 &bdy_control
 spec_bdy_width                      = 5,
 specified                           = .true.,.false.
 /

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /

   （3）连接由metgrid.exe创建的met_em文件

test@test:~/Build_WRF/WRF-4.1.2/test/em_real$ ln -sf ../../../WPS-4.1/met_em.d0*.2016-10* .

确保em_real中有以下文件：
在这里插入图片描述
   （4）运行real.exe

test@test:~/Build_WRF/WRF-4.1.2/test/em_real$ ./real.exe

得到以下文件：
wrfinput_d01
wrfbdy_d01
wrfinput_d02
   （5）运行wrf.exe

test@test:~/Build_WRF/WRF-4.1.2/test/em_real$ ./wrf.exe

可能是我的笔记本电脑的原因，大概运行了4个小时才出结果，得到以下文件：
wrfout_d01_2016-10-06_00:00:00
wrfout_d02_2016-10-06_00:00:00
在这里插入图片描述
   （6）查看输出结果

test@test:~/Build_WRF/WRF-4.1.2/test/em_real$ ncview wrfout_d02_2016-10-06_00:00:00

如下图，我选择的是3d vars=LH,总共模拟了48个小时，可以明显看出台风过境带来的变化。
在这里插入图片描述
结语：以上就是WRF官网嵌套实例的全部内容，如果有错误还请大家批评指正！

GISer凌
关注
  • 4
    点赞
  • 踩
  • 23
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
专栏目录
WRF-4.2.2-install-script-linux-64bit:这是一个脚本，它安装了所有库，软件，程序和地静数据，以运行“天气研究预报模型”（WRF-4.2.2）。 请分享和评论
03-31
这是一个脚本，用于安装所有库，软件，程序和地静数据，以运行Weather Research Forecast Model（WRF-4.2.2），并可以选择运行3DVAR和4DVAR观测数据。 安装 （确保将文件夹下载到您的主目录中）： git clone ...
处理WRF结果Process-WRF-master.zip
05-03
"Process_WRF-master.zip"文件很可能是包含WRF结果处理脚本和工具的项目，可能包含了从数据读取、转换、分析到可视化的全套流程。用户可以按照项目中的指导，结合自身需求，定制化处理WRF模型的输出结果。 总结，...
WRF示例运行（官方示例）
qq_41592222的博客
12-08 1356
【目的】该博客文用于备忘，有兴趣者可以照例实现，有疑问或需求在评论区留言即可；【内容】WRF的单层粗域实现；【参考】WRF实例运行系列（1）_钢筋火龙果的博客-CSDN博客_wrf运行实例【气象水文】超算服务器配置WRF环境_星空&逆行者的博客-CSDN博客_wrf安装环境变量配置【步骤】1、静态地理数据准备与模拟设置1.1 下载地理背景数据，下载地址见参考链接 1.2 进入wps，设置模拟域相关信息 1.3 运行geogrid.exe 实现地理背景数据插值 2、气象驱动数据提取2.1 示例
WRF实例运行系列（1）
Astronaut的博客
05-22 6169
## WRF实例运行1 ****单层嵌套模型**** step1：定义模型模拟域 下载静态地理数据并储存（以本人计算机实践为例） 下载地址：https://www2.mmm.ucar.edu/wrf/users/download/get_sources_wps_geog.html#specifi 储存地址： ​​​​​​​ ...
WRF模型的运行流程
最新发布
qq_44959558的博客
05-14 615
（1）。推荐大家使用WF Domain Wizard在线网站https://jiririchter.github.io/WRFDomainWizard/直接绘制自己的研究区域，且直接可以导出对应的namelist.wps文件（贼拉方便）。注意：设置parent_grid_ratio时建议大家使用官方推荐的值（3或5）。（2）。
WRF模式的运行（一）
weixin_52317507的博客
03-23 2313
WRF、中尺度数值模式
WRF模式运行及相关问题的解决
qq_46026699的博客
03-30 6726
一、稍微介绍一下WRF模式 WRF(Weather Research Forecast)模式系统是由许多美国研究部门及大学的科学家共同参与进行开发研究的新一代中尺度预报模式和同化系统。WRF模式系统的开发计划是在1997年由NCAR中小尺度气象处、NCEP的环境模拟中心、FSL的预报研究处和奥克拉荷马大学的风暴分析预报中心四部门联合发起建立的，并由国家自然科学基金和NOAA共同支持。现在，这项计划，得到了许多其他研究部门及大学的科学家共同参与进行开发研究。WRF模式系统具有可移植、易维护、可扩充、高效率、
WRF安装和运行的技术教程
热门推荐
ME1010的博客
04-02 1万+
想学习WRF？有这一篇文章就够了。 WRF和WPS安装保姆级教程，namelist中文入门教程，结合自己经验写成，非机翻，保证你能看懂。
超算入门笔记：大型机上如何运行WRF模式？一文总结（并行运算、Linux基础、作业调度、WRF运行）
weixin_43750300的博客
07-19 7229
接下来我们可以总结一下在超算中运行WRF，我们应当知道什么？首先是超算集群的基础构造与基本运行原理，让我们明白自己是在什么样的东西上运行程序。其次是Linux系统的命令，系统帮助我们让计算机实现各种功能。然后是一些在运行时我们需要了解到知识运行WRF时，我们需要安装压缩包与数据，这就需要我们学会远程登陆与数据传输。运行WRF时整个操作都是在Linux上进行，所以需要基础命令。配置WPS与WRF的namelist文件需要文本编辑，所以需要vi知识。httpshttpshttps。...
Python库 | wrf-python-1.1.3.tar.gz
04-19
资源分类：Python库 所属语言：Python 资源全名：wrf-python-1.1.3.tar.gz 资源来源：官方 安装方法：https://lanzao.blog.csdn.net/article/details/101784059
第三版WRF及WRF-Chem预处理以及后处理的Python脚本.zip
03-02
《基于Python的WRF与WRF-Chem预处理与后处理详解》 在现代气象学和环境科学研究中，Weather Research and Forecasting (WRF) 模型及其化学模块 WRF-Chem 是广泛应用的数值模拟系统。它能进行天气预报、气候变化分析...
wrf_python-1.3.2-cp310-cp310-win_amd64
08-19
wrf_python-1.3.2-cp310-cp310-win_amd64
wrf_python-1.3.2-cp37-cp37m-win_amd64
08-19
wrf_python-1.3.2-cp37-cp37m-win_amd64
WRF & WRF-chem 安装
07-11
WRF & WRF-Chem 安装 WRF（Weather Research and Forecasting）是一种气象预报模式，WRF-Chem 是其化学传输模块，用于研究大气污染物的传输和化学反应。以下是在 CentOS 6.5 系统中安装 WRF、WPS、WRF-Chem 及相应...
简单易用的WPS+WRF,easy-wrf-main.zip
05-03
在"easy-wrf-main.zip"这个压缩包中，我们找到了一个名为"easy-wrf-main"的文件夹，这通常包含了简化版的WPS和WRF配置，使得初学者或研究人员可以快速上手，避免复杂的初始设置过程。这个"easy-wrf-main"文件夹可能...
海浪模型的matlab代码-WRF-RCE-WTG:WRF-RCE-WTG
05-27
标题中的“海浪模型的matlab代码-WRF-RCE-WTG:WRF-RCE-WTG”指的是一个基于MATLAB实现的海浪模拟模型，该模型可能用于研究海洋环境中的波浪动力学。WRF（Weather Research and Forecasting）通常是一个气象学领域的...
WRF-Chem执行./real.exe时报错real: error opening wrfbdy for writing问题解决
qq_50207366的博客
01-04 408
WRFChem遇到错误及解决
WRF模式安装+详细运行教程
weixin_xiao5kou4chang6kai4的博客
12-13 4076
WRF模式在气象和相关学科领域有广泛的用途，但其安装对于许多人而言仍是一个难题
WRF实例运行（2）
Astronaut的博客
05-30 3556
## WRF实例运行1 ****单层嵌套模型**** step1：定义模型模拟域 下载静态地理数据并储存（以本人计算机实践为例） 下载地址：https://www2.mmm.ucar.edu/wrf/users/download/get_sources_wps_geog.html#specifi 储存地址： ...
python wrf_python-1.3.2-cp38
08-23
wrf-python-1.3.2-cp38是一个特定版本的Python软件包。WRF（Weather Research and Forecasting）是一种广泛用于大气科学和气象模拟研究的数值模式。而wrf-python是一个专门为WRF模型数据分析和可视化而设计的Python模块。 在WRF模式的输出数据中，包含了大量的气象参数和变量，如温度、湿度、风速、降水量等等。而wrf-python提供了一系列的工具和功能来处理和分析这些数据。它可以帮助用户读取WRF模式输出的NetCDF文件，并提供了各种方法来提取和操作其中的数据。例如，可以根据时间、空间坐标或特定变量来选择感兴趣的数据，计算统计信息，生成时间序列和空间图像等。 同时，wrf-python还提供了一系列用于可视化的工具。用户可以使用它绘制地图、等值线图、风向风速图和时序图等，以便更好地理解和展示WRF模式的模拟结果。这些功能使得wrf-python成为研究人员在大气科学领域中进行WRF模拟数据分析的重要工具之一。 wrf-python-1.3.2-cp38是适用于Python 3.8版本的wrf-python软件包的特定版本。这意味着它可以在Python 3.8的环境中运行，并提供对WRF模式数据处理和可视化的支持。如果你想在Python 3.8中使用wrf-python模块进行WRF数据分析，可以下载和安装这个特定版本，然后按照官方文档中的指南来使用它的功能。 总而言之，wrf-python-1.3.2-cp38是一个专门用于WRF模式数据分析和可视化的Python模块的特定版本，在大气科学和气象模拟研究中具有重要的作用。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值