basemap 报错

容与

想画中国范围的图，中国陆地以外的留白，但basemap每次都报错，gs文件如下：
'reinit'
'c'
'open D:\grads\shortclimate\shixi3\high\grid.ctl'
'open D:\grads\shortclimate\shixi3\high\grid1.ctl'
'set mpdset cnhimap'
'set map 1 3 1'
'draw map'
'set grads off'
'set grid off'
'set lon 73 135.5'
'set lat 17 55'

'define a=oacres(g,r.2,49,49,20,49,1)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'set gxout shaded'
'set csmooth on'
*'set black -0.15 0.15'
'd aa'
'basemap cnout 0 1'
'set gxout contour'
'set csmooth on'
'd aa'
*'D:\grads\shortclimate\shixi3\taizhan\map\river.gs'
*'D:\grads\shortclimate\shixi3\taizhan\map\china.gs'
*'D:\grads\shortclimate\shixi3\taizhan\map\hainan.gs'
*'D:\grads\shortclimate\shixi3\taizhan\map\taiwan.gs'
'D:\grads\shortclimate\shixi3\taizhan\map\southsea.gs'
;

                               
登录/注册后可看大图

                               
登录/注册后可看大图

river

配合basemap用的有好几个文件，你都有吗，路径都对吗？论坛里有cnbasemap，楼主自己搜索一下

容与

river 发表于 2014-11-16 10:27
配合basemap用的有好几个文件，你都有吗，路径都对吗？论坛里有cnbasemap，楼主自己搜索一下

谢了，顺便问一下，用了这个帖子（http://bbs.06climate.com/forum.php?mod=viewthread&tid=625&extra=page%3D1%26filter%3Dtypeid%26typeid%3D6）的cnbasemap后，有一个问题，我南海的地图盖掉了，这个怎么办

                               
登录/注册后可看大图

river

容与 发表于 2014-11-16 11:35
谢了，顺便问一下，用了这个帖子（http://bbs.06climate.com/forum.php?mod=viewthread&tid=625&extra=pa ...

你不是有southsea.gs吗，最后运行一下就出来了吧

容与

river 发表于 2014-11-16 11:44
你不是有southsea.gs吗，最后运行一下就出来了吧

那个gs文件用，但还是看不到，我又用了（http://bbs.06climate.com/forum.p ... typeid%26typeid%3D1）的南海地图，结果好一点，但是他的south_value.gs有个问题：
****************************************
if(ptype=contour|ptype= Contour)
'q contours'
rec=sublin(result,1)
ncol=subwrd(rec,5)
nlev=ncol
ii=1
levs=''
cols=''
while(ii<=ncol)
     rec=sublin(result,ii+1)
     col.ii=subwrd(rec,1)
     lev.ii=subwrd(rec,2)
   levs=levs%' '%lev.ii
   cols=cols%' '%col.ii   
ii=ii+1
endwhile
endif


作者说是自动选择shaded、contour、vector等信息，第二行的q 会在grads里报错，我删了它，并在出图那段改了一下，画阴影+等值图：
* 段1 尾: 在南海诸岛附图中标明经纬度。
* 段2 首: 在南海诸岛附图中显示与主图一致的物理量分布。
'set grid off'
'set grads off'
'set gxout 'ptype
*'set time aug2013'
*'set t 1'
*say 111111
*say ptype
if(ptype=Contour | ptype=Shaded | ptype=Fgrid | ptype=shade2b)
say levs
say cols
'set clevs 'levs
'set ccols 'cols
'set gxout shaded'
'set csmooth on'
'set black -0.15 0.15'
'set cint 0.1'
'd 'v1''
'set gxout contour'
'set cint 0.1'
'd 'v1''
endif


效果如下：

                               
登录/注册后可看大图

我想让子图的坐标系变成其他颜色，试着改'set xlopts 4 3 0.12'，但怎么改都是上面那样
* 段1 首: 在南海诸岛附图中标明经纬度。
'set lon 'lon1' 'lon2''
'set lat 'lat1' 'lat2''
'set cthick 1'
'set xlopts 2 3 0.12'
'set ylopts 2 3 0.12'
'set xlab on'
'set ylab on'
'set xlpos -20.0'
'set ylpos -20.0'
*'set ylevs 2 10 18'
*'set xlevs 110 120'
'set lon 105 123'
'set lat 0 25'
'set parea 'xleft' 'xright' 'ybottom' 'ytop''
* 段1 尾: 在南海诸岛附图中标明经纬度。

river

容与 发表于 2014-11-16 14:11
那个gs文件用，但还是看不到，我又用了（http://bbs.06climate.com/forum.php?mod=viewthread&tid=28492& ...

你说的第二行q contours在我这是不报错的，我用的grads2.1，你换个版本看看。还你说的那个命令是改变label的颜色，粗细，不是坐标轴的。

容与

river 发表于 2014-11-16 16:37
你说的第二行q contours在我这是不报错的，我用的grads2.1，你换个版本看看。还你说的那个命令是改变labe ...

谢谢river前辈这几天叨扰了（课太多了，回帖晚了），我改了一下程序就可以了（把cnbasemap 换成红色语句就可以了，子图色标不同在south_value里改一下就可以了（这个暂时没改，））
'reinit'
'c'
'open D:\grads\shortclimate\shixi3\high\grid.ctl'
'open D:\grads\shortclimate\shixi3\high\grid1.ctl'
'set grads off'
'set grid off'
'set lon 73 135.5'
'set lat 17 55'
'define a=oacres(g,r.2,49,49,20,49,1)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'set gxout shaded'
'set csmooth on'
'set black -0.15 0.15'
'd aa'
'set gxout contour'
'set ccolor 3'
'set csmooth on'
'd aa'
'basemap cnout 0 1'
*'cnbasemap aa 1'
'set mpdset cnhimap'
'set map 1 3 1'
'draw map'
*'D:\grads\shortclimate\shixi3\taizhan\map\river.gs'
'southsea_value -v1 aa -lab on'
*'province-basemap anhui aa 1'
*'D:\grads\shortclimate\shixi3\taizhan\map\china.gs'
*'D:\grads\shortclimate\shixi3\taizhan\map\hainan.gs'
*'D:\grads\shortclimate\shixi3\taizhan\map\taiwan.gs'
*'D:\grads\shortclimate\shixi3\taizhan\map\southsea.gs'
;

samvelly

{:eb502:}{:eb502:}

thornrone

菜鸟来观贴