c      ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c      cc     OutPut to Files                                      ccc
c      cc                   outpt.f                                ccc
c      ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
       subroutine outpt( qpai, temp, rhom )
c
        implicit real*8 ( a-h, o-z )
        dimension rhom( 1 ), temp( 1 ), qpai( 1 )
        common /param/ ndat, n
c
        open ( unit = 80, file = 'pair.dat' , status = 'unknown'
     &       , access = 'sequential', form = 'formatted'        )
        write (80,1100) 
        do 100 i = 1 , ndat
          t = dlog10( temp( i ) )
          r = dlog10( rhom( i ) )
c         -----------------------
          if (qpai( i ) .gt. 0 )
     &       then
               q = dlog10( qpai( i ) )
             else
               q = 9.9999999d-99
          endif
c         -----------------------
          write (80,1000) t, r, qpai( i ), q 
 100    continue
        return
 1000     format(3x,1p4e17.8)
 1100   format(/,3x,'    Log(T)       ','     Rho/Mu_e    ',
     &        '     Q_pair      ')
        end