55 lines
1.3 KiB
Fortran
55 lines
1.3 KiB
Fortran
SUBROUTINE PRDINI
|
|
c =================
|
|
c
|
|
c initialization of PRD
|
|
c
|
|
INCLUDE 'IMPLIC.FOR'
|
|
INCLUDE 'BASICS.FOR'
|
|
INCLUDE 'ATOMIC.FOR'
|
|
INCLUDE 'MODELQ.FOR'
|
|
c
|
|
ntrprd=0
|
|
do itr=1,ntrans
|
|
iprd(itr)=0
|
|
if(ifprd.gt.0.and.line(itr).and.indexp(itr).ne.0) then
|
|
ii=ilow(itr)
|
|
jj=iup(itr)
|
|
iat=iatm(ii)
|
|
c
|
|
c select Lyman alpha for PRD
|
|
c
|
|
if(iat.eq.iath.and.ii.eq.nfirst(ielh).and.
|
|
* fr0(itr).lt.2.5e15) then
|
|
ntrprd=ntrprd+1
|
|
iprd(itr)=ntrprd
|
|
itrtot(ntrprd)=itr
|
|
end if
|
|
c
|
|
c select Mg I resonance line for PRD
|
|
c
|
|
if(numat(iat).eq.12.and.iz(iel(ii)).eq.1.and.
|
|
* ii.eq.nfirst(iel(ii)).and.fr0(itr).lt.1.06e15) then
|
|
ntrprd=ntrprd+1
|
|
iprd(itr)=ntrprd
|
|
itrtot(ntrprd)=itr
|
|
end if
|
|
c
|
|
c select Mg II resonance lines for PRD
|
|
c
|
|
if(numat(iat).eq.12.and.iz(iel(ii)).eq.2.and.
|
|
* ii.eq.nfirst(iel(ii)).and.fr0(itr).lt.1.08e15) then
|
|
ntrprd=ntrprd+1
|
|
iprd(itr)=ntrprd
|
|
itrtot(ntrprd)=itr
|
|
end if
|
|
end if
|
|
end do
|
|
c
|
|
do itrp=1,ntrprd
|
|
do id=1,nd
|
|
pjbar(itrp,id)=0.
|
|
end do
|
|
end do
|
|
return
|
|
end
|