ADAS Subroutine cxphot
SUBROUTINE CXPHOT( IZ1 , TE , TR , V ,
& N , L , L1 , LP ,
& ISP , LT , LT1 , IS ,
& IRES , PREC , PION , PSTIM
& )
C
C
C-----------------------------------------------------------------------
C
C ****************** FORTRAN77 SUBROUTINE: CXPHOT *********************
C
C PURPOSE: CALCULATES PHOTO INTEGRALS USING GIIH BOUND-FREE
C GAUNT-FACTORS.
C
C CALLING PROGRAM: C6TBRC.
C
C INPUT : (I*4) IZ1 = ION CHARGE.
C INPUT : (R*8) TE = ELECTRON TEMPERATURE.
C UNITS: K
C INPUT : (R*8) TR = RADIATION TEMPERATURE.
C UNITS: K
C INPUT : (R*8) V = EFFECTIVE PRINCIPAL QUANTUM NUMBER OF BOUND
C ELECTRON.
C INPUT : (I*4) N = PRINCIPAL QUANTUM NUMBER OF BOUND ELECTRON.
C INPUT : (I*4) L = ORBITAL QUANTUM NUMBER OF BOUND ELECTRON.
C INPUT : (I*4) L1 = ORBITAL QUANTUM NUMBER OF FREE ELECTRON.
C INPUT : (I*4) LP = TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM
C NUMBER OF PARENT STATE.
C INPUT : (I*4) ISP = 2*SP+1 WHERE SP IS TOTAL SPIN OF PARENT
C STATE.
C INPUT : (I*4) LT = TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM
C NUMBER OF BOUND SYSTEM.
C INPUT : (I*4) LT1 = TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM
C NUMBER OF FREE SYSTEM.
C INPUT : (I*4) IS = 2*S+1 WHERE S IS TOTAL SPIN OF SYSTEM.
C INPUT : (I*4) IRES = LEVEL OF RESOLUTION.
C = 1 :
C = 2 : ABOVE LT1 SUM.
C = 3 : ABOVE LT SUM.
C = 4 : ABOVE S SUM.
C = 5 : UNRESOLVED GBF.
C
C OUTPUT: (R*8) PREC = RADIATIVE RECOMBINATION INTEGRAL.
C OUTPUT: (R*8) PION = PHOTOIONISATION INTEGRAL.
C OUTPUT: (R*8) PSTIM = STIMULATED RECOMBINATION INTEGRAL.
C
C
C PARAM : (I*4) MXT = 8.
C PARAM : (I*4) MXI = 4.
C PARAM : (R*8) P1 =
C PARAM : (R*8) P2 = 0.9
C
C (I*4) IB = FLAGS VALUE OF B.
C = 1 : B >= P2.
C = 2 : B < P2.
C (I*4) IB1 = FLAGS VALUE OF B1.
C = 1 : B1 >= P2.
C = 2 : B1 < P2.
C (I*4) I = LOOP INDEX.
C
C (R*8) Z1 = REAL VALUE = IZ1.
C (R*8) B = P1 * Z1**2 / ( V**2 * TE )
C (R*8) B1 = P1 * Z1**2 / ( V**2 * TR )
C (R*8) VVE =
C (R*8) GI = R8GIIH().
C (R*8) GI1 = R8GIIH().
C (R*8) GI2 = R8GIIH().
C (R*8) F1 =
C (R*8) F2 =
C (R*8) F3 =
C (R*8) F4 =
C (R*8) F5 =
C (R*8) F6 =
C (R*8) T =
C (R*8) T2 =
C (R*8) T4 =
C (R*8) T6 =
C (R*8) Q1 =
C (R*8) Q2 =
C (R*8) EX =
C (R*8) EX1 =
C (R*8) D1 =
C (R*8) D2 =
C (R*8) U2 =
C (R*8) U6 =
C
C (R*8) X1() =
C DIMENSION: MXT.
C (R*8) X2() =
C DIMENSION: MXT.
C (R*8) W1() =
C DIMENSION: MXT.
C (R*8) W2() =
C DIMENSION: MXT.
C
C ROUTINES:
C ROUTINE SOURCE BRIEF DESCRIPTION
C -------------------------------------------------------------
C R8GIIH ADAS RETURNS BOUND-FREE G-FACTORS.
C R8FEII ADAS
C
C AUTHOR: JONATHAN NASH (TESSELLA SUPPORT SERVICES PLC)
C K1/0/81
C JET EXT. 5183
C
C DATE: 05/11/93
C
C UNIX-IDL PORT:
C
C AUTHOR: WILLIAM OSBORN (TESSELLA SUPPORT SERVICES PLC)
C
C DATE: 22ND MAY 1996
C
C VERSION: 1.1 DATE: 22-05-96
C MODIFIED: WILLIAM OSBORN
C - FIRST VERSION. IBM VERSION NOT CHANGED
C
C-----------------------------------------------------------------------
C
C-----------------------------------------------------------------------
INTEGER IRES, IS, ISP, IZ1
INTEGER L, L1, LP, LT
INTEGER LT1, N
REAL*8 PION, PREC, PSTIM, TE
REAL*8 TR, V