- 积分
- 19
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-1-13
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
编写运行gs文件,想将每天的四个数据相加出一张图
'reinit'
'open /data1/trx/lq/model/naqpms/drydep/data/drydepd1.2006041506.ctl'
'set lat 0 50'
'set lon 100 150'
idate='20060415 20060416 20060417 20060418 20060419 20060420 20060421 20060422 20060423 20060424 20060425 20060426 20060427 20060428 20060429 20060430 20060501 20060502 20060503 20060504 20060505 20060506'
ii=1
t = 1
while(t <= 85 )
date=subwrd(idate,ii)
'c'
'set parea 1 10.0 0.5 7.5'
'set grads off'
'set gxout shaded'
'set t 't
'd sum(veldryNH3+veldryNO+veldryNO2+veldryHONO+veldryO3+veldryOH+veldryHO2+veldryH2O2+veldryCO+veldrySO2+veldryISOP+veldryTERP+veldryPM25+veldryPM10+veldryBC+veldryOC+veldrySOA1+veldrySOA2+veldrySOA3+veldrySOA4+veldrySOA5+veldrySOA6,t='t',t='t+3')'
t1=(t+3)/4
'cbarn 1 1 9 4'
'set string 1 c 7'
'set strsiz 0.26 0.24'
'draw string 5.5 7.8 'date''
'set strsiz 0.15'
'draw string 10 0.5 ug.m-3'
'printim mod.'date'.png white'
t=t+4
ii=ii+1
endwhile
;
出现下面这些错误
No hardcopy metafile open
All files closed; all defined objects released;
All GrADS attributes have been reinitialized
Data Request Warning: Request is completely outside file limits
Cannot contour grid - all undefined values
Cannot plot color bar: No shading information
DRAW error: Syntax is DRAW STRING x y string
哪里错了呢
|
|