大家帮我看看这个计算年平均水汽通量散度的程序哪里错了？

龙行天下

最近需要计算全国500hPa高度1962年~1991年年平均水汽通量散度，但是目前所有的这方面程序都是Grads的，自己对Grads刚刚入门，在论坛上找的代码，我下载了NECP的uwnd.mon.mean.nc，vwnd.mon.mean.nc，rhum.mon.mean.nc，air.mon.mean.nc数据，按照下载的代码修改了一下，但是一直运行不成功，请大家帮我指导一下，或者哪位好心人有Matlab代码也行的，不胜感激！'sdfopen E:\数据\NECP数据\uwnd.mon.mean.nc'
'sdfopen D:\m_d\vwnd.mon.mean.nc'
'sdfopen D:\m_d\rhum.mon.mean.nc'
'sdfopen D:\m_d\air.mon.mean.nc'
'set lon 70 140',
'set lat 15 55'
'set lev 500'
'set mpdset cnworld'
i=257

while(i<=528)
'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 E:\数据\NECP数据\500hPa.png white'
'c'

i=i+1
endwhile
;

下面是运行中的报错，

river

首先不说这个gs能不能达到你的目的，先说你说从论坛上找的代码，但是有水分享的代码里有中文路径了？你给举个例子最好贴个地址，我去学习一下。第二，即使是找的代码也应该有一定的基础才能去修改，才能最终变成自己的。但是我看你这个改的也不多，却出现好几个不应该有的问题，第一个是有一句多了一个 ,  ，第二个是居然出现了中文下面的分号。具体在哪自己去找。第三，对grads刚入门就应该先从简单的来，找一些论坛里的教程，把上面的例子先练习一下，搞懂最基本的东西，那修改起复杂的脚本也就能驾轻就熟了

龙行天下

river 发表于 2013-11-30 21:26

                               
登录/注册后可看大图

首先不说这个gs能不能达到你的目的，先说你说从论坛上找的代码，但是有水分享的代码里有中文路径了？你给举 ...

受教了，必须打好基础

xiaohezhou

来看看