站点资料的使用

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。该数据只有一个时次则相应的程序st

2、n.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',stat

3、us='unknown')open(99,file='d:commonrr.txt')do 100 i=1,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

4、).lt.0.0)thenwrite(*,*)lon(i),lat(i),r(i)endifenddo!write(*,*)i,zh(5000),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

5、 functionendvars3、用r.ctl生成r.map文件实际上，是在GrADS中运行可执行文件stnmap.exe，如何运行对于自带显示窗口的GrADS1.8版本的：进入命令窗口后，在提示符下输入“！stnmap”ga_>！stnmap在出现的信息提示处输入r.ctlEnter stn ctl filename: D:r.ctl4、生成格点文件因为GrADS只识别格点数据，所以需要把站点上的资料通过插值函数插值到某个格点文件上才能被识别。根据r.dat提供的信息来写格点文件grid.grd（二进制）要求：（1） grid.grd文件的精度要高于或等于r.dat的精度 （2） g

6、rid.grd文件的范围要大于或等于r.dat的范围 （3）grid.grd文件的每个点上均赋值1（范围与r.dat相同） （4）当grid.grd文件的范围大于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)=11

7、9.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)=1enddoenddowrite(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 prepar

8、ed for oacres functionendvars6、把r.grd插值到grid.grd的网格点上，并画图示范.gs文件如下：'reinit''open d:commongrid.ctl''OPEN d:commonr.ctl'*'enable print d:commonr.gmf''set grid off''set mpdset cnworld''define a1=oacres(g,r.2)''define b1=maskout(a1,g-0.5)''define bb1=smth9(b1)''set black -10000 0''set lon 119.8 121.4''set lat 30.7 32.1''