爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6999|回复: 0

[图形美化] GRADS画出来的图国境外的被屏蔽了,怎么样能显示出来

[复制链接]

新浪微博达人勋

发表于 2020-4-27 14:32:23 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
gs文件
                                                                                                                                                                                                                      p.1='a';p.2='b';p.3='c';p.4='d'
str.1='50km-Obs,DJF';str.2='50km-subex,DJF';str.3='50km-Topmkins,DJF';str.4='50km-WSM5,DJF'
'reinit'
'open /datadir2/home/xxk/Resolution/icbc1/ICBC50ALL.ctl'
'open /datadir2/home/xxk/Resolution/outputs1/subex50/ATMALL.ctl'
'open /datadir2/home/xxk/Resolution/outputs1/Topmkins50/ATMALL.ctl'
'open /datadir2/home/xxk/Resolution/outputs1/WSM550/ATMALL.ctl'
itmpa=2000
while(itmpa<=2010)
'set dfile 1'
'set map 1 1 2'
'set mpdraw on'
'set lon 70 150'
'set lat 15 60'
'set lev 850'
'set grid off'
'set grads off'
'define at'itmpa'=ave(t.1,time=00z01JAN'itmpa',time=18z31JAN'itmpa')'
'define au'itmpa'=ave(u.1,time=00z01JAN'itmpa',time=18z31JAN'itmpa')'
'define av'itmpa'=ave(v.1,time=00z01JAN'itmpa',time=18z31JAN'itmpa')'
itmpa=itmpa+1
endwhile
'define T1=(at2000+at2001+at2002+at2003+at2004+at2005+at2006+at2007+at2008+at2009+at2010)/11-273'
'define u1=(au2000+au2001+au2002+au2003+au2004+au2005+au2006+au2007+au2008+au2009+au2010)/11'
'define v1=(av2000+av2001+av2002+av2003+av2004+av2005+av2006+av2007+av2008+av2009+av2010)/11'

itmpb=2000
while(itmpb<=2010)
'set dfile 2'
'set map 1 1 2'
'set mpdraw on'
'set lon 70 150'
'set lat 15 60'
'set lev 850'
'set grid off'
'set grads off'
'define bt'itmpb'=ave(ta.2,time=00z01JAN'itmpb',time=18z31JAN'itmpb')'
'define bu'itmpb'=ave(ua.2,time=00z01JAN'itmpb',time=18z31JAN'itmpb')'
'define bv'itmpb'=ave(va.2,time=00z01JAN'itmpb',time=18z31JAN'itmpb')'
itmpb=itmpb+1
endwhile
'define T2=(bt2000+bt2001+bt2002+bt2003+bt2004+bt2005+bt2006+bt2007+bt2008+bt2009+bt2010)/11-273'
'define u2=(bu2000+bu2001+bu2002+bu2003+bu2004+bu2005+bu2006+bu2007+bu2008+bu2009+bu2010)/11'
'define v2=(bv2000+bv2001+bv2002+bv2003+bv2004+bv2005+bv2006+bv2007+bv2008+bv2009+bv2010)/11'

itmpc=2000
while(itmpc<=2010)
'set dfile 3'
'set map 1 1 2'
'set mpdraw on'
'set lon 70 150'
'set lat 15 60'
'set lev 850'
'set grid off'
'set grads off'
'define ct'itmpc'=ave(ta.3,time=00z01JAN'itmpc',time=18z31JAN'itmpc')'
'define cu'itmpc'=ave(ua.3,time=00z01JAN'itmpc',time=18z31JAN'itmpc')'
'define cv'itmpc'=ave(va.3,time=00z01JAN'itmpc',time=18z31JAN'itmpc')'
itmpc=itmpc+1
endwhile
'define T3=(ct2000+ct2001+ct2002+ct2003+ct2004+ct2005+ct2006+ct2007+ct2008+ct2009+ct2010)/11-273'
'define u3=(cu2000+cu2001+cu2002+cu2003+cu2004+cu2005+cu2006+cu2007+cu2008+cu2009+cu2010)/11'
'define v3=(cv2000+cv2001+cv2002+cv2003+cv2004+cv2005+cv2006+cv2007+cv2008+cv2009+cv2010)/11'

itmpd=2000
while(itmpd<=2010)
'set dfile 4'
'set map 1 1 2'
'set mpdraw on'
'set lon 70 150'
'set lat 15 60'
'set lev 850'
'set grid off'
'set grads off'
'define dt'itmpd'=ave(ta.4,time=00z01JAN'itmpd',time=18z31JAN'itmpd')'
'define du'itmpd'=ave(ua.4,time=00z01JAN'itmpd',time=18z31JAN'itmpd')'
'define dv'itmpd'=ave(va.4,time=00z01JAN'itmpd',time=18z31JAN'itmpd')'
itmpd=itmpd+1
endwhile
'define T4=(dt2000+dt2001+dt2002+dt2003+dt2004+dt2005+dt2006+dt2007+dt2008+dt2009+dt2010)/11-273'
'define u4=(du2000+du2001+du2002+du2003+du2004+du2005+du2006+du2007+du2008+du2009+du2010)/11'
'define v4=(dv2000+dv2001+dv2002+dv2003+dv2004+dv2005+dv2006+dv2007+dv2008+dv2009+dv2010)/11'

p=1
i=1
while(p<=4)
'set parea 1.0 11 1.0 8.0'
'subplot_x 1 4 'p''
'set map 1 1 2'
'set mpdraw on'
'set lon 70 150'
'set lat 15 60'
'set lev 850'
'set grid off'
'set grads off'
'set gxout shaded'
'run cbarn 1 0'
'set xlint 15'
'set ylint 10'
*'set xlopts 1 5 0.27'
*'set ylopts 1 5 0.27'
'set font 2'
'd T'p''
*'set gxoud contour'
'basemap cnout 0 15'
'river 15'
'xscs 1'
*'set string 1 tl 5'
*'draw string  1.6 7.9 'str.p''
p=p+1
i=i+1
endwhile
'printim tem1_uv.png'

画出来的图把国境外的屏蔽了,想要显示出来,不知道怎么设置,还有
'basemap cnout 0 15'
'river 15'
'xscs 1'
这三个没找到是设置什么的
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表