GrADS cnriver 黄河，长江，各省行政边界

ashleywu

请教各位，原先用set mpdset cnriver的时候出来的是只有长江黄河的地图。

现在的

今天重装了下Grads，set mpdset cnrive 出来的还有各省行政边界。

原先的

Grads 为了用cnbasemap.gs，重新拷了地图进去，是不是拷的时候有问题？出图的gs文件完全没有动过。
gs如下：
'reinit'
'open h:/GPCP/gpcp_v2.2_psg_swap.ctl'
'enable print h:\GPCP\strong5yr2.gmf'
'set vpage 0 11 0 8.5'
'set parea 1.0 10.0 1.0 8.0'
'set grads off'
'set grid off'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set xlint 10'
'set lat 18 54'
'SET lon 72 138'
'set ylint 10'
'set mpdset cnriver'
'set map 1 1 5'
'draw map'
'set ccolor 0'

'define a=precip(time=jun1985)+precip(time=jun1994)+precip(time=jun1997)+precip(time=jun1999)+precip(time=jun2002)'

'define c=precip(time=jul1985)+precip(time=jul1994)+precip(time=jul1997)+precip(time=jul1999)+precip(time=jul2002)'

'define e=precip(time=aug1985)+precip(time=aug1994)+precip(time=aug1997)+precip(time=aug1999)+precip(time=aug2002)'

'define g=ave(precip,time=jun1979,time=jun2010,1yr)'

'define i=ave(precip,time=jul1979,time=jul2010,1yr)'

'define k=ave(precip,time=aug1979,time=aug2010,1yr)'

'define m=(a/5+c/5+e/5-g-i-k)/3'
'set gxout shaded'
'set black -0.5 0.5'
'set csmooth on'
'd m'

'set gxout contour'
'set cthick 7'
'd m'
'cnbasemap m'
'run  H:\GPCP\scs_map.gs'
'set string 1 c 6 0'
'set strsiz 0.2 0.2'
'draw string 1.4 7.2 (a)'
'print'
'disable print'
'reinit'

kongfeng0824

感谢分享！！！！

jacci_zou

用同样的GS，两次出的图不一样么？，看GS里面你画了几次地图诶！

river

draw map 之前再加上一句set poli off 就可以了吧

ashleywu

jacci_zou 发表于 2013-7-3 18:08

                               
登录/注册后可看大图

用同样的GS，两次出的图不一样么？，看GS里面你画了几次地图诶！

是啊，gs没变，就是重装了grads，想再运行下gs文件，结果就不一样了，多了各省的行政边界。

ashleywu

river 发表于 2013-7-3 18:11

                               
登录/注册后可看大图

draw map 之前再加上一句set poli off 就可以了吧

set 以后中国边界，省边界，黄河都没有了。。。

mofangbao

打开cnbasemap里面设置地图的语句注释掉  因为你原来没有cnworld这个地图，所以cnbasemap的时候没出现，现在你重新安装的GrADS里面带了这个地图，所以就显示出来的了

ashleywu

mofangbao 发表于 2013-7-3 18:50

                               
登录/注册后可看大图

打开cnbasemap里面设置地图的语句注释掉  因为你原来没有cnworld这个地图，所以cnbasemap的时候没出现，现在 ...

大神就是大神，膜拜。cnbasemap.gs里
*     Draw the polygon
      'set line 'color
      cmd
    endif
  endwhile
'draw map'
注释掉*'draw map'就没有各省的行政边界了，唯一的不足是全国的边界有点不清晰了，断断续续的感觉。可有方法改进？

mofangbao

ashleywu 发表于 2013-7-3 19:28

                               
登录/注册后可看大图

大神就是大神，膜拜。cnbasemap.gs里
*     Draw the polygon
      'set line 'color

自己做个清晰点的底图就行了,用MeteoInfo可以方便的制作，论坛有教程，请搜索
还有 不要注释draw map  而是注释程序开始那个 set mpdset ...的一句话

风云翩跹

前来学习，感谢分享！