计算1957年的通过0摄氏度的活动积温，能不能看看哪里有问题

leilei3039

Program Main   
real Dat(1000,4),Tem(1000)
Integer yue1 , yue2 , ri1 , ri2 , ts , nian
Integer Day,Sum,i,j,k,n,m,n1,n2
Parameter (T=0)
n2 = 0

open(12,file='1.txt')
Open(300,file='D:\FORTRAN\Result.txt')
Write(300,*) '温度大于0℃积温'
Write(300,*) '年份  初月 日 终月 日  持续时间 活动积温'
112 read(12,*,end=113)
n2=n2+1
goto 112
113  rewind(12)
print*,'共',n2,'行'
Read(12,*) ((Dat(i,j),j=1,4),i=1,n2)
nian = Dat(1,1)
   !初日
   Do m=1,365
     if (Dat(m,4)>=T) Exit
   EndDo
   Do k=m,180
    Tem(k)=(Dat(k,4)+Dat(k+1,4)+Dat(k+2,4)+Dat(k+3,4)+Dat(k+4,4))/5
   Enddo
   Do k=180,m,-1
     If (Tem(k)<=T) Exit
   EndDo
   Do m=k,k+4
     If (Dat(m,4)>=T) Then
        yue1 = Dat(m,2)
        ri1  = Dat(m,3)
     Exit
     Else if ( (m==k+4).and.( Dat(m,4)<T ) ) Then
        yue1 = Dat(m+5,2)
        ri1 = Dat(m+5,3)
      Exit
     EndIf
   EndDo
  !终日
   Do k=180,361
    Tem(k)=(Dat(k,4)+Dat(k+1,4)+Dat(k+2,4)+Dat(k+3,4)+Dat(k+4,4))/5
   EndDo
   Do k=180,361
     If (Tem(k)<=T) Exit
   Enddo
   Do n=K,K+4
     If (Dat(n,4)<T) Then
         yue2 = Dat(n-1,2)
         ri2 = Dat(n-1,3)
     Exit
     EndIf
   Enddo
   !持续天数
   ts = n-m
   !积温
   Sum=0
   Do k=m,n
         If(Dat(k,4)>=T) Then
           Sum=Sum+Dat(k,4)
    EndIf
   Enddo
Write(300,*) nian , yue1 , ri1, yue2 , ri2 , ts , Sum
Close(12)
Close(300)
End

wml

请问楼主解决了吗？本人最近在计算活动积温