Hello WRFDA Community,
I am working on 4DVAR data assimilation using WRFDA version 4.6. I have successfully run real.exe with GFS data to generate the wrfinput_d01 and wrfbdy_d01 files for my domain.
For observational input, I am using the NCEP ADP Global Surface Observational Weather Data from this dataset for my case study.
Although obsproc runs without any errors, I am unable to obtain the expected output files like
obs_gts_2025-10-22_09:00:00.4DVAR and obs_gts_2025-10-22_12:00:00.4DVAR. This prevents me from proceeding with data assimilation.
I have attached my namelist.obsproc for reference in case it helps identify the issue.
Has anyone encountered this problem or can provide guidance on successfully generating obs_gts files with surface observations from this dataset using WRFDA 4.6? Any advice or examples of working configurations would be greatly appreciated.
iam atatching the log file also
Thank you in advance!
Paparao
&record1
obs_gts_filename = 'SURFACE_OBS:2025102212',
obs_err_filename = 'obserr.txt',
gts_from_mmm_archive = .false.,
/
&record2
time_window_min = '2025-10-22_09:00:00',
time_analysis = '2025-10-22_12:00:00',
time_window_max = '2025-10-22_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 = .false.,
Thining_SATOB = .false.,
Thining_SSMI = .false.,
Thining_QSCAT = .false.,
calc_psfc_from_qnh = .true.,
/
&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 = 3,
PHIC = 22.50,
XLONC = 82.0,
TRUELAT1= 0,
TRUELAT2= 0,
MOAD_CEN_LAT = 22.50,
STANDARD_LON = 82.0,
/
&record8
IDD = 1,
MAXNES = 1,
NESTIX = 425,
NESTJX = 409,
DIS = 9,
NUMC = 1,
NESTI = 1,
NESTJ = 1,
/
&record9
PREPBUFR_OUTPUT_FILENAME = '',
PREPBUFR_TABLE_FILENAME = '',
OUTPUT_OB_FORMAT = 2
use_for = '4DVAR',
num_slots_past = 3,
num_slots_ahead = 3,
write_synop = .true.,
write_metar = .true.,
/
&record10
/
I am working on 4DVAR data assimilation using WRFDA version 4.6. I have successfully run real.exe with GFS data to generate the wrfinput_d01 and wrfbdy_d01 files for my domain.
For observational input, I am using the NCEP ADP Global Surface Observational Weather Data from this dataset for my case study.
Although obsproc runs without any errors, I am unable to obtain the expected output files like
obs_gts_2025-10-22_09:00:00.4DVAR and obs_gts_2025-10-22_12:00:00.4DVAR. This prevents me from proceeding with data assimilation.
I have attached my namelist.obsproc for reference in case it helps identify the issue.
Has anyone encountered this problem or can provide guidance on successfully generating obs_gts files with surface observations from this dataset using WRFDA 4.6? Any advice or examples of working configurations would be greatly appreciated.
iam atatching the log file also
Thank you in advance!
Paparao
&record1
obs_gts_filename = 'SURFACE_OBS:2025102212',
obs_err_filename = 'obserr.txt',
gts_from_mmm_archive = .false.,
/
&record2
time_window_min = '2025-10-22_09:00:00',
time_analysis = '2025-10-22_12:00:00',
time_window_max = '2025-10-22_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 = .false.,
Thining_SATOB = .false.,
Thining_SSMI = .false.,
Thining_QSCAT = .false.,
calc_psfc_from_qnh = .true.,
/
&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 = 3,
PHIC = 22.50,
XLONC = 82.0,
TRUELAT1= 0,
TRUELAT2= 0,
MOAD_CEN_LAT = 22.50,
STANDARD_LON = 82.0,
/
&record8
IDD = 1,
MAXNES = 1,
NESTIX = 425,
NESTJX = 409,
DIS = 9,
NUMC = 1,
NESTI = 1,
NESTJ = 1,
/
&record9
PREPBUFR_OUTPUT_FILENAME = '',
PREPBUFR_TABLE_FILENAME = '',
OUTPUT_OB_FORMAT = 2
use_for = '4DVAR',
num_slots_past = 3,
num_slots_ahead = 3,
write_synop = .true.,
write_metar = .true.,
/
&record10
/