GrADS第6章變量和函數(shù)

1、第第6 6章章 變量和函數(shù)變量和函數(shù) 在在GrADS運(yùn)行環(huán)境中，參與操作運(yùn)行環(huán)境中，參與操作的數(shù)據(jù)對(duì)象可以是變量、函數(shù)和的數(shù)據(jù)對(duì)象可以是變量、函數(shù)和表達(dá)式形式。表達(dá)式形式。 6.1 6.1 變量名變量名 n完全的變量名形式完全的變量名形式n abbrev.file#（dimexpr， dimexpr，）n 其中，其中，abbrev：是數(shù)據(jù)描述文件中給：是數(shù)據(jù)描述文件中給出的變量名縮寫；出的變量名縮寫；file#：為包含此變量的：為包含此變量的已打開的文件序號(hào)。已打開的文件序號(hào)。 例如：ga-open d:model.ctl open d:model.le.ctl 上例中，上例中，model.c

2、tl文件的序列號(hào)為文件的序列號(hào)為1，model.le.ctl文件的序列號(hào)為文件的序列號(hào)為2。序列號(hào)。序列號(hào)1為缺為缺省情況，可以省略。省情況，可以省略。 括號(hào)內(nèi)的括號(hào)內(nèi)的dimexpr：是對(duì)當(dāng)前維數(shù)環(huán)境：是對(duì)當(dāng)前維數(shù)環(huán)境進(jìn)行的局域維數(shù)設(shè)置表達(dá)式，該表達(dá)式僅進(jìn)行的局域維數(shù)設(shè)置表達(dá)式，該表達(dá)式僅對(duì)該變量的維數(shù)環(huán)境進(jìn)行局域修正，因而對(duì)該變量的維數(shù)環(huán)境進(jìn)行局域修正，因而僅對(duì)固定的維有效。僅對(duì)固定的維有效。 設(shè)置形式可以采用以下兩種方式：設(shè)置形式可以采用以下兩種方式： 絕對(duì)維數(shù)表達(dá)式為：絕對(duì)維數(shù)表達(dá)式為： X|Y|Z|T|LON|LAT|LEV|TIME = value （value為絕對(duì)維數(shù)值）；為絕

3、對(duì)維數(shù)值）；相對(duì)維數(shù)表達(dá)式為：相對(duì)維數(shù)表達(dá)式為： X|Y|Z|T|LON|LAT|LEV|TIME +/- offset （offset為相對(duì)于當(dāng)前維數(shù)環(huán)境設(shè)置的偏差維為相對(duì)于當(dāng)前維數(shù)環(huán)境設(shè)置的偏差維數(shù)值）。數(shù)值）。 例如例如： z.3（lev=500）tv.1（time-12hr）rh 缺省的當(dāng)前文件中的變量缺省的當(dāng)前文件中的變量rhq.2（t-1,lev=850）2號(hào)文件中相對(duì)于當(dāng)前時(shí)刻的號(hào)文件中相對(duì)于當(dāng)前時(shí)刻的前一時(shí)刻，高度為前一時(shí)刻，高度為850hPa等等壓面上的變量壓面上的變量q表示文件表示文件3中高度為中高度為 500hPa 等壓面上的變量等壓面上的變量Z相對(duì)于當(dāng)前時(shí)刻之前相對(duì)于當(dāng)

4、前時(shí)刻之前12小時(shí)小時(shí)時(shí)刻的時(shí)刻的1號(hào)文件中的變量號(hào)文件中的變量tvndset model.le.dat noptions little_endian cray_32bit_ieeenUNDEF -2.56E33 nTITLE 5 Days of Sample Model OutputnXDEF 72 LINEAR 0.0 5.0 nYDEF 46 LINEAR -90.0 4.0 nZDEF 7 LEVELS 1000 850 700 500 300 200 100nTDEF 5 LINEAR 02JAN1987 1DY nvars 8 nps 0 99 Surface Pressurenu

5、 7 99 U Windsnv 7 99 V Windsnz 7 99 Geopotential Heightsnt 7 99 Temperaturenq 5 99 Specific Humiditynts 0 99 Surface Temperaturenp 0 99 PrecipitationnENDVARS 2. 舉例說(shuō)明舉例說(shuō)明ndset model.grbntitle Sample Model Data for lats4d Tutorialnundef 1e+20ndtype gribnindex model.gmpnxdef 72 linear 0.000000 5.000000

6、nydef 46 linear -90.000000 4.000000nzdef 7 levels 1000 850 700 500 300 200 100 ntdef 5 linear 0Z1jan1987 1dynvars 8nps 0 1, 1, 0, 0 Surface pressure hPanua 7 33,100 Eastward wind m/snva 7 34,100 Northward wind m/snzg 7 7,100 Geopotential height mnta 7 11,100 Air Temperature Knhus 7 51,100 Specific h

7、umidity kg/kgnts 0 11,105, 2 Surface (2m) air temperature Knpr 0 59, 1, 0, 0 Total precipitation rate kg/(m2*s)nendvarsga-open d:model.ctl open d:model.le.ctl set t 1 將時(shí)間維數(shù)固定在第將時(shí)間維數(shù)固定在第1個(gè)時(shí)次個(gè)時(shí)次 d z.2(lev=500) set z 3 設(shè)置垂直高度層次為第設(shè)置垂直高度層次為第3個(gè)層次個(gè)層次,即即700hPa d z(z+1)同時(shí)打開兩個(gè)文件同時(shí)打開兩個(gè)文件: model.ctl文件序列號(hào)為文件序列號(hào)為

8、1，model.le.ctl 文件序列號(hào)為文件序列號(hào)為2.顯示顯示2號(hào)文件中號(hào)文件中500hPa等壓面等壓面 上的變量上的變量z(絕對(duì)維數(shù)絕對(duì)維數(shù))顯示顯示1號(hào)文件中相對(duì)于第號(hào)文件中相對(duì)于第3個(gè)層次的后一個(gè)層次的后一個(gè)層次，即個(gè)層次，即500hPa面上的變量面上的變量z （相對(duì)（相對(duì)維數(shù)）維數(shù)）6.2 6.2 表達(dá)式表達(dá)式1 1、表達(dá)式的構(gòu)成、表達(dá)式的構(gòu)成n運(yùn)算符：運(yùn)算符： + + 加，加，- - 減，減，* * 乘，乘，/ / 除；除；n運(yùn)算域：運(yùn)算域： 變量變量, , 函數(shù)和常數(shù)；函數(shù)和常數(shù)；n括號(hào)：用于控制運(yùn)算的次序。括號(hào)：用于控制運(yùn)算的次序。2. 2. 運(yùn)算規(guī)則運(yùn)算規(guī)則n對(duì)相同網(wǎng)格點(diǎn)上

9、的不同變量進(jìn)行運(yùn)算時(shí)，對(duì)相同網(wǎng)格點(diǎn)上的不同變量進(jìn)行運(yùn)算時(shí)，只要有一個(gè)變量在某格點(diǎn)的值為缺測(cè)，則只要有一個(gè)變量在某格點(diǎn)的值為缺測(cè)，則該網(wǎng)格點(diǎn)的運(yùn)算結(jié)果為缺測(cè)值；當(dāng)被該網(wǎng)格點(diǎn)的運(yùn)算結(jié)果為缺測(cè)值；當(dāng)被0 0除除時(shí)的結(jié)果也設(shè)為缺測(cè)。時(shí)的結(jié)果也設(shè)為缺測(cè)。n對(duì)多個(gè)數(shù)據(jù)文件的數(shù)據(jù)做運(yùn)算時(shí)，兩種數(shù)對(duì)多個(gè)數(shù)據(jù)文件的數(shù)據(jù)做運(yùn)算時(shí)，兩種數(shù)據(jù)網(wǎng)格要一致，即運(yùn)算對(duì)象的維數(shù)變化范據(jù)網(wǎng)格要一致，即運(yùn)算對(duì)象的維數(shù)變化范圍要一致。圍要一致。n如某一格點(diǎn)的變動(dòng)維數(shù)多于其他格點(diǎn)，則如某一格點(diǎn)的變動(dòng)維數(shù)多于其他格點(diǎn)，則具有較少變動(dòng)維數(shù)的格點(diǎn)其維數(shù)環(huán)境將被具有較少變動(dòng)維數(shù)的格點(diǎn)其維數(shù)環(huán)境將被拓展，以便于運(yùn)算操作。拓展，以便于運(yùn)算操作。

10、例如：例如：z - z(t-1)t(lev=200)-t(lev=850) ave(z,t=1,t=5) z - ave(z,lon=0,lon=360,-b)變量變量Z的時(shí)間變化的時(shí)間變化200hPa與與850hPa等壓面上等壓面上的變量的變量t的變化的變化變量變量Z從第從第1時(shí)次到第時(shí)次到第5時(shí)次的平均時(shí)次的平均變量變量Z 的緯向偏差的緯向偏差6.3 6.3 函數(shù)函數(shù)1、調(diào)用方式、調(diào)用方式 采用直接引用方式，形式如：采用直接引用方式，形式如： 函數(shù)名函數(shù)名(參數(shù)參數(shù)1,參數(shù)參數(shù)2,)2. 2. 常用函數(shù)常用函數(shù)（1）ave函數(shù)函數(shù)n格式：格式： ave(expr,dexpr1,dexpr2

11、,tincr)n功能：通用的求平均函數(shù)。功能：通用的求平均函數(shù)。n說(shuō)明：說(shuō)明：expr是由是由dexpr1和和dexpr2定義的維數(shù)定義的維數(shù) 范圍內(nèi)求平均的量。范圍內(nèi)求平均的量。 例如：ga-open model.ctl set lev 500 d ave(z,t=1,t=5) 顯示顯示500hPa等壓面上從第等壓面上從第1時(shí)次到第時(shí)次到第5時(shí)次時(shí)次變量變量z的平均量的平均量set parea 1 8 1 6set grid offset grads offset xlint 20set xlopts 1 2 0.18set ylopts 1 2 0.18set lon 60 180set

12、lat 0 60set lev 200set t 1define a=ave(v,t=7,t=660,12)d aprintc例如：例如： 求緯向平均：求緯向平均： 如果如果ave(expr,lon=0,lon=360)在端點(diǎn)求在端點(diǎn)求 了兩次平均；了兩次平均； 如若使終端點(diǎn)取半數(shù)權(quán)重，可用：如若使終端點(diǎn)取半數(shù)權(quán)重，可用： ave(expr,lon=0,lon=360,-b)。 set xyrev onset lev 500set t 1set lat 0 90 set x 1define a=ave(v,lon=0,lon=360,-b)d aset cmark 0set ccolor 2d

13、 const(a,0)(2) const函數(shù)函數(shù)n格式：格式：const（expr,constant）n功能：設(shè)置部分網(wǎng)格點(diǎn)的值取為常數(shù)功能：設(shè)置部分網(wǎng)格點(diǎn)的值取為常數(shù)constantn說(shuō)明：所有非缺測(cè)格點(diǎn)處的說(shuō)明：所有非缺測(cè)格點(diǎn)處的expr值取為常數(shù)，值取為常數(shù)，flag為選項(xiàng)為選項(xiàng)，如果加上選項(xiàng)，如果加上選項(xiàng)-a，則所有網(wǎng)格點(diǎn)，則所有網(wǎng)格點(diǎn)值均設(shè)定為指定的常數(shù)，值均設(shè)定為指定的常數(shù)，如果加上選項(xiàng)如果加上選項(xiàng)-u，則，則只把只把缺測(cè)格點(diǎn)處缺測(cè)格點(diǎn)處的的expr值設(shè)定為常數(shù)。值設(shè)定為常數(shù)。n該函數(shù)對(duì)格點(diǎn)和臺(tái)站資料均適用。該函數(shù)對(duì)格點(diǎn)和臺(tái)站資料均適用。（3）aave函數(shù)函數(shù)n格式：格式： aav

14、e(expr,xdim1,xdim2,ydim1,ydim2)n功能：求面積平均。功能：求面積平均。n說(shuō)明：說(shuō)明：expr: 任何表達(dá)式任何表達(dá)式xdim1: X 維數(shù)方向的起始維數(shù)表達(dá)式維數(shù)方向的起始維數(shù)表達(dá)式xdim2: X 維數(shù)方向的終止維數(shù)表達(dá)式維數(shù)方向的終止維數(shù)表達(dá)式y(tǒng)dim1: Y 維數(shù)方向的起始維數(shù)表達(dá)式維數(shù)方向的起始維數(shù)表達(dá)式y(tǒng)dim2: Y 維數(shù)方向的終止維數(shù)表達(dá)式維數(shù)方向的終止維數(shù)表達(dá)式set lev 850set t 1 48set x 1set y 1define a=aave(v,lon=100,lon=120,lat=30,lat=40)d a（4）mag函數(shù)函數(shù)n

15、格式：格式：mag(aexpr,bexpr)n功能：計(jì)算功能：計(jì)算n說(shuō)明：說(shuō)明：aexpr表示表示a變量，變量，bexpr表示表示b變量，變量，該函數(shù)對(duì)格點(diǎn)和臺(tái)站資料均適合。該函數(shù)對(duì)格點(diǎn)和臺(tái)站資料均適合。22)exp()exp(rbraset lon 60 200set lat -20 80d mag(u.1,v.2)（5）hcurl函數(shù)函數(shù)n格式：格式：hcurl(uexpr,vexpr)n功能：計(jì)算功能：計(jì)算垂直渦度垂直渦度。n說(shuō)明：說(shuō)明： uexpr表示表示U風(fēng)速分量，風(fēng)速分量，vexpr表示表示V風(fēng)速分風(fēng)速分量，風(fēng)速單位用量，風(fēng)速單位用m/s。邊界上的渦度值設(shè)定。邊界上的渦度值設(shè)定為缺

16、測(cè)。為缺測(cè)。set lon 60 200set lat -20 80 d hcurl(u,v) （6）hdivg函數(shù)函數(shù)n格式：格式：hdivg(uexpr,vexpr)n功能：計(jì)算功能：計(jì)算水平散度水平散度。n說(shuō)明：說(shuō)明： uexpr表示表示U風(fēng)速分量，風(fēng)速分量，vexpr表示表示V風(fēng)速分風(fēng)速分量，風(fēng)速單位用量，風(fēng)速單位用m/s。set lon 60 200set lat -20 80 d hdivg (u,v) （7）skip函數(shù)函數(shù)n格式：格式：skip（expr,skipx,skipy）n功能：設(shè)定樣本的取樣密度。功能：設(shè)定樣本的取樣密度。n說(shuō)明：說(shuō)明：skipx,skipy 數(shù)值決定

17、數(shù)值決定X和和Y方向的取方向的取樣密度（取值樣密度（取值1可以省略不給）可以省略不給）*該函數(shù)主要用于對(duì)該函數(shù)主要用于對(duì)矢量場(chǎng)的稀疏化矢量場(chǎng)的稀疏化顯示。顯示。 d u.1;v.2d skip(u.1,2,2);skip(v.2,2,2)（8）cdiff函數(shù)n格式：格式：cdiff（expr,dim）n功能：中央差分函數(shù)n說(shuō)明：expr表示進(jìn)行差分運(yùn)算的量， dim表示進(jìn)行差分運(yùn)算的維數(shù)方向， 為X，Y，Z中的任一個(gè)字符。 邊界格點(diǎn)的運(yùn)算值設(shè)為缺測(cè)。邊界格點(diǎn)的運(yùn)算值設(shè)為缺測(cè)。n例如：計(jì)算溫度平流define dtx=cdiff(t,x)define dty=cdiff(t,y)define d

18、x=cdiff(lon,x)*3.1416/180define dy=cdiff(lat,y)*3.1416/180d -1*(u*dtx)/(cos(lat*3.1416/180)*dx)+v*dty/dy)/6.37e6（9 9）smth9smth9函數(shù)函數(shù)格式：格式：smth9(expr)smth9(expr)功能：功能：9 9點(diǎn)平滑函數(shù)點(diǎn)平滑函數(shù)說(shuō)明：說(shuō)明：exprexpr表示須平滑的表達(dá)式。表示須平滑的表達(dá)式。 該函數(shù)主要用于改善圖形的輸出質(zhì)量，該函數(shù)主要用于改善圖形的輸出質(zhì)量， 使其更加美觀。使其更加美觀。當(dāng)格點(diǎn)資料為當(dāng)格點(diǎn)資料為1 1維時(shí)，維時(shí)，9 9點(diǎn)平滑退化為點(diǎn)平滑退化為3

19、3點(diǎn)平滑。點(diǎn)平滑。d smth9(b1)（11）vint函數(shù)n格式：vint(psexpr,expr,top)n功能：質(zhì)量加權(quán)垂直積分函數(shù)n說(shuō)明：psexpr表示地面氣壓變量， expr表示被積變量， top表示積分上限(須為常數(shù)，hPa) 本函數(shù)要求定義水平維數(shù)環(huán)境。Vint( psexpr, expr, top)=f/g* sum(expr * Delta(level)f=100，g=9.8，sum是求和的意思是求和的意思 n例如：計(jì)算整層水汽通量fuq;fvq。open e:bsuq.ctlopen e:bspres.ctlset gxout fwriteset fwrite e:bsm

20、onzuq.grdi=1while(i=660)set t i set x 1 144set y 1 73set z 1define fuq=vint(v.2(z=1)/100,v,300)set t iset x 1 144set y 1 73set z 1d fuqi=i+1 endwhiledisable fwritereinit1998年夏季整層水汽通量矢量（箭頭）、水汽通量值（陰影）和水汽路徑（黑色曲線）站點(diǎn)畫圖站點(diǎn)畫圖（1）寫成）寫成站點(diǎn)格式站點(diǎn)格式grd數(shù)據(jù)；數(shù)據(jù)；（2）生成）生成map文件文件；（3）建立）建立格點(diǎn)文件格點(diǎn)文件（建立后可以繼續(xù)使用）；（建立后可以繼續(xù)使用）；（4

21、）將站點(diǎn)資料通過）將站點(diǎn)資料通過插值函數(shù)插值函數(shù)插值到格點(diǎn)文件插值到格點(diǎn)文件上，然后再畫圖。上，然后再畫圖。 文件的時(shí)間一定要與站點(diǎn)數(shù)據(jù)時(shí)間一致文件的時(shí)間一定要與站點(diǎn)數(shù)據(jù)時(shí)間一致。 范圍：范圍：lat：15 55 N lon：70 140 E 格距：格距：1*1ngrid.grd文件的每個(gè)點(diǎn)上均賦值為文件的每個(gè)點(diǎn)上均賦值為1格點(diǎn)文件grid.grd ：第二章數(shù)據(jù)處理中第二章數(shù)據(jù)處理中70-73頁(yè)寫出的格點(diǎn)數(shù)據(jù)頁(yè)寫出的格點(diǎn)數(shù)據(jù)grid.grd，對(duì)應(yīng)對(duì)應(yīng)ctl為為grid.ctl與該格點(diǎn)文件相對(duì)應(yīng)的數(shù)據(jù)描述文件grid.ctl如下：dset d:datagrid.grdundef -999.0ti

22、tle Sample GRID Dataxdef 71 linear 70 1ydef 41 linear 15 1zdef 1 linear 500 1tdef 1 linear jul1951 1movars 1g 0 99 grid data endvars有了這個(gè)格點(diǎn)文件及與其相對(duì)應(yīng)的數(shù)據(jù)描述文件后，就可以通過有了這個(gè)格點(diǎn)文件及與其相對(duì)應(yīng)的數(shù)據(jù)描述文件后，就可以通過GrADS插值函數(shù)將站點(diǎn)資料插值到該格點(diǎn)文件上然后使用了。插值函數(shù)將站點(diǎn)資料插值到該格點(diǎn)文件上然后使用了。 (12) oacres函數(shù)函數(shù)n格式：格式：oacres(gexpr, sexpr ) n功能：功能： cressm

23、an 插值函數(shù)插值函數(shù)n說(shuō)明：說(shuō)明： 常用于常用于 站點(diǎn)資料插值為格點(diǎn)資料。站點(diǎn)資料插值為格點(diǎn)資料。 gexpr 大于站點(diǎn)數(shù)據(jù)范圍的格點(diǎn)數(shù)據(jù)變量名大于站點(diǎn)數(shù)據(jù)范圍的格點(diǎn)數(shù)據(jù)變量名sexpr 站點(diǎn)數(shù)據(jù)變量名站點(diǎn)數(shù)據(jù)變量名radii 影響半徑，影響半徑， 默認(rèn)值為默認(rèn)值為: 10,7,4,2,1 n為避免插值出現(xiàn)虛假結(jié)果為避免插值出現(xiàn)虛假結(jié)果, , 參考網(wǎng)格的間參考網(wǎng)格的間隔與站點(diǎn)間距相近為好隔與站點(diǎn)間距相近為好將站點(diǎn)資料通過插值函數(shù)將站點(diǎn)資料通過插值函數(shù)插值插值到到某個(gè)格點(diǎn)文件某個(gè)格點(diǎn)文件上，上，然后再使用。然后再使用。與該格點(diǎn)文件相對(duì)應(yīng)的數(shù)據(jù)描述文件grid.ctl如下：dset d:data

24、grid.grdundef -999.0title Sample GRID Dataxdef 71 linear 70 1ydef 41 linear 15 1zdef 1 linear 500 1tdef 1 linear jul1951 1movars 1g 0 99 grid data endvarsdefine a=oacres(g,rain.2) 站點(diǎn)數(shù)據(jù)變量名站點(diǎn)數(shù)據(jù)變量名open d:datagrid.ctl open d:datarain7.ctl （13）maskout函數(shù)函數(shù)n格式：格式：maskout(expr,mask) n功能：功能： 標(biāo)記函數(shù)。當(dāng)標(biāo)記函數(shù)。當(dāng) mas

25、k所在的網(wǎng)格點(diǎn)值取為所在的網(wǎng)格點(diǎn)值取為負(fù)值時(shí)，對(duì)應(yīng)格點(diǎn)上的負(fù)值時(shí)，對(duì)應(yīng)格點(diǎn)上的 expr的值設(shè)為缺測(cè)值，不的值設(shè)為缺測(cè)值，不參與運(yùn)算或畫圖。參與運(yùn)算或畫圖。 即：只畫出即：只畫出mark大于大于0的的expr。n說(shuō)明：說(shuō)明： expr 需要處理的變量名需要處理的變量名mask 標(biāo)記變量名標(biāo)記變量名 常用于輸出常用于輸出 mask 代表的陸地區(qū)或海洋區(qū)；代表的陸地區(qū)或海洋區(qū)；也常用于我國(guó)站點(diǎn)資料繪圖。也常用于我國(guó)站點(diǎn)資料繪圖。open d:datagrid.ctlopen d:datar7.qht.ctlset grads offenable print d:datarain7.qht.gmfs

26、et mpdset hires cnworldset lon 73 135set lat 15 55define a=oacres(g(t=1),p.2,1.5) （“oacres”為插值函數(shù)）為插值函數(shù)）define a1=maskout(a,g-0.5) （“maskout”為標(biāo)記函數(shù)）為標(biāo)記函數(shù)）define aa=smth9(a1) （“smth9”9點(diǎn)平滑函數(shù)）點(diǎn)平滑函數(shù)）set grid offset gxout shadedset black -0.2976 0.2976d aaset gxout contourset cint 0.1d aa一般將格點(diǎn)數(shù)據(jù)作為第一個(gè)打開的文一般

27、將格點(diǎn)數(shù)據(jù)作為第一個(gè)打開的文件，在當(dāng)前件，在當(dāng)前缺省文件為格點(diǎn)數(shù)據(jù)缺省文件為格點(diǎn)數(shù)據(jù)時(shí)時(shí), 可可以用以用define命令將站點(diǎn)數(shù)據(jù)的客觀分命令將站點(diǎn)數(shù)據(jù)的客觀分析結(jié)果保存到變量中。析結(jié)果保存到變量中。printdisable printreinit;open e:datagrid.ctlopen e:datat.ctlset lon 70 140set lat 15 55define aa=oacres(g.1,t1.2,50,40,30,20,10,5,2)define aaa=maskout(aa,g.1-0.5)set gxout shadedd aaaset gxout contour

28、d aaaprintim e:datasample.png whitereinit；n地形高度資料：dxgd.ctl 變量名：hn700h Pa上，地形大于3000米以上的區(qū)域是沒有風(fēng)場(chǎng)的，所以可以利用maskout函數(shù)將該地區(qū)風(fēng)場(chǎng)不繪制出來(lái)reinitopen E:bsmonu.ctlopen E:bsmonv.ctlopen e:dxgd.ctlenable print E:uv.gmfset lev 700set t 1set lon 70 140set lat 10 60d maskout(v,3000-h.3(z=1,t=1);maskout(v.2,3000-h.3(z=1,t=1

29、)printdisable printreinit;dset e:dxgd.grdundef -9.99E+33xdef 144 linear 0 2.5ydef 73 linear -90 2.5zdef 1 linear 1 1tdef 1 linear 01jan1951 1movars 1h 1 99 hendvarsset t 1set lon 70 140set lat 10 60set gxout shadedset clevs 3000set ccols 0 1d h.3(z=1)set lev 700set t 1set lon 70 140set lat 10 60set

30、gxout vectord v;v.2printdisable printreinit;其他函數(shù)n數(shù)學(xué)函數(shù)數(shù)學(xué)函數(shù) 統(tǒng)計(jì)函數(shù)：統(tǒng)計(jì)函數(shù)： 網(wǎng)格處理函數(shù)n常數(shù) const (expr, const, | -u | -a )n剔除 maskout (expr, mask_expr )n跳點(diǎn)skip (expr, skip_x, skip_y )n加權(quán)9點(diǎn)平滑點(diǎn)平滑 smth9 ( expr ) 1-D 3點(diǎn)平滑點(diǎn)平滑n中央差分 cdiff ( expr, dim )物理量計(jì)算函數(shù)虛溫,相對(duì)濕度求比濕相對(duì)濕度求比濕-虛溫,相對(duì)濕度求溫度相對(duì)濕度求溫度-垂直渦度-中央差求水平散度-全風(fēng)速-6.4 6.

31、4 定義變量定義變量1、定義形式、定義形式define varname = expr 其中，其中，varname為新變量名，為新變量名，expr為表達(dá)為表達(dá)式。所定義的新變量可以用于隨后的表達(dá)式。所定義的新變量可以用于隨后的表達(dá)式中。式中。define命令命令只只適用于適用于格點(diǎn)數(shù)據(jù)格點(diǎn)數(shù)據(jù)2. 存儲(chǔ)形式存儲(chǔ)形式 新定義的變量新定義的變量varname只存在于內(nèi)存中，只存在于內(nèi)存中，所以建議不要定義過大的維數(shù)范圍。所以建議不要定義過大的維數(shù)范圍。 3.3.維數(shù)環(huán)境的設(shè)置維數(shù)環(huán)境的設(shè)置n用戶定義的變量可以有用戶定義的變量可以有04個(gè)變化的維數(shù)。個(gè)變化的維數(shù)。假設(shè)定義一個(gè)假設(shè)定義一個(gè)4維的變量，在

32、維的變量，在“define”命命令執(zhí)行后，記住改變維數(shù)環(huán)境，使之變化令執(zhí)行后，記住改變維數(shù)環(huán)境，使之變化的維數(shù)小于的維數(shù)小于4。n當(dāng)當(dāng)4個(gè)個(gè)維數(shù)同時(shí)變化時(shí)，維數(shù)同時(shí)變化時(shí)，“define”命令是命令是GrADS中所有命令中中所有命令中唯一唯一可行的命令。命令。 open e:uv.ctlset lon 60 150set lat 0 40set lev 850 200set t 1 48define temp=mag(u,v)set t 5set lev 200d tempcreinit定義后使用時(shí)改變其維數(shù)環(huán)境定義后使用時(shí)改變其維數(shù)環(huán)境定義數(shù)組維數(shù)環(huán)境（定義數(shù)組維數(shù)環(huán)境（4維）維）此時(shí)顯示的

33、變量此時(shí)顯示的變量temp，其時(shí)間和層次維是固定的，其時(shí)間和層次維是固定的定義數(shù)組定義數(shù)組n如果定義的變量具有某些固定的維數(shù)，隨如果定義的變量具有某些固定的維數(shù)，隨后使用這個(gè)變量，無(wú)論維數(shù)環(huán)境如何改變，后使用這個(gè)變量，無(wú)論維數(shù)環(huán)境如何改變，此變量在固定維數(shù)上保持不變。此變量在固定維數(shù)上保持不變。 set lon -180 0 set lat 0 90 set lev 500 set t 1 define zave = ave(z,t=1,t=5) X、Y維度不固定Z、t維度固定set t 5 set lev 200 d zave 顯示的是顯示的是500 h Pa，1-5時(shí)刻平均時(shí)刻平均的值的值

34、n當(dāng)定義的變量具有某些變化的維數(shù)環(huán)當(dāng)定義的變量具有某些變化的維數(shù)環(huán)境，然后再將這些維數(shù)固定起來(lái)，這境，然后再將這些維數(shù)固定起來(lái)，這時(shí)所得的變量將被固定在該維數(shù)上。時(shí)所得的變量將被固定在該維數(shù)上。set lon 0 180set lat -90 90 set lev 500 set t 5 define temp = z set lat 40 d temp Y維固定為40N，顯示的只能是40N 的值set lon 200set lat 40 d temp 超出定義范圍，超出定義范圍，顯示為缺測(cè)顯示為缺測(cè)n用戶也可以使用局地維數(shù)環(huán)境。用戶也可以使用局地維數(shù)環(huán)境。n如果局地環(huán)境是在定義變量時(shí)變化的如果局地環(huán)境是在定義變量時(shí)變化的維數(shù)維數(shù)環(huán)境內(nèi)環(huán)境內(nèi)，則，則顯示顯示局地維數(shù)環(huán)境下的局地維數(shù)環(huán)境下的變量值變量值。n如果這個(gè)維數(shù)在定義變量時(shí)是如果這個(gè)維數(shù)在定義變量時(shí)是固定的固定