王盤(pán)興,趙輝,任律,羅璇,郭棟
(南京信息工程大學(xué)大氣科學(xué)學(xué)院,江蘇南京 210044)
閉合氣壓系統(tǒng)中心位置指數(shù)的計(jì)算方案
王盤(pán)興,趙輝,任律,羅璇,郭棟
(南京信息工程大學(xué)大氣科學(xué)學(xué)院,江蘇南京 210044)
王盤(pán)興等在“閉合氣壓系統(tǒng)環(huán)流指數(shù)的定義及計(jì)算”一文中定義了閉合氣壓系統(tǒng)的面積指數(shù)S、強(qiáng)度指數(shù)P和中心位置指數(shù)(λc,φc),并給出了它們的計(jì)算方案。但其中(λc,φc)的計(jì)算方案對(duì)低緯系統(tǒng)的計(jì)算存在明顯誤差,誤差原因是原方案中的極點(diǎn)球面坐標(biāo)系不適合低緯系統(tǒng)的計(jì)算。本文建立了一套原點(diǎn)位于搜索區(qū)Ω中心的新球面坐標(biāo)系,通過(guò)坐標(biāo)轉(zhuǎn)換實(shí)現(xiàn)了低緯閉合氣壓系統(tǒng)中心位置指數(shù)(λc,φc)的準(zhǔn)確計(jì)算。對(duì)7月南亞高壓和1月蒙古高壓中心位置指數(shù)(λc,φc)的實(shí)際計(jì)算表明,它既消除了低緯系統(tǒng)(南亞高壓)(λc,φc)原計(jì)算方案的明顯誤差,又保持了與高緯系統(tǒng)(蒙古高壓)(λc,φc)原計(jì)算方案計(jì)算結(jié)果的一致。因此,本文給出了適于計(jì)算任意緯度閉合氣壓系統(tǒng)中心位置指數(shù)(λc,φc)的計(jì)算方案。
閉合氣壓系統(tǒng);中心位置指數(shù);計(jì)算方案
在大氣環(huán)流異常與短期氣候預(yù)測(cè)研究中,全年或季節(jié)出現(xiàn)于特定地理區(qū)域的高壓、低壓系統(tǒng)被稱為大氣活動(dòng)中心[1]。大氣活動(dòng)中心的性狀可以用環(huán)流指數(shù)簡(jiǎn)明扼要地描述,常用的環(huán)流指數(shù)有它們的面積、強(qiáng)度和中心位置等。以南亞高壓、蒙古高壓為例,文獻(xiàn)[2-3]列舉了它們的多種環(huán)流指數(shù)定義[4-9];可見(jiàn),對(duì)同一氣壓系統(tǒng),不同作者定義環(huán)流指數(shù)的方法存在明顯差別。這給研究不同大氣活動(dòng)中心的相互關(guān)系帶來(lái)一定困難。
文獻(xiàn)[10]按統(tǒng)一的方法定義閉合氣壓系統(tǒng)的環(huán)流指數(shù)。對(duì)季、月平均氣壓場(chǎng)中具有大氣活動(dòng)中心性質(zhì)的氣壓系統(tǒng),均可按該方法定義一組環(huán)流指數(shù),它包括面積指數(shù)S、強(qiáng)度指數(shù)P、中心位置指數(shù)(λc,φc)。對(duì)某閉合氣壓系統(tǒng)多年序列,在詳細(xì)分析歷史氣壓場(chǎng)的基礎(chǔ)上容易確定其搜索域Ω(該閉合氣壓系統(tǒng)主體部分所在區(qū)域)和特征等值線f0(界定閉合氣壓系統(tǒng)主體部分的等值線);據(jù)此可以確定t年的計(jì)算域D(t),它是Ω上由t年場(chǎng)中f0線圍成的區(qū)域。S(t)是單位半徑球面上D(t)的面積;P(t)是壓差(高度差)df(t)=f(t)-f0在D(t)上的積分;(λc(t),φc(t))是D(t)上df(t)場(chǎng)對(duì)應(yīng)重力場(chǎng)重心所在經(jīng)、緯度。(ˉλc,ˉφc)是該氣壓系統(tǒng)的氣候中心位置指數(shù),它是(λc(t),φc(t))以P(t)為權(quán)重的多年平均。
文獻(xiàn)[10]借助地理坐標(biāo)系(λ,φ)和一個(gè)球面曲線直角坐標(biāo)系oxy完成了北、南半球500hPa層極渦環(huán)流指數(shù)S、P、(λc,φc)的計(jì)算,坐標(biāo)系oxy的原點(diǎn)o在北(南)極點(diǎn),x、y軸是λ=0、π/2的經(jīng)線。S、P的計(jì)算可直接在(λ,φ)系中完成,(λc,φc)的計(jì)算則在oxy系中通過(guò)球面直角三角形計(jì)算公式[11]完成。中、高緯度多種閉合氣壓系統(tǒng)(500hPa極渦[10]、10hPa極地渦旋[12-15]、1 000hPa蒙古高壓[16-17]和阿留申低壓[18]等)(λc,φc)的計(jì)算結(jié)果驗(yàn)證了文獻(xiàn)[10]計(jì)算方案的正確。
但是,球面三角形計(jì)算公式只對(duì)三邊和小于π的球面三角形[11](注:指對(duì)單位半徑球面)精確成立,上述計(jì)算方案不能用于低緯系統(tǒng)(λc,φc)計(jì)算。為了解決這個(gè)問(wèn)題,本文引進(jìn)了一個(gè)原點(diǎn)設(shè)在λ-φ平面上搜索區(qū)Ω中心(λ0,φ0)上的新地理坐標(biāo)系(λ′,φ′),給出一個(gè)計(jì)算(λc,φc)的新方案,對(duì)7月南亞高壓、1月蒙古高壓(λc,φc)的計(jì)算結(jié)果,驗(yàn)證了它對(duì)低、高緯系統(tǒng)(λc,φc)的正確計(jì)算均適用。
(λc,φc)的計(jì)算方案涉及S、P的定義和計(jì)算,為敘述清楚,首先簡(jiǎn)要給出它們的定義和計(jì)算式。
對(duì)某季、月平均氣壓場(chǎng)(位勢(shì)高度場(chǎng))的多年(L年)序列
定義了t年面積、強(qiáng)度指數(shù),
式中:λ、θ為經(jīng)度、余緯,余緯θ=π/2-φ。D(t)是t年計(jì)算域,它是搜索域Ω上由t年f0線圍成的區(qū)域。圖1給出了2008年1月蒙古高壓和1998年7月南亞高壓的Ω、f0、D(t)三者關(guān)系。
計(jì)算在單位半徑球面上的地理坐標(biāo)系(λ,φ)中進(jìn)行。
NCEP/NCAR、ECMWF再分析資料[19-20]的要素場(chǎng)給出在λ-φ平面上均勻經(jīng)度、余緯格點(diǎn)網(wǎng)上。其基本參數(shù)為格距Δ λ=Δ θ=π/72=2.5°,緯圈、經(jīng)線等分?jǐn)?shù)m=144、n=72;故格點(diǎn)(λ,θ)與格點(diǎn)序數(shù)(i,j)關(guān)系為
f場(chǎng)的多年序列的離散形式為
S(t)、P(t)的差分計(jì)算式為
式中:dS(i,j)為(i,j)格點(diǎn)所在面元的元面積;dP(i,j,t)為t年D(t)上(i,j)格點(diǎn)所在面元上的元壓差(元位勢(shì)高度差)。根據(jù)文獻(xiàn)[10],
按文獻(xiàn)[10],t年的閉合氣壓系統(tǒng)中心位置指數(shù)(λc(t),φc(t))是D(t)上壓差場(chǎng)dP(λ,θ,t)=f(λ,θ,t)-f0對(duì)應(yīng)重力場(chǎng)的“重心”位置。以r記球面上任意點(diǎn)(λ,θ)的位置矢量,rc(t)是球面上矢端在重心(λc(t),θc(t))上的矢量,r、rc(t)分別是從球面上選定(λ0,θ0)點(diǎn)出發(fā)至(λ(t),θ(t))、(λc(t),θc(t))的與大圓弧重合的矢量。則rc(t)由下積分方程定義而該系統(tǒng)氣候中心位置指數(shù)(ˉλc,ˉθc)則由下式定義,
圖1 閉合氣壓系統(tǒng)Ω、f0、D(t)三者關(guān)系(粗實(shí)線為特征等高線f0,粗虛線矩形區(qū)為搜索區(qū)Ω,陰影為計(jì)算域D(t);·為中心位置) a.2008年1月1 000hPa上蒙古高壓(f0=22dagpm);b.1998年7月100hPa上南亞高壓(f0=1 678dagpm)Fig.1 Relationships ofΩ,f0and D(t)for the closed pressure system s(The thick solid line is the typical line f0;the rectangle is searching area Ω;the shading is computational domain D(t);·is the center position) a.Mongolia Highat 1000hPa in January2008(f0=22dagpm);b.South Asia H igh at100hPa in July1998(f0=1 678dagpm)
為了計(jì)算t年中高緯度閉合系統(tǒng)的rc(t),文獻(xiàn)[10]建立了坐標(biāo)系OXYZ和oxy(圖2)。OXYZ是原點(diǎn)在地心的三維直角坐標(biāo)系,其基向量i、j、k是矢端分別位于單位半徑球面上(λ,θ)為(0,π/2)、(π/2,π/2)、(λ,0)點(diǎn)的單位向量。oxy是原點(diǎn)在北極的球面曲線直角坐標(biāo)系,其基向量e1、e2是矢端分別位于單位半徑球面上(λ,θ)為(0,1)、(π/2,1)點(diǎn)的單位向量。均為標(biāo)準(zhǔn)正交基。
由圖2可知,球面上任意點(diǎn)q的地理坐標(biāo)(λ,θ)與(X,Y,Z)間的轉(zhuǎn)換關(guān)系為
而據(jù)球面直角三角形計(jì)算公式[11],可導(dǎo)出點(diǎn)q的地理坐標(biāo)(λ,θ)向球面直角坐標(biāo)(x,y)的坐標(biāo)轉(zhuǎn)換關(guān)系(稱為正變換),
據(jù)此可將場(chǎng)f(λ,θ,t)轉(zhuǎn)換為f(x,y,t),并在坐標(biāo)系oxy中求得t年的該閉合氣壓系統(tǒng)的中心位置,
圖2 地理坐標(biāo)(λ,θ)與單位半徑球直角坐標(biāo)系OXYZ、oxy三者關(guān)系(q是球面上位于(λ,φ)的任意點(diǎn),它在OXYZ中的坐標(biāo)為(X,Y,Z),在oxy中的坐標(biāo)為(x,y))Fig.2 Relationship of the geographic coordinate(λ,θ)and the rectangular coordinate system s OXY and oxy in an unit radius sphere(q is an arbitrary point at(λ,φ)on the spherical,corresponding to(X,Y,Z)in OXYZ and(x,y)in oxy)
而由(x,y)向(λ,θ)的坐標(biāo)變換關(guān)系(稱為逆變換)為
據(jù)此可將(xc(t),yc(t))轉(zhuǎn)換為(λc(t),θc(t)),求得逐年分析對(duì)象的中心位置指數(shù)序列
文獻(xiàn)[10]將(2)、(12)式給出的閉合系統(tǒng)強(qiáng)度指數(shù)及位置參數(shù)多年序列
視為一個(gè)力學(xué)系統(tǒng),定義了該系統(tǒng)的氣候中心位置指數(shù)
閉合氣壓系統(tǒng)的搜索域Ω一般為(或可處理為)λ-θ平面上的矩形域[λw~λe,θn~θs]。這里λw、λe為Ω的西、東界所在經(jīng)度,θn、θs為北、南界所在余緯(圖1)。為準(zhǔn)確計(jì)算低緯系統(tǒng)的(λc,φc),建立新地理坐標(biāo)系(λ′,θ′),其北極點(diǎn)o′取在Ω的中心點(diǎn)(λ0,θ0)上,
對(duì)南亞高壓(圖1b),λ0=70°E,θ0=60°(即φ=30°N)。經(jīng)線λ′=0從o′出發(fā)沿λ=λ0經(jīng)線指向南極,經(jīng)線λ′=π/2從o′點(diǎn)出發(fā)垂直于λ=λ0經(jīng)線指向東(即為與φ=π/2-θ0緯圈相切的大圓的向東部分);它們分別是新球面曲線直角坐標(biāo)系o′x′y′的x、y軸。簡(jiǎn)言之,新地理坐標(biāo)系(λ′,θ′)是原地理坐標(biāo)系(λ,θ)繞地軸逆時(shí)針旋轉(zhuǎn)λ0,再將北極點(diǎn)沿λ0經(jīng)線南移θ0角得到。對(duì)新地理坐標(biāo)系(λ′,θ′)建立相應(yīng)的o′x′y′和OX′Y′Z′坐標(biāo)系,它們的相互關(guān)系全同于(λ,θ)、oxy與OXYZ。
原點(diǎn)位于地心的直角坐標(biāo)系OXYZ、OX′Y′Z′的基向量間關(guān)系為
由此得球面上任意點(diǎn)q(λ,θ)的坐標(biāo)(X,Y,Z)、(X′,Y′,Z′)間的正、逆轉(zhuǎn)換關(guān)系為
利用正轉(zhuǎn)換關(guān)系(19)式,可將求(λc,φc)的問(wèn)題轉(zhuǎn)入新坐標(biāo)系中進(jìn)行。因?yàn)樾伦鴺?biāo)系的北極點(diǎn)(即o′)位于Ω中心,只要Ω不很大(保證計(jì)算涉及的球面直角三角形三邊和不超過(guò)π),新坐標(biāo)系中(λ′c,φ′c)的計(jì)算將是精確的。而(λ′c,φ′c)可經(jīng)逆轉(zhuǎn)換關(guān)系(9)、(20)式返還原坐標(biāo)系,最終求得低緯系統(tǒng)中心位置指數(shù)(λc,φc)。
因(8)式中(λ0,θ0)可在整個(gè)球面域上取值(當(dāng)然必須是Ω中心點(diǎn)),故低緯系統(tǒng)(λc,φc)計(jì)算方案理論上也適合于高緯系統(tǒng)(λc,φc)的計(jì)算,是求氣壓系統(tǒng)(λc,φc)的普適計(jì)算方案。
低緯系統(tǒng)(λc,φc)計(jì)算方案的流程可歸結(jié)為:
1)計(jì)算(Xij,Yij,Zij)Ω、(xij,yij)Ω(據(jù)(9)、(11)式);
2)計(jì)算(S(t))t、(dSij)Ω(據(jù)(5)、(6)式);
3)計(jì)算(P(t))t、(dPij)Ω,t(據(jù)(5)、(6)式);
4)(Xij,Yij,Zij)→(據(jù)(19)式);
最后,將換算為(λc(t),換算為
在資料、Ω、f0與文獻(xiàn)[2,16]相同情況下,用本文(λc,φc)計(jì)算方案分別計(jì)算了60a(1948—2007年)7月100hPa南亞高壓和60a(1949—2008年)1月1 000hPa蒙古高壓中心位置指數(shù),記為
相應(yīng)地,文獻(xiàn)[2,16]用高緯計(jì)算方案求得的南亞高壓、蒙古高壓中心位置指數(shù)記為
則t年兩方案計(jì)算結(jié)果之差(高緯方案減低緯方案)為
其中,εC(t)是^C(t)、C(t)間在球面上的距離,
上述統(tǒng)計(jì)量(λc、φc、ε λc、ε φc、εC)的60a平均值(以上標(biāo)“-”表示)、60a絕對(duì)平均值(以“”表示用于度量?jī)煞N方案計(jì)算同一系統(tǒng)中心位置的差異。
由表1可知,兩方案求得7月100hPa南亞高壓中心位置指數(shù)60a平均存在明顯差異。達(dá)2.63°(它相當(dāng)于中心偏差292km);它主要由緯向差異ε λc(t)引起達(dá)3.05°)。因?yàn)?月100hPa南亞高壓λc(t)、φc(t)的極差Δ λc、Δ φc分別為18.7°、6.7°,相對(duì)偏差近似為1/6、1/10,故用高緯計(jì)算方案求出的南壓高壓中心位置指數(shù)(λc(t),φc(t))只是中心位置的近似(圖3)。
表1 兩種方案7月100 hPa南亞高壓(λc,φc)的比較(1948—2007年)Table1 Comparisons of(λc,φc)by the two schemes for100hPa South Asia High in July from1948to2007
圖3 兩種方案1948—2007年7月南亞高壓中心位置指數(shù)標(biāo)準(zhǔn)化序列比較(實(shí)線為低緯計(jì)算方案,虛線為高緯計(jì)算方案) Fig.3 Standardized center local index of South Asia High in July of1948—2007by the two calculation schemes(The solid line is low latitude scheme,and the dash line is high latitude one)
由表2可知,兩方案求得之1月蒙古高壓中心位置指數(shù)60a平均相同,精度為0.1°。給出了個(gè)別年份月平均圖上兩方案中心位置偏差(總是正值)的多年平均,其值為0.14°(它相當(dāng)于15km的偏差);由知,εC(t)主要由ε λc(t)引起。因?yàn)?月1 000hPa蒙古高壓λc(t)、φc(t)的極差(極大值減極小值)Δ λc、Δ φc分別為12.3°、12.8°,兩方案的相對(duì)偏差在1%左右,可以忽略不計(jì)。因此,低緯系統(tǒng)(λc,φc)計(jì)算方案適于高緯系統(tǒng)。
驗(yàn)證表明,低緯氣壓系統(tǒng)(λc,φc)計(jì)算方案既適合于低緯氣壓系統(tǒng)(λc,φc)的計(jì)算,也適合于高緯氣壓系統(tǒng)(λc,φc)的計(jì)算,是求氣壓系統(tǒng)(λc,φc)普適計(jì)算方案。
綜上,本文分析了文獻(xiàn)[10]定義的閉合氣壓系統(tǒng)中心位置指數(shù)(λc,φc)計(jì)算方案在低緯氣壓系統(tǒng)計(jì)算中存在的誤差,通過(guò)建立一組原點(diǎn)位于搜索區(qū)Ω中心的新的球面坐標(biāo)系,實(shí)現(xiàn)了低緯閉合氣壓系統(tǒng)中心位置指數(shù)(λc,φc)的準(zhǔn)確計(jì)算。對(duì)7月南亞高壓和1月蒙古高壓中心位置指數(shù)(λc,φc)的實(shí)際計(jì)算表明,它既消除了低緯系統(tǒng)(南亞高壓)
(λc,φc)原計(jì)算方案的明顯誤差,又保持了與高緯系統(tǒng)(蒙古高壓)(λc,φc)原計(jì)算方案計(jì)算結(jié)果的一致。因此,本文給出了適于計(jì)算任意緯度閉合氣壓系統(tǒng)中心位置指數(shù)(λc,φc)的計(jì)算方案。
致謝:本文使用資料由國(guó)家自然科學(xué)基金委地球科學(xué)部南京大氣資料服務(wù)中心提供,謹(jǐn)致謝忱!
表2 兩種方案1月1 000 hPa蒙古高壓(λc,φc)的比較(1949—2008年)Table2 Comparisons of(λc,φc)by the two Schemes for1 000hPa M ongolian H igh in January from1949to2008
[1] 《大氣科學(xué)辭典》編委會(huì).大氣科學(xué)辭典[M].北京:氣象出版社,1994:93.
[2] 陳延聰,王盤(pán)興,周國(guó)華,等.夏季南亞高壓的一組環(huán)流指數(shù)及其初步分析[J].大氣科學(xué)學(xué)報(bào),2009,32(6):101-107.
[3] 孫曉娟,王盤(pán)興,智海,等.蒙古高壓若干環(huán)流指數(shù)及與我國(guó)冬溫異常相關(guān)的分析和比較[J].高原氣象,待發(fā)表.
[4] 施能,朱乾根.南半球澳大利亞、馬斯克林高壓氣候特征及其對(duì)我國(guó)東部夏季降水的影響[J].氣象科學(xué),1995,15(2):20-27.
[5] 崔錦,楊修群.馬斯克林高壓的變化及其與ENSO的關(guān)系[J].氣象科學(xué),2005,25(5):441-449.
[6] 顧思南,楊修群.北半球繞極渦的變異及其與我國(guó)氣候異常的關(guān)系[J].氣象科學(xué),2006,26(2):135-142.
[7] 朱乾根,施能,吳朝琿,等.近百年北半球冬季大氣活動(dòng)中心的長(zhǎng)期變化及其與中國(guó)氣候[J].氣象學(xué)報(bào),1997,55(6):750-758.
[8] 龔道溢,王紹武.西伯利亞高壓的長(zhǎng)期變化及全球變暖可能影響的研究[J].地理學(xué)報(bào),1999,54(2):125-133.
[9] 侯亞紅,楊修群,李剛.冬季西伯利亞高壓變化特征及其與中國(guó)氣溫的關(guān)系[J].氣象科技,2007,35(5):646-650.
[10] 王盤(pán)興,盧楚翰,管兆勇,等.閉合氣壓系統(tǒng)環(huán)流指數(shù)的定義及計(jì)算[J].南京氣象學(xué)院學(xué)報(bào),2007,30(6):730-735.
[11] 數(shù)學(xué)手冊(cè)編寫(xiě)組.數(shù)學(xué)手冊(cè)[M].北京:高等教育出版社,2005:52.
[12] 管樹(shù)軒,王盤(pán)興,麻巨慧,等.北半球10hPa極地渦旋環(huán)流指數(shù)定義及分析[J].高原氣象,2009,28(4):777-785.
[13] 麻巨慧,王盤(pán)興,郭棟,等.南半球10hPa極地渦旋環(huán)流的多尺度變化特征分析[J].高原氣象,2009,28(6):1299-1307.
[14] Li Liping,Guan Shuxuan,Wang Panxing,et al.Definition and analysis of the circulation indices of polar vortex at 10hPa in the Northern Hemisphere[J].Sciences in Cold and A rid Regions,2009,1(6):549-558.
[15] Li Liping,Ma Juhui,Wang Panxing,et al.Multi-scale variation characteristics of polar vortex at10hPa in the Southern Hemisphere[J].Sciences in Cold and A rid Regions,2010,2(1):73-81.
[16] 劉晴晴,王盤(pán)興,徐祥德,等.蒙古高壓一組環(huán)流指數(shù)與我國(guó)同期氣侯異常關(guān)系分析[J].熱帶氣象學(xué)報(bào),待發(fā)表.
[17] 麻巨慧,王盤(pán)興,李麗平,等.“0801南方雪災(zāi)”與同期蒙古高壓中期活動(dòng)關(guān)系的分析[J].大氣科學(xué)學(xué)報(bào),2009,32(5):652-660.
[18] W ang Panxing,W ang Julian X L,W ang Yukun,et al.Circulation indices of the Aleutian low pressure system and its definitions and relationships to N.H.climate anomalies[J].GRL,審稿中
[19] Kalnay E,Kanamitsu M,Kistler R,et al.The NCEP/NCAR40-year reanalysis project[J].Bull Amer Meteor Soc,1996,77(3):437-471.
[20] Simmons A J,Gibswson J K.The ERA-40project Plan[M]//ECMW F.ERA-40projects report series.Shinfield Park,Reading:ECMW F,2000:1-62.
Calculation Method of Center Position Index for Closed Pressure Systems
WANG Pan-xing,ZHAO Hui,REN Lü,LUO Xuan,GUO Dong
(School of Atmospheric Sciences,NU IST,Nanjing 210044,China)
Areaindex S,strength index P and center longitude/latitude position index(λc,φc)have been defined in the paper Definition and calculation of three circulation indices for closed pressure systems written by Panxing Wang,et al.But the calculated index(λc,φc)for the lower latitude system s exist obvious errors,which is induced by the unsuitable polar spherical coordinate.Thus,a new spherical coordinate with the original point locating in the center of searching area is established,in such a way that the(λc,φc)is accurately calculated.The practical calculations of(λc,φc)for South Asian High in July and Mongolian High in January show that the evident errors are eliminated and the calculated results are in accordance with the original one.Therefore,this method can be applied to random closed pressure system s.
closed pressure system;center location index;calculation scheme
P434
A
1674-7097(2010)05-0520-07
2010-02-06;改回日期:2010-09-10
國(guó)家自然科學(xué)基金資助項(xiàng)目(40633018);江蘇省研究生培養(yǎng)創(chuàng)新工程創(chuàng)新計(jì)劃項(xiàng)目(CX08B_016Z)
王盤(pán)興(1943—),男,江蘇常州人,教授,博士生導(dǎo)師,研究方向?yàn)榇髿猸h(huán)流異常、短期氣候預(yù)測(cè),wangpx@nuist.edu.cn.
王盤(pán)興,趙輝,任律,等.閉合氣壓系統(tǒng)中心位置指數(shù)的計(jì)算方案[J].大氣科學(xué)學(xué)報(bào),2010,33(5):520-526.Wang Pan-xing,Zhao Hui,Ren Lü,et al.Calculation method of centerposition index for closed pressure systems[J].Trans Atmos Sci,2010,33(5):520-526.
(責(zé)任編輯:倪東鴻)