GrADS站點(diǎn)資料的使用

1、.GrADS站點(diǎn)資料的使用臺(tái)站型離散資料畫圖的一般原則1、把臺(tái)站資料r.dat寫成二進(jìn)制文件r.grd一般而言，臺(tái)站資料都是文本格式的，須用專門的程序?qū)懗蓭в姓咎?hào)、經(jīng)度、緯度等的二進(jìn)制文件。例如：對(duì)某一時(shí)次的降水資料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。該數(shù)據(jù)只有一個(gè)時(shí)次則相應(yīng)的程序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數(shù)據(jù)的數(shù)據(jù)說明文件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文件實(shí)際上，是在GrADS中運(yùn)行可執(zhí)行文件stnmap.exe，如何運(yùn)行對(duì)于

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

6、d文件的范圍大于r.dat的范圍時(shí)，除了相同范圍的每個(gè)點(diǎn)賦值為1外，其余外圍的點(diǎn)也都賦值為1 （5）grid.grd文件的時(shí)間一定要與r.dat時(shí)間一致！生成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的網(wǎng)格點(diǎn)上，并畫圖示范.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