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

1、第6章變量和函數(shù)是在GrADS執(zhí)行環(huán)境中，操作所涉及的數(shù)據(jù)對(duì)象可以是變量、函數(shù)、式形式。 6.1變量名，完全的變量名形式abbrev.file#(dimexpr，dimexpr，)，其中，abbrev :是賦予數(shù)據(jù)描述文件的變量名的省略形式。 file# :包含此變量的打開文檔的序號(hào)。 例如，在ga-opend 3360 model.ctl opend 3360 model.le.ctl上的示例中，model.CTL文件的序列號(hào)是1，model.le.CTL文件的序列號(hào)是2，序列號(hào)1是缺省值括號(hào)中的dimexpr :當(dāng)前維環(huán)境的局部維設(shè)置公式，該公式僅對(duì)該變量的維環(huán)境進(jìn)行局部修改，因此僅對(duì)固定

2、維有效。 設(shè)定形式為x|y|z|t|lat|lev|time=value (value為絕對(duì)維度數(shù)值)； 相對(duì)維數(shù)表達(dá)式是x|y|z|t|lon|lev|time/- offset (offset是針對(duì)當(dāng)前維數(shù)環(huán)境設(shè)定的偏差維數(shù)值)。 例如，對(duì)于z.3(lev=500) tv.1(time-12hr) rh默認(rèn)的當(dāng)前文件的變量rh q.2(t-1、lev=850 )、第2文件當(dāng)前時(shí)刻的12小時(shí)前的時(shí)刻的第1文件中的變量tv， dset model.le.datoptionslittle _ endian Cray _ 32bit _ IEEE undef-2.56 e 33title 5大型機(jī)

3、模型72.0 zdef 7級(jí)別10008500 an 19871 dy var s8PS 099 surfacepressureu 799 uwindsv 799 vwindsz 799 geo端口初始化799臨時(shí)q 599規(guī)范099 surfacete dset model.grbtitlesamplemodeldataforlats4dtutorialundef1e 20 dtypegribindexmodel.gmpxdef 72線性電腦0.0000005.00000 ydef 46線性電腦000 50700500300200 TDE f 5線性0 z1Jan 19871 dy var

4、s8PS 0，1，0，0 surfacepressurehpaua 733，100 eastwardwindm/sva 734，100 northwardwindm/s 100 geopotentialheightmta 711，100 airtemperaturekhus 751，100 specifichumiditykg/kgts 011，105，2曲面(2m ) airtemperaturekpr 059 將0 totalprecipitationrate ga-opend : model.CTL opend 3360 model.le.CTL set 1時(shí)間維數(shù)固定為第一個(gè)時(shí)間維數(shù)d

5、 z.2(lev=500) set z 3的垂直高度電平設(shè)定為第三個(gè)電平， 即，設(shè)為7000的兩個(gè)文件： model.ctl文件的序列號(hào)為1，model.le.ctl文件的序列號(hào)為2 .顯示第二個(gè)文件的500hPa等壓面上的變量z (絕對(duì)維數(shù))，第一個(gè)文件的第三層次2 .在對(duì)相同網(wǎng)格點(diǎn)上的不同變量進(jìn)行運(yùn)算時(shí)，如果一個(gè)變量在某個(gè)網(wǎng)格點(diǎn)處的值是缺測(cè)，則將該網(wǎng)格點(diǎn)的運(yùn)算結(jié)果除以成為缺測(cè)值的運(yùn)算規(guī)則0時(shí)的結(jié)果也視為缺測(cè)。 在計(jì)算多個(gè)數(shù)據(jù)文件中的數(shù)據(jù)時(shí)，兩個(gè)數(shù)據(jù)網(wǎng)格必須匹配。 也就是說，運(yùn)算對(duì)象的維數(shù)的變化范圍必須一致。 當(dāng)一個(gè)網(wǎng)格點(diǎn)的變動(dòng)維數(shù)大于其他的網(wǎng)格點(diǎn)時(shí)，對(duì)于變動(dòng)維數(shù)少的網(wǎng)格點(diǎn)，擴(kuò)展維數(shù)環(huán)境以便

6、于運(yùn)算操作。 例如，z-z (t-1 ) t (lev=200 )-t (lev=850 ) ave (z，t=1，t=5)z- .調(diào)用方式采用直接引用方式，形式例如，函數(shù)名(殘奧儀表1、殘奧儀表2、)、2 .常用函數(shù)、(1)ave函數(shù)例如，ga-open model.ctlsetlev 500 Dave (z，t=1，t=5)表示500hPa等壓面上的從第一次到第五次的變量z的平均量，setp area 1816 setgridoffsetgradsoffs e a=ave(v，t=7，t=660，12 ) daprintc，例如求出緯向的平均： ave(expr，lon=00，lon=。

7、如果對(duì)終點(diǎn)取一半的權(quán)重，則可以使用ave(expr、lon=0、lon=360、-b )。setxyrevonsetlev 500 set t1set lat 090 setx1definea=ave (v、lon=0、lon=360、-b)dasetcmark0setcolor (2) const函數(shù)格式： const (退出恒定(constant )功能：將部分網(wǎng)格點(diǎn)的值設(shè)為常數(shù)的恒定說明：將所有非缺陷網(wǎng)格點(diǎn)的expr值設(shè)為常數(shù)，將flag設(shè)為可選項(xiàng)，加上選項(xiàng)-a，此函數(shù)同時(shí)適用于網(wǎng)格點(diǎn)和站數(shù)據(jù)。 (3)aave函數(shù)格式：求出aave(expr，xdim1，xdim2，ydim1，ydim

8、2)功能：面積平均。 說明： expr:任意式xdim1: X維方向的開始維數(shù)學(xué)式xdim 23360 x維方向的結(jié)束維數(shù)學(xué)式y(tǒng)dim 13360 y維方向的開始維數(shù)學(xué)式y(tǒng)dim 23360 y維方向的結(jié)束維數(shù)學(xué)式set lev 850 sett 148 setx1set y1definea=aavv lon=100、lon=120、lat=30、lat=40) d a、設(shè)定lon 60 200設(shè)定lat-2080 d mag (u.1、v.2 )、(5)hcurl函數(shù)格式： hcurl(uexpr、vexpr )功能邊界上的渦度值設(shè)定為缺失測(cè)量。 6)hdivg函數(shù)格式： hdivg(uex

9、pr，vexpr )功能：修正水平分散度。 說明： uexpr是u風(fēng)速成分，vexpr是v風(fēng)速成分，風(fēng)速單位是m/s。 set lon 60 200 set lat -20 80 d hdivg (u，v )、(7)skip函數(shù)格式： skip(expr，skipx，skipy )功能：設(shè)定樣本的采樣密度。 說明： skipx、skipy的數(shù)值決定x方向和y方向的采樣密度(取值1時(shí)可以省略) *此函數(shù)主要用于向量場(chǎng)的稀疏化顯示。 d u.1； v.2、德斯基普(u.1、2、2 )； skip (v.2，2，2 )，(8)cdiff函數(shù)格式： cdiff(expr，dim )功能：中央差分函數(shù)說

10、明： expr表示進(jìn)行差分運(yùn)算的量，dim表示進(jìn)行差分運(yùn)算的次元的方向。 邊界網(wǎng)格的運(yùn)算值被設(shè)定為預(yù)設(shè)值。 例如，通過校正溫度平流，設(shè)定定義dtx=c diff (t，x )定義DTY=c diff (t，y )定義dx=c diff (lon，x )。 y ) *3. 1416/180 d-1 * (u * dtx )/(cos (lat *3. 1416/180 ) * dx ) v * DTY/dy )此函數(shù)主要用于改善圖形的輸出品質(zhì)，使其更美觀。 格點(diǎn)數(shù)據(jù)為一維時(shí)，9點(diǎn)平滑化劣化為3點(diǎn)平滑化。 d smth9(b1)、(11)vint函數(shù)格式： vint(psexpr，expr，top

11、 )功能：質(zhì)量加權(quán)垂直積分函數(shù)說明： psexpr是地面氣壓變量，expr是vint(psexpr，expr，top )=f/g * sopene : bsuq.ctlopene : bsp RES.ctlsetgxoutfwritesetfwritee : bsmonzuq.grdi=1while (I=660 )設(shè)置1定義常見問題=vint SETI setx 1144 sety 173 setz1d fuqi=I 1生成endwhiledisaad (2) map文件(3)創(chuàng)建網(wǎng)格文檔(創(chuàng)建后也可以繼續(xù)使用) (4)通過插值函數(shù)將站點(diǎn)數(shù)據(jù)插值到網(wǎng)格文件中，并繪制文件的時(shí)間必須與站點(diǎn)的數(shù)據(jù)

12、時(shí)間一致。 范圍： lat:15 55 N lon:70 140 E格式：1*1 grid.grd文件中的每個(gè)點(diǎn)一個(gè)，格式文件grid.grd : 分配在第2章數(shù)據(jù)處理中寫入70-73頁(yè)的網(wǎng)格數(shù)據(jù)grid.grrd。與該網(wǎng)格文件對(duì)應(yīng)的數(shù)據(jù)描述文件grid.ctl， 作為dsetd 3360 datagrid.grd undef-999.0 titlesamplegriddataxdef 71線性化70 ydef 41線性化15 z de f 1線性化的f 1線性化jul 19511移動(dòng)1 g 099 gr 如果有該網(wǎng)格文件和與其對(duì)應(yīng)的數(shù)據(jù)描述文件，則可通過GrADS插值函數(shù)將站點(diǎn)數(shù)據(jù)插值到該網(wǎng)

13、格文件中來(lái)使用。 (12) oacres函數(shù)格式： oacres(gexpr，sexpr )功能： cressman插值函數(shù)說明：站點(diǎn)數(shù)據(jù)插值中經(jīng)常使用的網(wǎng)格數(shù)據(jù)。 為了避免gexpr大于站點(diǎn)數(shù)據(jù)范圍的網(wǎng)格數(shù)據(jù)變量名稱sexpr站點(diǎn)數(shù)據(jù)變量名稱radii影響半徑，缺省值為：10、7、4、2、1插值造成的虛假結(jié)果，參照網(wǎng)格的間隔最好接近站點(diǎn)間隔。 將站點(diǎn)數(shù)據(jù)插值到作為插值函數(shù)的網(wǎng)格文件后使用。 中的組合圖層性質(zhì)變更選項(xiàng)。 與此網(wǎng)格文件相對(duì)應(yīng)的數(shù)據(jù)描述文件grid.ctl是dsetd 3360 datagrid.grd undef-999.0 titlesamplegriddataxdef 71線性文件70 ydef