请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6953|回复: 12

WRF 3DVar相关问题请教

[复制链接]

新浪微博达人勋

发表于 2016-5-7 12:41:14 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 问题少女 于 2017-3-21 15:14 编辑

1、想请教一下各位前辈,运行完WRFVAR之后,我绘制了相关的cost function函数图像,为什么我的cost function画出来的图像并不是单调递减的呢?它是波动的,如下图所示。这种情况应该如何解决呢?
QQ截图20160507123835.png
具体情况如下:
当我设置calculate_cg_cost_fn=false,时,输出在rsl.out.0000中的cost function是这样的(cost_fn中仅输出最初的和最终的cost function):
**********************************************************************此行为标记行,以下是输出内容
Starting outer iteration :   1
Starting cost function:  1.850135433455994D+04, Gradient=  2.023651300694863D+03
For this outer iteration gradient target is:        2.023651300694863D+01
----------------------------------------------------------
Iter    Cost Function         Gradient             Step
  1      1.119339373800786D+04      1.637790701250011D+03      3.569068076092051D-03
  2      9.038978010503586D+03      4.844914329442357D+02      1.606359346078572D-03
  3      7.013372468989302D+03      6.790646189191600D+02      1.725888234321006D-02
  4      5.841891642782608D+03      3.201132391570301D+02      5.080927221582558D-03
  5      5.581638422784656D+03      2.193054304840096D+02      5.079475095252152D-03
  6      5.147189116033407D+03      2.853920602484674D+02      1.806634637462796D-02
  7      4.617793177194228D+03      2.225853243314326D+02      1.299950536907846D-02
  8      4.109406251576758D+03      1.662139369138801D+02      2.052254966630655D-02
  9      3.803526979580640D+03      1.432084967288243D+02      2.214344414551412D-02
10      3.545256201123293D+03      1.082956347688449D+02      2.518649273078374D-02
11      3.372375807480053D+03      9.753346091804764D+01      2.948178915308340D-02
12      3.228715870289940D+03      5.978758147628481D+01      3.020357825411410D-02
13      3.181905271889629D+03      7.028301623240058D+01      2.619100822168356D-02
14      3.110993957843571D+03      4.604783557093056D+01      2.871076379989686D-02
15      3.069703323556176D+03      3.603853572538664D+01      3.894602216747459D-02
16      3.040939399255598D+03      2.607822475116062D+01      4.429389362339538D-02
17      3.026036669200979D+03      2.593788271508789D+01      4.382680196726843D-02
18      3.012314913311424D+03      2.189480073276839D+01      4.079159047484099D-02
19      3.002920171279729D+03      1.403180221674693D+01      3.919519772221250D-02
----------------------------------------------------------

Inner iteration stopped after   19 iterations
**********************************************************************此行为标记行,以上是输出内容
cost function是单调递减的,但是当我设置calculate_cg_cost_fn=true,时,输出在cost_fn中的相关内容如下(rsl.out.0000中的也是这样):

**********************************************************************此行为标记行,以下是输出内容
Outer    EPS     Inner      J           Jb       Jo           Jc         Je         Jp         Js        jl
Iter             Iter
   1  0.100E-01     0  18501.354      0.000  18501.354      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01     1  17845.824     26.083  17819.741      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01     2  17228.158     47.189  17180.969      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01     3  17827.180    112.441  17714.739      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01     4  17691.344    192.664  17498.680      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01     5  17647.155    217.059  17430.096      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01     6  17805.551    275.774  17529.777      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01     7  18057.652    392.572  17665.080      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01     8  18098.367    567.794  17530.573      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01     9  18279.128    714.892  17564.236      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01    10  18016.271    876.562  17139.709      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01    11  17758.032   1012.987  16745.045      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01    12  17758.140   1152.544  16605.595      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01    13  17834.185   1205.404  16628.782      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01    14  17908.744   1297.951  16610.793      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01    15  18052.511   1361.117  16691.394      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01    16  18086.694   1411.629  16675.065      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01    17  18137.642   1441.146  16696.496      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01    18  18180.231   1471.664  16708.567      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01    19  18219.276   1495.110  16724.166      0.000      0.000      0.000      0.000      0.000
   1  0.100E-01    19  18219.276   1495.110  16724.166      0.000      0.000      0.000      0.000      0.000
**********************************************************************此行为标记行,以上是输出内容

并不是单调递减的?咨询一下各位前辈,这是为什么啊? 以及如何解决呢?
同化小白,希望大家给予指导,谢谢各位前辈了。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-7 17:51:58 | 显示全部楼层
你这个迭代一点都不稳定,说实话我估计你同化的效果应该是很糟糕
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-7 18:39:25 | 显示全部楼层
freekiller 发表于 2016-5-7 17:51
你这个迭代一点都不稳定,说实话我估计你同化的效果应该是很糟糕

首先感谢您的回复,其次是同化的效果也不是很糟糕,但是我想知道为什么会出现这样的现象,难道两个cost function不是一个东西吗?为什么打开一个开关(calculate_cg_cost_fn)会影响cost function的值呢?这是不是说明处理的有问题,还是怎么样呢?应该怎样改呢?麻烦您了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-7 18:47:09 | 显示全部楼层
freekiller 发表于 2016-5-7 17:51
你这个迭代一点都不稳定,说实话我估计你同化的效果应该是很糟糕

还有就是,我想知道这个到底是对的还是错的呢?打扰您了。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-7 19:50:02 | 显示全部楼层
问题少女 发表于 2016-5-7 18:39
首先感谢您的回复,其次是同化的效果也不是很糟糕,但是我想知道为什么会出现这样的现象,难道两个cost f ...

即使结果不糟糕,也很难让人信服啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-7 19:56:30 | 显示全部楼层
freekiller 发表于 2016-5-7 19:50
即使结果不糟糕,也很难让人信服啊

是这样的,我现在特别想知道的是这是不是说明我做的东西是错误的呢?还是怎么样?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-7 19:57:46 | 显示全部楼层
freekiller 发表于 2016-5-7 19:50
即使结果不糟糕,也很难让人信服啊

还是说同化的结果不好呢?打扰您了。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-7 20:05:24 | 显示全部楼层
问题少女 发表于 2016-5-7 19:56
是这样的,我现在特别想知道的是这是不是说明我做的东西是错误的呢?还是怎么样?

我不敢轻易说你做的东西对错,只是从常理来说有点不太合理,再检查一下你的整个数据
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-5-7 20:41:10 | 显示全部楼层
freekiller 发表于 2016-5-7 20:05
我不敢轻易说你做的东西对错,只是从常理来说有点不太合理,再检查一下你的整个数据

恩恩,好的,我再检查一下我的雷达数据,我只同化了雷达数据,不知道是不是因为这个。。。谢谢您的指导啦~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-5-10 18:28:15 | 显示全部楼层
谢谢楼主的分享,有机会可以了解一下
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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