- 积分
- 8270
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-12-27
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
想把NCEP ADP的高空资料同化到case的模拟中,修改了namelist.obsproc后,./obsproc.ex出了问题。
namelist.input 如下:
&time_control
run_days = 0,
run_hours = 30,
run_minutes = 0,
run_seconds = 0,
start_year = 2010, 2010,
start_month = 8, 8,
start_day = 29, 29,
start_hour = 0, 0,
start_minute = 00, 00,
start_second = 00, 00,
end_year = 2010, 2010,
end_month = 8, 8,
end_day = 30, 30,
end_hour = 6, 6,
end_minute = 00, 00,
end_second = 00, 00,
interval_seconds = 21600,
input_from_file = .true., .true.,
history_interval = 180, 60,
frames_per_outfile = 1000, 1000,
restart = .false.,
restart_interval = 5000,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
/
&domains
time_step = 180,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 100, 193,
e_sn = 79, 187,
e_vert = 28, 28,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 27000, 9000,
dy = 27000, 9000,
grid_id = 1, 2,
parent_id = 1, 1,
i_parent_start = 1, 16,
j_parent_start = 1, 9,
parent_grid_ratio = 1, 3,
parent_time_step_ratio = 1, 3,
feedback = 1,
smooth_option = 0,
/
&physics
mp_physics = 8, 8,
ra_lw_physics = 1, 1,
ra_sw_physics = 2, 2,
radt = 30, 30,
sf_sfclay_physics = 1, 1,
sf_surface_physics = 2, 2,
bl_pbl_physics = 1, 1,
bldt = 0, 0,
cu_physics = 3, 3,
cudt = 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
/
&fdda
/
&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0,
diff_6th_factor = 0.12, 0.12,
base_temp = 290.,
damp_opt = 0,
zdamp = 5000., 5000.,
dampcoef = 0.2, 0.2,
khdif = 0, 0,
kvdif = 0, 0,
non_hydrostatic = .true., .true.,
moist_adv_opt = 1, 1,
scalar_adv_opt = 1, 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,
nested = .false., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
-----------------------------------------------------------------------------------------------------------------------
namelist.obsproc 如下:
&record1
obs_gts_filename = '/home/mode/data/2010.08.29.ADP.Sounding/OBS_2010083000',
obs_err_filename = 'obserr.txt',
fg_format = 'WRF'
/
&record2
time_window_min = '2010-08-29_23:00:00',
time_analysis = '2010-08-30_00:00:00',
time_window_max = '2010-08-30_01: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 = .TRUE.,
domain_check_h = .true.,
Thining_SATOB = false,
Thining_SSMI = false,
Thining_QSCAT = false,
/
&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 = 28.0,
XLONC = 120.0,
TRUELAT1= 30.0,
TRUELAT2= 60.0,
MOAD_CEN_LAT = 28.0,
STANDARD_LON = 120.0,
/
&record8
IDD = 1,
MAXNES = 2,
NESTIX = 100, 193, 136, 181, 211,
NESTJX = 79, 187, 181, 196, 211,
DIS = 27, 9, 3.3, 1.1, 1.1,
NUMC = 1, 1, 2, 3, 4,
NESTI = 1, 16, 28, 35, 45,
NESTJ = 1, 9, 25, 65, 55,
/
&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.,
/
--------------------------------------------------------------------------------------------------------------------------------
报错信息如下:
READ NAMELIST FILE: namelist.obsproc
------------------
&RECORD1
OBS_GTS_FILENAME=/home/mode/data/2010.08.29.ADP.Sounding/OBS_2010083000
OBS_ERR_FILENAME=obserr.txt
FG_FORMAT=WRF
GTS_FROM_NCAR_ARCHIVE=F,
/
&RECORD2
TIME_WINDOW_MIN=2010-08-29_23:00:00
TIME_ANALYSIS=2010-08-30_00:00:00
TIME_WINDOW_MAX=2010-08-30_01:00:00
/
&RECORD3
MAX_NUMBER_OF_OBS= 400000,
FATAL_IF_EXCEED_MAX_OBS=T,
/
&RECORD4
QC_TEST_VERT_CONSISTENCY=T,
QC_TEST_CONVECTIVE_ADJ=T,
QC_TEST_ABOVE_LID=T,
REMOVE_ABOVE_LID=T,
DOMAIN_CHECK_H=T,
THINING_SATOB=F,
THINING_SSMI=F,
THINING_QSCAT=F,
/
&RECORD5
PRINT_GTS_READ=T,
PRINT_GPSPW_READ=T,
PRINT_RECOVERP=T,
PRINT_DUPLICATE_LOC=T,
PRINT_DUPLICATE_TIME=T,
PRINT_RECOVERH=T,
PRINT_QC_VERT=T,
PRINT_QC_CONV=T,
PRINT_QC_LID=T,
PRINT_UNCOMPLETE=T,
/
&RECORD6
PTOP= 5000.00000 ,
PS0= 100000.000 ,
TS0= 290.000000 ,
TLP= 50.0000000 ,
PIS0= 20000.0000 ,
TIS0= 215.000000 ,
BASE_PRES= 100000.000 ,
BASE_TEMP= 290.000000 ,
BASE_LAPSE= 50.0000000 ,
BASE_STRAT_TEMP= 215.000000 ,
BASE_TROPO_PRES= 20000.0000 ,
/
&RECORD7
IPROJ= 1,
PHIC= 28.0000000 ,
XLONC= 120.000000 ,
TRUELAT1= 30.0000000 ,
TRUELAT2= 60.0000000 ,
MOAD_CEN_LAT= 28.0000000 ,
STANDARD_LON= 120.000000 ,
/
&RECORD8
IDD= 1,
MAXNES= 2,
NESTIX= 100, 193, 136, 181, 211,,
5*0 ,
NESTJX= 79, 187, 181, 196, 211,,
5*0 ,
DIS= 27.0000000 , 9.00000000 , 3.29999995 , 2*1.10000002 , 5*0.00000000 ,,
NUMC= 2*1 , 2, 3, 4, 5*0 ,,
NESTI= 1, 16, 28, 35, 45,,
5*0 ,
NESTJ= 1, 9, 25, 65, 55,,
5*0 ,
/
&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=T,
WRITE_SHIP=T,
WRITE_METAR=T,
WRITE_BUOY=T,
WRITE_PILOT=T,
WRITE_SOUND=T,
WRITE_AMDAR=T,
WRITE_SATEM=T,
WRITE_SATOB=T,
WRITE_AIREP=T,
WRITE_GPSPW=T,
WRITE_GPSZTD=T,
WRITE_GPSREF=T,
WRITE_GPSEPH=T,
WRITE_SSMT1=T,
WRITE_SSMT2=T,
WRITE_SSMI=T,
WRITE_TOVS=T,
WRITE_QSCAT=T,
WRITE_PROFL=T,
WRITE_BOGUS=T,
WRITE_AIRS=T,
WRITE_TAMDAR=T,
/
&RECORD10
WIND_SD=F,
WIND_SD_SYNOP=F,
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,
/
Error reading NAMELIST record 10
Using default values for wind_sd
Error reading namelist file: namelist.obsproc
??? Should be STOPPed in PROC_NAMELIST
=== 3DVAR_OBSPROC is used for WRF 3DVAR ===
{number of slots = 1, length of the full-slot = 7200 sec.}
Setting up LAMBERT CONFORMAL map...
Computed cone factor: 0.715567
Computed pole i/j = 1.000 294.806
map_info: 3 1
code,lat1,lon1,dx,dlat,dlon,stdlon,truelat1,truelat2,hemi,cone,polei,polej,rsw,rebydx,knowni,knownj,init
3 28.0000000 120.000000 27000.0000 0.00000000 0.00000000 120.000000 30.0000000 60.0000000 1.00000000 0.715566695 1.00000000 294.805542 293.805542 235.925919 39.5000000 50.0000000 T
COARSE GRID CENTER IS AT X = 0.0 KM AND Y = -7932.8 KM. DIS(1)= 27.000
LLXY:lat, lon, (dot) xxi, yyj: 28.00000 120.00000 40.00000 50.50000
latlon_to_ij:lat, lon, (cross)xxi, yyj: 28.00000 120.00000 39.50000 50.00000
------------------------------------------------------------------------------
READ GTS OBSERVATIONS IN FILE /home/mode/data/2010.08.29.ADP.Sounding/OBS_2010083000
At line 116 of file module_date.f
Fortran runtime error: Bad value during integer read
Error termination. Backtrace:
#0 0x7fc288a29f1a
#1 0x7fc288a2aa75
#2 0x7fc288a2b22c
#3 0x7fc288af4f03
#4 0x7fc288af90b4
#5 0x55a121392a77
#6 0x55a12139e0d7
#7 0x55a12142b4d7
#8 0x55a121434587
#9 0x55a12138c982
#10 0x7fc2881473f0
#11 0x55a12138c9b9
#12 0xffffffffffffffff
Diagnostics in file obs_gts_read.diag
1527 modified the original lat = -90.00 to -89.5
--------------------------------------------------------------------------------------------------------------------
请问各位大神,是哪里设置出问题了么?
|
|