MHS资料同化没有数据

QMA

下载了2024年6月24日、25日、26日的MHS资料，使用WRFDA同化，但是不管是哪个时刻，00/06/12/18四个时刻都没有数据能够被应用。附上设置，求助各位老师！感谢！
namelist.input设置如下：
&wrfvar1
/
&wrfvar2
/
&wrfvar3
/
&wrfvar4
use_mhsobs = .true.
/
&wrfvar5
/
&wrfvar6
/
&wrfvar7
cv_options = 3,
/
&wrfvar8
/
&wrfvar9
/
&wrfvar10
/
&wrfvar11
check_rh = 1,
/
&wrfvar12
/
&wrfvar13
/
&wrfvar14
rtminit_nsensor = 1,
rtminit_platform = 1,
rtminit_satid = 18,
rtminit_sensor = 15,
thinning_mesh = 60.0,
thinning = .true.
qc_rad = .true.
write_iv_rad_ascii = .false.
write_oa_rad_ascii = .true.
rtm_option = 2
only_sea_rad = .false.
use_varbc = .true.
use_crtm_kmatrix = .true.
varbc_factor = 1.0,
varbc_nobsmin = 10,
use_blacklist_rad = .true.
crtm_coef_path="/home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs"
/
&wrfvar15
/
&wrfvar16
/
&wrfvar17
/
&wrfvar18
analysis_date="2024-06-24_00:00:00.0000",
/
&wrfvar19
/
&wrfvar20
/
&wrfvar21
time_window_min = "2024-06-23_21:00:00"
/
&wrfvar22
time_window_max = "2024-06-24_03:00:00"
/
&wrfvar23
/
&time_control
start_year                          = 2024,
start_month                         = 6,
start_day                           = 24,
start_hour                          = 00,
end_year                            = 2024,
end_month                           = 6,
end_day                             = 24,
end_hour                            = 06,
force_use_old_data=T,
/
&fdda
/
&domains
e_we   = 139,
e_sn   = 119,
e_vert = 51,
dx     = 15000,
dy     = 15000,
/
&dfi_control
/
&tc
/
&physics
mp_physics                          = 6,
ra_lw_physics                       = 1,
ra_sw_physics                       = 1,
radt                                = 15,
sf_sfclay_physics                   = 1,
sf_surface_physics                  = 2,
bl_pbl_physics                      = 1,
bldt                                = 0,
cu_physics                          = 1,
cudt                                = 5,
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,
sf_ocean_physics                    = 0,
/
&scm
/
&dynamics
hybrid_opt                          = 2,
w_damping                           = 0,
diff_opt                            = 1,
km_opt                              = 4,
diff_6th_opt                        = 0,
diff_6th_factor                     = 0.12,
base_temp                           = 290.,
damp_opt                            = 3,
zdamp                               = 5000.,
dampcoef                            = 0.2,
khdif                               = 0,
kvdif                               = 0,
non_hydrostatic                     = .true.,
moist_adv_opt                       = 1,
scalar_adv_opt                      = 1,
gwd_opt                             = 1,
/
&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
/
&grib2
/
&fire
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
&perturbation
/

QMA

rsl.out如下：taskid: 0 hostname: zl-virtual-machine
Namelist logging not found in namelist.input. Using registry defaults for variables in logging.

*** VARIATIONAL ANALYSIS ***
    WRFDA V4.0                                                                     

Ntasks in X            2 , ntasks in Y            2
*************************************
Parent domain
ids,ide,jds,jde            1         139           1         119
ims,ime,jms,jme           -4          76          -4          66
ips,ipe,jps,jpe            1          69           1          59
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,             258609156  bytes allocated
   Input data is acceptable to use: fg
hybrid_opt =  2
  
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS      1 IE     69 JS      1 JE     59
WRF NUMBER OF TILES =   1
Domain mapping info:
map_proj =     1
cen_lat   =    0.287530E+02
cen_lon   =    0.116906E+03
truelat1  =    0.274750E+02
truelat2  =    0.274750E+02
start_lat =    0.205811E+02
start_lon =    0.107026E+03
pole_lat  =    0.900000E+02
dsm       =    0.150000E+02
  
Set up observations (ob)
  
Using ASCII format observation input
  
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A >da_scan_obs_ascii.inc</a>  LINE:      68
Error    2 opening gts obs file ob.ascii
---------------------------------------------------------------
--------------------------- WARNING ---------------------------
WARNING FROM FILE: <A >da_read_obs_ascii.inc</a>  LINE:     106
Error    2 opening gts obs file ob.ascii
---------------------------------------------------------------
Using NCEP BUFR radiance 1b input
  
ACCoeff_ReadFile(Binary)(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/mhs_n18.SpcCoeff.bin;
ACCoeff RELEASE.VERSION:  1.04
   N_FOVS=90  N_CHANNELS=5
SpcCoeff_ReadFile(Binary)(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/mhs_n18.SpcCoeff.bin;
SpcCoeff RELEASE.VERSION:  8.01
   N_CHANNELS=5
Read_ODPS_Binary(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/mhs_n18.TauCoeff.bin;
   ODPS RELEASE.VERSION:  2.01  N_LAYERS=100  N_COMPONENTS=2  N_ABSORBERS=1  N_CHANNELS=5  N_COEFFS=10200
CloudCoeff_ReadFile(Binary)(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/CloudCoeff.bin;
CloudCoeff RELEASE.VERSION:  3.04  N_FREQUENCIES(MW)=  31  N_FREQUENCIES(IR)=  61  N_RADII(MW)=10  N_RADII(IR)=10  N_TEMPERATURES= 5  N_DENSITIES= 3  N_LEGENDRE_TERMS=38  N_PHASE_ELEMENTS= 1
AerosolCoeff_ReadFile(Binary)(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/AerosolCoeff.bin;
AerosolCoeff RELEASE.VERSION:  3.04  N_WAVELENGTHS=  61  N_RADII= 36  N_TYPES= 8  N_RH= 36  N_LEGENDRE_TERMS=37  N_PHASE_ELEMENTS= 1
FitCoeff_ReadFile(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 2
FitCoeff_ReadFile(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 6
FitCoeff_ReadFile(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 3, 7, 2
FitCoeff_ReadFile(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 10, 4, 3
FitCoeff_ReadFile(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 8
FitCoeff_ReadFile(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
FitCoeff RELEASE.VERSION : 1.5; DIMENSIONS= 3, 6, 2
MWwaterLUT_ReadFile(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
MWwaterLUT RELEASE.VERSION:  1.01
   N_ANGLES=15  N_FREQUENCIES=18  N_TEMPERATURES=12  N_WIND_SPEEDS=16
MWwaterCoeff_ReadFile(INFORMATION) : FILE: /home/zl/Build_WRF/WRFDA/CRTM/crtm_coeffs/FASTEM5.MWwater.EmisCoeff.bin;
MWwaterCoeff RELEASE.VERSION: 1.5
Reading radiance 1b data from mhs.bufr
Bufr file date is         2024           6          24           0 mhs                 
num_tovs_file num_tovs_global num_tovs_local num_tovs_used num_tovs_thinned
   2157939         0         0         0         0
Observation summary
   ob time  1

----------------------------------------------------------
[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
  
VARBC: Reading VARBC.in file
VARBC: no matching for platform/satid/sensor   1  15   3
VARBC: no matching for platform/satid/sensor   1  16   3
VARBC: no matching for platform/satid/sensor   1  18   3
VARBC: no matching for platform/satid/sensor   1  19   3
VARBC: no matching for platform/satid/sensor  23   1  40
VARBC: no matching for platform/satid/sensor  23   2  40
VARBC: no matching for platform/satid/sensor  10   1   3
VARBC: no matching for platform/satid/sensor  10   2   3
VARBC: no matching for platform/satid/sensor   1  15   4
VARBC: no matching for platform/satid/sensor   1  16   4
VARBC: no matching for platform/satid/sensor   1  17   4
VARBC: no matching for platform/satid/sensor   1  19  15
VARBC: no matching for platform/satid/sensor  23   1  41
VARBC: no matching for platform/satid/sensor  23   2  41
VARBC: no matching for platform/satid/sensor  10   1  15
VARBC: no matching for platform/satid/sensor  10   2  15
VARBC: no matching for platform/satid/sensor  17   0  19
VARBC: no matching for platform/satid/sensor   2  16  10
VARBC: no matching for platform/satid/sensor   9   2   3
VARBC: no matching for platform/satid/sensor   9   2  11
VARBC: no matching for platform/satid/sensor  19   0   0
VARBC: no matching for platform/satid/sensor   1  11   0
VARBC: no matching for platform/satid/sensor   1  12   0
VARBC: no matching for platform/satid/sensor   1  14   0
VARBC: no matching for platform/satid/sensor  19   0   1
VARBC: no matching for platform/satid/sensor   1  11   1
VARBC: no matching for platform/satid/sensor   1  12   1
VARBC: no matching for platform/satid/sensor   1  14   1
VARBC: no matching for platform/satid/sensor  10   1  16
VARBC: no matching for platform/satid/sensor  10   2  16
VARBC: no matching for platform/satid/sensor  12   2  21
VARBC: no matching for platform/satid/sensor  12   3  21
VARBC: no matching for platform/satid/sensor  29   1  63
VARBC: no matching for platform/satid/sensor   4  13  22
VARBC: no matching for platform/satid/sensor   4  14  22
VARBC: no matching for platform/satid/sensor   4  15  22
Calculate innovation vector(iv)
  
VARBC: Cold-starting noaa-18-mhs    3 -->  0.00
VARBC: Cold-starting noaa-18-mhs    4 -->  0.00
VARBC: Cold-starting noaa-18-mhs    5 -->  0.00
VARBC: Applying bias correction for noaa-18-mhs
VARBC: Estimate Hessian for preconditioning
VARBC:     0 active observations for noaa-18-mhs channel    3
VARBC:     0 active observations for noaa-18-mhs channel    4
VARBC:     0 active observations for noaa-18-mhs channel    5
Minimize cost function using CG method
  
  
Diagnostics
   Final cost function J       =         0.00
  
   Total number of obs.        =        0
   Final value of J            =         0.00000
   Final value of Jo           =         0.00000
   Final value of Jd           =         0.00000
   Final value of Jb           =         0.00000
   Final value of Jc           =         0.00000
   Final value of Je           =         0.00000
   Final value of Jp           =         0.00000
   Final value of Jl           =         0.00000
   Final value of Jm           =         0.00000
   Jb factor used              =         1.00000
   Je factor used              =         1.00000
   VarBC factor used           =         1.00000
  
   Total number of radiances    =        0
   Cost function for radiances  =         0.00000
  
Writing radiance OMA ascii file
  
VARBC: Updating bias parameters
VARBC: Not enough data to keep statistics for noaa-18-mhs    3
VARBC: Not enough data to keep statistics for noaa-18-mhs    4
VARBC: Not enough data to keep statistics for noaa-18-mhs    5
VARBC: Writing information in VARBC.out file
*** WRF-Var completed successfully ***

QMA

顶顶顶

任性小王子

QMA 发表于 2025-3-20 19:45
rsl.out如下：taskid: 0 hostname: zl-virtual-machine
Namelist logging not found in namelist.input.  ...

mhs.bufr (including MHS data from NOAA-18/19 and METOP-1 and -2)     参照这个数据格式gdas1.t00z.1bmhs.tm00.bufr_d

mhs.bufr
我没做过 MHS的同化。

QMA

任性小王子 发表于 2025-3-24 09:44
mhs.bufr (including MHS data from NOAA-18/19 and METOP-1 and -2)     参照这个数据格式gdas1.t00z.1b ...

好的，感谢！

任性小王子

QMA 发表于 2025-3-24 11:49
好的，感谢！

不客气哟，没帮到你。

QMA

顶一顶，有做过相关的大佬可以帮帮忙吗？