WPS处理最近的FNL数据时，met_em_d01*土壤温度出现异常值

hzh_1211 · 发表于 2017-9-23 07:53:40

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 hzh_1211 于 2017-9-23 07:53 编辑

我在用WPS3.7.1处理数据时，发现met_em.d01*文件中的soil temperature 出现异常，有些点的温度非常低，接近绝对零度，如图所示，我查看了原始资料，发现并没有这些异常值，因此肯定是在ungrid或者metgrid的过程中出现了错误，ungrid中我是用的是Vtable.GFS,并在官网中（http://www2.mmm.ucar.edu/wrf/users/download/free_data.html）更新了Vtable，metgrid中是用的是METGRID.TBL.ARW，由于我不知道怎么打开ungrid产生的WPS格式的文件，不能确定是这两个过程中的哪个产生的这一问题。
为了验证这一问题的普遍性，我选用了FNL1°和0.25°两种数据产品，在ucar官网中用户自定义下载的，所用数据的日期是2017.09.08-2017.09.18，区域为西太平洋，发现都存在这一问题，但是都能通过real.exe过程，在wrf.exe时会出现以下错误提示：

Flerchinger USEd in NEW version. Iterations= 10
forrtl: severe (174): SIGSEGV, segmentation fault occurred

复制代码

当我将以上异常点全部改为285K的海洋温度后，在运行了几十个步长后，依然出现以下提示：

d01 2017-09-08_01:53:24 MAX AT i,j,k: 6 123 3 vert_cfl,w,d(eta)= 41.78369 -167.2362 1.2999952E-02
Flerchinger USEd in NEW version. Iterations= 10
forrtl: severe (174): SIGSEGV, segmentation fault occurred

复制代码

注意：时间步长我是按照dx*6=162s, 出现错误后进一步减小步长到60s，但错误依然存在。
根据论坛http://forum.wrfforum.com/viewtopic.php?f=6&t=2531&start=10中的解释，是由于土壤水分出现负值导致的，但是我检查了wrfinput_d01,土壤水分没有负值，但是有的点等于0，不知道这会不会引起模式不稳定。
后来我将研究区域缩小到没有保护这些异常值的中国南海区域，但是模式在模拟了7小时后会出现十几处报错：

<blockquote>d01 2017-09-09_07:45:18 2 points exceeded cfl=2 in domain d01 at time 2017-09-09_07:45:18 hours

复制代码

但是能够继续运行，之后模式稳定，正常运行到结束。

因此通过人为的编辑met_em_d0*文件中的ST并不能解决问题，这个应该是新的数据与模式不匹配的问题。
我用CFSv2 2016年的数据做Sarika台风模拟时并没有出现这一情况，并且在做2017年的模拟时，照搬了Sarika案例的物理参数化方案，证明我的处理过程并没有问题。
总结：
FNL 1°和0.25° met_em_d01*出现ST异常值（接近绝对零度）--> 模式迅速出错-->
人为编辑剔除异常值-->模式运行几十步后出错-->
减小模拟区域，使区域中不包含异常点--> 模式运行7h后提示超出CFL判据，但模式继续运行直到稳定，运行成功问题
met_em_d0*文件中出现的ST异常值是如何产生的，如何才能消除，为何人为消除后模式还是不稳定？

附录

资料

NCEP GDAS Final Analysis ds083.3
NCEP Final Analysis (GFS-FNL) ds083.2

WPS设置

&share
wrf_core = 'ARW',
max_dom = 2,
start_date = '2017-09-08_00:00:00','2017-09-08_00:00:00',
end_date = '2017-09-18_00:00:00','2017-09-08_00:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/

&geogrid
parent_id       = 1, 1,
parent_grid_ratio = 1, 3,
i_parent_start = 1,  20,
j_parent_start = 1,  20,
e_we          =  245, 550,
e_sn             =  204,  307,
geog_data_res    = '2m','2m',
dx =    27000,
dy =    27000,
map_proj = 'lambert',
ref_lat =  30,
ref_lon =  130,
truelat1  =  30.0,
truelat2  =  60.0,
stand_lon =  130,
geog_data_path = '/THL6/home/henry/zhhan/coawst/geog'
opt_geogrid_tbl_path = '../geogrid'
/