2 模擬實(shí)驗(yàn)方案設(shè)計(jì)
　　WRF 模式是由美國(guó)NCAR、NOAA 聯(lián)合美國(guó)多所大學(xué)和研究機(jī)構(gòu)開(kāi)發(fā)的完全可壓縮、非靜力中尺度模式，模式的水平網(wǎng)格采用ArakawaC 格式。模式包含了陸面過(guò)程、行星邊界層物理過(guò)程、次網(wǎng)格湍流擴(kuò)散過(guò)程、大氣和地球表層輻射、微物理過(guò)程和積云對(duì)流等物理過(guò)程。
　　CALMET 模式是美國(guó)環(huán)保局（Environmental ProtectionAgency）開(kāi)發(fā)的復(fù)雜地形動(dòng)力診斷模式，主要用來(lái)對(duì)氣象數(shù)據(jù)包括MM5 等氣象數(shù)值模式及常規(guī)地面、探空觀測(cè)等，進(jìn)行地形動(dòng)力學(xué)、傾斜流和熱動(dòng)力學(xué)等診斷。以發(fā)散最小化原理求解三維風(fēng)場(chǎng)，根據(jù)湍流參數(shù)化方法，計(jì)算湍流尺度參數(shù)。

圖2 CALMET模擬區(qū)域地形分布及測(cè)風(fēng)塔分布

　　根據(jù)Liu 和Yocke 的理論，CALMET 模式假設(shè)地形強(qiáng)迫產(chǎn)生的垂直速度w 可由公式（2.1-2.2）求得。地形強(qiáng)迫對(duì)水平風(fēng)場(chǎng)的影響，體現(xiàn)在為使三維散度小于某一指定的臨界值為目的進(jìn)行的水平風(fēng)場(chǎng)調(diào)整上。

　　V 是模擬網(wǎng)格平均風(fēng)，ht 是地形高度，k 是穩(wěn)定度相關(guān)的衰減系數(shù)，z 是垂直坐標(biāo)，N 是布倫特—維塞拉頻率，|V|是模擬網(wǎng)格平均風(fēng)速，w 為笛卡爾坐標(biāo)系物理垂直風(fēng)分量。斜坡氣流的速度采用經(jīng)驗(yàn)的方法：

　　其中e S 是斜坡氣流的平衡風(fēng)速， e L 是平衡尺度。障礙物阻擋的熱力和動(dòng)力效應(yīng)用局地弗勞德數(shù)Fr 衡量：

　　　其中△ h1 是障礙物的有效高度。如果局地弗勞德數(shù)小于等于臨界弗勞德數(shù)切網(wǎng)格點(diǎn)風(fēng)速有上坡分量，則調(diào)整風(fēng)向與地形的切線一致，風(fēng)速不變；反之，不調(diào)整。
　　本實(shí)驗(yàn)中首先采用WRF 模式進(jìn)行模擬，模式積分計(jì)算36 小時(shí)，根據(jù)全球環(huán)流模式再分析資料（NCEP）進(jìn)行降尺度，然后利用calwrf 模塊對(duì)WRF 模擬的結(jié)果的后24 小時(shí)進(jìn)行分析處理，并形成CALMET 可用的格式，最后使用CALMET模式進(jìn)行動(dòng)力診斷模擬分析。兩個(gè)模式結(jié)果輸出時(shí)間間隔均為1 小時(shí)。