1、高拴柱高拴柱國家氣象中心國家氣象中心2008年年12月月7日日Grid Analysis and Display Systemreal ps(72,46,5), z(72,46,7,5),u(72,46,7,5),v(72,46,7,5),q(72,46,5,5).read(*)ps,z,u,v,q.open (8,file=example.dat,form=unformatted,access=direct,recl=72*46*4)irec = 0do it=1,5 irec=irec+1 write(8,rec=irec) (ps(j,i,it),j=1,72),i=1,46) do k

2、=1,7 irec=irec+1 write(8,rec=irec+1) (z(j,i,k,it),j=1,72),i=1,46) enddo .u,v,q.enddoopen (9,file=example.ctl)write(9, (a) dset example.datwrite(9, (a) title 5 days of model outputwrite(9, (a) undef -2.56e33write(9, (a) xdef 72 linear 100.0 1.0write(9, (a) ydef 46 linear 10.0 1.0write(9, (a) zdef 7 l

3、evels 1000,850,700,500,300,200,100write(9, (a) tdef 5 linear 02jan2008 1dywrite(9, (a) vars 5write(9, (a) ps 0 99 Surface Pressurewrite(9, (a) z 7 99 Geopotential Heightswrite(9, (a) u 7 99 U windswrite(9, (a) v 7 99 V windswrite(9, (a) q 5 99 Specific Humiditywrite(9, (a) endvarsclose(9)data zlev/1

4、000,925,850,700,500,400,300,250,200,150,100/tim=0.0；ilev=11；nflag=0open(2,file=filename.dat,form=unformatted,access=direct,recl=73*4)irecord=0do j=1,irecirecord=irecord+1write(2,rec=irecord)abc/sta(j)(1:5),lat(j),lon(j),tim,ilev,nflag,+ (zlev(i),z(j,i),t(j,i),td(j,i),d(j,i),f(j,i),i=1,11)enddoilev=0

5、irecord=irecord+1write(2,rec=irecord)abcdefgh,lat(1),lon(1),tim,ilev,nflag,+ (zlev(i),i=1,11),(z(1,i),i=1,11),(t(1,i),i=1,11),+ (td(1,i),i=1,11),(d(1,i),i=1,11),(f(1,i),i=1,11)close(2) open(3,file=filename.ctl,form=formatted,status=unknown) write(3,(a)dset /filename.dat write(3,(a)dtype station writ

6、e(3,(a)stnmap sounding.map write(3,(a)undef -999.0 write(3,(a)title sounding write(3,(a)tdef 1 linear 02JAN1987 1DY write(3,(a)vars 5 write(3,(a)zz 11 99 /staa/ height1 write(3,(a)tt 11 99 /staa/ temper1 write(3,(a)t_td 11 99 /staa/ t_td1 write(3,(a)uu 11 99 /staa/ u wind1 write(3,(a)vv 11 99 /staa/

7、 v wind1 write(3,(a)endvarsclose(3)ga打開打開ctl文件文件命令：命令：open filename1.ctlga- open strike024_041714Scanning description file: strike024_041714Data file strike024_041714.dat is open as file 1LON set to 100 150LAT set to 0 40LEV set to 1000 1000Time values set: 2005:4:17:0 2005:4:17:0 命令：命令：open filenam

8、e2.ctlga- open strike024_041720Scanning description file: strike024_041720Data file strike024_041720.dat is open as file 2ga- close 1Close Error: Only last file may be closedga- close 2File 2 has been closed查詢命令查詢命令命令：命令：query filesga- q filesFile 1 : TC0801-StrkPro024 2008041714 Descriptor: strike0

9、24_041714.ctl Binary: strike024_041714.datFile 2 : TC0801-StrkPro024 2008041720 Descriptor: strike024_041720.ctl Binary: strike024_041720.dat命令：命令：query file nga- q file 2File 2 : TC0801-StrkPro024 2008041720 Descriptor: strike024_041720.ctl Binary: strike024_041720.dat Type = Gridded Xsize = 101 Ys

10、ize = 81 Zsize = 2 Tsize = 1 Number of Variables = 2 strp 0 0 strike probability strn 0 0 strike number查詢命令查詢命令命令：命令：query ctlinfoga- q ctlinfodset strike024_041714.dattitle TC0801-StrkPro024 2008041714undef 9999xdef 101 linear 100 0.5ydef 81 linear 0 0.5zdef 2 levels 1000 500tdef 1 linear 00Z19JUL2

11、005 60mnvars 2strp 0 0 strike probabilitystrn 0 0 strike number命令：命令：query dimsga- q dimsDefault file number is: 1X is varying Lon = 100 to 150 X = 1 to 101Y is varying Lat = 0 to 40 Y = 1 to 81Z is fixed Lev = 1000 Z = 1T is fixed Time = 00Z19JUL2005 T = 11.等值線等值線 命令：命令： display expressga- d strpCo

12、ntouring: 0 to 100 interval 10標量場分析標量場分析 1. 等值線等值線 設置設置 set gxout contour （缺?。ㄈ笔。?2. 陰影圖陰影圖 設置設置 set gxout shaded設置：設置： set * #維數參數設置維數參數設置Set lon 100 120 Set lat 10 25Set lev 1000Set t 1等值線設置等值線設置set cint 15 指定等值線間隔指定等值線間隔set ccolor 3 指定等值線顏色指定等值線顏色 0：黑； 1：白；2：紅；3：綠；4：藍；5：；6：；7：黃；8：桔黃； 9：紫；10：黃綠；1

13、1：中藍；12：深藍；13：；14：深紫；15：灰色set cstyle 2 指定等值線線形指定等值線線形 1：實線；2：長斷線；3：短斷線；4：一長一短斷線； 5：點線；6：點劃線；7：點點劃線 set cthick 10 設定線粗細設定線粗細(1-10)set cmin 10 不畫低于不畫低于10（含）以下的等值線（含）以下的等值線 set cmax 90 不畫超過不畫超過90（不含）以上的等值線（不含）以上的等值線 set clevs 495 523 534 556 560 564 572 584 588set ccols 1 2 3 4 5 當指定顏色少于等值線條數時，最后一種顏色適用

14、于超出的所有等值線。 等值線標注設置等值線標注設置set clab %.1f 等值線標記方式等值線標記方式 set clab on（等值線標數值）/off（不標記）/forded（強制標記）/%gK（在標記后加上字符K）/%g%（加上字符%）/%.2f(保留2位小數)/ %03.1f（整數部分保留3為，若不足3位，前部用0補齊）set clskip 2 每隔一條等值線標記數值每隔一條等值線標記數值set clopts 3 0.1 0.2 等值線標記的顏色等值線標記的顏色 粗細粗細 1.陰影圖陰影圖 命令：命令： display expressga- cga- set gxout shadedg

15、a- d strpContouring: 0 to 100 interval 10ga- cbarn 1 0ga- cbarn 1 1標量場分析標量場分析 1. 等值線等值線 設置設置 set gxout contour （缺?。ㄈ笔。?2. 陰影圖陰影圖 設置設置 set gxout shaded矢量場分析矢量場分析display u;v set gxout barb 風標表示的風矢量風標表示的風矢量 set gxout vector 以箭頭方式表示矢量場以箭頭方式表示矢量場 set gxout stream 以流線方式表示矢量場以流線方式表示矢量場 set gxout vector 以箭

16、頭方式表示矢量場以箭頭方式表示矢量場set arrlab on|off 設置是否標注矢量圖題注（箭頭）設置是否標注矢量圖題注（箭頭）set arrowhead size 矢量圖題注箭頭尺寸（英寸）矢量圖題注箭頭尺寸（英寸） set arrscl 把把值矢量用長度為值矢量用長度為的箭頭表示的箭頭表示set gxout barb 風標表示的風矢量風標表示的風矢量也可以設置風向桿的顏色、粗細、大小也可以設置風向桿的顏色、粗細、大小set gxout stream 以流線方式表示矢量場以流線方式表示矢量場set strmden value 設置流線密度設置流線密度1-10,缺省為缺省為5draw title 2008.09.26.08 TC0815 寫圖標題寫圖標題 set string 1 c 2 45 設置字串的顏色、位置、粗細、角度設置字串的顏色、位置、粗細、角度set strsiz 0.3 0.6 字水平大小、垂直大小字水平大小、垂直大小draw string 8 3 first group 在（在（8,3）位置寫字符串）位置寫字符串firs