- 积分
- 1236
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-11-18
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
程序目的是想在显示地形情况下,再将风速填色显示。单独画可以实现。合在一块,地形就看不到了。请问这个要怎么做。
'reinit'
'open d:/dz.ctl'
'sdfopen d:/vwnd.mon.mean.nc'
'sdfopen d:/omega.mon.mean.nc'
'set grads off'
'set grid off'
'set lon 100'
'set lat 15 55'
'set lev 1000 100'
'set csmooth on'
'set xlopts 1 4 0.14'
'set ylopts 1 4 0.14'
'set dfile 1'
'define zz=dz'
'set gxout shaded'
'set cint 8000'
'set rgb 18 125 125 125'
'set csmooth on'
'set rbcols 1'
'set cmax 0'
'set xlpos -10'
'set ylpos -10'
'd smth9(zz)'
'set dfile 2'
'set dfile 3'
'set lat 35 55'
'set lon 80'
'set lev 1000 100'
'set gxout shaded'
'define uu=ave((ave(vwnd.2,t=51,t=708,12)+ave(vwnd.2,t=52,t=708,12)+ave(vwnd.2,t=53,t=708,12))/3,lon=80,lon=120)'
'define ww=-100*(ave((ave(omega.3,t=51,t=710,12)+ave(omega.3,t=52,t=710,12)+ave(omega.3,t=53,t=710,12))/3,lon=80,lon=120))'
'd mag(uu,ww)'
'set gxout stream'
'set ccolor 1'
'set cthick 6'
'set xlpos -10'
'set ylpos -10'
'd uu;ww'
'run e:/axis.gs -type L -position o -start 1000 -end 10 -interval 100'
'run e:/axis.gs -type b -position o -start 80 -end 120 -interval 20 -suffix `3.'
|
|