爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5457|回复: 2

[分享资料] gs文件一次性画出1958年5-9月,逐候500hPa环流场

[复制链接]

新浪微博达人勋

发表于 2012-3-8 09:41:23 | 显示全部楼层 |阅读模式

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

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

x
想用gs文件一次性画出1958年5-9月,逐候500hPa环流场,这是我的gs文件,请大家看看有什么问题,可以输出但是一个文件有N个图叠加,不知道原因,请大家指点!
'reinit'

iyear=1958
while(iyear<=1958)
'sdfopen D:\lw2011\huanliu\hgt.1958.nc'
'set lat 0 90'
'set lon 40 180'
*'set t 1'
'set lev 500'
'set xlopts 1  1  0.17'
'set ylopts 1  1  0.17'
'set grid off'
'set grads off'
'set map 1 1 10'
'set mpdset hires'
'set mpdset cnworld'
mon.5=may
mon.6=jun
mon.7=jul
mon.8=aug
mon.9=sep

*-----------------------------------------------
a=5
while(a<=9)
  say mon.a'  '
  b=1
  while(b<=5)
  t1=(b-1)*5+1
  t2=(b-1)*5+5
  'define zz=ave(hgt,time='t1''mon.a''iyear',time='t2''mon.a''iyear')'
  'd zz/10.0'
pull dummy
'enable print D:\lw2011\huanliu\'t1''mon.a'-'t2''mon.a''iyear'.500.gmf'
'print'
'disable print'

  b=b+1
  endwhile
  t1=26
  t2=31
  
  if(a=6);t2=30;endif
  if(a=9);t2=30;endif
  'define zz=ave(hgt,time='t1''mon.a''iyear',time='t2''mon.a''iyear')'
'set cint 2'
'set gxout contour'
'set cthick 6'
'set cterp on'
'set csmooth on'
'd zz/10.0'
pull dummy
'enable print D:\lw2011\huanliu\'t1''mon.a'-'t2''mon.a''iyear'.500.gmf'
'print'
'disable print'
a=a+1
endwhile

'close 1'
iyear=iyear+1
endwhile
'reinit'

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 19710
发表于 2012-3-8 09:46:58 | 显示全部楼层
画完一张图后面清除掉就行啦,加个'c'~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-3-8 10:11:33 | 显示全部楼层
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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