利用Grads画站点图（contour、shaded、grfill等）

1. 利用Fortran程序将数据输出为grd格式。参考程序如下：

parameter(num=160)   （站点数）

character*8 sta(num)    (站名名数组，可任意）

dimension xlon(num),ylat(num),rc(num)   (经、纬度数组及其上对应的数值)

open(30,file='cor.grd',form='unformatted')  (工作站)

open(30,file='cor.grd',form='binary')       (微机)

do 100 lev=1,20        (共输出20个时间上的观测或结果)

tim=0.0

ilev=1

nflag=1

do 10 i=1,num

sta(i)=char(i)

write(30) sta(i),ylat(i),xlon(i),tim,ilev,nflag

write(30) rc(i)

10 continue

   ilev=0

write(30)sta(num),ylat(num),xlon(num),tim,ilev,nflag  (每个时次的结束)

100 continue

close(30)

2. ctl文件 (创建与步骤1中输出的站点数据配对的station.ctl文件)： 

dset  cor.grd

format  sequential

dtype  station

stnmap  cor.map

undef  -9.99e33

tdef  20  linear jun1958 1mo

vars  1

r  0  99  correlation

endvars

3. 为插值函数准备格点数据grid.grd及对应的grid.ctl

dset  grid.grd

undef  -99.0

xdef  71  linear  70   1

ydef  41  linear  15   1

zdef  1   linear  1    1

tdef  20  linear  jun1958  1mo

vars  1

g   0  99  grid data prepared for oacres function

endvars

注意：两个ctl（station.ctl和grid.ctl）文件中的时间要严格一致。

4. 利用命令！stnmap -i  station.ctl创建站点映射文件（如cor.map）运行画图文件gs，利用插值函数oacres实现插值运算，作图。