上文中已经生成了.map文件,要想生成等值线,还需要生成grid.dat文件。生成grid的Fortran代码比较简单如下所示:
描述Grid的二进制文件为:
program grid
real gtd(100,100)
open(10,file='D:/YK/grid.dat',form='binary')
do i=1,100
do j=1,100
gtd(i,j)=0
end do
end do
write(10)((gtd(i,j),i=1,100),j=1,100)
end
要生成等值线,还需要输入Grad命令,或写成如下的gs文件:
最后运行Grads命令:grads -cl D:/YK/YK.gs即可得到等值线图。