很多時候，在利用WRF模式三維變分系統(tǒng)（3DVar）同化觀測資料需要用到適用研究區(qū)域的背景誤差協(xié)方差矩陣（簡稱B矩陣）。WRF 3DVar可以從CV3、CV5和CV6任選一種，CV3是由NCEP依據(jù)“NMC”方法（Parrish and Derber 1992）統(tǒng)計而得，適用區(qū)域或全球范圍。根據(jù)實際經(jīng)驗，區(qū)域模式中CV5的同化效果優(yōu)于CV3，現(xiàn)將CV5的統(tǒng)計過程記錄如下：

一、生成同一時刻不同預報時效的預報文件

一般是統(tǒng)計一個月的預報文件，預報時間間隔為12h?？梢宰孕袑懸粋€shell腳本實現(xiàn)。

二、運行$WRFDA/var/scripts/gen_be/gen_be_wrapper.ksh腳本后生成B矩陣

WRF 3DVar將產(chǎn)生B矩陣的可執(zhí)行程序封裝在ksh腳本gen_be_wrapper.ksh中，“天河二號”中指定ksh解釋器的路徑為 /BIGDATA1/app/ast/20180418-gcc-4.8.5/bin/ksh。在運行這個腳本之前，你需要在腳本中設(shè)置好相關(guān)的環(huán)境變量。

注意1：在stage0階段（從模式預報文件中計算集合擾動），只需指定一個core即可。

注意2：stage4階段耗時最長，可以將多個作業(yè)同時提交到一個node上，提高計算效率。