- 积分
- 64
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-1-9
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2015-4-19 14:05:09
|
显示全部楼层
画水平场,已知dz.ctl和grd文件,怎样添加地形?这是我的gs文件,帮忙看一下怎么加,谢谢了!'reinit'
'clear'
'set grads off'
'open d:/data/fnl_20130711_06_00.ctl'
'open d:/data/dz.ctl'
'set mpdset cnworld'
*'set mpdset cnscs'
'set map 0 1 1'
'set mproj scaled'
'set lat 20 70'
'set lon 40 145'
'set lev 500'
'set xlopts 1 4 0.2'
'set ylopts 1 4 0.2'
'set ylint 10'
'set xlint 10'
'set xlopts 1 6'
'set ylopts 1 6'
'set grid on'
'set map 3 3 2'
'set csmooth on'
'set cthick 5'
'enable print d:/data/201307090000500.gif'
'set parea 2. 10.0 1. 8.0'
'set grid on'
'set mpdraw on'
'set dignum 1'
'set digsize 0.15'
'set cint 20'
*'set clevs 10 20 30 40'
*'set ccols 7 12 8 2 6'
'set gxout contour'
'set ccolor 2'
'd HGTprs'
'q w2xy 104.06 30.67'
x=subwrd(result,3)
y=subwrd(result,6)
'draw string 'x' 'y' Chengdu'
'q w2xy 104.06 30.67'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3 'x' 'y' 0.15'
'set gxout barb'
'set digsize 0.05'
'd skip(UGRDprs*2.5,4);skip(VGRDprs*2.5,4)'
'draw title 2013:07:11:06 500hPa height [gpm] & wind [m/s]'
'draw xlab Longitude(deg)'
'draw ylab Latitude(deg)'
'print'
'disable print' |
|