求助！grads提取ERA-40中slp.nc文件指定年份月份区域并转为grd

may592331632

各位大神，我用的ERA-40nc资料，想提取slp在1961-2005年6-8月指定区域的资料并转化成grd文件，这是我的程序，图片是运行时出的问题，求各路大神帮帮忙看下怎么改、、
'reinit'
i=1961
while(i<=2005)
'set gxout fwrite'
'set fwrite D:\NHMM\ERA40-2.5\surface\slp.grd'
'sdfopen D:\NHMM\ERA40-2.5\surface\slp.nc'
j=i-1960
if（j=4|j=8|j=12|j=16|j=20|j=24|j=28|j=32|j=36|j=40|j=44)
*这里是想判断闰年
  to=153
else
  to=152
endif

   k=to+91
   while(to<=k)
     'set t 'to
        'set lev 0 0'
        'set lon 105 120'
        'set lat 12.5 35'
*该经纬度为我所需要的区域
        'd msl'
      to=to+1
   endwhile
'disablefwrite'
i=i+1
endwhile