用gs画螺旋度，总是提示错误non-numeric args to numeric operation,求大神帮助看看

yuejun0026

我用在论坛上下的画螺旋度的gs文件，可是一直提示出现问题undefine warning:name not found;non-numeric args to numeric operation;error occurred on line 51;in file f:\x\luoxuandu1.gs。知道这个上面说的是错误在51行，不过我看不出来51行有什么错误啊，求大神帮忙看看。下面是我的gs文件：
*水平螺旋度的gs文件:*  
   'reinit'
   'clear'
    prefix.1='u' ; prefix.2='v'
    pi=3.1415926
   'open f:\x\fnl.ctl'
    levstr=' 850 700 600 500 400 300 '
   'set lon 80 110'
   'set lat 30 50'
   'set time 00Z15jul2010'
   'define cr=ugrdprs'
**************************************************************************************************
    i=1
    while(i<=9)
    j=1
    while(j<=13)
    uave.i.j=0.0
    vave.i.j=0.0
    k=1
    while(k<=5)
   'set lev 'subwrd(levstr,k)
   'undefine fld'
   'define fld=ugrdprs.1'
   'q defval fld '32+j' '48+i
    nwval=subwrd(result,3)
    uave.i.j=uave.i.j+nwval           这就是51行的位置
   'undefine fld'
   'define fld=vgrdprs.2'
   'q defval fld '32+j' '48+i
    nwval=subwrd(result,3)
    vave.i.j=vave.i.j+nwval
*   say ' nwval='nwval ' vave.'i'.'j'='vave.i.j
    k=k+1
    endwhile
    uave.i.j=uave.i.j/5
    vave.i.j=vave.i.j/5
    if(vave.i.j!=0)
    qq.i.j=math_atan(math_abs(vave.i.j/uave.i.j))
    else
    if(uave.i.j>0)
    qq.i.j=0
    endif
    if(uave.i.j<0)
    qq.i.j=pi
    endif
    if(uave.i.j=0)
    say ' chaos has arrised!'
    endif
    endif
    if(uave.i.j>=0 & vave.i.j>=0);qq.i.j=qq.i.j-40.0    ;endif;
    if(uave.i.j<0 & vave.i.j>=0) ;qq.i.j=180-40.0-qq.i.j;endif;
    if(uave.i.j<=0 & vave.i.j<=0);qq.i.j=180-40.0+qq.i.j;endif;
    if(uave.i.j>0 & vave.i.j<=0) ;qq.i.j=270-40.0+qq.i.j;endif;
*    if(uave.i.j>=0 & vave.i.j>=0);qq.i.j=qq.i.j-40.0*3.1415926/180.0      ;endif;
*    if(uave.i.j<0 & vave.i.j>=0) ;qq.i.j=(180-40.0)*3.1415926/180.0-qq.i.j;endif;
*    if(uave.i.j<=0 & vave.i.j<=0) ;qq.i.j=(180-40.0)*3.1415926/180.0+qq.i.j;endif;
*    if(uave.i.j>0 & vave.i.j<=0)  ;qq.i.j=(270-40.0)*3.1415926/180.0+qq.i.j;endif;
    uave.i.j=math_sqrt(uave.i.j*uave.i.j+vave.i.j*vave.i.j)*math_cos(qq.i.j*3.1415926/180.0)*0.75
    vave.i.j=math_sqrt(uave.i.j*uave.i.j+vave.i.j*vave.i.j)*math_sin(qq.i.j*3.1415926/180.0)*0.75
*    uave.i.j=math_sqrt(uave.i.j*uave.i.j+vave.i.j*vave.i.j)*math_cos(qq.i.j)*0.75
*    vave.i.j=math_sqrt(uave.i.j*uave.i.j+vave.i.j*vave.i.j)*math_sin(qq.i.j)*0.75

    j=j+1
    endwhile
    i=i+1
    endwhile
*************************************************************************************************************************************

*************************************************************************************************************************************
    i=1
    while(i<=9)
    j=1
    while(j<=13)
    hd.i.j=0.
    k=1
    while(k<=5)
   'set lev 'subwrd(levstr,k)
   'undefine fld'
   'define fld=ugrdprs.1'
   'q defval fld '32+j' '48+i
    u0=subwrd(result,3)
   'undefine fld'
   'define fld=vgrdprs.2'
   'q defval fld '32+j' '48+i
    v0=subwrd(result,3)
   'set lev 'subwrd(levstr,k+1)
   'undefine fld'
   'define fld=ugrdprs.1'
   'q defval fld '32+j' '48+i
    u1=subwrd(result,3)
   'undefine fld'
   'define fld=vgrdprs.2'
   'q defval fld '32+j' '48+i
    v1=subwrd(result,3)
    hd.i.j=hd.i.j+(u1-uave.i.j)*(v0-vave.i.j)-(v1-vave.i.j)*(u0-uave.i.j)
    k=k+1
    endwhile
    'set defval cr '32+j' '48+i'  'hd.i.j
*  say 'set defval sucess' ; pull c
*    'd     'hd.i.j
*    'print '
    j=j+1
    endwhile
    i=i+1
    endwhile
   'enable print f:\x\luoxuandu.gmf'
   'set cint 40'
   'd cr'
   'print '
   'disable print'

*   'reinit'
*    say 'u0='u0' u1='u1
*    say 'v0='v0' v1='v1
*    say ' i='i' j='j
*    say ' hd.'i'.'j'='hd.i.j

yuejun0026

求大神帮忙看看啊！！！

yuejun0026

为啥没有大神帮忙看看啊，急需啊，毕业设计就在眼前啊！！！

yuejun0026

为什么没有人看啊？求大神指导一下下啊！！！！

爱公爵

没大神来帮忙吗

1417920146

需要根据你的格点的分辨率或尺度调整13,9这些数字，不能直接套用，除非你是2.5*2.5格距的

成都蒋妹妹

我也是相同问题，你解决了么？

yuejun0026

成都蒋妹妹 发表于 2013-12-24 15:19

                               
登录/注册后可看大图

我也是相同问题，你解决了么？

额  忘了哦  这是好久之前的了   后面一直没弄了  不好意思了

成都蒋妹妹

yuejun0026 发表于 2013-12-25 14:42

                               
登录/注册后可看大图

额  忘了哦  这是好久之前的了   后面一直没弄了  不好意思了

你周围有人搞过这个没？我还是没有解决呃。。。。55555

成都蒋妹妹

1417920146 发表于 2013-4-2 12:39

                               
登录/注册后可看大图

需要根据你的格点的分辨率或尺度调整13,9这些数字，不能直接套用，除非你是2.5*2.5格距的

你好，我也改正了，但是对于程序中那个32+j,48+i不知道怎么改，改后仍然还是超出。。。不知问题在哪里，可否赐教