cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c                                                              c
c         HAFT,RAFFELT,& WEISS                                 c 
c                 Plasma process fitting formulae              c 
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c
       implicit real*8 (a-z)       
      
            
cccccccccccccccccccc parameter ccccccccccccccccccccccccccccccccc
 
      open (unit=80,file='init.par',status='unknown') 
      read(80,10) tinit
      read(80,10) tfin
      read(80,10) tstep
      read(80,10) rinit
      read(80,10) rfin
      read(80,10) rstep
      read(80,10) n

      close(80)
 10   format(8x,f10.4)

cccccccccccccccccccccc  main  cccccccccccccccccccccccccccccccc
      
      

      do 2000 t=tinit,tfin,tstep
      do 1000 r=rinit,rfin,rstep
c     
        
      call formulae ( t, r, qn, q, n ) 
 100  format(3x,1p4e17.8)

      open (unit=80,file='plasma.dat',status='unknown')
 
      write(80,100) t, r, qn

 1000  continue
 2000 continue
 
      end