gen_be出错

一两

请教一下各位前辈，在测试官网案例，运行gen_bc希望生成cv5选项的BE时，出现如下报错从什么方面去排查呢？
gen_be_stage0_wrf.2008020400.log文件中：
---------------------------- FATAL ERROR -----------------------Fatal error in file:  gen_be_stage1.b  LINE:     145
Could not open /share/home/wangwk02/lyf/Build_WRF/WRFDA/var/scripts/gen_be/gen_be5_cv5/working/stage0/pert.2008020500.e001

屏幕上显示的是：WRFVAR_DIR is /share/home/wangwk02/lyf/Build_WRF/WRFDA Unversioned directory
这是为什么呢，这个路径是真实存在的呀，反复核对了也是正确的路径，pert.2008020500.e001的大小是0kb，应该是没有生成

一两

顶一下

Nozigaka46

运行脚本贴出来看看

一两

Nozigaka46 发表于 2024-5-14 21:31
运行脚本贴出来看看

#! /bin/ksh
#-----------------------------------------------------------------------
# Script gen_be_wrapper.ksh
#
# Purpose: Calculates background error statistics for WRF-Var.
#-----------------------------------------------------------------------

#[1] Define job by overriding default environment variables:

export WRFVAR_DIR=/share/home/wangwk02/lyf/Build_WRF/WRFDA

export RUN_GEN_BE_STAGE0=true
export RUN_GEN_BE_STAGE1=true
export RUN_GEN_BE_STAGE2=true
export RUN_GEN_BE_STAGE2A=true
export RUN_GEN_BE_STAGE3=true
export RUN_GEN_BE_STAGE4=true
export RUN_GEN_BE_DIAGS=true
export RUN_GEN_BE_DIAGS_READ=true
export RUN_GEN_BE_MULTICOV=true


export NL_CV_OPTIONS=5       # Specify control variable options
                             # NL_CV_OPTIONS = 5: default wind control variables (psi and chi_u)
                             #                 7: u/v wind control variables

export BIN_TYPE=5            # How data is binned for calculating statistics
                             # BIN_TYPE = 5 (default, recommended): Average over all horizontal points;
                             #                                        i.e. only one bin per vertical level
                             #            0: No binning
                             #            1: X-direction mean
                             #            2: Bin by latitude and height
                             #            3: Bin by latitude and vertical level
                             #            4: Fixed number of horizontal bins and by vertical level
                             #            6: Average over all points (only 1 bin)

export START_DATE=2008020612 # the first perturbation valid date
export END_DATE=2008020700   # the last perturbation valid date
export NUM_LEVELS=39         # = bottom_top = e_vert - 1
#export DATA_ON_LEVELS=.true. # "False if fields projected onto modes."
#export ALLOW_MISSING_DATES=.false.  # Set to true if you want GEN_BE to attempt to continue with gaps in data (missing dates)

export BE_METHOD=NMC
#Example of changes required for "be_method=ENS":
#export BE_METHOD=ENS
#export NE=2 # 30

export FC_DIR=/share/home/wangwk02/lyf/testdata/gen_be_forecasts_200802/gen_be_fc   # where wrf forecasts are
export RUN_DIR=/share/home/wangwk02/lyf/gen_be${BIN_TYPE}_cv${NL_CV_OPTIONS}
export DOMAIN=01             # For nested domains, set to the appropriate domain number
export FCST_RANGE1=24        # Longer forecast time for the NMC method (i.e. for 24-12 NMC, FCST_RANGE1=24, for 36-24 NMC, FCST_RANGE1=36)
export FCST_RANGE2=12        # Shorter forecast time for the NMC method (i.e. for 24-12 NMC, FCST_RANGE2=12, for 36-24 NMC, FCST_RANGE2=24)
export INTERVAL=12           # The interval between your forecast initial times
export STRIDE=1              # STRIDE=1 calculates correlation for every model grid point.
                             # STRIDE=2 calculates correlation every 2nd model gridpoint. 3 means every 3rd grid point, etc.
export NOCOLONS=false        # Optional, for use if you set nocolons=true in your WRF namelists
export USE_RFi=true                # use recursive filters?
#[2] Run gen_be:
if ${USE_RFi}; then
   ${WRFVAR_DIR}/var/scripts/gen_be/gen_be.ksh
else                          # loop over wavelet filter lengths:
   export DO_NORMALIZE=.false.        # normalize before wavelet transform?
   NEW_SUF=
   export RUN_DIR=${RUN_DIR}.
   for L in 7;do
      export WAVELET_NBAND=$L
      for N in C;do          # possible WAVELET_NAME values: B C D V
         export WAVELET_NAME=$N
         if [[ $WAVELET_NAME == B ]];then
            export ISTRT=18      
            export IINC=1
            export IFIN=$ISTRT
         elif [[ $WAVELET_NAME == C ]];then
            export ISTRT=30
            export IINC=6
            export IFIN=30
         elif [[ $WAVELET_NAME == D ]];then
            export ISTRT=6
            export IINC=2
            export IFIN=20
         elif [[ $WAVELET_NAME == V ]];then
            export ISTRT=24      
            export IINC=1
            export IFIN=$ISTRT
         fi
         for I in `seq ${ISTRT} ${IINC} ${IFIN}`; do
            export WAVELET_FILT_LEN=$I
            OLD_SUF=${NEW_SUF}
            NEW_SUF=${WAVELET_NBAND}${WAVELET_NAME}${WAVELET_FILT_LEN}n
            export RUN_DIR=${RUN_DIR%${OLD_SUF}}${NEW_SUF}
            ${WRFVAR_DIR}/var/scripts/gen_be/gen_be.ksh
         done
      done
   done
fi

一两

Nozigaka46 发表于 2024-5-14 21:31
运行脚本贴出来看看

前辈，在窗口的报错是这样的，在log中显示错误是：
Could not open /share/home/wangwk02/lyf/Build_WRF/WRFDA/var/scripts/gen_be/gen_be5_cv5/working/stage0/pert.2008020612.e001
我看了一下，pert.2008020612.e001确实没有生成，大小是0kb，不知道是为什么，能不能请前辈指点一二

Nozigaka46

一两 发表于 2024-5-15 09:51
前辈，在窗口的报错是这样的，在log中显示错误是：
Could not open /share/home/wangwk02/lyf/Build_WRF ...

你确定你的这些个20080203**/文件夹里面，对应所需要的wrfout文件都有吗？

一两

Nozigaka46 发表于 2024-5-15 12:28
你确定你的这些个20080203**/文件夹里面，对应所需要的wrfout文件都有吗？

2008020300文件夹下有wrfout_d01_2008-02-03_12_00_00  和  wrfout_d01_2008-02-04_00_00_00，就是官网下载的测试数据，我也有看到有些帖子说，要包含T,T+12,T+24，那我是把其他文件夹的文件直接复制过来吗

一两

Nozigaka46 发表于 2024-5-15 12:28
你确定你的这些个20080203**/文件夹里面，对应所需要的wrfout文件都有吗？

前辈，我刚才尝试了把开始预报时刻的wrfout文件也放进去了，还是不行，依然是相同的报错

Nozigaka46

一两 发表于 2024-5-15 12:38
前辈，我刚才尝试了把开始预报时刻的wrfout文件也放进去了，还是不行，依然是相同的报错

是一个文件夹有两个文件，分别是T+12和T+24。我建议你ls一下，因为我这里看不到你全部文件夹里面的信息

一两

Nozigaka46 发表于 2024-5-16 11:53
是一个文件夹有两个文件，分别是T+12和T+24。我建议你ls一下，因为我这里看不到你全部文件夹里面的信息

文件夹下是这样的