用Grads画等值线（五）-------等值线边界控制

      上文中声称了带底图的等值线，但此等值线是越界的。要控制等值线边界，在Grads中要使用maskout函数。要使用此函数，需要设置边界内外不同的值。因此对grid.dat文件进行了如下修改：

      program grid real gtd(26,24) real mm(26,2) real yCol(24) real xCol(26) xCol(1)=121.9 xCol(2)=121.95 xCol(3)=122.0 xCol(4)=122.05 xCol(5)=122.1 xCol(6)=122.15 xCol(7)=122.2 xCol(8)=122.25 xCol(9)=122.3 xCol(10)=122.35 xCol(11)=122.4 xCol(12)=122.45 xCol(13)=122.5 xCol(14)=122.55 xCol(15)=122.6 xCol(16)=122.65 xCol(17)=122.7 xCol(18)=122.75 xCol(19)=122.8 xCol(20)=122.85 xCol(21)=122.9 xCol(22)=122.95 xCol(23)=123 xCol(24)=123.05 xCol(25)=123.1 xCol(26)=123.15 yCol(1)=39.9 yCol(2)=39.95 yCol(3)=40.0 yCol(4)=40.05 yCol(5)=40.1 yCol(6)=40.15 yCol(7)=40.2 yCol(8)=40.25 yCol(9)=40.3 yCol(10)=40.35 yCol(11)=40.4 yCol(12)=40.45 yCol(13)=40.5 yCol(14)=40.55 yCol(15)=40.6 yCol(16)=40.65 yCol(17)=40.7 yCol(18)=40.75 yCol(19)=40.8 yCol(20)=40.85 yCol(21)=40.9 yCol(22)=40.95 yCol(23)=50.0 yCol(24)=50.05 mm(1,1)=0 mm(1,2)=0 mm(2,1)=0 mm(2,2)=0 mm(3,1)=40.05 mm(3,2)=40.18 mm(4,1)=40.06 mm(4,2)=40.20 mm(5,1)=40.05 mm(5,2)=40.31 mm(6,1)=40.04 mm(6,2)=40.33 mm(7,1)=40.01 mm(7,2)=40.41 mm(8,1)=39.96 mm(8,2)=40.87 mm(9,1)=39.98 mm(9,2)=40.88 mm(10,1)=39.97 mm(10,2)=40.91 mm(11,1)=39.96 mm(11,2)=40.93 mm(12,1)=39.95 mm(12,2)=40.94 mm(13,1)=39.94 mm(13,2)=40.93 mm(14,1)=40.06 mm(14,2)=40.76 mm(15,1)=40.06 mm(15,2)=40.59 mm(16,1)=40.08 mm(16,2)=40.60 mm(17,1)=40.08 mm(17,2)=40.60 mm(18,1)=40.11 mm(18,2)=40.58 mm(19,1)=40.11 mm(19,2)=40.56 mm(20,1)=41.15 mm(20,2)=40.53 mm(21,1)=40.33 mm(21,2)=40.53 mm(22,1)=40.42 mm(22,2)=40.51 mm(23,1)=0 mm(23,2)=0 mm(24,1)=0 mm(24,2)=0 mm(25,1)=0 mm(25,2)=0 mm(26,1)=0 mm(26,2)=0 ! redefine the grids numbers when necessary open(10,file='D:/YK/grid.dat',form='binary') do i=1,26 do j=1,24 if(yCol(j)<=mm(i,1) .or. yCol(j)>=mm(i,2)) then gtd(i,j)=1 print *,gtd(i,j) else gtd(i,j)=3 print *,gtd(i,j) endif end do end do write(10)((gtd(i,j),i=1,26),j=1,24) end

 

       最后生成的gs文件如下所示：

         'open D:/YK/grid.ctl 1' 'open D:/YK/rain.ctl 2' 'set mpdset D:/Grads/dat/cn3000' '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 yy=maskout(xx,grainval.1-2)' 'enable print D:/YK/YK.gmf' 'set ccolor rainbow' 'set gxout shaded' 'd yy' 'set stid on' 'set ccolor 8' 'set digsize 0.10' 'd srainval.2' 'run D:/YK/YKBoard.gs' 'printim D:/YK/YK5.gif white'

 

        

 

 

           但是这种控制等值线的方式并不可取，首先是人为控制太复杂，再者效果也不理想。好像MeteoInfo能生成用于控制等值线边界的

    maskout所需文件，但我还没有试过。

          还有很多地方要弄，因为项目原因了，所以要告别Grads一段时间，希望以后有时间能继续研究。