[金币求助]WRFDA运行出错

lee309

生成be.dat成功，生成observation文件成功，然后运行da_wrfvar.exe，出错
Namelist logging not found in namelist.input. Using registry defaults for variables in logging.
Namelist logging not found in namelist.input. Using registry defaults for variables in logging.
***  VARIATIONAL ANALYSIS ***
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain             1 ,                 209693392  bytes allocated
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS      1 IE     59 JS      1 JE     60
WRF NUMBER OF TILES =   1
Set up observations (ob)
  
Using ASCII format observation input
  
Observation summary
   ob time  1
      synop                116 global,     116 local
      metar                 28 global,      28 local

Set up background errors for regional application for cv_options =   5

  

   Using the averaged regression coefficients for unbalanced part
  
   WRF-Var dry control variables are: psi, chi_u, t_u and ps_u
   Humidity control variable is rh
  
Vertical truncation for psi    =   0(  99.00%)

Vertical truncation for chi_u  =   0(  99.00%)

Vertical truncation for t_u    =   0(  99.00%)

Vertical truncation for rh     =   0(  99.00%)

  
   Scaling: var, len, ds:   0.100000E+01   0.100000E+01   0.300000E+05
   Scaling: var, len, ds:   0.100000E+01   0.100000E+01   0.300000E+05
   Scaling: var, len, ds:   0.100000E+01   0.100000E+01   0.300000E+05
   Scaling: var, len, ds:   0.100000E+01   0.100000E+01   0.300000E+05
   Scaling: var, len, ds:   0.100000E+01   0.100000E+01   0.300000E+05
Calculate innovation vector(iv)
  
Minimize cost function using CG method
  
da_transform_vtovv_adj.inc: be%v1%mz=0
da_transform_vtovv_adj.inc: be%v2%mz=0
da_transform_vtovv_adj.inc: be%v3%mz=0
da_transform_vtovv_adj.inc: be%v4%mz=0
Starting outer iteration :   1
Starting cost function:  4.254539167082756D+02, Gradient=                    NaN
For this outer iteration gradient target is:                          NaN
----------------------------------------------------------
Iter    Cost Function         Gradient             Step
----------------------------------------------------------

Inner iteration stopped after    1 iterations

da_transform_vtovv.inc: be%v1%mz=0
da_transform_vtovv.inc: be%v2%mz=0
da_transform_vtovv.inc: be%v3%mz=0
da_transform_vtovv.inc: be%v4%mz=0
da_transform_vtovv_adj.inc: be%v1%mz=0
da_transform_vtovv_adj.inc: be%v2%mz=0
da_transform_vtovv_adj.inc: be%v3%mz=0
da_transform_vtovv_adj.inc: be%v4%mz=0
Final:   1 iter, J=                   NaN, g=                   NaN
----------------------------------------------------------
da_transform_vtovv.inc: be%v1%mz=0
da_transform_vtovv.inc: be%v2%mz=0
da_transform_vtovv.inc: be%v3%mz=0
da_transform_vtovv.inc: be%v4%mz=0
  
Diagnostics
   Final cost function J       =          NaN
  
   Total number of obs.        =      384
   Final value of J            =             NaN
   Final value of Jo           =             NaN
   Final value of Jb           =         0.00000
   Final value of Jc           =         0.00000
   Final value of Je           =         0.00000
   Final value of Jp           =         0.00000
   Final value of Jl           =         0.00000
   Final J / total num_obs     =             NaN
   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        ---------------------------- WARNING --------------------------
WARNING FROM FILE: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_read_omb_tmp.html">da_read_omb_tmp.inc</a>  LINE:     425
---------------------------------------------------------------
---------------------------- FATAL ERROR -----------------------
----------------------------------------------------------------
wrf_abort
      =         1.00000
   VarBC factor used           =         1.00000
  
read error on unit:  26
  
Fatal error in file: <A HREF="http://www.mmm.ucar.edu/people/wrfhelp/wrfvar/code/trunk/da_read_omb_tmp.html">da_read_omb_tmp.inc</a>  LINE:     412
Got unknown obs_type string:                   1 on unit  26
1. 开头的Namelist logging not found in namelist.input. Using registry defaults for variables in logging.不知道是什么意思？
2. 有很多数据都是0很奇怪，如psi, chi_u, t_u and ps_u等，这是怎么回事，是因为观测数据中没有这些变量么？
3. 最后的Got unknown obs_type string又是什么意思，貌似重点就在这，哪位做过同化懂这方面的指点下？
虽然observation文件生成了，但是运行obsproc.exe也有点error和警告
READ NAMELIST FILE: namelist.obsproc
------------------
&RECORD1
OBS_GTS_FILENAME = /usr/WRF/test/WRFDA/TestData/SURFACE_OBS:2013030112                             ,
OBS_ERR_FILENAME = obserr.txt                                                                      ,
FG_FORMAT = MM5,
GTS_FROM_MMM_ARCHIVE =  F
/
&RECORD2
TIME_WINDOW_MIN = 2013-03-01_11:00:00,
TIME_ANALYSIS = 2013-03-01_12:00:00,
TIME_WINDOW_MAX = 2013-03-01_13:00:00
/
&RECORD3
MAX_NUMBER_OF_OBS =       400000,
FATAL_IF_EXCEED_MAX_OBS =  T
/
&RECORD4
QC_TEST_VERT_CONSISTENCY =  T,
QC_TEST_CONVECTIVE_ADJ =  T,
QC_TEST_ABOVE_LID =  T,
REMOVE_ABOVE_LID =  T,
DOMAIN_CHECK_H =  T,
THINING_SATOB =  F,
THINING_SSMI =  F,
THINING_QSCAT =  F,
CALC_PSFC_FROM_QNH =  F
/
&RECORD5
PRINT_GTS_READ =  T,
PRINT_GPSPW_READ =  T,
PRINT_RECOVERP =  T,
PRINT_DUPLICATE_LOC =  T,
PRINT_DUPLICATE_TIME =  T,
PRINT_RECOVERH =  T,
PRINT_QC_VERT =  T,
PRINT_QC_CONV =  T,
PRINT_QC_LID =  T,
PRINT_UNCOMPLETE =  T
/
&RECORD6
PTOP =    1000.000    ,
PS0 =    100000.0    ,
TS0 =    290.0000    ,
TLP =    50.00000    ,
PIS0 =    20000.00    ,
TIS0 =    215.0000    ,
BASE_PRES =    100000.0    ,
BASE_TEMP =    290.0000    ,
BASE_LAPSE =    50.00000    ,
BASE_STRAT_TEMP =    215.0000    ,
BASE_TROPO_PRES =    20000.00   
/
&RECORD7
IPROJ =            1,
PHIC =    40.00001    ,
XLONC =    122.0000    ,
TRUELAT1 =    30.00000    ,
TRUELAT2 =    60.00000    ,
MOAD_CEN_LAT =    28.00000    ,
STANDARD_LON =    122.0000   
/
&RECORD8
IDD =            1,
MAXNES =            1,
NESTIX =           60,
                   200,
                   136,
                   181,
                   211,
                     0,
                     0,
                     0,
                     0,
                     0,
NESTJX =           90,
                   200,
                   181,
                   196,
                   211,
                     0,
                     0,
                     0,
                     0,
                     0,
DIS =    60.00000    ,
          10.00000    ,
          3.300000    ,
          1.100000    ,
          1.100000    ,
          0.000000    ,
          0.000000    ,
          0.000000    ,
          0.000000    ,
          0.000000    ,
NUMC =            1,
                   1,
                   2,
                   3,
                   4,
                   0,
                   0,
                   0,
                   0,
                   0,
NESTI =            1,
                   40,
                   28,
                   35,
                   45,
                    0,
                    0,
                    0,
                    0,
                    0,
NESTJ =            1,
                   60,
                   25,
                   65,
                   55,
                    0,
                    0,
                    0,
                    0,
                    0
/
&RECORD9
PREPBUFR_OUTPUT_FILENAME = prepbufr_output_filename                                                        ,
PREPBUFR_TABLE_FILENAME = prepbufr_table_filename                                                         ,
OUTPUT_OB_FORMAT =            2,
USE_FOR = 3DVAR,
NUM_SLOTS_PAST =            3,
NUM_SLOTS_AHEAD =            3,
WRITE_SYNOP =  T,
WRITE_SHIP =  T,
WRITE_METAR =  T,
WRITE_BUOY =  T,
WRITE_PILOT =  T,
WRITE_SOUND =  T,
WRITE_AMDAR =  T,
WRITE_SATEM =  T,
WRITE_SATOB =  T,
WRITE_AIREP =  T,
WRITE_GPSPW =  T,
WRITE_GPSZTD =  T,
WRITE_GPSREF =  T,
WRITE_GPSEPH =  T,
WRITE_SSMT1 =  T,
WRITE_SSMT2 =  T,
WRITE_SSMI =  T,
WRITE_TOVS =  T,
WRITE_QSCAT =  T,
WRITE_PROFL =  T,
WRITE_BOGUS =  T,
WRITE_AIRS =  T,
WRITE_TAMDAR =  T
/
&RECORD10
WIND_SD =  F,
WIND_SD_SYNOP =  F,
WIND_SD_SHIPS =  F,
WIND_SD_METAR =  F,
WIND_SD_BUOY =  F,
WIND_SD_SOUND =  F,
WIND_SD_QSCAT =  F,
WIND_SD_PILOT =  F,
WIND_SD_AIREP =  F,
WIND_SD_GEOAMV =  F,
WIND_SD_TAMDAR =  F,
WIND_SD_PROFILER =  F
/
Error reading NAMELIST record 10
Using default values for wind_sd
COARSE GRID CENTER IS AT X =      0.0 KM AND Y =  -6616.7 KM.   DIS(1)=    60.000
------------------------------------------------------------------------------
READ GTS OBSERVATIONS IN FILE /usr/WRF/test/WRFDA/TestData/SURFACE_OBS:2013030112

Diagnostics in file obs_gts_read.diag


   304  modified the original lat =  -90.00 to -89.5


  1274  modified the original lat =  -90.00 to -89.5

Have reached the end of observation file.

------------------------------------------------------------------------------
GTS OBSERVATIONS READ:

SYNOP reports: 27673
SHIPS reports:  1455
BUOYS reports: 16207
BOGUS reports:     0
METAR reports: 46057
PILOT reports:     0
SOUND reports:     0
AMDAR reports:     0
SATEM reports:     0
SATOB reports:     0
GPSPW reports:     0
GPSZD reports:     0
GPSRF reports:     0
GPSEP reports:     0
AIREP reports:     0
TAMDAR reports:     0
SSMT1 reports:     0
SSMT2 reports:     0
SSMI  reports:     0
TOVS  reports:     0
QSCAT reports:     0
PROFL reports:     0
AIRST reports:     0
OTHER reports:     0
Total reports: 91392

Number of observations read:             91392
Number of empty observations:               35
Number of observations out of domain:    89717
Number of observations for ingestion:     1640


AIRCRAFT DATA: Total=     0  Above cut_height=     0
------------------------------------------------------------------------------
PRESSURE RECOVERED FROM HEIGHTS:

Diagnostics in file obs_recover_pressure.diag

------------------------------------------------------------------------------
CHECK PRESENCE OF HEIGHT OR/AND PRESSURE

Diagnostics in file obs_check_pressure.diag

------------------------------------------------------------------------------
REMOVE DUPLICATE STATIONS BY LOCATION:

Diagnostics in file obs_duplicate_loc.diag

Found   155 location duplicate stations that have been merged.

------------------------------------------------------------------------------
REMOVE DUPLICATE STATIONS BY TIME:

Diagnostics in file obs_duplicate_time.diag_2013-03-01_12:00:00

Found   594 time duplicate stations that have been removed.

------------------------------------------------------------------------------
DIAGNOSTICS: U, V, RH, QV etc...:

------------------------------------------------------------------------------
HEIGHT RECOVERED FROM P, T, Q,..:

Diagnostics in file obs_recover_height.diag

------------------------------------------------------------------------------
CHECK PRESENCE OF HEIGHT OR/AND PRESSURE

Diagnostics in file obs_check_height.diag

------------------------------------------------------------------------------

<AFWA> OBSERVATIONAL ERRORS PER TYPE:

HEIGHT SENSOR ERRORS:  RAOBS, PIBALS, AIREPS, TAMDARS, SURFACE SHIP, BUOY, DMSP - A RETRIEVAL,
DMSP - A RETRIEVAL,
NOAA - CLEAR PATH RETRIEVAL,
NOAA - A RETRIEVAL, PROFL, AIRSRET,

UV SENSOR ERRORS:  RAOBS, PIBALS, AIREPS, TAMDARS, US LOW LEVEL WINDS, PROFL, BUOY,

SPD SENSOR ERRORS:  RAOBS, PIBALS, AIREPS, TAMDARS, US LOW LEVEL WINDS, PROFL, BUOY,

DIR SENSOR ERRORS:  RAOBS, PIBALS, TAMDARS, AIREPS, US LOW LEVEL WINDS, PROFL, BUOY,

TEMP SENSOR ERRORS:  RAOBS, PIBALS, AIREPS, TAMDARS, SURFACE LAND, SURFACE SHIP, BUOY, DMSP - A RETRIEVAL, DMSP - A RETRIEVAL, NOAA - A RETRIEVAL, AIRSRET,

RH SENSOR ERRORS ( % ):  RAOBS, AIREPS, TAMDARS, PIBALS, DMSP - A RETRIEVAL, NOAA - A RETRIEVAL, SURFACE LAND, AIRSRET,

PRESSURE SENSOR ERRORS:  RAOBS, PIBALS, AIREPS, TAMDARS, SURFACE LAND, SURFACE SHIP, DMSP - A RETRIEVAL, NOAA - A RETRIEVAL, AIRSRET,

Number of processed stations:              891 =    867 levels.
Number of processed surface stations:      891 =    891 surface levels.
Number of processed upper-air stations:      0 =    -24 upper-air levels.

------------------------------------------------------------------------------
VERTICAL CONSISTENCY TEST QC:

Diagnostics in file obs_qc1.diag

------------------------------------------------------------------------------
CONVECTIVE ADJUSTEMENT TEST QC:

Diagnostics in file obs_qc1.diag

------------------------------------------------------------------------------
CHECK OBS HEIGHT...:

Diagnostics in file obs_qc2.diag

------------------------------------------------------------------------------
LOOK FOR UNCOMPLETE DATA:

Diagnostics in file obs_uncomplete.diag
Remove      1 surface stations.
Remove      0 upper-air levels.

------------------------------------------------------------------------------
INGESTED OBSERVATION AFTER CHECKS:

                 READ   EMPTY OUTSIDE LOCDUPL TIMDUPL UNCOMPL INGESTD
SYNOP reports: 27673       0   26400     150     358       0     765
SHIPS reports:  1455       0    1424       5       2       0      24
BUOYS reports: 16207      35   16170       0       1       1       0
BOGUS reports:     0       0       0       0       0       0       0
METAR reports: 46057       0   45723       0     233       0     101
PILOT reports:     0       0       0       0       0       0       0
SOUND reports:     0       0       0       0       0       0       0
AMDAR reports:     0       0       0       0       0       0       0
SATEM reports:     0       0       0       0       0       0       0
SATOB reports:     0       0       0       0       0       0       0
AIREP reports:     0       0       0       0       0       0       0
GPSPW reports:     0       0       0       0       0       0       0
GPSRF reports:     0       0       0       0       0       0       0
GPSEP reports:     0       0       0       0       0       0       0
SSMT1 reports:     0       0       0       0       0       0       0
SSMT2 reports:     0       0       0       0       0       0       0
SSMI  reports:     0       0       0       0       0       0       0
TOVS  reports:     0       0       0       0       0       0       0
QSCAT reports:     0       0       0       0       0       0       0
RPOFL reports:     0       0       0       0       0       0       0
AIRS  reports:     0       0       0       0       0       0       0
TAMDAR reports:     0       0       0       0       0       0       0
OTHER reports:     0       0       0       0       0       0       0
Total reports: 91392      35   89717     155     594       1     890

------------------------------------------------------------------------------
REDUCE QC FROM 7 TO 2 DIGITS:



slot= 1  time_min, time_fg, time_max:  2013-03-01_11:00:00  2013-03-01_12:00:00  2013-03-01_12:59:59
------------------------------------------------------------------------------
OBSERVATIONS FOR OUTPUT SLOT 01

SYNOP reports:   765
SHIPS reports:    21
BUOY  reports:     0
BUGUS reports:     0
METAR reports:   100
PILOT reports:     0
SOUND reports:     0
AMDAR reports:     0
SATEM reports:     0
SATOB reports:     0
AIREP reports:     0
TAMDAR reports:     0
GPSPW reports:     0
GPSZD reports:     0
GPSRF reports:     0
GPSEP reports:     0
SSMT1 reports:     0
SSMT2 reports:     0
SSMI  reports:     0
TOVS  reports:     0
QSCAT reports:     0
PROFL reports:     0
AIRST reports:     0
OTHER reports:     0
Total reports:   886 =    886 single +      0 multi-level reports.

------------------------------------------------------------------------------
Write 3DVAR 3DVAR GTS observations in file obs_gts_2013-03-01_12:00:00.3DVAR (wrfvar version 3.0)

Wrote    1751 lines of data in file: obs_gts_2013-03-01_12:00:00.3DVAR

No SSMI observations available.

99999

又是那隻貓

问题是什么

andrewsoong

恕我浅薄，没看懂，帮你顶~

lee309

andrewsoong 发表于 2014-12-30 19:50
恕我浅薄，没看懂，帮你顶~

谢版主，确实挺复杂的，说不清楚

小虹

1." 开头的Namelist logging not found in namelist.input. Using registry defaults for variables in logging."

这些没关系，大意就是说你的namelist里面有的选项没有指明，所以用默认的。

2. "有很多数据都是0很奇怪，如psi, chi_u, t_u and ps_u等，这是怎么回事，是因为观测数据中没有这些变量么？"

我也不是很懂，但是我看了我的结果，是有数据的。例如：
Vertical truncation for psi    =  14(  99.00%)

Vertical truncation for chi_u  =  19(  99.00%)

Vertical truncation for t_u    =  28(  99.00%)

Vertical truncation for rh     =  25(  99.00%)

3. ”最后的Got unknown obs_type string又是什么意思，貌似重点就在这，哪位做过同化懂这方面的指点下？"

我猜这个可能是你的观测文件写的不对，你再重新检查一下的观测文件。

lee309

小虹 发表于 2014-12-31 11:26
1." 开头的Namelist logging not found in namelist.input. Using registry defaults for variables in log ...

嗯，我正在重新生成观测文件，只是namelist.obsproc中的有些设置不知道什么意思。
我从2013-03-01-12预报到2013-03-02-00,单层嵌套，网格为60X61,精度dx和dy都是30km，evert=30分析的时刻是初始时刻2013-03-01-12（分析2013-03-02-00总是出错，暂时不管它）
你帮我看下我的namelist.obsproc有没有问题吧，谢谢！！！！！！！！！！！
&record1
obs_gts_filename = '/usr/WRF/test/WRFDA/TestData/SURFACE_OBS:2013030112',
obs_err_filename = 'obserr.txt',
/

&record2
time_window_min  = '2013-03-01_11:00:00',
time_analysis    = '2013-03-01_12:00:00',
time_window_max  = '2013-03-01_13:00:00',
/

&record3
max_number_of_obs        = 400000,
fatal_if_exceed_max_obs  = .TRUE.,
/

&record4
qc_test_vert_consistency = .TRUE.,
qc_test_convective_adj   = .TRUE.,
qc_test_above_lid        = .TRUE.,
remove_above_lid         = .TRUE.,
domain_check_h           = .true.,
Thining_SATOB            = false,
Thining_SSMI             = false,
Thining_QSCAT            = false,
/

&record5
print_gts_read           = .TRUE.,
print_gpspw_read         = .TRUE.,
print_recoverp           = .TRUE.,
print_duplicate_loc      = .TRUE.,
print_duplicate_time     = .TRUE.,
print_recoverh           = .TRUE.,
print_qc_vert            = .TRUE.,
print_qc_conv            = .TRUE.,
print_qc_lid             = .TRUE.,
print_uncomplete         = .TRUE.,
/

&record6
ptop =  1000.0,
base_pres       = 100000.0,
base_temp       = 290.0,
base_lapse      = 50.0,
base_strat_temp = 215.0,
base_tropo_pres = 20000.0
/

&record7
IPROJ = 1,
PHIC  = 28.0,
XLONC = 122.0,
TRUELAT1= 30.0,
TRUELAT2= 60.0,
MOAD_CEN_LAT = 28.0,
STANDARD_LON = 122.0,
/

&record8
IDD    =   1,
MAXNES =   1,
NESTIX =  61,  200,  136,  181,  211,
NESTJX =  60,  200,  181,  196,  211,
DIS    =  30,  10.,  3.3,  1.1,  1.1,
NUMC   =    1,    1,   2,     3,    4,
NESTI  =    1,   40,  28,    35,   45,
NESTJ  =    1,   60,  25,    65,   55,
/

&record9
PREPBUFR_OUTPUT_FILENAME = 'prepbufr_output_filename',
PREPBUFR_TABLE_FILENAME = 'prepbufr_table_filename',
OUTPUT_OB_FORMAT = 2
use_for          = '3DVAR',
num_slots_past   = 1,
num_slots_ahead  = 1,
write_synop = .true.,
write_ship  = .true.,
write_metar = .true.,
write_buoy  = .true.,
write_pilot = .true.,
write_sound = .true.,
write_amdar = .true.,
write_satem = .true.,
write_satob = .true.,
write_airep = .true.,
write_gpspw = .true.,
write_gpsztd= .true.,
write_gpsref= .true.,
write_gpseph= .true.,
write_ssmt1 = .true.,
write_ssmt2 = .true.,
write_ssmi  = .true.,
write_tovs  = .true.,
write_qscat = .true.,
write_profl = .true.,
write_bogus = .true.,
write_airs  = .true.,
/

lee309

又是那隻貓 发表于 2014-12-30 19:46
问题是什么

1. 开头的Namelist logging not found in namelist.input. Using registry defaults for variables in logging.不知道是什么意思？
2. 有很多数据都是0很奇怪，如psi, chi_u, t_u and ps_u等，这是怎么回事，是因为观测数据中没有这些变量么？
3. 最后的Got unknown obs_type string又是什么意思，貌似重点就在这，哪位做过同化懂这方面的指点下？
虽然observation文件生成了，但是运行obsproc.exe也有点error和警告

夭夭坂田淘

想问一下，WRFDA是什么。。。。

lee309

夭夭坂田淘 发表于 2015-8-9 21:20
想问一下，WRFDA是什么。。。。

WRF同化系统

夭夭坂田淘

lee309 发表于 2015-10-10 21:50
WRF同化系统

噢~~好的，谢谢啦