请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4460|回复: 3

请问plotskew.gs参数中关于外框的是哪一个

[复制链接]

新浪微博达人勋

发表于 2015-4-24 10:46:03 | 显示全部楼层 |阅读模式
GrADS
系统平台: win7
问题截图: -
问题概况: 参照传说中的谁--GrADS中斜温图的绘制(Skew-T)帖子试图while批量出图,除了第一张有外框外,其余均无,想请教这个参数是哪一个。或是应该怎么处理?
我看过提问的智慧: 看过
自己思考时长(天): 2

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

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

x
参照“@传说中的谁 ”的帖子http://bbs.06climate.com/forum.p ... =%D0%B1%CE%C2%CD%BC试图一次性批量出图,但找不到plotskew.gs里面的外框参数,所以只有第一张有外框,其它无。求指点~
QQ截图20150424102702.png 1.png

gs如下:
'reinit'
'open E:\FN\data\fnl.ctl'
'set lon 100 125'
'set lat 15 45'
'set lev 1000 100'
'set parea 1 10.5 1.0 8.0'
'set mpdset cnworld'
'set grid off'
'set grads off'
'set frame on'
*'set display color white'
*'c'

nb=4
tt=1
while(tt<=nb)
'set t 'tt''
'set grads off'
'set grid off'
*********计算整层可降水量*************
'set lat 28.33'
'set lon 109.72'
'set lev 1000 100'

'define t=tmpprs-273.16'
'define RH=rhprs'
'define td=t-((14.55+0.114*t)*(1-0.01*RH) + pow((2.5+0.007*t)*(1-0.01*RH),3) + (15.9+0.37*T)*pow((1-0.01*RH),14))'
*计算露点温度
'define uv=mag(vgrdprs,ugrdprs)'
*计算风速
'define pi=3.14159'
if (ugrdprs>=0.)
if (vgrdprs>=0.)
'define dir=180+atan2(ugrdprs,vgrdprs)*180/pi'
endif
endif
if (vgrdprs<=0.)
if (ugrdprs>=0.)
'define dir=270+atan2(ugrdprs,vgrdprs)*180/pi'
endif
endif
if (ugrdprs<=0.)
if (vgrdprs>=0.)
'define dir=180-atan2(ugrdprs,vgrdprs)*180/pi'
endif
endif
if (vgrdprs<0.)
if (ugrdprs<=0.)
'define dir=270-atan2(ugrdprs,vgrdprs)*180/pi'
endif
endif
*if(ugrdprs > 0.);'define dir=90-180/pi*atan2(vgrdprs,ugrdprs)';endif
*if(ugrdprs < 0.);'define dir=270-180/pi*atan2(vgrdprs,ugrdprs)';endif
rc=plotskew(t,td,uv,dir)

'q dims'
rec=sublin(result,5)
recs=subwrd(rec,6)
' draw title 'recs''
'printim E:\FN\result\tlog-'tt'.png white x1200 y1000'
'print'
'c'
tt=tt+1
endwhile

*'disable print'


function plotskew(sndtemp,snddewp,sndspd,snddir)

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

新浪微博达人勋

发表于 2015-4-24 10:59:27 | 显示全部楼层
不会啊,我原来也批量画过的
  1. 'reinit'
  2. 'sdfopen Z:\ncep\uwnd.2009.nc'
  3. 'sdfopen Z:\ncep\vwnd.2009.nc'
  4. 'sdfopen Z:\ncep\air.2009.nc'
  5. 'sdfopen Z:\ncep\rhum.2009.nc'
  6. tt=1
  7. while(tt<=1460)
  8. 'set t   'tt
  9. 'set grads off'
  10. 'set grid off'
  11. 'set csmooth on'
  12. 'set lon 112.5'
  13. 'set lat 22.5'
  14. 'set lev 1000 200'
  15. 'define t=air.3-273.15'
  16. 'define rH=rhum.4'
  17. 'define td=t-((14.55+0.114*t)*(1-0.01*RH) + pow((2.5+0.007*t)*(1-0.01*RH),3) + (15.9+0.37*T)*pow((1-0.01*RH),14))'
  18. 'define uv=mag(uwnd,vwnd.2)'
  19. 'define pi=3.14159'
  20. if (uwnd>=0.)
  21. if (vwnd.2>=0.)
  22. 'define dir=180+atan2(uwnd,vwnd.2)*180/pi'
  23. endif
  24. endif
  25. if (vwnd.2<=0.)
  26. if (uwnd>=0.)
  27. 'define dir=270+atan2(uwnd,vwnd.2)*180/pi'
  28. endif
  29. endif
  30. if (uwnd<=0.)
  31. if (vwnd.2>=0.)
  32. 'define dir=180-atan2(uwnd,vwnd.2)*180/pi'
  33. endif
  34. endif
  35. if (vwnd.2<0.)
  36. if (uwnd<=0.)
  37. 'define dir=270-atan2(uwnd,vwnd.2)*180/pi'
  38. endif
  39. endif

  40. rc=plotskew(t,td,uv,dir)

  41. 'q time'
  42. str=subwrd(result,3)
  43. year=substr(str,9,4)
  44. month=substr(str,6,3)
  45. day=substr(str,4,2)
  46. hour=substr(str,1,2)
  47. if(month=JAN)
  48.     month=01  
  49. endif
  50. if(month=FEB)
  51.     month=02  
  52. endif
  53. if(month=MAR)
  54.     month=03  
  55. endif
  56. if(month=APR)
  57.     month=04  
  58. endif

  59. if(month=MAY)
  60.     month=05  
  61. endif
  62. if(month=JUN)
  63.     month=06  
  64. endif
  65. if(month=JUL)
  66.     month=07  
  67. endif
  68. if(month=AUG)
  69.     month=08  
  70. endif
  71. if(month=SEP)
  72.     month=09
  73. endif
  74. if(month=OCT)
  75.     month=10  
  76. endif
  77. if(month=NOV)
  78.     month=11  
  79. endif
  80. if(month=DEC)
  81.     month=12  
  82. endif
  83. say year '-' month '-'day' 'hour':00 UTC'
  84. 'draw title 'year'-' month'-'day' 'hour':00 UTC '
  85. 'printim E:\NCEP\tlogp\'year''month''day''hour'.png black'
  86. 'c'
  87. tt=tt+1
  88. endwhile


  89. function plotskew(sndtemp,snddewp,sndspd,snddir)
复制代码
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-4-24 11:49:31 | 显示全部楼层
传说中的谁 发表于 2015-4-24 10:59
不会啊,我原来也批量画过的

还是老样子,可能我的哪个地方被不小心动着了。我再查查看。不过又学到一招关于画时间的好方法。赞!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-5-19 07:04:22 | 显示全部楼层
看两位高人在谈论天书
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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