用Grads画等值线（三）-----等值线的生成

     上文中已经生成了.map文件，要想生成等值线，还需要生成grid.dat文件。生成grid的Fortran代码比较简单如下所示：

       program grid real gtd(26,24) open(10,file='D:/YK/grid.dat',form='binary') do i=1,26 do j=1,24 gtd(i,j)=0 end do end do write(10)((gtd(i,j),i=1,26),j=1,24) end

 

       描述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文件：

     'open D:/YK/grid.ctl 1' 'open D:/YK/rain.ctl 2' 'set lon 121.900 123.10' 'set lat 39.900 41.00' 'set t 1' 'define xx=oacres(grainval.1,srainval.2,20,10,5)' 'define'enable print D:/YK/YK.gmf' 'set ccolor rainbow' 'set gxout shaded' 'd xx' 'set stid on' 'set ccolor 8' 'set digsize 0.10' 'd srainval.2' 'printim D:/YK/YK5.gif white'

 

     最后运行Grads命令：grads -cl D:/YK/YK.gs即可得到等值线图。