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

clare

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


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)

传说中的谁

不会啊，我原来也批量画过的

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. 
    
41. rc=plotskew(t,td,uv,dir)
    
42. 
    
43. 'q time'
    
44. str=subwrd(result,3)
    
45. year=substr(str,9,4)
    
46. month=substr(str,6,3)
    
47. day=substr(str,4,2)
    
48. hour=substr(str,1,2)
    
49. if(month=JAN)
    
50.     month=01  
    
51. endif
    
52. if(month=FEB)
    
53.     month=02  
    
54. endif
    
55. if(month=MAR)
    
56.     month=03  
    
57. endif
    
58. if(month=APR)
    
59.     month=04  
    
60. endif
    
61. 
    
62. if(month=MAY)
    
63.     month=05  
    
64. endif
    
65. if(month=JUN)
    
66.     month=06  
    
67. endif
    
68. if(month=JUL)
    
69.     month=07  
    
70. endif
    
71. if(month=AUG)
    
72.     month=08  
    
73. endif
    
74. if(month=SEP)
    
75.     month=09
    
76. endif
    
77. if(month=OCT)
    
78.     month=10  
    
79. endif
    
80. if(month=NOV)
    
81.     month=11  
    
82. endif
    
83. if(month=DEC)
    
84.     month=12  
    
85. endif
    
86. say year '-' month '-'day' 'hour':00 UTC'
    
87. 'draw title 'year'-' month'-'day' 'hour':00 UTC '
    
88. 'printim E:\NCEP\tlogp\'year''month''day''hour'.png black'
    
89. 'c'
    
90. tt=tt+1
    
91. endwhile
    
92. 
    
93. 
    
94. function plotskew(sndtemp,snddewp,sndspd,snddir)

复制代码

clare

传说中的谁 发表于 2015-4-24 10:59
不会啊，我原来也批量画过的

还是老样子，可能我的哪个地方被不小心动着了。我再查查看。不过又学到一招关于画时间的好方法。赞！

nzy718d

看两位高人在谈论天书