SpectraRust/tlusty/extracted/dietot.f
2026-03-19 14:05:33 +08:00

29 lines
743 B
Fortran

SUBROUTINE DIETOT
C =================
C
C modification of the photoionization cross-section
C for taking into account dielectronic recombination
C
INCLUDE 'IMPLIC.FOR'
INCLUDE 'BASICS.FOR'
INCLUDE 'ATOMIC.FOR'
INCLUDE 'MODELQ.FOR'
C
do ion=1,nion
i=nfirst(ion)
ia=numat(iatm(i))
io=iz(ion)
do id=1,nd
t=temp(id)
xpx=dens(id)/wmm(id)/ytot(id)
call dielrc(ia,io,t,xpx,dirt,sig0)
diesig(ion,id)=sig0
if(id.eq.1.or.id.eq.35.or.id.eq.nd) then
write(99,699) ion,ia,io,id,i,nnext(ion),dirt,sig0
end if
end do
end do
699 format(6i5,1p2e12.4)
return
end