画多种形状的点以及台风路径的gs

‘reinit’
　　’open E:/eddy-statistic/ssh/test-ssh-2008.ctl’
　　t1=’29’
　　’set t ‘t1
　　’q time’
　　utc = subwrd(result,3)
　　’enable print E:/eddy-statistic/PRC/2008/2008-ssh-argo-weekly-‘utc’mark.gmf’
　　xsta=’120′
　　xend=’160 ‘
　　ysta=’15’
　　yend=’30 ‘
　　xsta1=’0′
　　xend1=’360 ‘
　　ysta1=’-90 ‘
　　yend1=’90’
　　’set x 1′
　　’set y 1′
　　’set xlopts 1 6 0.12′
　　’set ylopts 1 6 0.12′
　　’set grads off’
　　’set grid off’
　　
　　’set xlevs 121 123 125 127 129 131 133 135 137 139 141 143 145 147 149 151 153 155 157 159′
　　’set ylevs 11 13 15 17 19 21 23 25 27 29′
　　’set xlab on’
　　’set lon ‘xsta1’ ‘xend1
　　’set lat ‘ysta1’ ‘yend1
　　’set gxout shaded’
　　’define f=2.0*sin(lat*3.1416/180)*7.292e-5′
　　’define dhx = cdiff(ssh,x)/(cdiff(lon,x)*3.1416/180*cos(lat*3.1416/180))’
　　’define dhy = cdiff(ssh,y) /(cdiff(lat,y)*3.1416/180)’
　　’define wo=1e6*hdivg(dhx*9.8/100.0/f,dhy*9.8/100.0/f)/6.37e6′
　　’set clevs -7 -5 -3 -1 0 1 3 5 7 ‘
　　’rgbsetnew’
　　’set ccols 77 76 74 72 0 0 64 65 66 67 68 ‘
　　’set ccols 77 76 74 72 71 62 64 65 66 67 68 ‘
　　’set lon ‘xsta’ ‘xend
　　’set lat ‘ysta’ ‘yend
　　’d smth9(wo) ‘
　　’set gxout contour’
　　
　　’set cthick 4′
　　’set clevs -30 -15 -9 -6 -3 0 3 6 9 15 30 ‘
　　’set ccols 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1′
　　’d ssh’
　　’basemap l 1 1′
　　
　　cnt = 100
　　lll =1000
　　i = 1
　　*lat=’8.4000 8.0000 8.1000 8.0000 7.9000 7.8000 8.5000 9.1000 9.9000 10.5000 11.4000 11.8000 12.2000 12.5000 13.0000 13.5000 13.9000 14.1000 14.5000 14.8000 15.2000 15.5000 16.1000 16.8000 17.6000 18.4000 19.2000 20.0000 20.8000 21.8000 23.0000 24.4000 26.1000 27.1000 29.3000 ‘
　　*lon=’149.5000 148.2000 146.6000 145.7000 144.5000 144.0000 143.3000 142.4000 141.5000 140.4000 138.8000 137.8000 137.1000 136.7000 136.4000 135.9000 135.4000 134.9000 134.4000 133.8000 133.2000 132.9000 132.4000 131.8000 131.3000 131.0000 130.9000 131.0000 131.8000 132.6000 133.8000 135.5000 137.6000 139.1000 140.3000 ‘
　　*************************************draw track*************************************
　　col = -1
　　’set line 1 5′
　　’set line 1 1 6′
　　while (i<=cnt)
　　’q w2xy ‘%subwrd(lon,i)%’ ‘%subwrd(lat,i)
　　x = subwrd(result,3)
　　y = subwrd(result,6)
　　if (i>1)
　　’draw line ‘xold’ ‘yold’ ‘x’ ‘y
　　endif
　　xold = x
　　yold = y
　　i = i + 1
　　endwhile
　　i = 1
　　while (i　　 if (i>lll); col=2; endif;
　　 ‘q w2xy ‘%subwrd(lon,i)%’ ‘%subwrd(lat,i)
　　 x = subwrd(result,3)
　　 y = subwrd(result,6)
　　 if (i　　’set cmark 1 1′
　　’draw mark 5 ‘x’ ‘y’ 0.1′
　　’draw string ‘x’ ‘y’ ‘i
　　endif
　　 i = i + 1
　　endwhile
　　
　　************************************************another type of mark*******************************************
　　
　　lat2=’17.3890 17.4450 17.3860 16.6230 19.1260 17.5620 16.6090 16.4540 ‘
　　lon2=’124.6540 124.4540 130.4050 129.1810 125.4050 131.8590 132.9470 140.4480 ‘
　　
　　i = 1
　　while (i　　 if (i>lll); col=2; endif;
　　 ‘q w2xy ‘%subwrd(lon2,i)%’ ‘%subwrd(lat2,i)
　　 x = subwrd(result,3)
　　 y = subwrd(result,6)
　　 if (i　　’set cmark 1 1′
　　’draw mark 9 ‘x’ ‘y’ 0.1′
　　*’draw string ‘x’ ‘y’ ‘i
　　endif
　　 i = i + 1
　　endwhile
　　’cbarn 1 0 4.9 0.7′
　　’draw title ‘utc
　　’print’
　　’disable print’
　　’reinit’

转载自：https://blog.csdn.net/chht/article/details/5256020