program resistivewall
use bmad
implicit none
real(rp) Z0/376.7/,E0,I_ave,Qv,sigma, circumference/768.428/
real(rp) b1, s1, b2, s2
real(rp) I_A/17.045E3/ !Amps
real(rp) gamma
real(rp) Rate
real(rp) constant
real(rp) omega
integer lun
namelist/input/E0,I_ave,Qv,sigma,b1,s1,b2,s2

 lun = lunget()
 open(unit=lun,file='input_params.dat')
 read(lun,nml = input)
 write(6,nml=input)
 close(unit=lun)

 gamma = E0/m_electron
 constant = (2*twopi/Z0/c_light)*(1./twopi) &
        *(c_light/gamma) * (I_ave/I_A) * (1/twopi/Qv) &
        * sqrt(Z0 * c_light**2/2/twopi/sigma) &
        * sqrt(circumference/(1-(Qv-int(Qv))))

 Rate = constant * (s1/b1**3 + s2/b2**3)

 print '(a,es12.4,a,es12.4)', 'Z0 = ', Z0,' gamma = ',gamma
 print '(a,es12.4,a)',' Rate = ', rate,' /second'
 print '(a,es12.4,a)',' Rate = ', rate * (circumference/c_light),' /turn'
 print '(a,es12.4,a)',' Growth time = ', 1/(rate*circumference/c_light),' turns'
 omega = twopi*c_light/circumference
 print '(a,es12.4)',' Omega_rev = ',omega
 print '(a,es12.4)','Skin depth at revolution frequency = ',sqrt(2*c_light/Z0/sigma/omega)
end program