相关系数求取问题

hpa666

这是我写的两个subroutine
!线性趋势检验
subroutine ls(n,x,y,yz)
integer::n,i
real::x(n),y(n),yx(n),yz(n),b2,b1,s1,s2,s3,s4,q,u,f,as3
s1=0.0
s2=0.0
s3=0.0
s4=0.0
u=0.0
q=0.0
do i=1,n
s1=s1+x(i)*y(i)
s2=s2+x(i)
s3=s3+y(i)
s4=s4+x(i)**2
enddo
b2=(n*s1-s2*s3)/(n*s4-s2**2)
b1=(s3-b2*s2)/n
do i=1,n
yx(i)=b1+b2*x(i)
yz(i)=y(i)-yx(i)
enddo

as3=s3/n
   do i=1,n
  q=q+(y(i)-yx(i))**2
  u=u+(yx(i)-as3)**2
enddo
f=(u/1)/(q/(n-1-1))
print*,'podu',b2
print*,'jianyan f',f
end subroutine

!相关系数
subroutine r2(n,x,y,corr)
implicit none
integer::n,i
real::x(n),y(n),s1,s2,s3,s4,s5,corr
s1=0.0
s2=0.0
s3=0.0
s4=0.0
s5=0.0
do i=1,n
s1=s1+x(i)*y(i)
s2=s2+x(i)
s3=s3+y(i)
s4=s4+x(i)**2
s5=s5+y(i)**2
enddo
corr=(n*s1-s2*s3)/sqrt((n*s4-s2**2)*(n*s5-s3**2))
print *, corr
end subroutine
end module

这是主程序的代码
  call indexit(ymd(1,1),12,31)
  if(v<365)then
  call ls((soy-1),ytemminave(2:soy),soya(2:soy),yz)

  call r2((soy-1),soya(2:soy),ytemminave(2:soy),corr)

   else
   call ls(soy,ytemminave,soya,yz)

   call r2(soy,soya,ytemminave,corr)

  end if

然后就是我这个筛选的意思就是：如果这个数据中的第一年是一整年的话我就把这个年份数组从头至尾都要算上，如果第一年不是一整年，我就要从第二年开始算起
问题就出在标红的这一行，编译的时候没有问题，但是运行时会出现这个错误，如果去掉就会正确，加上就会错误
也出现了有时能得出结果，有时不能得出结果的情况