- 积分
- 606
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-4-17
- 最后登录
- 1970-1-1

|
发表于 2016-12-19 20:34:39
|
显示全部楼层
为什么垂直剖面图画不出来呢?
求大神解答
gs脚本
'reinit'
'open E:\ac\4-22tongji\chang\H\f-decadal-eof1\WD-LEV5\COR\58-12\wd-t.ctl'
*'enable print E:\ac\4-22tongji\chang\H\f-decadal-eof1\WD-LEV5\COR\58-12\xg58-eof1u.gmf'
'set parea 1.5 10.0 1.4 7'
'set xlopts 1 6 0.3'
'set ylopts 1 6 0.3'
'set lon 0'
*'set lon 0'没有也行
'set lat 20 80'
'set lev 1000 150'
'set ylevs 1000 850 700 500 400 300 250 200 150'
*'set xyrev on' *xy轴交换
'set xlint 10'
*'set ylint 10'
'set mpdraw on'
'set poli off'
'set map 1 1 1'
'set mpdset mres'
*'set grid on 2 15'
'set grid off'
**********************
'set gxout shaded'
*'set cterp on'
'set csmooth on'
'set rgb 20 32 96 255'
'set rgb 21 32 159 255'
'set rgb 22 32 191 255'
'set rgb 23 0 207 255'
'set rgb 24 42 255 255'
'set rgb 25 85 255 255'
'set rgb 26 127 255 255'
'set rgb 27 170 255 255'
'set rgb 90 255 255 255'
'set rgb 28 255 255 84'
'set rgb 29 255 240 0'
'set rgb 30 255 191 0'
'set rgb 31 255 168 0'
'set rgb 32 255 138 0'
'set rgb 33 255 112 0'
'set rgb 34 255 77 0'
'set rgb 35 255 0 0'
'set rbcols auto'
*99%,0.336;95%,0.259——58
'set clevs -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 0.6 0.7'
'set ccols 20 21 22 23 24 25 27 90 28 30 31 32 33 34 35'
*'set clopts 1 0.1 0.2'
'd uwnd*1.5'
'run cbar_interp 1 1.2 0 5.6 0.5'
*'cbar 1.0 0'
'a=uwnd*1.5'
'myp=maskout(a,abs(a)-0.26)'
'set gxout shp'
'set shp -pt shppt'
'set shpopts -1 3 0.05'
'd myp'
'draw shp shppt'
*99%,0.336;95%,0.259
'set strsiz 0.3 0.3'
'set string 1 Br 8 0'
'set font 0'
'draw string 10. 7.1 cor(AC,H)'
'set strsiz 0.3 0.3'
'set string 1 Bl 8 0'
'set font 0'
'draw string 1.5 7.1 80-140E'
'set font 0'
'draw ylab Level'
'gxprint E:\ac\4-22tongji\chang\H\f-decadal-eof1\WD-LEV5\COR\58-12\xg58-eof1u--fade.eps white'
'print'
'disable print'
;
提示错误:附图
|
-
|