使用 GIS4WRF 模拟 2018 年欧洲热浪(以阿姆斯特丹为中心)
本博客结合GIS4WRF官网提供教程,总结使用 GIS4WRF 模拟 2018 年欧洲热浪(以阿姆斯特丹为中心)的详细步骤。
案例介绍
GIS4WRF官网提供教程-The 2018 European heat wave
本教程演示如何使用 GIS4WRF 插件在 QGIS 中快速配置并运行一个真实的天气模拟案例,模拟时间为:
- 2018年7月15日 12:00 至 18:00(6小时)
- 中心位置:荷兰阿姆斯特丹
预计总耗时:20 分钟以内(包括模拟)
模拟运行时间:< 5 分钟(在普通双核/8GB 内存电脑上)
✅ 开始前的准备工作清单:
- 已正确 安装并配置 GIS4WRF 插件
QGIS和GIS4WRF 插件的安装可参见另一博客-【QGIS】软件下载安装及GIS4WRF插件使用 - 已启用 WPS/WRF 集成,并下载了 预编译的 WRF/WPS 二进制文件
可参见另一博客-【QGIS】软件下载安装及GIS4WRF插件使用 - 若使用 MPI 版本,进程数不超过 6
- 已启用 与 NCAR 数据档案的集成
使用 GIS4WRF 模拟总流程
🧭 全流程共 6 步:
- 创建项目
- 定义模拟区域(Domain)
- 下载输入数据(地理+气象)
- 选择输入数据
- 配置并运行 WPS & WRF
- 查看模拟结果
创建项目
1、打开 QGIS,点击菜单:插件 > GIS4WRF
2、在:Simulation > General中点击 Create a GIS4WRF Project
3、新建一个文件夹命名为:2018_07_15_Amsterdam
4、选择该文件夹作为项目目录
🚨 GIS4WRF 会自动保存所有配置到项目文件中,无需手动保存。
定义模拟区域(Domain)
1、在菜单 Simulation > Domain 中设置:
参数 | 值 |
---|---|
坐标系统 | Lambert Conformal |
True Latitude 1 | 3.5 |
True Latitude 2 | 7 |
Standard Longitude | 4 |
网格间距 | 3000 米 |
中心经度 | 4.8952 |
中心纬度 | 52.3702 |
水平网格数量 | 30 |
垂直网格数量 | 30 |
子域与父域比例 | 3 |
边缘填充 | 上/右/下/左:10 |
2、启用高级配置
3、启用子域嵌套(Parenting)
可在图层中查看域边界并“缩放至组”以聚焦区域
GIS4WRF 从内到外定义域。这使您可以完全控制最内层域的定位。
下载输入数据(地理+气象)
需要两类数据:
- 地理数据(Geo):用于 WPS 地理处理
- 气象数据(Met):用于初始边界条件
⚠ 注意:必须等地理数据和气象数据都下载完成后,才能进入下一步。
🔷 地理数据
菜单:Dataset > Geo
点击:Select Mandatory Fields in Lowest Resolution
然后点击:Download Selected Datasets
静态地理数据下载完成后,再进行气象数据的下载。
🔶 气象数据
菜单:Dataset > Met
设置如下:
参数 | 值 |
---|---|
数据集 | ds083.3 |
产品 | Analysis |
起始时间 | 2018年7月15日 12:00 |
结束时间 | 2018年7月15日 18:00 |
点击“Subset”,在图层面板中选择 Domain 2
点击“Set from Active Layer”
然后点击“Download”开始下载
选择输入数据
菜单:Simulation > Data
在 Domain 1 和 Domain 2 输入框中输入 lowres
选择数据时间段 2018-07-15 12:00 - 2018-07-15 18:00
点击:Use Dataset Selection from List
GIS4WRF 会自动填充 WPS 和 WRF 所需配置(namelists)
配置并运行 WPS & WRF
🌀 运行 WPS
WPS 包括三个程序:
- Geogrid
- Ungrib
- Metgrid
菜单:Simulation > Run
依次点击这三个按钮,每个步骤完成后会提示成功。
🌬️ 运行 WRF
WRF 运行包括两个步骤:
- Real
- WRF
虽然 GIS4WRF 会自动填充大部分配置,但物理参数需手动设置。
✅ 设置 namelist.input 文件
点击 Open configuration,将以下内容粘贴进去并点击保存:
&time_control
start_year = 2018, 2018
start_month = 7, 7
start_day = 15, 15
start_hour = 12, 12
end_year = 2018, 2018
end_month = 7, 7
end_day = 15, 15
end_hour = 18, 18
interval_seconds = 21600
input_from_file = .true., .true.
history_interval = 10, 10
frames_per_outfile = 100, 100
restart = .false.
restart_interval = 7200
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
start_minute = 0, 0
start_second = 0, 0
end_minute = 0, 0
end_second = 0, 0
nocolons = .true.
/
&domains
time_step = 40
time_step_fract_num = 0
time_step_fract_den = 1
max_dom = 2
e_we = 31, 31
e_sn = 31, 31
e_vert = 33, 33
p_top_requested = 5000
num_metgrid_levels = 32
num_metgrid_soil_levels = 4
dx = 9000.0, 3000.0
dy = 9000.0, 3000.0
grid_id = 1, 2
parent_id = 1, 1
i_parent_start = 1, 11
j_parent_start = 1, 11
parent_grid_ratio = 1, 3
parent_time_step_ratio = 1, 3
feedback = 1
smooth_option = 0
/
&physics
physics_suite = 'CONUS'
mp_physics = 0, 0
cu_physics = 0, 0
radt = 9, 3
bldt = 0, 0
cudt = 0, 0
icloud = 0
num_land_cat = 21
sf_urban_physics = 1, 1
/
&dynamics
hybrid_opt = 2
w_damping = 0
diff_opt = 1, 1
km_opt = 4, 4
diff_6th_opt = 0, 0, 0
diff_6th_factor = 0.12, 0.12
base_temp = 290.0
damp_opt = 3
zdamp = 5000.0, 5000.0
dampcoef = 0.2, 0.2
khdif = 0, 0
kvdif = 0, 0
/
&bdy_control
spec_bdy_width = 5
specified = .true.
/
&namelist_quilt
/
完成后,依次运行 Real 和 WRF。
查看模拟结果
点击菜单:Visualize Output
打开文件:wrfout_d01_2018-07-15_12_00_00
可以:
- 切换显示变量(如温度、风速等)
- 切换时间步长
- 双击图层自定义颜色方案(在“Symbology”中)