ADAS Subroutine aphotdw
SUBROUTINE APHOTDW(B,B1,V,N,L,L1,LP,ISP,LT,LT1,IS,PREC,PION, &PSTIM,IRESOL,ndgnt,gaunt,energy,ie) IMPLICIT REAL*8(A-H,O-Z) C --------------------------------------------------------------------- C C VERSION OF PHOTDW FOR USE BY ADASRRC WHICH USES AGIIDW C ******************* H.P.SUMMERS, JET 30 JUNE 1992 *************** C C PURPOSE: CALCULATE PHOTO INTEGRALS USING GIIDW BOUND-FREE C GAUNT-FACTORS C C SAME AS RECOM.FORT(PHOTO5) BUT CALLS GIIDW C ******************* H.P.SUMMERS, JET 19 AUG. 1984**************** C INPUT C B=1.5789D5*Z**2/(V**2*TE) C B1=1.5789D5*Z**2/(V**2*TR) C WHERE C TE=ELECTRON TEMPERATURE (K) C TR=RADIATION TEMPERATURE (K) C Z=BOUND STATE ION CHARGE +1 C (THUS Z**2/V***2 IS THE IONISATION POTENTIAL (RYD)) C V=EFFECTIVE PRINCIPAL QUANTUM NUMBER OF BOUND ELECTRON C N=PRINCIPAL QUANTUM NUMBER OF BOUND ELECTRON C L=ORBITAL QUANTUM NUMBER OF BOUND ELECTRON C L1=ORBITAL QUANTUM NUMBER OF FREE ELECTRON C ISP=2*SP+1 WHERE SP IS TOTAL SPIN OF PARENT STATE C LP=TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM NUMBER OF PARENT STATE C LT=TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM NUMBER OF BOUND SYSTEM C LT1=TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM NUMBER OF FREE SYSTEM C IS=2*S+1 WHERE S IS TOTAL SPIN OF SYSTEM C ndgnt = max number of Gaunt factors allowed C OUTPUT C PREC=RADIATIVE RECOMBINATION INTEGRAL C PION=PHOTOIONISATION INTEGRAL C PSTIM=STIMULATED RECOMBINATION INTEGRAL C WHERE C IRESOL=1 FOR ((LP,SP)N L LT S,(LP,SP)L1 LT1 S) C =2 FOR ((LP,SP)N L LT S,(LP,SP)L1 S) =ABOVE LT1 SUM C =3 FOR ((LP,SP)N L S,(LP,SP)L1 S) = ABOVE LT SUM C =4 FOR ((LP,SP)N L,(LP,SP)L1) = ABOVE S SUM C =5 FOR NO L RESOLUTION USING GBF C gaunt() = Bound-free Gaunt factor at energy C energy() = v**2*e C where e = (free electron energy)/z**2 (ryd) C v = effective principal quantum number C of bound electron C ie = number of Gaunt/energy pairs C C C UNIX-IDL PORT: C C AUTHOR: WILLIAM OSBORN (TESSELLA SUPPORT SERVICES PLC) C C DATE: 4TH JULY 1996 C C VERSION: 1.1 DATE: 04-07-96 C MODIFIED: WILLIAM OSBORN C - FIRST VERSION. C C VERSION: 1.2 DATE: 19-12-01 C MODIFIED: Martin O'MULLANE C - Removed junk from > column 72. C C VERSION: 1.3 DATE: 02-02-05 C MODIFIED: Paul Bryans C - Returns Gaunt factor, associated vve and number C of Gaunt/vve pairs C C VERSION: 1.4 DATE: 16-05-07 C MODIFIED: Allan Whiteford C - Modified comments as part of subroutine documentation C procedure. C CC ---------------------------------------------------------------------- C INTEGER IE, IRESOL, IS, ISP INTEGER L, L1, LP, LT INTEGER LT1, N, NDGNT REAL*8 B, B1, ENERGY(NDGNT) REAL*8 GAUNT(NDGNT), PION, PREC REAL*8 PSTIM, V