c     ccccccccccccccccccccccccccccccccccccccccccccccccccccc
c     c     a quadratic interpolation                     c
c     c            main.f                                 c
c     c                             Copyright  Naoki Itoh c  
c     ccccccccccccccccccccccccccccccccccccccccccccccccccccc   

      
      real*8  y(15000), qpl(15000),qmi(15000) ,t,dam,ta,tb
      integer bb,cc,i    

cccccccccccccccccc  parameter  cccccccccccccccccccccccccccc      
      

            
       open(unit=1,file='init.par',status='unknown')  

   10  format(10x,e6.3)  

       read(1,10) ta
       read(1,10) tb
       read(1,10) dam
       close(1)
       write(*,*)'log T log rho/mue   Qpair(+)   Qpair(-)'       
       
       t=ta
       

       do 1000 i=1,15000 


       call tinit(   t,   y,   qpl,   qmi, bb, cc   )

       call spline( t,   y,   qpl,   qmi,pl,bb,cc,dam)
       
        t= t + 0.05d0 
       
       if (sngl(t).gt.sngl(tb)) goto 1200

 1000  continue
 1200  continue
        end