Running gen_be 制作背景场错误信息
the B0 shown in the WRFDA flowchart
WRFDA/var/scripts/gen_be/gen_be_wrapper.ksh 是用到的主要问价
- cp /WRFDA/var/scripts/gen_be/gen_be_wrapper.ksh . 将gen_be_wrapper.ksh复制到当前工作文件夹下
- vi gen_be_wrapper.ksh
export WRFVAR_DIR=/WRFDA
export NL_CV_OPTIONS=5 # 指定控制变量选项 有CV3,CV5,CV6,CV7
export BIN_TYPE=5 # 如何对数据进行分组以计算统计数据
export START_DATE=2008020300 # 初始扰动有效时间
export END_DATE=2008020912 # 最后扰动有效时间
export NUM_LEVELS=40 # e_vert -1
export BE_METHOD=NMC
export FCST_RANGE=12 # No meaning for WRFDA
export FC_DIR=/fc #WRF 预报文件所在位置
export RUN_DIR=`pwd`/gen_be #GEN_BE的运行和输出文件的位置
export DOMAIN=01 # 对于嵌套域,设置为适当的域号
export FCST_RANGE1=24 # NMC方法的第二次预测时间 (i.e. for 24-12 NMC, FCST_RANGE1=24, for 36-12 NMC, FCST_RANGE1=36)
export FCST_RANGE2=12 # No meaning for WRFDA
export INTERVAL=12 # 预测初始时间之间的间隔
export STRIDE=1 # STRIDE=1 计算每个模型网格点的相关性。STRIDE=2 每 2 个模型网格点计算相关性。3 表示每 3 个网格点,依此类推。
export NOCOLONS=false # 可选,用于在 WRF 名称列表中设置 nocolons=true 时使用
CV_OPTIONS
CV3
- CV3是NCEP背景误差协方差。
- 每个变量的方差及其二阶导数的方差用于估计其水平尺度。例如,可以根据涡度和流函数的方差估计流函数的水平尺度。
- 垂直尺度是用每个变量的垂直相关性来估计的。将构建一个表来涵盖变量的垂直刻度范围。然后使用该表查找垂直网格单位的比例。过滤器轮廓和垂直相关性在本地拟合。表中最佳拟合的尺度被指定为每个纬度在该垂直水平上的变量比例。请注意,垂直尺度是局部定义的,因此不包括垂直方向上更远的负相关。
- 从理论上讲,**CV3 BE是一个通用的背景错误统计文件,可用于任何情况。**在您自己的情况下使用 CV3 非常简单。要在您的情况下使用 CV3 BE 文件,请在工作目录中的 namelist.input 中的 &wrfvar7 中设置 cv_options=3,并使用位于 WRFDA/var/run/be.dat.cv3 中的 be.dat。
- 由于 CV3 是通用背景误差选项,因此有必要调整每个实验的默认背景误差尺度长度和方差值。这些可以在运行时使用下面描述的一系列名称列表变量进行控制。
&wrfvar7
cv_options = 3,
as1 = 0.25, 1.00, 1.50, #流函数
as2 = 0.25, 1.00, 1.50, #不平衡的速度势
as3 = 0.25, 1.00, 1.50, #不平衡的温度
as4 = 0.25, 1.00, 1.50, #假相对湿度
as5 = 0.25, 1.00, 1.50, #不平衡的地面气压
/
#第一列:方差比例因子 第二列:水平长度比例因子 第三列:垂直长度比例因子
#对于多个外部循环,必须填充每个向量的接下来 3 个元素。因此,对于 2 个外部循环 (max_ext_its=2),要使用默认值,名称列表应如下所示:
&wrfvar7
cv_options = 3,
as1 = 0.25, 1.00, 1.50, 0.25, 1.00, 1.50,
as2 = 0.25, 1.00, 1.50, 0.25, 1.00, 1.50,
as3 = 0.25, 1.00, 1.50, 0.25, 1.00, 1.50,
as4 = 0.25, 1.00, 1.50, 0.25, 1.00, 1.50,
as5 = 0.25, 1.00, 1.50, 0.25, 1.00, 1.50,
/
CV6
CV6 是 WRFDA 中的多变量后台错误统计选项。通过设置 namelist 变量 cv_options=6 来激活它。此选项在分析控制变量的平衡部分的定义中引入了六个额外的相关系数。因此,通过这种实现,水汽分析是多变量的,因为温度和风可能导致水汽增加,反之亦然。更多详细信息可在以下位置查看:http://www2.mmm.ucar.edu/wrf/users/wrfda/Docs/WRFDA_updated_for_cv6.pdf
&wrfvar7
cv_options = 6,
psi_chi_factor = 1, #控制流函数在定义速度势平衡部分时的贡献的参数,默认是1,>1贡献增加,<1贡献减小
psi_t_factor = 1, #控制流函数在定义温度平衡部分时的贡献的参数,默认是1,>1贡献增加,<1贡献减小
psi_ps_factor = 1, #控制流函数在定义表面压力平衡部分的贡献的参数,默认是1,>1贡献增加,<1贡献减小
psi_rh_factor = 1, #控制流函数在定义水汽平衡部分的贡献的参数,默认是1,>1贡献增加,<1贡献减小
chi_u_t_factor = 1, #控制速度势不平衡部分在定义温度平衡部分时的贡献的参数,默认是1,>1贡献增加,<1贡献减小
chi_u_ps_factor = 1, #控制速度势不平衡部分在定义表面压力平衡部分时的贡献的参数,默认是1,>1贡献增加,<1贡献减小
chi_u_rh_factor = 1, #控制速度势不平衡部分在定义水汽平衡部分时的贡献的参数,默认是1,>1贡献增加,<1贡献减小
t_u_rh_factor = 1, #控制温度不平衡部分对水汽平衡部分的贡献的参数,默认是1,>1贡献增加,<1贡献减小
ps_u_rh_factor = 1, #控制表面压力不平衡部分在定义水汽平衡部分时的贡献的参数,默认是1,>1贡献增加,<1贡献减小
WRFDA 的 CV6 后台错误统计信息是通过执行脚本 gen_be/wrapper_gen_mbe.ksh 生成的,该脚本在SCRIPTS下。该过程的其余部分与正常运行gen_be实用程序时相同。成功的运行将在工作文件夹目录中创建一个 be.dat 文件。