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