- 积分
- 7756
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-10-28
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
帖子里引用的代码,做水汽通量,运行时总提示错误,画不出风场,请问是什么原因呢?谢谢!
'reinit'
'sdfopen D:\ncep\day\uwnd.2012.nc'
'sdfopen D:\ncep\day\vwnd.2012.nc'
'sdfopen D:\ncep\day\rhum.2012.nc'
'sdfopen D:\ncep\day\air.2012.nc'
'set lon 70 140'
'set lat 15 55'
'set lev 700'
i=89
while(i<=89)
'set t 'i
'set grads off'
'set grid off'
*'set gxout shaded'
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(air.4-273.15))/(air.4-29.65)))'
'define qs=(0.62197*es/(prs-0.378*es))'
'define q=qs*rhum.3/100'
'd mag(q*uwnd.1/g,q*vwnd.2/g)*1000'
'set gxout vector'
'd uwnd.1/g;vwnd.2/g;mag(uwnd.1,vwnd.2)'
*'run cbarn'
'printim D:\ncep\rh200808'i'-700hPa.png white'
'print'
*'c'
i=i+1
endwhile
'disable print'
;
错误提示如下:
ga-> run d:\ncep\sqtl.gs
No hardcopy metafile open
All files closed; all defined objects released
All GrADS attributes have been reinitialized
Syntax Error: Expected operator or ')'
Error ocurred at column 9
DISPLAY error: Invalid expression
Expression = vwnd.2/g;mag(uwnd.1,vwnd.2)
Hardcopy output file is closed
|
-
|