求解WRFDA同化使用obsproc.exe提取有效观测的问题

SAchinb

自行制作了littleR格式的逐小时synop观测数据，在使用obsproc.exe过滤时，设置了如下的namelist.obsproc（见第一个代码块），但只是提取了同化时刻（2017-11-01_12:00:00）的观测数据，有很多其他数据被过滤掉了，（部分过滤信息见第二个代码块），请问如何设置namelist来保留观测时刻临近时间窗（±3h）的观测值，或者请大家帮忙解惑是否因为我对OBSPROC.exe功能的理解有误。

namelist.obsproc

1. &record1
   
2. obs_gts_filename = 'obs_2017-11-01_12',
   
3. fg_format        = 'WRF',
   
4. obs_err_filename = 'obserr.txt',
   
5. gts_from_mmm_archive = .true.,
   
6. ! first_guess_file = 'fg',
   
7. /
   
8. 
   
9. &record2
   
10. time_window_min  = '2017-11-01_09:00:00',
    
11. time_analysis    = '2017-11-01_12:00:00',
    
12. time_window_max  = '2017-11-01_15:00:00',
    
13. /
    
14. 
    
15. &record3
    
16. max_number_of_obs        = 400000,
    
17. fatal_if_exceed_max_obs  = .TRUE.,
    
18. /
    
19. 
    
20. &record4
    
21. qc_test_vert_consistency = .TRUE.,
    
22. qc_test_convective_adj   = .TRUE.,
    
23. qc_test_above_lid        = .TRUE.,
    
24. remove_above_lid         = .false.,
    
25. domain_check_h           = .true.,
    
26. Thining_SATOB            = .false.,
    
27. Thining_SSMI             = .false.,
    
28. Thining_QSCAT            = .false.,
    
29. calc_psfc_from_qnh       = .true.,
    
30. /
    
31. 
    
32. &record5
    
33. print_gts_read           = .TRUE.,
    
34. print_gpspw_read         = .false.,
    
35. print_recoverp           = .TRUE.,
    
36. print_duplicate_loc      = .TRUE.,
    
37. print_duplicate_time     = .TRUE.,
    
38. print_recoverh           = .false.,
    
39. print_qc_vert            = .false.,
    
40. print_qc_conv            = .false.,
    
41. print_qc_lid             = .false.,
    
42. print_uncomplete         = .false.,
    
43. /
    
44. 
    
45. &record6
    
46. ptop             = 1000.0,
    
47. base_pres          = 100000.0,
    
48. base_temp          = 290.0,
    
49. base_lapse         = 50.0,
    
50. base_strat_temp    = 215.0,
    
51. base_tropo_pres    = 20000.0,
    
52. /
    
53. 
    
54. &record7
    
55. IPROJ          = 1,
    
56. PHIC           = 30.02,
    
57. XLONC          = 112.19,
    
58. TRUELAT1      = 30.0,
    
59. TRUELAT2      = 60.0,
    
60. MOAD_CEN_LAT    = 30.02,
    
61. STANDARD_LON    = 112.19,
    
62. /
    
63. 
    
64. &record8
    
65. IDD    =   1,
    
66. MAXNES =   2,
    
67. NESTIX =  120,  106,
    
68. NESTJX =   96,  82,
    
69. DIS    =   27,  9,
    
70. NUMC   =    1,    1,
    
71. NESTI  =    1,   43,
    
72. NESTJ  =    1,   35,
    
73. /
    
74. 
    
75. &record9
    
76. PREPBUFR_OUTPUT_FILENAME = 'prepbufr_output',
    
77. PREPBUFR_TABLE_FILENAME = 'prepbufr_table_filename',
    
78. OUTPUT_OB_FORMAT = 2,                ! 1,PREPBUFR OBS file; 2,ASCII OBS file; 3,Both PREPBUFR and ASCII OBS file.
    
79. use_for          = '3DVAR',         ! 3DVAR, FGAT, 4DVAR
    
80. num_slots_past   = 3,
    
81. num_slots_ahead  = 3,
    
82. write_synop = .true.,
    
83. write_ship  = .false.,
    
84. write_metar = .false.,
    
85. write_buoy  = .false.,
    
86. write_pilot = .false.,
    
87. write_sound = .false.,
    
88. write_amdar = .false.,
    
89. write_satem = .false.,
    
90. write_satob = .false.,
    
91. write_airep = .false.,
    
92. write_gpspw = .false.,
    
93. write_gpsztd= .false.,
    
94. write_gpsref= .false.,
    
95. write_gpseph= .false.,
    
96. write_ssmt1 = .false.,
    
97. write_ssmt2 = .false.,
    
98. write_ssmi  = .false.,
    
99. write_tovs  = .false.,
    
100. write_qscat = .false.,
     
101. write_profl = .false.,
     
102. write_bogus = .false.,
     
103. write_airs  = .false.,
     
104. /
     
105. 
     
106. &record10
     
107. wind_sd        = .true.
     
108. WIND_SD_SYNOP   = .true.
     
109. /

复制代码

过滤信息（部分，提取自obs_duplicate_time.diag_2017-11-01_12_00_00）,提示有大量观测点重复被剔除

1. Found duplicated stations:
   
2. Station 1 name and ID =   FM-12 SYNOP                                30.170  108.560          F
   
3. Station 2 name and ID =   FM-12 SYNOP                                30.170  108.560          F
   
4. Analysis  time = 2017-11-01_12:00:00
   
5. Station 1 time = 2017-11-01_09:00:00 = ta -  10800s
   
6. Station 2 time = 2017-11-01_10:00:00 = ta -   7200s
   
7. Keep station 2 and reject station 1.
   
8. 
   
9. Found duplicated stations:
   
10. Station 1 name and ID =   FM-12 SYNOP                                30.170  108.560          F
    
11. Station 2 name and ID =   FM-12 SYNOP                                30.170  108.560          F
    
12. Analysis  time = 2017-11-01_12:00:00
    
13. Station 1 time = 2017-11-01_10:00:00 = ta -   7200s
    
14. Station 2 time = 2017-11-01_11:00:00 = ta -   3600s
    
15. Keep station 2 and reject station 1.
    
16. 
    
17. Found duplicated stations:
    
18. Station 1 name and ID =   FM-12 SYNOP                                30.170  108.560          F
    
19. Station 2 name and ID =   FM-12 SYNOP                                30.170  108.560          F
    
20. Analysis  time = 2017-11-01_12:00:00
    
21. Station 1 time = 2017-11-01_11:00:00 = ta -   3600s
    
22. Station 2 time = 2017-11-01_12:00:00 = ta +      0s
    
23. Keep station 2 and reject station 1.
    
24. 
    
25. Found duplicated stations:
    
26. Station 1 name and ID =   FM-12 SYNOP                                30.170  108.560          F
    
27. Station 2 name and ID =   FM-12 SYNOP                                30.170  108.560          F
    
28. Analysis  time = 2017-11-01_12:00:00
    
29. Station 1 time = 2017-11-01_12:00:00 = ta +      0s
    
30. Station 2 time = 2017-11-01_13:00:00 = ta +   3600s
    
31. Keep station 1 and reject station 2.

复制代码

一两

请问您解决了吗

任性小王子

楼主，是不是范围的设置问题呐，我跑过案例，没出现过错误。