爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 11554|回复: 18

在wrfda 中利用ETKF更新集合成员有问题

[复制链接]
发表于 2013-10-27 22:03:51 | 显示全部楼层 |阅读模式

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

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

x
请问大家,有没有用过WRFDA中的ETKF更新集合成员呢?
我用测试数据发现通过ETKF更新后的集合成员和更新前的集合成员是不一样的,主要表现在更新后的集合扰动空间上看起来更为的平滑,但是用我自己的数据发现更新后的集合成员1和集合平均一样,而其它几个集合成员并没有发生变化,请问这是怎么回事啊?
我感觉问题主要出现在./gen_be_etkf.exe,我对gen_be_etkf_nl.nl里面的一些参数不是很懂,还请大侠们指教一下! 不甚感激!谢谢谢谢!
下面是我的gen_be_etkf_nl.nl
  &gen_be_etkf_nl
    num_members = 10,
    nv = 7,
    cv = 'U', 'V', 'W', 'PH', 'T', 'QVAPOR', 'MU',
    naccumt1 = 0,
    naccumt2 = 0,
    nstartaccum1 = 1,
    nstartaccum2 = 1,
    nout = 1,
    tainflatinput = 1,
    rhoinput = 1,
    infl_fac_file = 'inflation_factor.dat',
    infl_let_file = 'nflation_letkf.dat',
    eigen_val_file = 'eigen_value.dat',
    inno2_val_file = 'innovation_value.dat',
    proj2_val_file = 'projection_value.dat',
    infl_fac_TRNK = .false.,
    infl_fac_WG03 = .false.,
    infl_fac_WG07 = .true.,
    infl_fac_BOWL = .false.,
    letkf_flg=.false.,
    rand_filt = .false.,
    rnd_seed = 2008012712,
    rnd_nobs = 5000
    etkf_erro_max = 20.
    etkf_erro_min = .00001
    etkf_inno_max = 20.
    etkf_inno_min = .00001
    etkf_erro_flg = .true.
    etkf_inno_flg = .true.
    etkf_wrfda = .false.
/

密码修改失败请联系微信:mofangbao
发表于 2013-11-4 00:13:22 | 显示全部楼层
我也在用这个,可是我在运行gen_be_etkf.exe时只生成ob.etkf.filt.e* etkf_output根本没更新。。。而且这个怎么不输出rsl.error和rsl.out文件呢?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2013-11-4 11:31:02 | 显示全部楼层

我是在前台运行,让它输出在屏幕上的,不知道你的有没有提示什么错误?
密码修改失败请联系微信:mofangbao
发表于 2013-11-5 11:42:41 | 显示全部楼层
小虹 发表于 2013-11-4 11:31
我是在前台运行,让它输出在屏幕上的,不知道你的有没有提示什么错误?

我提交到银河机上,没有错误提示,但是程序退出。我也是是前台运行。。。
密码修改失败请联系微信:mofangbao
发表于 2013-11-6 12:34:08 | 显示全部楼层
简单乐活 发表于 2013-11-5 11:42
我提交到银河机上,没有错误提示,但是程序退出。我也是是前台运行。。。

错误如下:forrt1:SEVERE(174):SIGSEGV,segmentation fault occurred
密码修改失败请联系微信:mofangbao
发表于 2013-11-6 12:34:42 | 显示全部楼层
小虹 发表于 2013-11-4 11:31
我是在前台运行,让它输出在屏幕上的,不知道你的有没有提示什么错误?

forrt1:SEVERE(174):SIGSEGV,segmentation fault occurred
密码修改失败请联系微信:mofangbao
发表于 2014-4-13 23:07:26 | 显示全部楼层
LZ我想问一下,你先前的集合成员是否是基于观测资料上加扰动之后预报6小时所得到的。用的是哪种观测资料?新手也在做这个  不胜感激
密码修改失败请联系微信:mofangbao
发表于 2014-4-20 00:24:28 | 显示全部楼层
小虹 发表于 2013-11-4 11:31
我是在前台运行,让它输出在屏幕上的,不知道你的有没有提示什么错误?

上次问题解决了,属于电脑编译器的问题。。。
密码修改失败请联系微信:mofangbao
发表于 2014-5-5 00:35:27 | 显示全部楼层
简单乐活 发表于 2013-11-6 12:34
forrt1:SEVERE(174):SIGSEGV,segmentation fault occurred

求真相。同样的问题,费了很大劲,没搞明白。

编译器的哪些问题?能否给个详细的介绍,也好让大家都学习下,多谢了。
密码修改失败请联系微信:mofangbao
发表于 2014-5-23 07:51:52 | 显示全部楼层
wudi1201 发表于 2014-5-5 00:35
求真相。同样的问题,费了很大劲,没搞明白。

编译器的哪些问题?能否给个详细的介绍,也好让大家都 ...

你是用的什么编译器,我开始用的internel后面改为PGI就好了
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表