WRFDA-4Dvar数据同化过程出错，寻求帮助！

足无际

大家好！有做过WRFDA-4Dvar的同学不？最近模拟个降水过程，为提高模拟精度，试图采用WRFDA-4Dvar方法。但是不知为什么同化这一步始终过不去不过去，错误提示如下。在运行结果后面附有namelist.obsproc和WRFDA中的namelist.input文件，希望做过这方面工作的同学给帮忙看看，谢谢啦！
*** VARIATIONAL ANALYSIS ***
    WRFDA V3.9.1

DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,             545050944  bytes allocated
hybrid_opt = -1

Timing for processing lateral boundary for domain        1:    0.02738 elapsed seconds
num_fgat_time is:          4

Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS      1 IE     99 JS      1 JE    117
WRF NUMBER OF TILES =   1
Domain mapping info:
map_proj =     1
cen_lat   =    0.305500E+02
cen_lon   =    0.119570E+03
truelat1  =    0.305000E+02
truelat2  =    0.305000E+02
start_lat =    0.143034E+02
start_lon =    0.106459E+03
pole_lat  =    0.900000E+02
dsm       =    0.300000E+02

Set up observations (ob)

Using ASCII format observation input

Observation summary
   ob time  1
   ob time  2
   ob time  3
   ob time  4
      airep                 38 global,      38 local

----------------------------------------------------------
[3.0] Set up background errors (be) for cv_option:  3
----------------------------------------------------------

    DA_Setup_Background_Errors: 3DVAR dry control variables are:
    psi, chi_u, t_u and psfc
    DA_Setup_Background_Errors: 3DVAR humidity control variable is q/qsg

Number of vertical level for stats =         42
Number of latitude           nlath =         96

Number of vertical level for WRFVar=    29


Serial code ==> Write CV3 BE to fort.503
ids, ide, jds, jde, kds, kde, ku, samp:    1   99    1  117    1   29   22        0.58191
--- WARNING: traj_opt is zero, but num_traj is not zero; setting num_traj to zero.
--- NOTE: sst_update is 0, setting io_form_auxinput4 = 0 and auxinput4_interval = 0 for all domains
--- NOTE: grid_fdda is 0 for domain      1, setting gfdda interval and ending time to 0 for that domain.
--- NOTE: both grid_sfdda and pxlsm_soil_nudge are 0 for domain      1, setting sgfdda interval and ending time to 0 for that domain.
--- NOTE: obs_nudge_opt is 0 for domain      1, setting obs nudging interval and ending time to 0 for that domain.
--- NOTE: bl_pbl_physics /= 4, implies mfshconv must be 0, resetting
Need MYNN PBL for icloud_bl = 1, resetting to 0
*************************************
No physics suite selected.
Physics options will be used directly from the namelist.
*************************************
--- NOTE: RRTMG radiation is in use, setting:  levsiz=59, alevsiz=12, no_src_types=6
--- NOTE: num_soil_layers has been set to      4
WRF V3.9.1 MODEL
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,            1290211200  bytes allocated
  med_initialdata_input: calling input_input
Timing for processing wrfinput file (stream 0) for domain        1:    0.17828 elapsed seconds
Max map factor in domain 1 =  1.04. Scale the dt in the model accordingly.
INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          33  CATEGORIES           2  SEASONS WATER CATEGORY =           17  SNOW CATEGORY =           15
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS
LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND          20  CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND          19  CATEGORIES
open ETAMPNEW_DATA.expanded_rain, in fer_hires
  xtraj linked list is initialized
  ad_forcing linked list is initialized
  tl_pert linked list is initialized
In dolph, deltat =   360.00 taus =    3600.00
In dolph: M,N =       30      61
In dolph: THETAS (STOP-BAND EDGE) =      0.628
In dolph: R,DB =      0.000  -160.466
In dolph: TIME, W =   0.000000      0.5612023E-01
In dolph: TIME, W =   1.000000      0.5558921E-01
In dolph: TIME, W =   2.000000      0.5402456E-01
In dolph: TIME, W =   3.000000      0.5150921E-01
In dolph: TIME, W =   4.000000      0.4817360E-01
In dolph: TIME, W =   5.000000      0.4418507E-01
In dolph: TIME, W =   6.000000      0.3973453E-01
In dolph: TIME, W =   7.000000      0.3502203E-01
In dolph: TIME, W =   8.000000      0.3024247E-01
In dolph: TIME, W =   9.000000      0.2557311E-01
In dolph: TIME, W =  10.000000      0.2116353E-01
In dolph: TIME, W =  11.000000      0.1712911E-01
In dolph: TIME, W =  12.000000      0.1354809E-01
In dolph: TIME, W =  13.000000      0.1046196E-01
In dolph: TIME, W =  14.000000      0.7878942E-02
In dolph: TIME, W =  15.000000      0.5779498E-02
In dolph: TIME, W =  16.000000      0.4123162E-02
In dolph: TIME, W =  17.000000      0.2855761E-02
In dolph: TIME, W =  18.000000      0.1916250E-02
In dolph: TIME, W =  19.000000      0.1242581E-02
In dolph: TIME, W =  20.000000      0.7762616E-03
In dolph: TIME, W =  21.000000      0.4654370E-03
In dolph: TIME, W =  22.000000      0.2665784E-03
In dolph: TIME, W =  23.000000      0.1449647E-03
In dolph: TIME, W =  24.000000      0.7425032E-04
In dolph: TIME, W =  25.000000      0.3543180E-04
In dolph: TIME, W =  26.000000      0.1550860E-04
In dolph: TIME, W =  27.000000      0.6080120E-05
In dolph: TIME, W =  28.000000      0.2051719E-05
In dolph: TIME, W =  29.000000      0.5513138E-06
In dolph: TIME, W =  30.000000      0.9622389E-07
In dolph: SUM OF WEIGHTS W2 =     1.0000
Max map factor in domain 1 =  1.04. Scale the dt in the model accordingly.
INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          33  CATEGORIES           2  SEASONS WATER CATEGORY =           17  SNOW CATEGORY =           15
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS
LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND          20  CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND          19  CATEGORIES
open ETAMPNEW_DATA.expanded_rain, in fer_hires
Timing for processing lateral boundary for domain        1:    0.03450 elapsed seconds
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS      1 IE    100 JS      1 JE    118
WRF NUMBER OF TILES =   1
WARN5: NaN temperature; I,J,L,P=          54          38          15                       NaN
Timing for main: time 2013-10-06_06:06:00 on domain   1:   10.78650 elapsed seconds
d01 2013-10-06_06:06:00           47  points exceeded cfl=2 in domain d01 at time 2013-10-06_06:06:00 hours
d01 2013-10-06_06:06:00  MAX AT i,j,k:           54          38          24  vert_cfl,w,d(eta)=   5.1724315072496951       -342.28425589034572        2.0671606411806435E-002
d01 2013-10-06_06:06:00         1158  points exceeded cfl=2 in domain d01 at time 2013-10-06_06:06:00 hours
d01 2013-10-06_06:06:00  MAX AT i,j,k:           48          37          11  vert_cfl,w,d(eta)=   215.76577410438631       -1517837.6864491212        4.8249425144130394E-002
d01 2013-10-06_06:06:00        12957  points exceeded cfl=2 in domain d01 at time 2013-10-06_06:06:00 hours
d01 2013-10-06_06:06:00  MAX AT i,j,k:           54          48          29  vert_cfl,w,d(eta)=   35216.454761867979        78717.843906799011        9.2214660239857173E-003
WARN1: Water saturation T<180K; I,J,L,TC,P=          96          30           4  -94.725712598460689        83.231162758703491
WARN2: Too many iterations in function DEPOSIT;  I,J,L,TC,SSAT,DWV=          11          38          16  -261.94635734700609        2.7822079571199292        7.4258191042679018E-003
Timing for main: time 2013-10-06_06:12:00 on domain   1:    1.66117 elapsed seconds

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0  0x7f2bb6c3b2da in ???
#1  0x7f2bb6c3a503 in ???
#2  0x7f2bb605ef1f in ???
#3  0x7f2bb9d0db5e in ???
#4  0x7f2bb9d0fa4f in ???
#5  0x7f2bb9d22cf0 in ???
#6  0x7f2bb9d25288 in ???
#7  0x7f2bb94252fc in ???
#8  0x7f2bb9571e81 in ???
#9  0x7f2bb91a305c in ???
#10  0x7f2bb8f996fc in ???
#11  0x7f2bb81f12c7 in ???
#12  0x7f2bb74335ef in ???
#13  0x7f2bb74e7768 in ???
#14  0x7f2bb7468489 in ???
#15  0x7f2bb747685b in ???
#16  0x7f2bb74768b2 in ???
#17  0x7f2bb743244f in ???
#18  0x7f2bb6041b96 in ???
#19  0x7f2bb7432539 in ???
#20  0xffffffffffffffff in ???
Segmentation fault (core dumped)



相关namelist文件如下：
（1）obsproc的namelist文件设置如下：
&record1
obs_gts_filename = './OBS:2013100606',
obs_err_filename = 'obserr.txt',
gts_from_mmm_archive = .true.,
/

&record2
time_window_min  = '2013-10-06_03:00:00',
time_analysis    = '2013-10-06_06:00:00',
time_window_max  = '2013-10-06_09: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         = .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 = 1,
PHIC  = 30.55,
XLONC = 119.57,
TRUELAT1= 30.5,
TRUELAT2= 30.5,
MOAD_CEN_LAT = 30.55,
STANDARD_LON = 119.57
/

&record8
IDD    =   1,
MAXNES =   3,
NESTIX =  118,  106,  94,  181,  211,
NESTJX =  100,  88,  76,  196,  211,
DIS    =  9,  3,  1,  1.1,  1.1,
NUMC   =    1,    1,   2,     3,    4,
NESTI  =    1,   36,  32,    35,   45,
NESTJ  =    1,   42,  38,    65,   55,
/

&record9
PREPBUFR_OUTPUT_FILENAME = 'prepbufr_output_filename',
PREPBUFR_TABLE_FILENAME = 'prepbufr_table_filename',
OUTPUT_OB_FORMAT = 2
use_for          = '4DVAR',
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.,
/
（2）WRFDA的namelist文件设置如下：

&wrfvar1
var4d=true,
var4d_lbc=false,
var4d_bin=3600,
print_detail_outerloop=false,
print_detail_grad=false,
/
&wrfvar2
/
&wrfvar3
ob_format=2,
/
&wrfvar4
use_synopobs=true,
use_shipsobs=true,
use_metarobs=true,
use_soundobs=true,
use_pilotobs=true,
use_airepobs=true,
use_geoamvobs=true,
use_polaramvobs=true,
use_bogusobs=true,
use_buoyobs=true,
use_profilerobs=true,
use_satemobs=true,
use_gpspwobs=true,
use_gpsrefobs=true,
use_qscatobs=true,
use_rainobs=false,
/
&wrfvar5
check_max_iv=true,
/
&wrfvar6
max_ext_its=5,
ntmax=10,
orthonorm_gradient=true,
/
&wrfvar7
cv_options=3,
/
&wrfvar8
/
&wrfvar9
/
&wrfvar10
test_transforms=false,
test_gradient=false,
/
&wrfvar11
cv_options_hum=1,
check_rh=0,
sfc_assi_options=1,
calculate_cg_cost_fn=false,
/
&wrfvar12
/
&wrfvar13
/
&wrfvar14
/
&wrfvar15
/
&wrfvar16
/
&wrfvar17
/
&wrfvar18
analysis_date="2013-10-06_06:00:00.0000",
/
&wrfvar19
/
&wrfvar20
/
&wrfvar21
time_window_min="2013-10-06_06:00:00.0000",
/
&wrfvar22
time_window_max="2013-10-06_09:00:00.0000",
/
&time_control
run_hours=6,
start_year=2013,
start_month=10,
start_day=06,
start_hour=06,
end_year=2013,
end_month=10,
end_day=06,
end_hour=09,
interval_seconds=21600,
debug_level=0,
/
&domains
time_step=360,
e_we=100,
e_sn=118,
e_vert=30,
p_top_requested=5000,
interp_type=1,
t_extrap_type=1,
dx=30000,
dy=30000,
smooth_option=0,
nproc_x=0,
hypsometric_opt=2,
/
&fdda
/
&dfi_control
/
&tc
/
&physics
mp_physics=5,
ra_lw_physics=1,
ra_sw_physics=4,
radt=9,
sf_sfclay_physics=2,
sf_surface_physics=2,
bl_pbl_physics=2,
cu_physics=1,
cudt=5,
num_soil_layers=4,
mp_zero_out=2,
co2tf=0,
/
&scm
/
&dynamics
w_damping=0,
diff_opt=1,
km_opt=4,
dampcoef=0.2,
time_step_sound=4,
base_temp=290.0,
use_baseparam_fr_nml=true,
/
&bdy_control
specified=true,
real_data_init_type=3,
/
&grib2
/
&namelist_quilt
/
&perturbation
trajectory_io=true,
enable_identity=false,
jcdfi_use=false,
jcdfi_diag=1,
jcdfi_penalty=1000.0,
/

ljchen1989

30km用360s积分步长，1:12啊~
是不是积分溢出了

足无际

非常感谢您的提示，dx=30000, time_step=6*da(km)，既time_step=180。我设置之后程序会执行，但是不知道为什么运行到最后一个时间点会突然跳出！！！！结果显示如下，如果可以加个QQ好友可以吗？我QQ：820465176。麻烦您啦！

Timing for main: time 2013-10-06_11:53:20 on domain   1:    1.30991 elapsed seconds
Timing for main: time 2013-10-06_11:55:00 on domain   1:    1.09922 elapsed seconds
Timing for main: time 2013-10-06_11:56:40 on domain   1:    1.14620 elapsed seconds
Timing for main: time 2013-10-06_11:58:20 on domain   1:    0.99070 elapsed seconds
Timing for main: time 2013-10-06_12:00:00 on domain   1:    1.17788 elapsed seconds
Calculate innovation vector(iv)


Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0  0x7f1f5803b2da in ???
#1  0x7f1f5803a503 in ???
#2  0x7f1f5745ef1f in ???
#3  0x7f1f588c9f0c in ???
#4  0x7f1f58dbda44 in ???
#5  0x7f1f58dbfc9c in ???
#6  0x7f1f58dce240 in ???
#7  0x7f1f5886857d in ???
#8  0x7f1f5887685b in ???
#9  0x7f1f588768b2 in ???
#10  0x7f1f5883244f in ???
#11  0x7f1f57441b96 in ???
#12  0x7f1f58832539 in ???
#13  0xffffffffffffffff in ???
Segmentation fault (core dumped)

足无际

ljchen1989 发表于 2019-2-28 21:43
30km用360s积分步长，1:12啊~
是不是积分溢出了

非常感谢您的提示，dx=30000, time_step=6*da(km)，既time_step=180。我设置之后程序会执行，但是不知道为什么运行到最后一个时间点会突然跳出！！！！结果显示如下，如果可以加个QQ好友可以吗？我QQ：820465176。麻烦您啦！

Timing for main: time 2013-10-06_11:53:20 on domain   1:    1.30991 elapsed seconds
Timing for main: time 2013-10-06_11:55:00 on domain   1:    1.09922 elapsed seconds
Timing for main: time 2013-10-06_11:56:40 on domain   1:    1.14620 elapsed seconds
Timing for main: time 2013-10-06_11:58:20 on domain   1:    0.99070 elapsed seconds
Timing for main: time 2013-10-06_12:00:00 on domain   1:    1.17788 elapsed seconds
Calculate innovation vector(iv)


Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0  0x7f1f5803b2da in ???
#1  0x7f1f5803a503 in ???
#2  0x7f1f5745ef1f in ???
#3  0x7f1f588c9f0c in ???
#4  0x7f1f58dbda44 in ???
#5  0x7f1f58dbfc9c in ???
#6  0x7f1f58dce240 in ???
#7  0x7f1f5886857d in ???
#8  0x7f1f5887685b in ???
#9  0x7f1f588768b2 in ???
#10  0x7f1f5883244f in ???
#11  0x7f1f57441b96 in ???
#12  0x7f1f58832539 in ???
#13  0xffffffffffffffff in ???
Segmentation fault (core dumped)

SAchinb

有没有试过将积分步长减小到120s试试呢？有什么样的效果

足无际

SAchinb 发表于 2019-3-16 20:07
有没有试过将积分步长减小到120s试试呢？有什么样的效果

谢谢您，不知为何，将e_vert设为31就可以了！这个如何解释呢？？？

tyfcwgl

slots 设置为3 为什么 ob time  只有四个？ 而且只有airep data 吗？

SAchinb

足无际 发表于 2019-3-17 23:03
谢谢您，不知为何，将e_vert设为31就可以了！这个如何解释呢？？？

e_vert是垂直方向的层数，因为wrfinput_d01里面的层数就是你指的31了，这是由wps里面namelist.wps的设置来决定的，也就是说同化的namelist.input应该和WPS的namelist.wps中的e_vert保持一致。我也在做同化，欢迎加我QQ 602798957

任性小王子

是时间设置问题吧，obsproc的namelist文件obs_gts_filename = './OBS:2013100606', &record2 time_window_min  = '2013-10-06_03:00:00', time_analysis    = '2013-10-06_06:00:00', WRFDA的namelist文件 &wrfvar18 analysis_date="2013-10-06_06:00:00.0000", &wrfvar21
time_window_min="2013-10-06_06:00:00.0000"。所以你的 WRFDA的namelist文件 &wrfvar21 time_window_min  = '2013-10-06_03:00:00',不应该是time_window_min="2013-10-06_06:00:00.0000"