WRF 3DVar相关问题请教

问题少女

1、想请教一下各位前辈，运行完WRFVAR之后，我绘制了相关的cost function函数图像，为什么我的cost function画出来的图像并不是单调递减的呢？它是波动的，如下图所示。这种情况应该如何解决呢？

具体情况如下：
当我设置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
**********************************************************************此行为标记行，以上是输出内容

并不是单调递减的？咨询一下各位前辈，这是为什么啊？ 以及如何解决呢？
同化小白，希望大家给予指导，谢谢各位前辈了。

freekiller

你这个迭代一点都不稳定，说实话我估计你同化的效果应该是很糟糕

问题少女

freekiller 发表于 2016-5-7 17:51
你这个迭代一点都不稳定，说实话我估计你同化的效果应该是很糟糕

首先感谢您的回复，其次是同化的效果也不是很糟糕，但是我想知道为什么会出现这样的现象，难道两个cost function不是一个东西吗？为什么打开一个开关（calculate_cg_cost_fn）会影响cost function的值呢？这是不是说明处理的有问题，还是怎么样呢？应该怎样改呢？麻烦您了！

问题少女

freekiller 发表于 2016-5-7 17:51
你这个迭代一点都不稳定，说实话我估计你同化的效果应该是很糟糕

还有就是，我想知道这个到底是对的还是错的呢？打扰您了。

freekiller

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

即使结果不糟糕，也很难让人信服啊

问题少女

freekiller 发表于 2016-5-7 19:50
即使结果不糟糕，也很难让人信服啊

是这样的，我现在特别想知道的是这是不是说明我做的东西是错误的呢？还是怎么样？

问题少女

freekiller 发表于 2016-5-7 19:50
即使结果不糟糕，也很难让人信服啊

还是说同化的结果不好呢？打扰您了。

freekiller

问题少女 发表于 2016-5-7 19:56
是这样的，我现在特别想知道的是这是不是说明我做的东西是错误的呢？还是怎么样？

我不敢轻易说你做的东西对错，只是从常理来说有点不太合理，再检查一下你的整个数据

问题少女

freekiller 发表于 2016-5-7 20:05
我不敢轻易说你做的东西对错，只是从常理来说有点不太合理，再检查一下你的整个数据

恩恩，好的，我再检查一下我的雷达数据，我只同化了雷达数据，不知道是不是因为这个。。。谢谢您的指导啦~

zhchlue

谢谢楼主的分享，有机会可以了解一下