- 积分
- 1711
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-3-19
- 最后登录
- 1970-1-1
|
发表于 2018-6-30 20:22:13
|
显示全部楼层
本帖最后由 gingili 于 2018-6-30 20:23 编辑
我想画的是大于45m/s的风场阴影处,但是在rainbow里设置只能是set clevs,画出来不是只有大于45m/s的。
gs:'reinit'
'sdfopen E:\20160503\nc\uwnd.2016.nc'
'sdfopen E:\20160503\nc\vwnd.2016.nc'
'set mpdset cnworld'
'set poli on'
pages.491=q1
pages.492=q2
pages.493=q3
pages.494=q4
name.491=a
name.492=b
name.493=c
name.494=d
'set rgb 20 68 0 252'
'set rgb 21 0 239 252'
'set rgb 22 0 252 51'
'set rgb 23 252 222 0'
tt=491
while(tt<=494)
'run D:\opengrads\Contents\Resources\Scripts\page.gs 'pages.tt''
*'run E:\20160503\4colors.gs'
'set parea 0.5 10 1.0 8.0'
'set lat 10 70'
'set lon 80 160'
'set lev 200'
'set grid off'
'set grads off'
'set t 'tt
'define uv=mag(uwnd.1,vwnd.2)'
'set lat 30 55'
'set lon 100 140'
'set cmin 45'
'set gxout shaded'
'run D:\opengrads\Contents\Resources\Scripts\rainbow.gs'
'set strsiz 0.3 0.3'
'd uv'
'define div=hdivg(uwnd.1,vwnd.2)'
'set gxout contour'
'set csmooth on'
'set cint 0.5'
'set cthick 6'
'set ccolor 4'
'd div*1e5'
'set csmooth on'
'set gxout barb'
*'set ccolor 3'
'set cthick 5'
'd uwnd.1;vwnd.2'
tt=tt+1
endwhile
'run D:\opengrads\Contents\Resources\Scripts\plotcbar.gs D:\opengrads\Contents\Resources\Scripts\rainbow.gs 1 0.8 0 5.5 0.1'
'printim E:\20160503\nc\1\200div.png white'
;
E:\20160503\nc\1\200div.png |
-
|