论文需要fortran出数据，问题不知道怎么解决，求助...

佩佩123

本人是学生党，写论文需要出数据，刚开始接触fortran ，程序显示没有错误，但是运行到最后，就出现了这样。求助各位大神。

程序本身是这样的
!变荷载和温度下多元复合地基

implicit none
real rw,re,rs,H,Th1,Ew,Es        !尺寸参数
real a,b,T
real kh,kw,ks              !渗透性参数
real s,n,Y,C1,C2,Th                                     !压缩模量
real F,F1,F2,F3,M
real Fa,Tm1,Tm2,Tm3,Tm
real Pi                       !定义圆周率pai
real U                        !平均固结度
integer mm,mmE                !mm--循环变量，mmE--循环中止数
integer NN                    !为ut累积控制变量
real temp,utemp
!定义数组
real ut
dimension ut(0:100)

!打开数据文件
open(16,file='Input.dat',status='old')   
!读入数据
read(16,*) rw,re,rs,H,Th1
read(16,*) a,b,T
read(16,*) kh,kw,ks
read(16,*) Ew,Es

write(*,*) rw, re, rs,H,Th1
write(*,*) a,b,T
write(*,*) kh,kw,ks
write(*,*) Ew,Es
pause

!参数计算部分
Pi=3.141592654
n=re/rw
s=rs/rw
Y=Ew/Es
F1=(log(n/s)+(kh/ks)*log(s)-0.75)*(n*n/(n*n-1))
F2=(1-kh/ks)*(1-(s*s)/(4*n*n))*(s*s)/(n*n-1)
F3=(kh/ks)*(1/(n*n-1))*(1-1/(4*n*n))
F=F1+F2+F3


Th=0.0001

10 IF (Th.LE.100) THEN
      temp=0.0
          utemp=0.0
          U=0.0
          mmE=100

  IF (Th.LE.Th1) THEN
          DO mm=1,mmE
              ! temp=0.0
                                          
                   M=(2.0*mm-1.0)*Pi/2.0
                   C1=1/(a*T+b)
                   C2=(a*T+b)
Tm1=(4*kw*(n*n-1+Y)*C2)/(kh*n*n*(n*n-1))
Tm2=(H*H*C2)/(M*M)
Tm3=(F*re*re*kw)/(2.0*n*n*kh)
Tm=Tm1/(Tm2+Tm3)                !就是tao m

temp=1.0/(Tm*Th1)*2.0/(M*M)*(1-exp(-Tm*Th))
utemp=utemp+temp
          ENDDO
U=Th/Th1-Utemp       !原

         ELSEIF (Th.GE.Th1)        THEN
            DO mm=1,mmE
                   M=(2.0*mm-1.0)*Pi/2.0
Tm1=(4*kw*(n*n-1+Y)*C2)/(kh*n*n*(n*n-1))
Tm2=(H*H*C2)/(M*M)
Tm3=(F*re*re*kw)/(2.0*n*n*kh)
Tm=Tm1/(Tm2+Tm3)                !就是tao m
         temp=1.0/(Tm*Th1)*(2.0/(M*M))*(exp(Tm*(Th1-Th))-exp(-Tm*Th))
                 Utemp=Utemp+temp
ENDDO
               U=1.0-Utemp
ENDIF      
          !打开输出文件
          open(16,file='Output.dat',status='Old')
          write(16,*) Th,U

          IF (Th.LT.0.001) THEN                                       
                        Th=Th+0.0001
                        goto 10                                 
          ELSEIF(Th.LT.0.01) THEN
             Th=Th+0.001
                        goto 10         
          ELSEIF(Th.LT.0.1) THEN
                        Th=Th+0.01
                        goto 10         
          ELSEIF(Th.LT.1.0) THEN
                        Th=Th+0.1
                        goto 10       

          ELSEIF(Th.LT.10) THEN
                        Th=Th+1
                        goto 10          
          ELSE
                    Th=Th+10
          ENDIF
          GOTO 10
   ENDIF

END     ! 程序结束符号