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