关于画图时条件语句的使用以及数据问题的困惑

李文帅

'reinit''sdfopen D:\Data\ThirdGroup\nc\omega\w_2017.05.01_15.nc' 'sdfopen D:\Data\ThirdGroup\nc\wind\four_times_daily\uwnd.2017.05.01_15.nc''sdfopen D:\Data\ThirdGroup\nc\wind\four_times_daily\vwnd.2017.05.01_15.nc''set parea 0.5 10.5 0.9 8''set grid off''set grads off''set zlog on''set csmooth on'lon1=115.25lon2=120.25lat1=31.0lat2=36.0'define alfa=atan2('lat2-lat1','lon2-lon1')''set x 1''set y 1''set lev 1000 300''set t 22'lon=lon1'collect 1 free''collect 2 free'while(lon<=lon2)lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)'collect 1 gr2stn(-w.1*50,'lon','lat')''collect 2 gr2stn(u.2*cos(alfa)+v.3*sin(alfa),'lon','lat')'lon=lon+0.25endwhile'define lws = -w.1*50''define sonny = u.2*cos(alfa)+v.3*sin(alfa)''define isaac = mag(lws,sonny)'*'set lat 34 39''set lon 115.3 120.3' 'set lev 1000 280' 'set ylevs 925 850 700 600 500 400 350 300''set xlabs 32.5N,115.3E|33.2N,116.3E|33.9N,117.3E|34.6N,118.3E|35.3N,119.3E|36N,120.3E'if isaac < 100    'set gxout vector'    'd skip(coll2gr(2,-u),1,1);coll2gr(1,-u)'else'set gxout barb''d skip(coll2gr(2,-u),1,1);coll2gr(1,-u)'endif'set gxout shaded''d isaac''draw title May 6th 14:00BJT''printim D:\Research_Room\Data\Output\meteorology\practice.png white';

李文帅

发出来好像格式不整齐了，我重新发一下：
'reinit'
'sdfopen D:\Data\ThirdGroup\nc\omega\w_2017.05.01_15.nc'
'sdfopen D:\Data\ThirdGroup\nc\wind\four_times_daily\uwnd.2017.05.01_15.nc'
'sdfopen D:\Data\ThirdGroup\nc\wind\four_times_daily\vwnd.2017.05.01_15.nc'
'set parea 0.5 10.5 0.9 8'
'set grid off'
'set grads off'
'set zlog on'
'set csmooth on'
lon1=115.25
lon2=120.25
lat1=31.0
lat2=36.0
'define alfa=atan2('lat2-lat1','lon2-lon1')'
'set x 1'
'set y 1'
'set lev 1000 300'
'set t 22'
lon=lon1
'collect 1 free'
'collect 2 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(-w.1*50,'lon','lat')'
'collect 2 gr2stn(u.2*cos(alfa)+v.3*sin(alfa),'lon','lat')'
lon=lon+0.25
endwhile

'define lws = -w.1*50'
'define sonny = u.2*cos(alfa)+v.3*sin(alfa)'
'define isaac = mag(lws,sonny)'
*'set lat 34 39'
'set lon 115.3 120.3'
'set lev 1000 280'
'set ylevs 925 850 700 600 500 400 350 300'
'set xlabs 32.5N,115.3E|33.2N,116.3E|33.9N,117.3E|34.6N,118.3E|35.3N,119.3E|36N,120.3E'

if isaac < 100
    'set gxout vector'
    'd skip(coll2gr(2,-u),1,1);coll2gr(1,-u)'
else
'set gxout barb'
'd skip(coll2gr(2,-u),1,1);coll2gr(1,-u)'
endif
'set gxout shaded'
'd isaac'

'draw title May 6th 14:00BJT'
'printim D:\Research_Room\Data\Output\meteorology\practice.png white'
;

river

首先这个if isaac < 100 应该写成if(isaac < 100)
再有必须先画阴影图，不然会覆盖其他图形，你这个不知什么原因阴影图根本没画出来······
你先改了试试看吧

李文帅

river 发表于 2017-8-22 07:07
首先这个if isaac < 100 应该写成if(isaac < 100)
再有必须先画阴影图，不然会覆盖其他图形，你这个不知什 ...

好的，谢谢！！！

李文帅

river 发表于 2017-8-22 07:07
首先这个if isaac < 100 应该写成if(isaac < 100)
再有必须先画阴影图，不然会覆盖其他图形，你这个不知什 ...

好的，谢谢！！！