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)_