登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 夏朗的芒果 于 2016-5-5 20:42 编辑
WRF 中 Bogus 算法过程1. 在初始场 TC 中心的 300 km 半径范围内移除相关涡度、散度与地转涡度,达到移除背景场中的 TC 环流的目的。 2. 通过观测得到的最大风速半径以及某一个 WRF 自带的风廓线,在观测 TC 位置处移入一个轴对称的兰金涡旋。
引言WRF ARW中提供了一个简单的 TC Bogusing 方案。能移除已经存在的热带气旋,并且可以重新植入一个兰金涡旋。程序的输入文件为单时刻单区域的 metgrid 数据,并且需要在namelist.input 中添加一些定义参数来描述 bogus TC 的位置与强度。输出文件同样为 metgrid 类似数据文件。该方案现在只能作用于等压面数据。在运行 tc.exe 程序后,用户需要重命名输出文件,这样 real.exe 才能读取更改后的 metgrid 文件。
Namelist 中参数含义TC 方案的 namelist 信息位于一个可选的 namelist 记录 &tc 后。一次只能作用于一个区域。嵌套网格用户应该将通常气象场水平插值到精细网格区域中。或者,用户可以分别对不同区域的的 metgrid 文件逐一运行tc.exe, 但是不建议此操作。 insert_bogus_storm | 逻辑值,true or false | remove_storm | 逻辑值,true or false | num_storm | 整数,bogus TC 数目,现阶段只能为1 | latc_loc | 实数,bogus 纬度 ( + 北纬,- 南纬 ) | lonc_loc | 实数,bogus 经度 ( + 东经,- 西经 ) | vmax_meters_per_second | 实数,观测持续风速最大值 (m/s) | rmax | 实数,涡旋中心至最大风速出现处半径 (m) | vmax_ratio | 实数,模式兰金涡旋的尺度大小 |
注意: 如果 insert_bogus_storm 设置为 true, 那么remove_storm 应该设置成 false. 如果remove_storm 设置为 true, 那么insert_bogus_storm 设置成 fasle。
对于 45 km 分辨率区域,vmax_ratio 设置成 0.75;对于 15km 分辨率区域,设置成0.90。vmax_ratio 是一个控制比例大小的参数。因为当网格区域较粗时,不可能整个网格元都能达到观测最大风速值。
例如,假设观测 TC 中心位于 25°N, 75°W,最大持续风速达到 120 节,最大风速出现在距离中心 90 km 处,那么对于 45 km 的网格来说, namelist.input 应该设置成
&tc insert_bogus_storm = .true. remove_storm = .false. latc_loc = 25.0 lonc_loc = -75.0 vmax_meters_per_second = 61.7 rmax = 90000.0 vmax_ratio = 0.75 /
tc.exe 程序tc.exe 自动随 wrf 编译而成,但是只能进行串行运行。为节约时间,可使用串行与无嵌套选项进行编译。
运行 tc.exe1) 正常运行所有 WPS 程序 ( geogrid, ungrib 和 metgrid ) 2) 链接 met 文件 3) 编辑 namelist.input 文件,添加 tc.exe 相关参数。同时编辑模拟时间为单时刻 4) 运行tc.exe 5) 重命名输出文件 auxinput1_d01_ 为 met_em.d01. 6) 编辑 namelist.input 文件中的模拟时刻,改为所计划的时间,运行 real.exe
注意,tc.exe 只能对 d01 的初始时刻进行操作,如果有多重网格,可在 namelist.input 中添加如下设置: &time_control input_from_file = .true. , .true., .true., fine_input_stream = 0, 2, 2, io_form_auxinput2 = 2 /
fine_input_stream 参数设置为 2 允许嵌套网格的 3D 气象场从 d01 插值得到,而静态场(如地形,土地利用)和地表数据任由嵌套网格的 wrfinput 提供;
|