Grads：绘制风流畅

1、画风分布图和风流场

folder = 50mAvgMonthPic
i=0
res='E:\Pic\chenxia\drawPic\name.txt'
while(i<17)
    'reinit'
    i=i+1
    text=read(res)
    ln=sublin(text,2)
     filename=subwrd(ln,1)
     titlestr=subwrd(ln,2)
    'open E:\Pic\chenxia\drawPic\ctl\'filename'.ctl'
    'set grads off'
    'set grid off'
    'run E:\Pic\chenxia\drawPic\Colors_Wind.gs'
    'set lon 112.8 115.2'
    'set lat 37.5 39.5'
    'set t 1'
    'set ylopts 1 5 0.13'
    'set xlopts 1 5 0.13'
    'set xlint 0.5'
    'set ylint 0.5'
    'set gxout shaded'
    'set csmooth on'
    'd fengsu'
    'cbarn 0.8 1'
    'set gxout vector'
    'set arrscl 0.2'
    'set ccolor 1'
    'd -fengsu*sin(fengxiang*3.1415926/180)*2.5;-fengsu*cos(fengxiang*3.1415926/180)*2.5'
    'draw title 'titlestr''filename'_Avg_uv50m(m/s)'
    'set line 1 1 5'
    'draw shp E:\Pic\chenxia\shp\pingshan.shp'
    'printim E:\Pic\chenxia\drawPic\'folder'\'filename'.png png x1000 y800 white'
endwhile
'quit'
;

name.txt 为了循环加题目，格式如下：

chun 
xia 
qiu 
dong 
nian 
01 (1979_2019_month)_
02 (1979_2019_month)_
03 (1979_2019_month)_
04 (1979_2019_month)_
05 (1979_2019_month)_
06 (1979_2019_month)_
07 (1979_2019_month)_
08 (1979_2019_month)_
09 (1979_2018_month)_
10 (1979_2018_month)_
11 (1979_2018_month)_
12 (1979_2018_month)_