GrADS站點資料的使用

1、.GrADS站點資料的使用臺站型離散資料畫圖的一般原則1、把臺站資料r.dat寫成二進制文件r.grd一般而言，臺站資料都是文本格式的，須用專門的程序寫成帶有站號、經度、緯度等的二進制文件。例如：對某一時次的降水資料r.dat有如下形式：L lonL lat Precipitation119.830.60119.830.6160119.830.6320119.830.6480119.830.6640119.830.680119.830.6960119.830.7120119.830.7280119.830.7440119.830.760119.830.7760。該數據只有一個時次則相應的程序s

2、tn.f如下：parameter(n=10201) integer r1(n)real lat(n),lon(n),r(n)character*8 zh(n)open(1,file=d:common1.txt,form=formatted,status=old)do i=1,nread(1,*)lon(i),lat(i),r1(i)r(i)=r1(i)enddodo i=1,nzh(i)=char(i)enddo!open(9,file=d:commonr.grd,form=binary,status=unknown)open(99,file=d:commonrr.txt)do 100 i=1

3、,ntim=0.0nlev=1nflag=1write(9)zh(i),lat(i),lon(i),tim,nlev,nflagwrite(9)r(i)write(99,200)lon(i),lat(i),r(i)200format(f8.3,f8.3,f6.2)100 continuenlev=0write(9)zh(i-1),lat(i-1),lon(i-1),tim,nlev,nflagclose(9)!do i=1,160do i=1,nif(r(i).lt.0.0)thenwrite(*,*)lon(i),lat(i),r(i)endifenddo!write(*,*)i,zh(50

4、00),lat(5000),lon(5000),tim,nlev,nflag,r(5000)pause!enddoend2、編寫r.grd數據的數據說明文件r.ctldset d:commonr.grddtype stationstnmap d:commonr.mapundef -999.0title station datatdef 1 linear jan1998 1movars 1r 0 99 grid data prepared for oacres functionendvars3、用r.ctl生成r.map文件實際上，是在GrADS中運行可執(zhí)行文件stnmap.exe，如何運行對于

5、自帶顯示窗口的GrADS1.8版本的：進入命令窗口后，在提示符下輸入“！stnmap”ga_！stnmap在出現的信息提示處輸入r.ctlEnter stn ctl filename: D:r.ctl4、生成格點文件因為GrADS只識別格點數據，所以需要把站點上的資料通過插值函數插值到某個格點文件上才能被識別。根據r.dat提供的信息來寫格點文件grid.grd（二進制）要求：（1） grid.grd文件的精度要高于或等于r.dat的精度 （2） grid.grd文件的范圍要大于或等于r.dat的范圍 （3）grid.grd文件的每個點上均賦值1（范圍與r.dat相同） （4）當grid.gr

6、d文件的范圍大于r.dat的范圍時，除了相同范圍的每個點賦值為1外，其余外圍的點也都賦值為1 （5）grid.grd文件的時間一定要與r.dat時間一致！生成grid.grd文件的程序grid.for!1parameter(n=31)real lat(n),lon(n)real s(n,n)open(1,file=d:commongrid.grd,form=binary)lat(1)=30.0lon(1)=119.0do i=1,n-1lat(i+1)=lat(i)+0.1lon(i+1)=lon(i)+0.1enddo do i=1,ndo j=1,n s(i,j)=1enddoenddow

7、rite(1)send5、編寫grid.grd的.ctl文件：grid.ctldset d:commongrid1.grdundef -999.0title Sample GRIB Dataxdef 71 linear 70 1ydef 41 linear 15 1zdef 1 linear 1000 1tdef 1 linear jan1998 1movars 1g 0 99 grid data prepared for oacres functionendvars6、把r.grd插值到grid.grd的網格點上，并畫圖示范.gs文件如下：reinitopen d:commongrid.ctlOPEN d:commonr.ctl*enable print d:commonr.gmfset grid offset mpdset cnworlddefine a1=oacres(g,r.2)define b1=maskout(a1,g-0.5)define bb1=smth9(b1)set black -10