WRFDA数据同化运行gen_be_wrapper.ksh出错

齐天笑猴 · 发表于 2018-1-1 20:56:33

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

运行gen_be_wrapper.ksh时，在$RUN_DIR directory生成了be.dat文件。但是出现
cd: /home/sxw/program/WRFDA/var/scripts/gen_be/gen_be5_cv5/working: [No such file or directory]
运行代码及过程如下：

[sxw@localhost gen_be]$ ./gen_be_wrapper.ksh
Mon Jan 1 18:04:39 EST 2018 Start
WRFVAR_DIR is /home/sxw/program/WRFDA exported
RUN_DIR is /home/sxw/program/WRFDA/var/scripts/gen_be/gen_be5_cv5
---------------------------------------------------------------
Run Stage 0: Calculate ensemble perturbations from model forecasts.
---------------------------------------------------------------
Beginning CPU time: Mon Jan 1 18:04:39 EST 2018
gen_be_stage0_wrf: Calculating standard perturbation fields valid at time 2015010200
2015010100 /home/sxw/program/WRFDA/var/scripts/gen_be/forcasts/2015010100/wrfout_d02_2015-01-02_00:00:00 /home/sxw/program/WRFDA/var/scrip ts/gen_be/forcasts/2015010112/wrfout_d02_2015-01-02_00:00:00
gen_be_stage0_wrf: Calculating standard perturbation fields valid at time 2015010212
2015010112 /home/sxw/program/WRFDA/var/scripts/gen_be/forcasts/2015010112/wrfout_d02_2015-01-02_12:00:00 /home/sxw/program/WRFDA/var/scrip ts/gen_be/forcasts/2015010200/wrfout_d02_2015-01-02_12:00:00
Ending CPU time: Mon Jan 1 18:04:52 EST 2018
---------------------------------------------------------------
Run Stage 1: Read standard fields, and remove time/ensemble/area mean.
---------------------------------------------------------------
Beginning CPU time: Mon Jan 1 18:04:52 EST 2018
Ending CPU time: Mon Jan 1 18:04:53 EST 2018
---------------------------------------------------------------
Run Stage 2: Calculate regression coefficients.
---------------------------------------------------------------
Beginning CPU time: Mon Jan 1 18:04:53 EST 2018
Ending CPU time: Mon Jan 1 18:04:54 EST 2018
---------------------------------------------------------------
Run Stage 2a: Calculate control variable fields.
---------------------------------------------------------------
Beginning CPU time: Mon Jan 1 18:04:54 EST 2018
Ending CPU time: Mon Jan 1 18:04:54 EST 2018
---------------------------------------------------------------
Run Stage 3: Read 3D control variable fields, and calculate vertical covariances.
---------------------------------------------------------------
Beginning CPU time: Mon Jan 1 18:04:54 EST 2018
Ending CPU time: Mon Jan 1 18:04:57 EST 2018
---------------------------------------------------------------
Run Stage 4: Calculate horizontal covariances (regional lengthscales).
---------------------------------------------------------------
Beginning CPU time: Mon Jan 1 18:04:57 EST 2018
Ending CPU time: Mon Jan 1 18:17:44 EST 2018
Ending CPU time: Mon Jan 1 18:17:44 EST 2018
/home/sxw/program/WRFDA/var/scripts/gen_be/gen_be_cov3d.ksh[10]: cd: /home/sxw/program/WRFDA/var/scripts/gen_be/gen_be5_cv5/working: [No s uch file or directory]
/home/sxw/program/WRFDA/var/scripts/gen_be/gen_be_cov3d.ksh[10]: cd: /home/sxw/program/WRFDA/var/scripts/gen_be/gen_be5_cv5/working: [No s uch file or directory]
/home/sxw/program/WRFDA/var/scripts/gen_be/gen_be_cov2d.ksh[10]: cd: /home/sxw/program/WRFDA/var/scripts/gen_be/gen_be5_cv5/working: [No s uch file or directory]
---------------------------------------------------------------
Run gen_be_cov2d.
---------------------------------------------------------------
Mon Jan 1 18:17:48 EST 2018 Finished

复制代码

我的gen_be_wrapper.ksh脚本如下设置

export WRFVAR_DIR=/home/sxw/program/WRFDA
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=2015010200 # the first perturbation valid date
export END_DATE=2015010212 # the last perturbation valid date
export NUM_LEVELS=56 # = bottom_top = e_vert - 1
#export DATA_ON_LEVELS=.true. # "False if fields projected onto modes."
export BE_METHOD=NMC
export FCST_RANGE=12
#Example of changes required for "be_method=ENS":
#export BE_METHOD=ENS
#export NE=2 # 30
export FC_DIR=/home/sxw/program/WRFDA/var/scripts/gen_be/forcasts # where wrf forecasts are
export RUN_DIR=`pwd`/gen_be${BIN_TYPE}_cv${NL_CV_OPTIONS}
export DOMAIN=02
export FCST_RANGE1=24
export FCST_RANGE2=12
export INTERVAL=12
export STRIDE=1
export USE_RFi=true
#[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

复制代码

不知道哪儿有什么问题？还请各位大神指教！非常感谢啦！

杨智灵 · 发表于 2021-9-22 10:03:58

请问大神，怎么解决的这个问题啊

Soaring · 发表于 2022-4-13 11:52:32

应该是你的.bashrc文件环境变量里面指定了workdir，导致这里生成的目录和文件自动放到那边的workdir里了。我这两天也碰到了这个问题，把.bashrc里面指定工作目录workdir注释掉就可以了

任性小王子 · 发表于 2024-4-7 11:09:35

非常感谢楼主分享经验，避免踩坑

		自动登录	找回密码
密码			立即注册