很多时候，在利用WRF模式三维变分系统（3DVar）同化观测资料需要用到适用研究区域的背景误差协方差矩阵（简称B矩阵）。WRF 3DVar可以从CV3、CV5和CV6任选一种，CV3是由NCEP依据“NMC”方法（Parrish and Derber 1992）统计而得，适用区域或全球范围。根据实际经验，区域模式中CV5的同化效果优于CV3，现将CV5的统计过程记录如下：
一、生成同一时刻不同预报时效的预报文件
一般是统计一个月的预报文件，预报时间间隔为12h。可以自行写一个shell脚本实现。
二、运行$WRFDA/var/scripts/gen_be/gen_be_wrapper.ksh脚本后生成B矩阵
WRF 3DVar将产生B矩阵的可执行程序封装在ksh脚本gen_be_wrapper.ksh中，“天河二号”中指定ksh解释器的路径为 /BIGDATA1/app/ast/20180418-gcc-4.8.5/bin/ksh。在运行这个脚本之前，你需要在脚本中设置好相关的环境变量。
注意1：在stage0阶段（从模式预报文件中计算集合扰动），只需指定一个core即可。
注意2：stage4阶段耗时最长，可以将多个作业同时提交到一个node上，提高计算效率。
