同化NCEP ADP高空资料出错

小其其格

想把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
--------------------------------------------------------------------------------------------------------------------

请问各位大神，是哪里设置出问题了么？