c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc c cc output routine ccc c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc c subroutine output( rho , tem ,fityp , qb ,gamma ,nn ) implicit real*8 (a-h,o-z) integer fityp real*8 logt, logrho, logqb common /nika2/ a , z , ic dimension qb(nn) , tem(nn), rho(nn), fityp(nn),gamma(nn) c open ( unit = 50, file = 'brems.dat' , status = 'unknown' & , access = 'sequential', form = 'formatted') c write(50,1010) a , z write(50,1020) do 100 i = 1, ic logt = log10(tem(i)) logrho = log10(rho(i)) c ---------------------------------- if ( qb( i ) .le. 0.d0 ) & then logqb = -9.999999d99 else logqb = log10( qb(i) ) endif c ---------------------------------- vmue = a/z tes = logqb-logrho-log10( vmue ) write(50,1000) logt, logrho, qb( i ), & logqb, fityp( i ) 100 continue 1000 format(5x,1p2e12.3,3x,1p2e15.6,3x,i2) 1010 format(/,5x,'A = ',f3.0,' Z = ',f3.0) 1020 format(5x,' Log T ','Log rho/mu_e' & ,' Q_Bremss ',' Log Q_bremss ','Fitting type' & ) close(50) return end