WRFDA3.2 with radiance data assimilation

Issues with running 3DVAR.

WRFDA3.2 with radiance data assimilation

Postby nalssi » Thu Jun 24, 2010 11:30 pm

Dear all,

It seems that wrfda v3.2 need to many iterations(about 2 times) than v3.1.1 using radiance data assimilation.

Does anyone has any advise?

Thanks
Cho
nalssi
 
Posts: 4
Joined: Thu Jun 24, 2010 3:32 am

Re: WRFDA3.2 with radiance data assimilation

Postby nalssi » Thu Jul 22, 2010 5:37 am

Pasted below is rsl.out.0000 using WRFDA v3.1.1 and v3.2, respectively.


<<<<v3.1.1>>>>

Code: Select all
Reading radiance 1b data from amsua.bufr
 Bufr file date is          2010            7           22            0
 amsua.bufr         
num_tovs_file num_tovs_global num_tovs_local num_tovs_used num_tovs_thinned
    206616      3550       106        67        39
Allocating space for radiance innov structure  2  noaa-16-amsua                  67
Reading radiance 1b data from mhs.bufr
 Bufr file date is          2010            7           22            0
 mhs.bufr           
num_tovs_file num_tovs_global num_tovs_local num_tovs_used num_tovs_thinned
   1103496      2249         0         0         0
Observation summary
   ob time  1
      sound                146 global,       2 local
      synop                898 global,      45 local
      pilot                  5 global,       3 local
      satem                  0 global,       0 local
      geoamv              3908 global,      66 local
      polaramv               0 global,       0 local
      airep               2852 global,       6 local
      gpspw                  0 global,       0 local
      gpsrf                  0 global,       0 local
      metar                709 global,       0 local
      ships                 37 global,       0 local
      ssmi_rv                0 global,       0 local
      ssmi_tb                0 global,       0 local
      ssmt1                  0 global,       0 local
      ssmt2                  0 global,       0 local
      qscat               1638 global,       0 local
      profiler              39 global,       0 local
      buoy                  84 global,       7 local
      bogus                 22 global,       7 local
      pseudo                 0 global,       0 local
      radar                 95 global,       0 local
      radiance            5799 global,      67 local
      airs retrieval         0 global,       0 local
      sonde_sfc            146 global,       2 local
      mtgirs                 0 global,       0 local
      tamdar                 0 global,       0 local

........................

Starting outer iteration :   1
Starting cost function:  2.58140359D+04, Gradient=  8.47581871D+02
For this outer iteration gradient target is:        8.47581871D+00
----------------------------------------------------------
Iter    Cost Function         Gradient             Step
  1      2.46294729D+04      8.36996086D+02      3.29780388D-03
  2      2.23860852D+04      6.40969059D+02      6.40453325D-03
  3      2.16445601D+04      6.18703873D+02      3.60978826D-03
  4      2.05313186D+04      7.03261624D+02      5.81639233D-03
  5      1.94980886D+04      5.07129262D+02      4.17823774D-03
  6      1.86135392D+04      4.33513132D+02      6.87883289D-03
  7      1.76216533D+04      3.72283229D+02      1.05557041D-02
  8      1.68820230D+04      3.30876292D+02      1.06732767D-02
  9      1.62125931D+04      2.82048471D+02      1.22293603D-02
 10      1.57963615D+04      2.53941204D+02      1.04644795D-02
 11      1.53535156D+04      2.15033529D+02      1.37346065D-02
 12      1.51217578D+04      2.15327129D+02      1.00242521D-02
 13      1.48882884D+04      1.65985156D+02      1.00707676D-02
 14      1.46742887D+04      1.66359324D+02      1.55347592D-02
 15      1.45562635D+04      1.42889324D+02      8.52924392D-03
 16      1.44469071D+04      1.23776945D+02      1.07121005D-02
 17      1.43391639D+04      1.08072859D+02      1.40651031D-02
 18      1.42475323D+04      8.45253405D+01      1.56891150D-02
 19      1.41775452D+04      7.38073752D+01      1.96341213D-02
 20      1.41341000D+04      7.08051501D+01      1.52661659D-02
 21      1.41257052D+04      1.18039493D+02      7.46534554D-03
 22      1.40909989D+04      5.66093183D+01      3.33520660D-03
 23      1.40868904D+04      6.82004431D+01      1.53896717D-02
 24      1.40502539D+04      5.80497710D+01      2.90538690D-03
 25      1.40408706D+04      3.65499687D+01      1.31788008D-02
 26      1.40244280D+04      3.20639290D+01      2.33178635D-02
 27      1.40160126D+04      3.02748303D+01      2.04986733D-02
 28      1.40020163D+04      4.37634007D+01      1.06267463D-02
 29      1.40062312D+04      2.45970031D+01      4.53786878D-03
 30      1.39985408D+04      1.98141281D+01      2.07820297D-02
 31      1.39945740D+04      1.70886002D+01      2.07409816D-02
 32      1.39917848D+04      1.50157515D+01      1.90399813D-02
 33      1.39895914D+04      1.21395127D+01      1.94636877D-02
 34      1.39878326D+04      1.05577734D+01      2.38837890D-02
 35      1.39866541D+04      8.89699782D+00      2.08168126D-02
 36      1.39860487D+04      8.56326439D+00      2.29535770D-02
 37      1.39836424D+04      1.52670813D+01      9.56170927D-03
----------------------------------------------------------
 
Inner iteration stopped after   38 iterations
 
Final:  38 iter, J= 1.39830474D+04, g= 6.42276995D+00
----------------------------------------------------------
 
Diagnostics
   Final cost function J       =     13983.05
 
   Total number of obs.        =    43596
   Final value of J            =     13983.04745
   Final value of Jo           =     11032.83018
   Final value of Jb           =      2939.06137
   Final value of Jc           =         0.00000
   Final value of Je           =         0.00000
   Final value of Jp           =        11.15590
   Final J / total num_obs     =         0.32074
   Jb factor used(1)           =         1.00000
   Jb factor used(2)           =         1.00000
   Jb factor used(3)           =         1.00000
   Jb factor used(4)           =         1.00000
   Jb factor used(5)           =         1.00000
   Jb factor used              =         1.00000
   Je factor used              =         1.00000
   VarBC factor used           =         1.00000

   Total number of radiances    =     5794
   Cost function for radiances  =      1489.13865



=====================================================================================
<<<<v3.2>>>>

Code: Select all
Reading radiance 1b data from amsua.bufr
 Bufr file date is          2010            7           22            0
 amsua.bufr         
num_tovs_file num_tovs_global num_tovs_local num_tovs_used num_tovs_thinned
    206616      3550        17        11         6
Allocating space for radiance innov structure  2  noaa-16-amsua                  11
Reading radiance 1b data from mhs.bufr
 Bufr file date is          2010            7           22            0
 mhs.bufr           
num_tovs_file num_tovs_global num_tovs_local num_tovs_used num_tovs_thinned
   1103496      2249         0         0         0
Observation summary
   ob time  1
      sound                146 global,       0 local
      synop                898 global,      37 local
      pilot                  5 global,       2 local
      satem                  0 global,       0 local
      geoamv              3908 global,      47 local
      polaramv               0 global,       0 local
      airep               2852 global,       4 local
      gpspw                  0 global,       0 local
      gpsrf                  0 global,       0 local
      metar                709 global,       0 local
      ships                 37 global,       0 local
      ssmi_rv                0 global,       0 local
      ssmi_tb                0 global,       0 local
      ssmt1                  0 global,       0 local
      ssmt2                  0 global,       0 local
      qscat               1638 global,       0 local
      profiler              39 global,       0 local
      buoy                  84 global,       1 local
      bogus                 22 global,       1 local
      pseudo                 0 global,       0 local
      radar                 95 global,       0 local
      radiance            5799 global,      11 local
      airs retrieval         0 global,       0 local
      sonde_sfc            146 global,       0 local
      mtgirs                 0 global,       0 local
      tamdar                 0 global,       0 local

............
Starting outer iteration :   1
Starting cost function:  2.58140359D+04, Gradient=  8.50030529D+02
For this outer iteration gradient target is:        8.50030529D+00
----------------------------------------------------------
Iter    Cost Function         Gradient             Step
  1      2.46160795D+04      8.43803251D+02      3.31590391D-03
  2      2.23292179D+04      6.48364320D+02      6.42373321D-03
  3      2.15672561D+04      6.28253605D+02      3.62514150D-03
  4      2.04186333D+04      7.13703279D+02      5.82019401D-03
  5      1.93569827D+04      5.13833041D+02      4.16846563D-03
  6      1.84495624D+04      4.39153347D+02      6.87376066D-03
  7      1.74324020D+04      3.77131807D+02      1.05484110D-02
  8      1.66721232D+04      3.35877111D+02      1.06909567D-02
  9      1.59813357D+04      2.86769742D+02      1.22465475D-02
 10      1.55499707D+04      2.58820934D+02      1.04907899D-02
 11      1.50898086D+04      2.19374148D+02      1.37385868D-02
 12      1.48474654D+04      2.20638786D+02      1.00714010D-02
 13      1.46015376D+04      1.70992480D+02      1.01035487D-02
 14      1.43730235D+04      1.72291705D+02      1.56310791D-02
 15      1.42455793D+04      1.49182472D+02      8.58660390D-03
 16      1.41253972D+04      1.30412171D+02      1.08002551D-02
 17      1.40049360D+04      1.14756022D+02      1.41659257D-02
 18      1.39008404D+04      9.07990414D+01      1.58070657D-02
 19      1.38187829D+04      8.08382705D+01      1.99606827D-02
 20      1.37651482D+04      8.06390894D+01      1.55855264D-02
 21      1.37564178D+04      1.31156047D+02      6.66526520D-03
 22      1.37099876D+04      6.54251596D+01      3.82124374D-03
 23      1.36982750D+04      7.92751197D+01      1.61812887D-02
 24      1.36530450D+04      7.21361848D+01      3.06977973D-03
 25      1.36361103D+04      4.64491632D+01      1.28961826D-02
 26      1.36084917D+04      4.33678490D+01      2.47301481D-02
 27      1.35883805D+04      4.18728708D+01      2.21075048D-02
 28      1.35662969D+04      4.37036583D+01      1.94090206D-02
 29      1.35730956D+04      5.64264258D+01      6.51578995D-03
 30      1.35487593D+04      3.49277794D+01      9.41248872D-03
 31      1.35358224D+04      3.34658440D+01      2.27008233D-02
 32      1.35237660D+04      3.31730618D+01      2.13695973D-02
 33      1.35117292D+04      3.07339415D+01      2.18930937D-02
 34      1.34987563D+04      3.08661201D+01      2.74851715D-02
 35      1.34871749D+04      3.02802666D+01      2.39944872D-02
 36      1.34780254D+04      3.73130306D+01      2.59752162D-02
 37      1.34670510D+04      5.03727073D+01      4.23698857D-03
 38      1.34638136D+04      2.91308279D+01      6.95067161D-03
 39      1.34550129D+04      2.97987693D+01      2.83492277D-02
 40      1.34305359D+04      5.71687668D+01      9.21113636D-03
 41      1.34403444D+04      2.84956592D+01      5.22790248D-03
 42      1.34265795D+04      2.59950262D+01      2.99116721D-02
 43      1.34185911D+04      2.82804428D+01      2.86734711D-02
 44      1.34074555D+04      3.74852882D+01      1.39948878D-02
 45      1.34068479D+04      2.43067536D+01      7.97570950D-03
 46      1.33984250D+04      2.63956399D+01      2.57400091D-02
 47      1.33948303D+04      2.87683316D+01      9.42085006D-03
 48      1.33891159D+04      2.10811847D+01      1.34953168D-02
 49      1.33876504D+04      2.42764848D+01      2.17103869D-02
 50      1.33790444D+04      2.67486797D+01      5.95813656D-03
 51      1.33805004D+04      2.33335059D+01      7.74424016D-03
 52      1.33777208D+04      1.98981880D+01      8.44681432D-03
 53      1.33767008D+04      1.96219381D+01      1.46131185D-02
 54      1.33703129D+04      2.34601120D+01      7.62583154D-03
 55      1.33729680D+04      2.46408707D+01      3.57949526D-03
 56      1.33708767D+04      1.38662791D+01      6.11464258D-03
 57      1.33680822D+04      1.26735369D+01      2.74151323D-02
 58      1.33637721D+04      1.33338843D+01      2.49339630D-02
 59      1.33708455D+04      1.93251282D+01      5.30162838D-03
 60      1.33635023D+04      1.18746142D+01      8.63517231D-03
 61      1.33623796D+04      1.05878156D+01      2.22981764D-02
 62      1.33606351D+04      1.00455505D+01      2.85836465D-02
 63      1.33596594D+04      9.65188235D+00      2.89487237D-02
 64      1.33566297D+04      1.17949901D+01      2.60771525D-02
 65      1.33583849D+04      1.08496857D+01      1.27879198D-02
 66      1.33557836D+04      9.11526375D+00      1.73504293D-02
 67      1.33550385D+04      8.54651960D+00      2.72917638D-02
 68      1.33543930D+04      9.36905655D+00      2.50298659D-02
 69      1.33508739D+04      1.76360693D+01      1.92491724D-02
 70      1.33530313D+04      9.59294022D+00      2.09568807D-03
 71      1.33532592D+04      9.94355466D+00      1.76247578D-02
 72      1.33499423D+04      1.06304391D+01      1.05213543D-02
 73      1.33520619D+04      1.07907130D+01      5.06045910D-03
 74      1.33502837D+04      6.97110119D+00      1.53398023D-02
----------------------------------------------------------
 
Inner iteration stopped after   74 iterations
 
Final:  74 iter, J= 1.32871026D+04, g= 6.97110119D+00
----------------------------------------------------------
 
Diagnostics
   Final cost function J       =     13287.10
 
   Total number of obs.        =    43596
   Final value of J            =     13287.10264
   Final value of Jo           =     10656.04632
   Final value of Jb           =      2615.94789
   Final value of Jc           =         0.00000
   Final value of Je           =         0.00000
   Final value of Jp           =        15.10843
   Final J / total num_obs     =         0.30478
   Jb factor used(1)           =         1.00000        1.00000        1.00000        1.00000        1.00000
        1.00000        1.00000        1.00000        1.00000        1.00000
   Jb factor used(2)           =         1.00000        1.00000        1.00000        1.00000        1.00000
        1.00000        1.00000        1.00000        1.00000        1.00000
   Jb factor used(3)           =         1.00000        1.00000        1.00000        1.00000        1.00000
        1.00000        1.00000        1.00000        1.00000        1.00000
   Jb factor used(4)           =         1.00000        1.00000        1.00000        1.00000        1.00000
        1.00000        1.00000        1.00000        1.00000        1.00000
   Jb factor used(5)           =         1.00000        1.00000        1.00000        1.00000        1.00000
        1.00000        1.00000        1.00000        1.00000        1.00000
   Jb factor used              =         1.00000
   Je factor used              =         1.00000
   VarBC factor used           =         1.00000
 
   Total number of radiances    =     5794
   Cost function for radiances  =      1172.40133
nalssi
 
Posts: 4
Joined: Thu Jun 24, 2010 3:32 am

Re: WRFDA3.2 with radiance data assimilation

Postby hnlim » Mon Aug 02, 2010 7:32 pm

Hi

Does setting "use_varbc" further increase the number of inner iterations required significantly in order to fufil the condition set for solution to be achieved?


Agnes
hnlim
 
Posts: 74
Joined: Wed Apr 16, 2008 11:51 am

Re: WRFDA3.2 with radiance data assimilation

Postby hnlim » Tue Aug 03, 2010 12:46 pm

Hi

In addition which forward model did you use? RTTOV or CRTM

Agnes
hnlim
 
Posts: 74
Joined: Wed Apr 16, 2008 11:51 am

Re: WRFDA3.2 with radiance data assimilation

Postby nalssi » Tue Aug 10, 2010 6:04 am

Hi

I used RTTOV as the radiative transfer model in two WRFDA versions.
Also I set "use_varbc" in both WRFDA namelist files, but WRFDA 3.2 increases iterations only.

Thank you for your reply.

Cho.
nalssi
 
Posts: 4
Joined: Thu Jun 24, 2010 3:32 am

Re: WRFDA3.2 with radiance data assimilation

Postby sravan » Fri Sep 03, 2010 5:54 am

Hi
I am using wrfda 3.2. For the radinace data assimilation I have used the BUFR and CRTM libraries to get the executables. I am using amsua and amsub radiances from ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/ website. I have linked the following files
fg, ob.ascii, da_wrfvar.exe, be.dat, LANDUSE.TBL and NCEP bufr files amsua and amsub.
As I give the command ./da_wrfvar.exe following is coming.
---------------------------- WARNING --------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_read_obs_bufr.html">da_read_obs_bufr.inc</a> LINE: 169
Error 209 opening PREPBUFR obs file ob.bufr
---------------------------------------------------------------
da_read_obs_bufr: num_report, num_outside_all, num_outside_time, num_thinned: 0 0 0 0
0: DEALLOCATE: memory at (nil) not allocated

Kindly help me.

Sravan
sravan
 
Posts: 3
Joined: Fri Sep 03, 2010 5:30 am

Re: WRFDA3.2 with radiance data assimilation

Postby hnlim » Fri Sep 03, 2010 6:46 am

What did you set for ob_format in wrfvar namelist 3?
hnlim
 
Posts: 74
Joined: Wed Apr 16, 2008 11:51 am

Re: WRFDA3.2 with radiance data assimilation

Postby cycle_13 » Tue Dec 14, 2010 5:35 am

Hi ,do you have met the errors like as following?

**************BUFR ARCHIVE LIBRARY ABORT*****************
BUFRLIB: OPENBF - ERROR READING INPUT FILE CONNECTED TO UNIT 26 WHEN CHECKING FOR 'BUFR' IN FIRST 4 BYTES OF RECORD
**************BUFR ARCHIVE LIBRARY ABORT*****************

I use pgi compiler to run radiance testcase on Linux with wrfda 3.2 and crtm 1.2, the WRFDA website said there have some tips to solve the problems, but it says it is only take effect on GFORTRAN and IFORT compiler while without anything about PGI compiler. So I want to ask you for some help, can you give me some hints? Thank you very much.
cycle_13
 
Posts: 27
Joined: Mon Jun 21, 2010 8:48 pm


Return to Runtime Problems

Who is online

Users browsing this forum: No registered users and 1 guest