在p坐标系下画剖面图会有如下几条红色语句语句:
lon=lon1
'collect 1 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(U*cos(alfa)+V*sin(alfa),'lon','lat')'
lon=lon+0.02
endwhile
*'set lon 'lon1' 'lon2''
'set lon 109.2 109.8'
*'set lon 'lon1' 'lon2''
'set lev 0.500000 10.00000'
*'set ylab on'
'set gxout shaded'
*'set clevs -200 -150 -100 -50 0 50 100 150 200'
'set clevs -120 -90 -60 -30 0 30 60 90 120'
'f:\vortexline\output.gs'
*'set clevs 3 6 9 12 15 18 21 24 27 30'
*'F:\output20160721(1).gs'
*'set clevs 0 5 10 15 20 25 30 35 40 45 50'
*'F:\output20150121.gs'
*'d coll2gr(4,u)*100'
'd coll2gr(1,-u)*100'
但如果改成Z坐标系下,我把上面的-u改成了u,即‘d coll2gr(1,u)’但是会出现上面的warning,经过摸索,发现这里的2nd arg指的是u,因此把u去掉警告会消失。其实,警告消失后出来的图跟带警告的图一样没变化,即有警告并不是说出错。单用
‘d coll2gr(1)’更舒服点,毕竟有警告让人不爽。