- 积分
- 221
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-9-12
- 最后登录
- 1970-1-1
|
发表于 2025-2-20 09:01:43
|
显示全部楼层
还有想请教楼主一个WRFDA的问题。就是我在做地面观测数据前处理的时候,下载了区域范围内的数据处理成littleR格式的,执行./obsproc.exe后显示所有数据在区域外,但是数据显示都读入进去了,检查了namelist.obsproc和namelist.wps的区域设置是相同的,想问下楼主可能是哪里出问题了呢?
namelist.wps如下:
&share
wrf_core = 'ARW'
max_dom = 2
start_date = '2023-07-09_00:00:00','2023-07-09_00:00:00'
end_date = '2023-07-19_00:00:00','2023-07-19_00:00:00'
interval_seconds = 21600
io_form_geogrid = 2
debug_level = 0
/
&geogrid
parent_id = 1, 1
parent_grid_ratio = 1, 3
i_parent_start = 1, 160
j_parent_start = 1, 130
e_we = 304, 304
e_sn = 304, 304
geog_data_res = 'default','default'
dx = 9000
dy = 9000
map_proj = 'lambert'
ref_lat = 38
ref_lon = 89
truelat1 = 39
truelat2 = 41
stand_lon = 89
geog_data_path = '/lza/g1/xbrcnwp/MODEL/WPS_GEOG'
/
&ungrib
out_format = 'WPS'
prefix = 'FILE'
/
&metgrid
fg_name = 'FILE','GFSFILE'
io_form_metgrid = 2
/
namelist.obsproc如下:
&record1
obs_gts_filename = 'fort.2023070906',
obs_err_filename = 'obserr.txt',
FG_FORMAT = 'LITTLE_R',
gts_from_mmm_archive = .false.,
/
&record2
time_window_min = '2023-07-09_05:45:00',
time_analysis = '2023-07-09_06:00:00',
time_window_max = '2023-07-09_06:15: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 = .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 = 5000.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 = 38,
XLONC = 89,
TRUELAT1= 39,
TRUELAT2= 41,
MOAD_CEN_LAT = 38,
STANDARD_LON = 89,
/
&record8
IDD = 1,
MAXNES = 1,
NESTIX = 304,
NESTJX = 304,
DIS = 9,
NUMC = 1,
NESTI = 1,
NESTJ = 1,
/
&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 = .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.,
/
&RECORD10
WIND_SD = T,
WIND_SD_SYNOP = T,
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
/
|
|