- 积分
- 1235
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-11-18
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 酒国 于 2015-11-2 11:12 编辑
程序是求几个时次水汽通量和水平风场的平均值。加粗字体为调用ave函数,却只显示前面一个时次的水汽和风场。
运行的错误提示是
ga-> D:\fnl\shan\stream.gs
All files closed; all defined objects released;
All GrADS attributes have been reinitialized
Close Error: Missing file number
Close Error: Missing file number
GXPRINT error: missing output file name
DISABLE error: Invalid keyword
ga->
'reinit'
'open D:\fnl\shan\fnl_20121024_00_00_c.ctl'
'set grid off'
'set grads off'
'set lon 107 112'
'set lat 17 20.5'
'set lev 850'
'set t 14'
'set mpdset hires'
'set mpdset cnworld'
'set map 1 1 8'
'set cthick 6'
'set ccolor 1'
'set csmooth on'
'set xlopts 1 6 0.16'
'set ylopts 1 6 0.16'
'set xlint 0.5'
'set ylint 0.5'
'set gxout shaded'
'set rgb 21 234 233 244'
'set rgb 22 196 195 206'
'set rgb 23 158 157 168'
'set rgb 24 120 119 130'
'set rgb 25 82 81 92'
'set rgb 26 44 43 54'
'set clevs 0.013 0.014 0.016 0.019 0.022'
'set ccols 0 22 23 24 25 26'
'define e=(6.112*exp((17.67*(tmpprs-273.16))/(tmpprs-29.65)))'
'define qs=(0.622*e/(850-0.378*e))'
'define q=qs*rhprs/100'
'set cmin 0.014'
'define q1=ave(q,t=14,t=18)'
'd q1'
'run cbar_interp 1 1 0.4 9.9 4'
'close'
'set lon 107 112'
'set lat 17 20.5'
'set lev 850'
'set t 14'
'set gxout contour'
'set cthick 6'
'set csmooth on'
'set ccolor 1'
'set clab forced'
'set xlopts 1 6 0.16'
'set ylopts 1 6 0.16'
'set xlint 2.0'
'set ylint 2.0'
'set clopts 1 0.1 0.15'
'set cmin 12'
'define w1=sqrt(UGRDprs*UGRDprs+VGRDprs*VGRDprs)'
'define w0=ave(w1,t=14,t=18)'
'd w0'
'close'
'gxprint'
'printim D:\fnl\hai\2twtegstream.png white'
'disable gxprint'
|
|