色标问题（已解决）

风儿飘飘

'reinit'
'open F:/REGanalysis/Regressdata/OLR/regolr.ctl'
'open F:/REGanalysis/Regressdata/OLR/Ffc.ctl'
'enable print F:/REGanalysis/Regressdata/OLR/REGolr.gmf'
'set z 1'
'set t 1'
'set parea 0.75 10. 0.5 7.5'
'set grads off'
'set grid off'
'set mpdset mres'
*'set mpdset cnworld'
'set gxout shaded'
'set csmooth on'
'set annot 1 8'
'set xlopts 1 5 0.15 '
'set ylopts 1 5 0.15 '
'set xlint 10'
'set ylint 10'
'd regxs'
'cbarn 1 1 10.25 4. '
lon1=60
lat1=0
while(lon1<=160&lat1<=60)
'set lon 'lon1''
'set lat 'lat1''
if(f.2>=4.10)
'q w2xy 'lon1' 'lat1''
x1=subwrd(result,3)
y1=subwrd(result,6)
*say 'x1='x1
*say 'y1='y1
'set cmax -100'
'd f.2'
'draw mark 3 'x1' 'y1' 3'
else
break
endif
lon1=lon1+1
lat1=lat1+1
endwhile
*'set xlpos -20'   
*'set ylpos -20'
*'set black -0.3 0.3'
*'set rbcols 4'
*'set ccolor 15'
'draw title \\Regression of the AOD left coefficient onto OLR '
'printim F:/REGanalysis/Regressdata/OLR/REGolr.png X800 Y600 white'
'print'
'disable print'
;

望晴檐

谢谢楼主的知识普及了~~

shirly

mofangbao 发表于 2013-2-28 22:46

                               
登录/注册后可看大图

不加enable print的话，直接print 就是ps格式的矢量图了，加了的话就是gmf格式的矢量图，可以通过grads.exe ...

版主你好，你说的这句话“不加enable print的话，直接print 就是ps格式的矢量图了”，我不太理解，希望不要鄙视我，我已经看过GrADS Documentation Index了，但是还是不太明白，它的例子是这样的：gxps -cr -i mytest.mf -o mytest.ps 。但是按它的写gs，运行会出错。我是这样做的：之前enable print了，后面也print生成了一个a0.gmf的文件，接着再转化成PS格式：'!gxps -o -i d:/lw/dat/pet/regr/tu/a0.gmf -c d:/lw/dat/pet/regr/tu/a1.ps'，最后在我规定的文件夹里形成了一个2M多的gmf文件和一个6M多的PS文件，麻烦你看一眼我这样写对不对。谢谢清风了

mofangbao

shirly 发表于 2014-3-20 23:41

                               
登录/注册后可看大图

版主你好，你说的这句话“不加enable print的话，直接print 就是ps格式的矢量图了”，我不太理解，希望不 ...

你贴完整的出来看看   还有打开你的os文件看看

Nuist_jyc

{:eb513:}{:eb513:}

shirly

mofangbao 发表于 2014-3-21 00:10

                               
登录/注册后可看大图

你贴完整的出来看看   还有打开你的os文件看看

首先感谢版主回复~！上面这个gs比较长，分区域在一张上画了六个图~可以跳阅~
因为是第一次接触这个，正在学习怎么打开PS格式图，所以暂时不能查看，只能麻烦清风大侠帮我看一下这种生成PS图的方法对不对。感谢

'reinit'
'open d:/lw/dat/pet/china.1950.2009.tmp.x.ctl'
'open d:/lw/dat/pet/china.1950.2009.cld.x.ctl'
'open d:/lw/dat/inter/sta2grdx_grid.dswf.ctl'
'open d:/lw/dat/pet/china.1950.2009.tmn.x.ctl'
'open d:/lw/dat/pet/china.1950.2009.tmx.x.ctl'
'open d:/lw/dat/pet/china.1950.2009.dtr.x.ctl'
'open d:/lw/dat/pet/regr/c_time.ctl'
'enable print d:/lw/dat/pet/regr/tu/a0.gmf'
'define_colors'
'set grads off'
'set grid off'
'set mpdset hires'

'set x 1 '
'set y 1 '
'set t 1 720'
'define a=aave(tim.7,lon=72.75,lon=135.75,lat=16.75,lat=55.75)'

'set parea 0.35 3.55 4.8 8.0'
'set x 1 127'
'set y 1 79'
'set z 1'
'set t 1'
'set gxout shaded'
'define regr2=tregr(a,tmp.1,t=1,t=720)'
'set gxout shaded'
'set clevs    -0.05 0  0.1 0.2 0.3 0.4 0.5 0.6'
'set ccols   44 43  21 22 23 24 25 27 28 '
'cnbasemap regr2*120'
'set strsiz 0.15'
'draw string 0.5 7.3 a)'
'cbarn 0.4 0 2.1 4.7'

'set parea 4.0 7.2 4.8 8.0'
'set x 1 127'
'set y 1 79'
'set z 1'
'set t 1'
'set gxout shaded'
'define regr3=tregr(a,cld.2,t=1,t=720)'
'set gxout shaded'
'set clevs  -1.5 -1.2 -0.9 -0.6 -0.3 0 0.3 0.6 0.9 1.2 1.5 '
'set ccols  48 47 46 45 44 42  22  24 25 26 27 28 '
'cnbasemap regr3*120'
'set strsiz 0.15'
'draw string 4.2 7.5 b)'
'cbarn 0.4 0 5.65 4.7'

'set parea 7.6 10.8 4.8 8.0'
'set x 1 127'
'set y 1 79'
'set z 1'
'set t 1'
'set gxout shaded'
'define regr1=tregr(a,dswf.3,t=1,t=720)'
'set gxout shaded'
'set clevs -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 4 5'
'set ccols 45 44 43 42 0 21 22 23 24 25 26 27 28'
'cnbasemap regr1*120'
'set strsiz 0.15'
'draw string 7.8 7.5 c)'
'cbarn 0.4 0 9.2 4.7'

'set parea 0.35 3.55 1.4 4.6'
'set x 1 127'
'set y 1 79'
'set z 1'
'set t 1'
'set gxout shaded'
'define regr4=tregr(a,tmn.4,t=1,t=720)'
'set gxout shaded'
'set clevs    -0.05 0  0.1 0.2 0.3 0.4 0.5 0.6'
'set ccols   44 43  21 22 23 24 25 27 28 '
'cnbasemap regr4*120'
'set strsiz 0.15'
'draw string 0.5 4.1 d)'
'cbarn 0.4 0 2.1 1.3'

'set parea 4.0 7.2 1.4 4.6'
'set x 1 127'
'set y 1 79'
'set z 1'
'set t 1'
'set gxout shaded'
'define regr5=tregr(a,tmx.5,t=1,t=720)'
'set gxout shaded'
'set clevs    -0.05 0  0.1 0.2 0.3 0.4 0.5 0.6'
'set ccols   44 43  21 22 23 24 25 27 28 '
'cnbasemap regr5*120'
'set strsiz 0.15'
'draw string 4.2 4.1 e)'
'cbarn 0.4 0 5.65 1.3'

'set parea 7.6 10.8 1.4 4.6'
'set x 1 127'
'set y 1 79'
'set z 1'
'set t 1'
'set gxout shaded'
'define regr6=tregr(a,dtr.6,t=1,t=720)'
'set gxout shaded'
'set clevs    -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0  0.05 0.1 0.15'
'set ccols 48 47 46 45 44 43 42  22 23 24 25'
'cnbasemap regr6*120'
'set strsiz 0.15'
'draw string 7.8 4.1 f)'
'cbarn 0.4 0 9.2 1.3'
'print d:/lw/dat/pet/regr/tu/a0.gmf'
'disable print'

'!gxps -o -i d:/lw/dat/pet/regr/tu/a0.gmf -c d:/lw/dat/pet/regr/tu/a1.ps'

'reinit'

*'!gxeps -i d:/lw/dat/pet/regr/tu/a0.gmf -c -o d:/lw/dat/pet/regr/tu/a1.ps'
***这种方法也能出现一个ps格式的文件

台风下的宁静

支持一下楼主