求助：用WRFDA做卫星辐射率偏差订正run出现错误

平凡的幸福

请问论坛里面的大神：我用WRFDA做卫星辐射率偏差订正的时候，如果用变分偏差订正，输入文件varbc.in，提交./da_wrfvar.exe，可以同化成功。但是如果用离线统计的偏差订正方法，生成noaa-18-amsua.bcor文件，放到workiing/biascorr目录下，运行./da_wrfvar.exe出现如下错误，请问是什么原因：
1.   My rsl.error.0000 files:
taskid: 0 hostname: fb04
module_io_quilt_old.F        2931 T
Namelist logging not found in namelist.input. Using registry defaults for varia
bles in logging.
Ntasks in X            2, ntasks in Y            4
*************************************
Parent domain
ids,ide,jds,jde            1         283           1         229
ims,ime,jms,jme           -4         148          -4          64
ips,ipe,jps,jpe            1         141           1          57
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1,              905633648 bytes allocated
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS      1 IE    141 JS      1 JE     57
WRF NUMBER OF TILES =   1
forrtl: severe (153): allocatable array or pointer is not allocated
Image              PC                Routine            Line        Source            
da_wrfvar.exe      00000000016703A6  Unknown               Unknown  Unknown
da_wrfvar.exe      00000000010679B1  Unknown               Unknown  Unknown
da_wrfvar.exe      00000000010508DD  Unknown               Unknown  Unknown
da_wrfvar.exe      0000000000A0E38C  Unknown               Unknown  Unknown
da_wrfvar.exe      000000000041B498  Unknown               Unknown  Unknown
da_wrfvar.exe      0000000000417C12  Unknown               Unknown  Unknown
da_wrfvar.exe      000000000040F86D  Unknown               Unknown  Unknown
da_wrfvar.exe      000000000040590E  Unknown               Unknown  Unknown
libc.so.6          000000347341ECDD  Unknown               Unknown  Unknown
da_wrfvar.exe      00000000004057A9  Unknown               Unknown  Unknown


2. rsl.out.0000 files:
......Reading radiance 1b data from amsua.bufr
Bufr file date is         2016           7          18          12
amsua               
num_tovs_file num_tovs_global num_tovs_local num_tovs_used num_tovs_thinned
    400491      1385       155        50       105
Allocating space for radiance innov structure  1  noaa-18-amsua                  50
Observation summary
   ob time  1
      geoamv               121 global,      27 local
      radiance            1385 global,      50 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=    50

Calculate innovation vector(iv)

Performing bias correction for noaa-18-amsua

3. noaa-18-amsua.bcor file:
   15    30    18     4
    1        88  202.42   14.42   -0.09    3.00
    2        86  174.53   11.45   -0.08    3.34
    3        96  237.33    8.67   -0.24    2.37
    4        95  264.78    2.11   -0.10    0.45
    5        99  258.37    2.42    0.02    0.26
    6        98  241.57    2.49    0.01    0.14
    7        96  228.80    1.97   -0.02    0.12
    8       101  218.83    1.89   -0.01    0.21
    9        96  210.77    2.38    0.05    0.14
   10        98  215.78    1.90   -0.06    0.17
   11        93  223.95    1.49    0.07    0.20
   12        95  233.35    1.43   -0.09    0.48
   13        92  244.85    1.56    0.04    0.60
   14        96  255.63    1.61    0.07    0.80
   15        88  251.35   11.76   -0.03    2.80
    1    -0.02946    -0.02609    -0.02415    -0.08686   506.99227
    2    -0.01513    -0.03698    -0.82611    -0.03595   690.68902
    3    -0.01823    -0.02923    -0.58268     0.01185   584.73767
    4    -0.00163    -0.00246    -0.03675    -0.00391    46.53369
    5    -0.00174    -0.00052     0.03221    -0.00939    11.46845
    6     0.00009    -0.00009     0.00046    -0.00436    -1.05511
    7     0.00056    -0.00012    -0.02901    -0.00022     2.71288
    8     0.00148    -0.00011    -0.04089    -0.00354    -1.90358
    9     0.00190     0.00016    -0.04777     0.00285    -6.03233
   10     0.00088     0.00028    -0.00569     0.00019    -9.21375
   11    -0.00059    -0.00040     0.00542     0.00233     6.82853
   12    -0.00001    -0.00073    -0.05427     0.01823    21.03778
   13     0.00174    -0.00108    -0.12829     0.02770    28.10876
   14    -0.00023    -0.00026    -0.04260     0.01237    14.34453
   15    -0.03899    -0.03490    -0.01514    -0.20166   669.89459
RELATIVE SCAN BIASES
           1      2      3      4      5      6      7      8      9     10     11     12     13     14     15     16     17     18     19     20     21     22     23     24     25     26     27     28     29     30
    1   0.00   0.00 -11.30   0.01   1.93  -1.71  -3.55   2.70   0.87   4.10   0.84  -0.09  -0.35   0.80  -0.66   0.66  -1.93  -4.54   5.74   1.67  -1.02  -2.19  -7.83  -7.20  -5.23  -9.00  -9.09   0.00   0.00   0.00
    2   0.00   0.00   0.00 -11.53   0.55  -4.61  -6.41  -1.53  -2.12   1.49  -1.33  -2.07  -0.38  -1.69  -0.90   0.90  -2.42  -3.65   5.32   2.24   0.32  -4.11  -8.76  -7.48  -7.86 -10.39 -11.16   0.00   0.00   0.00
    3   0.00   0.00  -8.04  -4.18  -1.20  -2.75  -3.11   1.71  -1.06   0.85  -0.33  -2.05  -0.66  -1.74  -1.04   1.04  -3.28  -2.95   1.27   2.30  -1.66  -3.62  -6.45  -4.91  -5.24  -6.36  -8.01   0.00   0.00   0.00
    4   0.00   0.00  -3.17  -2.57  -1.92  -1.43  -1.40  -1.26  -0.48  -0.33  -0.51  -0.42  -0.29  -0.47  -0.05   0.05  -0.61  -0.70   0.07  -0.04  -0.37  -1.30  -1.39  -1.64  -1.68  -2.26  -2.85   0.00   0.00   0.00
    5   0.00   0.00  -1.42  -1.04  -0.92  -0.59  -0.21  -0.59  -0.07  -0.24  -0.04   0.06  -0.01   0.05   0.05  -0.05   0.12  -0.11   0.21   0.01   0.09  -0.73  -0.33  -0.38  -0.74  -1.10  -1.40   0.00   0.00   0.00
    6   0.00   0.00  -0.87  -0.72  -0.69  -0.51  -0.38  -0.48  -0.18  -0.08  -0.09   0.01   0.09   0.03   0.05  -0.05   0.07   0.04  -0.08   0.05   0.08  -0.01  -0.07   0.03  -0.09  -0.33  -0.59   0.00   0.00   0.00
    7   0.00   0.00  -0.43  -0.34  -0.25  -0.23  -0.08  -0.31  -0.09  -0.10  -0.06  -0.11  -0.01   0.18  -0.13   0.13   0.13  -0.02   0.13   0.22  -0.06   0.31   0.10   0.26   0.18   0.11  -0.19   0.00   0.00   0.00
    8   0.00   0.00  -0.81  -0.64  -0.47  -0.49  -0.27  -0.45  -0.31  -0.13  -0.29  -0.08  -0.09   0.08   0.07  -0.07  -0.59  -0.19   0.02   0.56   0.10   0.24   0.19   0.32   0.11  -0.14  -0.60   0.00   0.00   0.00
    9   0.00   0.00  -0.30  -0.57  -0.38  -0.36  -0.34  -0.33  -0.45  -0.37  -0.32  -0.23  -0.20  -0.09  -0.05   0.05  -0.02   0.00   0.01  -0.02  -0.04  -0.07  -0.37  -0.36  -0.46  -0.53   0.00   0.00   0.00   0.00
   10   0.00   0.00   0.32   0.21  -0.01   0.10   0.07   0.07  -0.15   0.13  -0.01   0.05   0.06  -0.26   0.20  -0.20   0.32   0.42  -0.18   0.04   0.17  -0.22  -0.09  -0.14  -0.15  -0.31  -0.43   0.00   0.00   0.00
   11   0.00   0.00   0.00   0.68   0.18   0.55   0.30   0.36   0.16   0.42   0.13   0.29   0.20  -0.04   0.26  -0.26  -0.13  -0.05  -0.02  -0.20   0.02  -0.36   0.25  -0.38  -0.42  -0.25  -0.26   0.00   0.00   0.00
   12   0.00   0.00   1.95   1.12   0.58   0.77   0.44   0.86   1.00   0.45   0.54  -0.08   0.11  -0.14   0.16  -0.16  -0.16  -0.80  -0.66   0.81   0.29   0.82   0.09  -0.14   0.01  -0.18  -0.12   0.00   0.00   0.00
   13   0.00   0.00   1.64  -0.09  -0.43   0.14   0.08   0.20   0.64  -0.60   0.64  -0.75  -0.41  -1.57   0.34  -0.34  -1.34  -0.43  -1.24   0.43  -0.76   1.27  -0.07   0.16   0.38  -0.48  -0.36   0.00   0.00   0.00
   14   0.00   0.00   0.93  -0.35  -0.77   0.75   0.48   1.67   0.79  -0.11   0.72  -0.61  -0.43  -0.85  -0.59   0.59   0.53   0.89  -0.18   2.06   0.35   2.39  -0.40   0.05   1.38   0.89   1.32   0.00   0.00   0.00
   15   0.00   0.00  -4.37   1.02  -0.21   1.94   0.50   1.13   2.90   4.72   1.39   2.41  -0.45   1.15   1.11  -1.11  -1.72   0.37  -1.10   5.00  -1.13  -0.99  -5.84  -0.70   0.27  -2.16  -5.74   0.00   0.00   0.00
4.namelist.input file:
&wrfvar1
print_detail_grad=.true.,
/
&wrfvar2
/
&wrfvar3
ob_format =2 !1:bufr 2:ascii
/
&wrfvar4
use_amsuaobs  =.true.
use_amsubobs  =.false.
use_hirs3obs  =.false.
use_hirs4obs  =.false.
use_mhsobs    =.false.
use_airsobs   =.false.
use_eos_amsuaobs  =.false.
/
&wrfvar5
  check_max_iv=.false.,
  max_error_t=6.0,
  max_error_uv=5.0,
  max_error_pw=2.0,
  max_error_ref=5.0,
  max_error_q=8.0,
  max_error_p=5.0,
  max_error_thickness=5.0,
  max_error_rv=5.0,
  max_error_rf=5.0,
  max_error_rain=5.0,
/
&wrfvar6
!!!!ntmax=0
/
&wrfvar7
cv_options=3,
/
&wrfvar8
/
&wrfvar9
/
&wrfvar10
/
&wrfvar11
/
&wrfvar12
/
&wrfvar13
/
&wrfvar14
rtm_option=2,
RTMINIT_NSENSOR = 1                          !11 !!!!!!!!!!!!!!6 AMSUA; 3 AMSUB; 3 MHS; 1 AIRS; 1 SSMIS;3 HIRS3;3 HIRS4;
RTMINIT_PLATFORM =1,            !1, 1, 1, 10, 1, 1,  1, 10,   1,  1, 1,
RTMINIT_SATID=18,           !15,16,18, 2,  16,17, 18, 2,   16, 17,18,
RTMINIT_SENSOR=3,           !3, 3, 3, 3,  4, 4,  15,15,   0,  0,  0,

qc_rad=true,
only_sea_rad=true,
use_varbc                           = .false.  
freeze_varbc                        = .false.
read_biascoef=true  
biascorr=true       
write_iv_rad_ascii=true,
write_oa_rad_ascii=true,
! biasprep=true
thinning=true,
thinning_mesh=120.0
/
&wrfvar15
/
&wrfvar16
/
&wrfvar17
/
&wrfvar18
analysis_date="2016-07-18_12:00:00.0000",
/
&wrfvar19
/
&wrfvar20
/
&wrfvar21
time_window_min="2016-07-18_09:00:00.0000",
/
&wrfvar22
time_window_max="2016-07-18_15:00:00.0000",
/
&wrfvar23
/
&time_control
start_year=2016,
start_month=07,
start_day=18,
start_hour=12,
end_year=2016,
end_month=07,
end_day=18,
end_hour=12,
/
&fdda
/
&domains
e_we=283,
e_sn=229,
e_vert=51,
dx=9000.00,
dy=9000.00,
/
&dfi_control
/
&tc
/
&physics
mp_physics=8,
ra_lw_physics=4,
ra_sw_physics=4,
radt=20,
sf_sfclay_physics=2,
sf_surface_physics=2,
sf_urban_physics=0,
bl_pbl_physics=2,
cu_physics=3,
cudt=0,
num_soil_layers=4,
mp_zero_out=2,
co2tf=0,
num_land_cat=20,
/
&scm
/
&dynamics
/
&bdy_control
/
&grib2
/
&fire
/
&namelist_quilt
/

平凡的幸福

一楼自己坐，希望哪位做过相关工作的一起讨论，谢谢

平凡的幸福

问题已解决，WRFDA还有是bug的

我不喜欢吃生姜

平凡的幸福 发表于 2017-3-14 10:10
问题已解决，WRFDA还有是bug的

请问楼主是如何解决的

钟小继

请问楼主我同化的是卫星辐射率资料，资料是6小时间隔的，我模拟时24小时的，那我应该把5个观测数据文件全部放进去同化吗？那文件名字怎么编啊？还是说只放初始时刻的文件就可以了？

任性小王子

平凡的幸福 发表于 2017-3-14 10:10
问题已解决，WRFDA还有是bug的

可以分享一下经验吗？楼主

任性小王子

钟小继 发表于 2017-12-18 19:57
请问楼主我同化的是卫星辐射率资料，资料是6小时间隔的，我模拟时24小时的，那我应该把5个观测数据文件全部 ...

一般情况下是用初始时刻的文件