Error when filtering obs from littleR file with obsproc.exe

Problems with data in 3DVAR.

Error when filtering obs from littleR file with obsproc.exe

Postby jiaojiu » Sun Mar 24, 2019 11:10 pm

I am filtering synop observations from littleR, and set my namelist.obsproc
namelist.obsproc.txt.
Code: Select all
&record1
 ! obs_gts_filename = 'obsHbEnv',
 obs_gts_filename = 'obsWhMet',
 fg_format        = 'WRF',
 obs_err_filename = 'obserr.txt',
 gts_from_mmm_archive = .true.,
 ! first_guess_file = 'fg',
/

&record2
 time_window_min  = '2017-11-01_09:00:00',
 time_analysis    = '2017-11-01_12:00:00',
 time_window_max  = '2017-11-01_15: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         = .false.,
 domain_check_h           = .true.,
 Thining_SATOB            = .false.,
 Thining_SSMI             = .false.,
 Thining_QSCAT            = .false.,
 calc_psfc_from_qnh       = .true.,
/

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

&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           = 30.02,
 XLONC          = 112.19,
 TRUELAT1      = 30.0,
 TRUELAT2      = 60.0,
 MOAD_CEN_LAT    = 30.02,
 STANDARD_LON    = 112.19,
/

&record8
 IDD    =   1,
 MAXNES =   2,
 NESTIX =  120,  106,
 NESTJX =   96,  82,
 DIS    =   27,  9,
 NUMC   =    1,    1,
 NESTI  =    1,   43,
 NESTJ  =    1,   35,
 /

&record9
 PREPBUFR_OUTPUT_FILENAME = 'prepbufr_output',
 PREPBUFR_TABLE_FILENAME = 'prepbufr_table_filename',
 OUTPUT_OB_FORMAT = 2,                ! 1,PREPBUFR OBS file; 2,ASCII OBS file; 3,Both PREPBUFR and ASCII OBS file.
 use_for          = '3DVAR',         ! 3DVAR, FGAT, 4DVAR
 num_slots_past   = 3,
 num_slots_ahead  = 3,
 write_synop = .true.,
 write_ship  = .false.,
 write_metar = .false.,
 write_buoy  = .false.,
 write_pilot = .false.,
 write_sound = .false.,
 write_amdar = .false.,
 write_satem = .false.,
 write_satob = .false.,
 write_airep = .false.,
 write_gpspw = .false.,
 write_gpsztd= .false.,
 write_gpsref= .false.,
 write_gpseph= .false.,
 write_ssmt1 = .false.,
 write_ssmt2 = .false.,
 write_ssmi  = .false.,
 write_tovs  = .false.,
 write_qscat = .false.,
 write_profl = .false.,
 write_bogus = .false.,
 write_airs  = .false.,
 /

 &record10
 wind_sd        = .true.
 WIND_SD_SYNOP   = .true.
 /

But I found that only observations at the analysis time(2017-11-01_12:00:00) was processed, not all the observations in the time window of 6 hours around the analysis time. This puzzled me.
As show in the file(Attached),
Code: Select all
Found duplicated stations:
Station 1 name and ID =   FM-12 SYNOP                                30.170  108.560          F
Station 2 name and ID =   FM-12 SYNOP                                30.170  108.560          F
Analysis  time = 2017-11-01_12:00:00
Station 1 time = 2017-11-01_09:00:00 = ta -  10800s
Station 2 time = 2017-11-01_10:00:00 = ta -   7200s
Keep station 2 and reject station 1.

Found duplicated stations:
Station 1 name and ID =   FM-12 SYNOP                                30.170  108.560          F
Station 2 name and ID =   FM-12 SYNOP                                30.170  108.560          F
Analysis  time = 2017-11-01_12:00:00
Station 1 time = 2017-11-01_10:00:00 = ta -   7200s
Station 2 time = 2017-11-01_11:00:00 = ta -   3600s
Keep station 2 and reject station 1.

Found duplicated stations:
Station 1 name and ID =   FM-12 SYNOP                                30.170  108.560          F
Station 2 name and ID =   FM-12 SYNOP                                30.170  108.560          F
Analysis  time = 2017-11-01_12:00:00
Station 1 time = 2017-11-01_11:00:00 = ta -   3600s
Station 2 time = 2017-11-01_12:00:00 = ta +      0s
Keep station 2 and reject station 1.
A large amount of observations are considerd as duplicated and removed.
SO anybody could point out my error in the namelist.obsproc, or the understanding of obsproc.exe program.
Thanks in advance!
jiaojiu
 
Posts: 3
Joined: Sat Aug 05, 2017 10:00 pm

Return to Data Issues

Who is online

Users browsing this forum: No registered users and 3 guests

cron