资源描述：
《GrADS站点资料的使用.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、GrADS站点资料的使用台站型离散资料画图的一般原则1、把台站资料r.dat写成二进制文件r.grd一般而言，台站资料都是文本格式的，须用专门的程序写成带有站号、经度、纬度等的二进制文件。例如：对某一时次的降水资料r.dat有如下形式：LlonLlatPrecipitation119.830.60119.830.6160119.830.6320119.830.6480119.830.6640119.830.680119.830.6960119.830.7120119.830.7280119.830.7440119.830.760119.830.7760。。。。。。该数据
2、只有一个时次则相应的程序stn.f如下：parameter(n=10201)integerr1(n)reallat(n),lon(n),r(n)character*8zh(n)open(1,file='d:common1.txt',form='formatted',status='old')doi=1,nread(1,*)lon(i),lat(i),r1(i)r(i)=r1(i)enddodoi=1,nzh(i)=char(i)enddo!!!!!!!!!!!!!!!!!!!!!!!!!!!!!open(9,file='d:commonr.grd',form='
3、binary',status='unknown')open(99,file='d:commonrr.txt')do100i=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)100continuenlev=0write(9)zh(i-1),lat(i-1),lon(i-1),tim,nlev,nflagclose(9)!doi=1,160doi=1,ni
4、f(r(i).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.ctldsetd:commonr.grddtypestationstnmapd:commonr.mapundef-999.0titlestationdatatdef1linearjan19981movars1r099griddatapreparedforoac
5、resfunctionendvars3、用r.ctl生成r.map文件实际上，是在GrADS中运行可执行文件stnmap.exe，如何运行对于自带显示窗口的GrADS1.8版本的：进入命令窗口后，在提示符下输入“！stnmap”ga_>！stnmap在出现的信息提示处输入r.ctlEnterstnctlfilename:D:r.ctl4、生成格点文件因为GrADS只识别格点数据，所以需要把站点上的资料通过插值函数插值到某个格点文件上才能被识别。根据r.dat提供的信息来写格点文件grid.grd（二进制）要求：（1）grid.grd文件的精度要高于或等于r.dat的精
6、度（2）grid.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)reallat(n),lon(n)reals(n,n)open(1,file='d:commongrid.grd',form='binary')lat(1)=30.0lo
7、n(1)=119.0doi=1,n-1lat(i+1)=lat(i)+0.1lon(i+1)=lon(i)+0.1enddodoi=1,ndoj=1,ns(i,j)=1enddoenddowrite(1)send5、编写grid.grd的.ctl文件：grid.ctldsetd:commongrid1.grdundef-999.0titleSampleGRIBDataxdef71linear701ydef41linear151zdef1linear10001tdef1linearjan19981movars1g099griddatap
