爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3658|回复: 0

[秀图] grads绘图警告提示,该怎么解决?

[复制链接]

新浪微博达人勋

发表于 2014-6-21 23:37:46 | 显示全部楼层 |阅读模式

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

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

x
绘图出错信息:gs文件如下:请帮忙纠错!
*light white to dark gray   
'set rgb 80 255 255 255'
'set rgb 81 216 216 216'
'set rgb 82 191 191 191'
'set rgb 83 165 165 165'
'set rgb 84 145 145 145'
'set rgb 85 128 128 128'   
'set rgb 86 100 100 100'  
'set rgb 87 88 88 88'
'set rgb 88 77 77 77'
'set rgb 89 65 65 65'

'reinit'
'open D:\2012tshan\newpic\20130617\20130617.ctl'
'enable print D:\2012tshan\newpic\20130617\11\pm41N.gmf'

t=1
while(t<=17)
'c'


'set t 't
'open D:\2012tshan\newpic\20130617\11\topo10.ctl'

'query time'
  res=subwrd(result,3)
  year=substr(res,11,2)
  month=substr(res,6,3)
  day=substr(res,4,2)
  hour=substr(res,1,2)
  if(month='JAN'); month='01'; mtype='big'; endif
  if(month='FEB'); month='02'; mtype='feb'; endif
  if(month='MAR'); month='03'; mtype='big'; endif
  if(month='APR'); month='04'; mtype='sma'; endif
  if(month='MAY'); month='05'; mtype='big'; endif
  if(month='JUN'); month='06'; mtype='sma'; endif
  if(month='JUL'); month='07'; mtype='big'; endif
  if(month='AUG'); month='08'; mtype='big'; endif
  if(month='SEP'); month='09'; mtype='sma'; endif
  if(month='OCT'); month='10'; mtype='big'; endif
  if(month='NOV'); month='11'; mtype='sma'; endif
  if(month='DEC'); month='12'; mtype='big'; endif
  endmonth=FALSE
  if(hour<16)
    hour=hour+8
  else
    hour=hour+8-24
    if((mtype='big' & day=31) | (mtype='sma' & day=30))
      month=month+1
      if(month<10); month='0'%month; endif
      if(month=13); month='01'; endif
      endmonth=TRUE
    endif
    if(mtype='feb' & day=28)
      month='03'
      endmonth=TRUE
    endif
    if(endmonth=FALSE)
      day=day+1
      if(day<10); day='0'%day; endif
    else
      day='01'
    endif
  endif
  if(hour<10); hour='0'%hour; endif
'set ccolor 4 '
  'set strsiz 0.3 0.8 '
  'draw title 'year'-'month'-'day' 'hour':00'

'set parea 1.0 10.0 1.2 8.0'
'set lev 1000 100'
'set lat 30 50'
'set lon 60 100'
'define vor=hcurl(ugrdprs,vgrdprs)*1e5'

************************thetase***********************
'define tt=TMPprs'
'define es=6.11*exp(17.26*(tt-273.16)/(tt-35.86))'
'define qs=0.622*es/(lev-0.378*es)'
'define qq=RHprs/100.*qs'
'define prs=lev*100'
'define el=log((qq/0.62197)*(prs/100.)/(1.+(qq/0.62197)))'
'define td=273.15+(243.5*el-440.8)/(19.48-el)'
'define tc=tt-1.21*(tt-td)'
'define e=qq/(qq+0.622)*prs'
'define cpm=1004*(1+0.86*qq)'
'define rm=287.04*(1+0.608*qq)'
'define pc=prs*pow(tc/tt,cpm/rm)'
'define st=tc*pow(100000/(pc-e),287.04/1004)'
'define lvc=2500.8*1000-2274.4*(tt-273.15)'
'define qse=st*exp(qq*lvc/1004/tc)'


'set ylab on'
'set xlab on'

'set zlog on'
'set lat 41'
'set lon 75 85'
'set xlopts 1 6 0.16'
'set ylopts 1 6 0.16'
'set mpdset hires'
'set grads off'
'set grid off'
'set lev 1000 100'
'set gxout shaded'
'd qse'


'set gxout contour'
'set cthick 6'
*'set cstyle 4'
'set ccolor 1'
*'set cmax 0'
'd vvelprs*100'

'set cstyle 4'
'set cthick 9'
'set ccolor 1'
'd qse'

'set gxout vector'
'set ccolor 1'
'set arrscl 0.2 20'
'd ugrdprs;vgrdprs'

'set dfile 2'
'set ylab off'
'set xlab off'
'set x 1'
'set y 1'
'set vrange  0  17000'
'set rgb 16 144 144 144'
'set gxout linefill'
'set lfcols 16 0'
*'set lat 37 45'
'set lon 75 85'
'set z 1'
'set t 1'
*'define cc=ave(tp,lon=70,lon=80)'
*'set lon 77'
'define cc=ave(tp,lat=40.5,lat=41.5)'
'set lat 41'
'define a0=tp-tp'
'd cc;a0'

'close 2'

'print'
'c'
t=t+1
endwhile
'disable print'

return

'reinit'



绘图出错信息!.JPG
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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