文章目录
前言
为了更好的理解WRF运行过程中的各种参数设置,就找了一篇资料比较好找、分辨率没那么高的个例,自己来实际临摹一下这个模拟的过程。详细的记录了模拟的整个过程,包括WRFdomain设置模拟范围,WRF前处理、模拟流程,namelist.wps和namelist.input设置参数详情,ncdump查看wrfout结果的头文件,最后和原文对比500hPa形势可以看出模拟效果较好。(硕士论文比期刊文献描述的要详细一些)
参考文献:黄玉蓉. 2013年“菲特”台风的数值模拟研究[D].上海海洋大学,2015.
一、下载数据
NCEP再分析资料FNL资料:为模式提供三维大气初始条件和随时间变化的边界条件。https://rda.ucar.edu/datasets/ds083.2/index.html#!access.
观测、模拟跨度:2013年10月5日20时-7日14时
使用python下载FNL资料
登录后按所需时间选Grib1或Grib2(2007年以后)。点击web file listing和complete file list。选择需要的时次。
(也可以手动下载)点击python下载得到.py文件,在共享文件夹下运行,输入密码后可直接等待.
二、WPS前处理
0.安装WRFdomain
参考:https://blog.csdn.net/weixin_42487488/article/details/114035586
官网的下载网址: https://www.esrl.noaa.gov/gsd/wrfportal/DomainWizard.html
mkdir WRFDomainWizard/
cd WRFDomainWizard/
unzip WRFDomainWizard.zip
设置环境变量
chmod +x run_DomainWizard #使运行程序生效
run_DomainWizard #运行
######### set envirment for WRFDomainWizard #####
export PATH=/home/tian/WRF/WRFDomainWizard:$PATH
1.区域设置——namelist.wps
论文提供信息:
区域中心外侧网格(125°E,25°N),水平方向500×500,水平格距9km(先计算外层嵌套模拟)
内侧网格(118.74°E,22.3°N),水平方向721×721,水平格距3km
模拟垂直分层41层,层顶10hPa
外层网格时间步长90s,内侧网格时间步长30s
WRFDomain设置范围
namelist.wps【&geogrid】和namelist.input【&domains】在文件夹里面生成
namelist.wps参数详解:
原文如下:
&share
wrf_core = 'ARW',
max_dom = 2,
start_date = '2013-10-05_00:00:00','2013-10-05_00:00:00',
end_date = '2013-10-07_18:00:00','2013-10-07_18:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/
&geogrid
parent_id = 1, 1,
parent_grid_ratio = 1, 3,
i_parent_start = 1, 16,
j_parent_start = 1, 32,
e_we = 150, 217,
e_sn = 150, 217,
geog_data_res = 'default','default'