NCL两次计算同一个程序结果不同

dannybear

请教一下，各位
  我用NCL计算同一个程序，为什么计算出来的结果会不同，我第一运行的是一个二元回归的脚本，接着，我运行了一元回归的脚本，我发现两个的结果竟然是一样的，用的都是相同的变量，于是我把两个脚本写在了一起，就是在一个脚本中包含了一元回归，也包含了2元回归，最后计算出来的一元回归的结果和二元回归的结果就不一样了，请问大家这个是怎么回事
我把脚本贴出来：
if(.not.ismissing(to1(lat,lon)))  then
           Y11 = yb(:)
         
           X11  = new ( (/3,nn/), "float" )  
       X11(0,:) = 1.0                           ; constant term on rhs
       X11(1,:) = to1(lat,lon)
          X11(2,:) = x1b(:)
           X11@_FillValue=9.99e+08
           beta1 = reg_multlin (Y11,X11,False)
           y2b(shix,iii,lat,lon) = beta1(0)+beta1(1)*to1(lat,lon)+beta1(2)*x10  ；；这是二元回归

           delete(X11)
           delete(Y11)
           delete(beta1)


            rc1   = regline (yb,x1b)

                temp=rc1@yintercept+rc1*x10  ；；这是一元回归
               
          y2b1(shix,iii,lat,lon) = temp
          
          
            delete(rc1)       

         else if(ismissing(to1(lat,lon)))  then

    mo(shix,iii,lat,lon)= mo@_FillValue
        obs(shix,iii,lat,lon)= obs@_FillValue
         
           end if
       end if          


当我把一元回归和二元回归的脚本分开写的时候，计算出来的结果一模一样，请问这个是为什么？

藕茭白空心菜

  rc1   = regline (yb,x1b)
temp=rc1@yintercept+rc1*x10  ；
和
rc = regline_stats(year, xavg)
一样嘛
请问~