- 积分
- 2473
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-3-27
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 东风急流 于 2016-10-22 10:18 编辑
参照清风大神的一页多图统一色标轴,画了张图,但是色标上数字有点小问题,最后一个数字后面有个单引号。。。感觉不那么完美。
这个玩意是从哪里出来的呢
我的脚本如下,请大家指教
'reinit'
'open d:\tu3\kz6.ctl' ;*利用此前的ARWpost输出文件读取地形数据
'open d:\tu3\cp6.ctl' ;*利用此前的ARWpost输出文件读取地形数据
*颜色设置(用清风调色板克隆)
'set rgb 20 160 0 200'
'set rgb 21 130 0 220'
'set rgb 22 30 60 255'
'set rgb 23 0 160 255'
'set rgb 24 0 200 200'
'set rgb 25 0 210 140'
'set rgb 26 0 220 0'
'set rgb 27 160 230 50'
'set rgb 28 230 220 50'
'set rgb 29 230 175 45'
'set rgb 30 240 130 40'
'set rgb 31 250 60 60'
'set rgb 32 240 0 130'
'run D:\tu3\lujin\13colors.gs';*为后续统一色标轴做准备
*画上图
'set vpage 0.5 5.5 1 8'
'set dfile 1'
'set grads off' ;*关闭grads角标 以便编辑个性的下标
'set grid off'
'set mproj latlon'
'set poli on' ;*画出国界省界线
'set mpdset hires' ;*设置高分辨率地图
'set xlint 5' ;*定义横坐标的标记间隔
'set ylint 5' ;*定义纵坐标的标记间隔
'set xlopts 1 4 0.15' ;*横坐标 颜色 线宽 大小
'set ylopts 1 4 0.15' ;*纵坐标 颜色 线宽 大小
'set gxout shaded'
'set lon 110 125'
'set lat 15 25'
'set t 7'
'd SST-273.15'
*****************************************************************
*控制试验路径
i=1
while(i<=13) ;*文件内容n行之内
read_file=read('d:/tu3/kz6.txt')
read_code=sublin(read_file,1) ;*读取文件打开代码
if(read_code>0);
break;
endif ;*读取正常为0,其他要么出错,要么结束
xbefore=x0
ybefore=y0
read_line=sublin(read_file,2) ;*读取文件行内容
lon0=subwrd(read_line,1) ;*取得第1部分 经度
lat0=subwrd(read_line,2) ;*取得第2部分 纬度
date0=subwrd(read_line,3) ;*取得第3部分 时间
'q w2xy ' lon0 ' ' lat0 ;*将经纬度坐标转换为GRADS坐标
x0=subwrd(result,3) ;*取得转换后的X坐标
y0=subwrd(result,6) ;*取得转换后的Y坐标
'draw mark 5 'x0' 'y0' 0.1' ;*作出各个坐标点 24
if(i>1) ;*第一个点作为起始点
'draw line 'xbefore' 'ybefore' 'x0' 'y0'' ;
endif;
i=i+1
Endwhile
***************************************************************
*画下图
'set vpage 5.5 10.5 1 8'
*'set parea 6 11 1 8'
'set dfile 2'
'set grads off' ;*关闭grads角标 以便编辑个性的下标
'set grid off'
'set mproj latlon'
'set poli on' ;*画出国界省界线
'set mpdset hires' ;*设置高分辨率地图
'set xlint 5' ;*定义横坐标的标记间隔
'set ylint 5' ;*定义纵坐标的标记间隔
'set xlopts 1 4 0.15' ;*横坐标 颜色 线宽 大小
'set ylopts 1 4 0.15' ;*纵坐标 颜色 线宽 大小
'set gxout shaded'
'set lon 110 125'
'set lat 15 25'
'set t 7'
'd SST.2-273.15' ;
*****************************************************************
*耦合试验路径
i=1
while(i<=13) ;*文件内容n行之内
read_file=read('d:/tu3/cp6.txt')
read_code=sublin(read_file,1) ;*读取文件打开代码
if(read_code>0);
break;
endif ;*读取正常为0,其他要么出错,要么结束
xbefore=x0
ybefore=y0
read_line=sublin(read_file,2) ;*读取文件行内容
lon0=subwrd(read_line,1) ;*取得第1部分 经度
lat0=subwrd(read_line,2) ;*取得第2部分 纬度
date0=subwrd(read_line,3) ;*取得第3部分 时间
'q w2xy ' lon0 ' ' lat0 ;*将经纬度坐标转换为GRADS坐标
x0=subwrd(result,3) ;*取得转换后的X坐标
y0=subwrd(result,6) ;*取得转换后的Y坐标
'draw mark 5 'x0' 'y0' 0.1' ;*作出各个坐标点 24
if(i>1) ;*第一个点作为起始点
'draw line 'xbefore' 'ybefore' 'x0' 'y0'' ;
endif;
i=i+1
Endwhile
***************************************************************
*'xcbar 9.5 9.7 1 1.2 -fw 0.1 -fh 0.15 -ft 3 -fs 1 -fo 0 -line off'
'run plotcbar D:\tu3\lujin\13colors.gs 1 0.9 0 5.5 1.8'
'printim d:/tu3/lujin/sst_lujin.png white' ;*图形输出
;
|
|