關於WRFDA的obsproc.exe運行時的錯

max533

我從ncep網站上下載已經轉換成little_R格式的觀測資料，想利用obsproc.exe轉換成obs_gts_xxxx-xx-xx_xx:xx:xx.3DVAR
發現它出現如下列的錯誤訊息，但不知道從何修正，想各位高手幫忙一起找出原因，解決問題。

我的編譯器是 gcc + gfortran

At line 123 of file module_date.f
Fortran runtime error: Bad value during integer read

obsproc.log內容如下：

1. At line 123 of file module_date.f
   
2. Fortran runtime error: Bad value during integer read
   
3. 
   
4. READ NAMELIST FILE: namelist.obsproc
   
5. ------------------
   
6. &RECORD1
   
7. OBS_GTS_FILENAME="obs_2009080706                                                                  ",
   
8. OBS_ERR_FILENAME="obserr.txt                                                                      ",
   
9. FG_FORMAT="WRF",
   
10. /
    
11. &RECORD2
    
12. TIME_WINDOW_MIN="2009-08-07_05:00:00",
    
13. TIME_ANALYSIS="2009-08-07_06:00:00",
    
14. TIME_WINDOW_MAX="2009-08-07_07:00:00",
    
15. /
    
16. &RECORD3
    
17. MAX_NUMBER_OF_OBS=     400000,
    
18. FATAL_IF_EXCEED_MAX_OBS=T,
    
19. /
    
20. &RECORD4
    
21. QC_TEST_VERT_CONSISTENCY=T,
    
22. QC_TEST_CONVECTIVE_ADJ=T,
    
23. QC_TEST_ABOVE_LID=T,
    
24. REMOVE_ABOVE_LID=T,
    
25. DOMAIN_CHECK_H=T,
    
26. THINING_SATOB=F,
    
27. THINING_SSMI=F,
    
28. THINING_QSCAT=F,
    
29. /
    
30. &RECORD5
    
31. PRINT_GTS_READ=T,
    
32. PRINT_GPSPW_READ=T,
    
33. PRINT_RECOVERP=T,
    
34. PRINT_DUPLICATE_LOC=T,
    
35. PRINT_DUPLICATE_TIME=T,
    
36. PRINT_RECOVERH=T,
    
37. PRINT_QC_VERT=T,
    
38. PRINT_QC_CONV=T,
    
39. PRINT_QC_LID=T,
    
40. PRINT_UNCOMPLETE=T,
    
41. /
    
42. &RECORD6
    
43. PTOP=  1000.0000    ,
    
44. PS0=  100000.00    ,
    
45. TS0=  290.00000    ,
    
46. TLP=  50.000000    ,
    
47. PIS0=  20000.000    ,
    
48. TIS0=  215.00000    ,
    
49. BASE_PRES=  100000.00    ,
    
50. BASE_TEMP=  290.00000    ,
    
51. BASE_LAPSE=  50.000000    ,
    
52. BASE_STRAT_TEMP=  215.00000    ,
    
53. BASE_TROPO_PRES=  20000.000    ,
    
54. /
    
55. &RECORD7
    
56. IPROJ=          1,
    
57. PHIC=  40.000011    ,
    
58. XLONC= -95.000000    ,
    
59. TRUELAT1=  30.000000    ,
    
60. TRUELAT2=  60.000000    ,
    
61. MOAD_CEN_LAT=  40.000011    ,
    
62. STANDARD_LON= -95.000000    ,
    
63. /
    
64. &RECORD8
    
65. IDD=          1,
    
66. MAXNES=          1,
    
67. NESTIX=         60,        200,        136,        181,        211,
    
68.   5*0          ,
    
69. NESTJX=         90,        200,        181,        196,        211,
    
70.   5*0          ,
    
71. DIS=  60.000000    ,  10.000000    ,  3.3000000    , 2*1.1000000      , 5*0.0000000      ,
    
72. 
    
73. NUMC= 2*1          ,          2,          3,          4, 5*0          ,
    
74. 
    
75. NESTI=          1,         40,         28,         35,         45,
    
76.   5*0          ,
    
77. NESTJ=          1,         60,         25,         65,         55,
    
78.   5*0          ,
    
79. /
    
80. &RECORD9
    
81. PREPBUFR_OUTPUT_FILENAME="prepbufr_output_filename                                                        ",
    
82. PREPBUFR_TABLE_FILENAME="prepbufr_table_filename                                                         ",
    
83. OUTPUT_OB_FORMAT=          2,
    
84. USE_FOR="3DVAR",
    
85. NUM_SLOTS_PAST=          3,
    
86. NUM_SLOTS_AHEAD=          3,
    
87. WRITE_SYNOP=T,
    
88. WRITE_SHIP=T,
    
89. WRITE_METAR=T,
    
90. WRITE_BUOY=T,
    
91. WRITE_PILOT=T,
    
92. WRITE_SOUND=T,
    
93. WRITE_AMDAR=T,
    
94. WRITE_SATEM=T,
    
95. WRITE_SATOB=T,
    
96. WRITE_AIREP=T,
    
97. WRITE_GPSPW=T,
    
98. WRITE_GPSZTD=T,
    
99. WRITE_GPSREF=T,
    
100. WRITE_GPSEPH=T,
     
101. WRITE_SSMT1=T,
     
102. WRITE_SSMT2=T,
     
103. WRITE_SSMI=T,
     
104. WRITE_TOVS=T,
     
105. WRITE_QSCAT=T,
     
106. WRITE_PROFL=T,
     
107. WRITE_BOGUS=T,
     
108. WRITE_AIRS=T,
     
109. WRITE_TAMDAR=T,
     
110. /
     
111. COARSE GRID CENTER IS AT X =      0.0 KM AND Y =  -6616.7 KM.   DIS(1)=    60.000
     
112. ------------------------------------------------------------------------------
     
113. READ GTS OBSERVATIONS IN FILE obs_2009080706
     
114. 
     
115. Diagnostics in file obs_gts_read.diag

复制代码

先在這邊謝謝大家。

max533

後來自行解決，解決方法是去上WRFDA官方網站，在V3.6.1版時有個obsproc.tar的已知錯誤和修正，把它下載下來並重新compile即可。

白兔糖0217

max533 发表于 2015-7-1 17:58
後來自行解決，解決方法是去上WRFDA官方網站，在V3.6.1版時有個obsproc.tar的已知錯誤和修正，把它下載下來 ...

楼主点赞！！！对了，3.5的有需要这样做吗？我写成了prepbufr, 应该可以直接同化的吧？

白兔糖0217

max533 发表于 2015-7-1 17:58
後來自行解決，解決方法是去上WRFDA官方網站，在V3.6.1版時有個obsproc.tar的已知錯誤和修正，把它下載下來 ...

我用的3.7,也有这样的问题。。。看官网说，obsproc的一些问题解决了啊。嗨，甩个连接呢，我看看是不是找对地方了

lucky_girl

楼主，你知道怎么把资料转换为LITTLE_R格式吗？

白兔糖0217

有木有人跟我一样，到现在像样的工作还没找到的

Sunnystone24

楼主请问只是用下的V3.6.1的官网修改包修改重新编译了就可以了吗，我的3.4.1版本 也出现了同样的问题。。。

Sunnystone24

刚才下载了也是楼主说的再v3.6.1版本处的那个修补的tar，然后重新编译后解决了 读 little_R报错的问题。。。

任性小王子

白兔糖0217 发表于 2015-12-21 20:46
有木有人跟我一样，到现在像样的工作还没找到的

慢慢找吧，也在找工作呢