program temca
    implicit none
    integer n,n1,i,j
    real pi,u0
    real(8) a,i0
    ! nΪn1Ϊʱ
    parameter(n=3,n1=30,a=200,i0=1)
    real(8) p(2*n-1),u(n1)
    real t(n1)
    !call getdata(t,p,n1,n)
    pi=3.1415926
    u0=4.e-7*pi
   
   open(111,file='data.txt')

   do i=1,2*n-1
   read(111,*)p(i)
   
   enddo
   print*,'**********************************'
   print*,p
!   pause
   
   do j=1,n1
   read(111,*)t(j)
   end do
   close(111)
   print*,t
!   do i=1,n1
!       t(i)=t(i)
!   enddo  
!   print*,'**********************************'
!   print*,t
   pause
    print*,4
    !print*,p,t
    call Tem_calculate(p,n,t,u,n1,a,i0)
    print*,4
    call output(t,u,n1)
    end program temca
!********************************************