ADAS Subroutine c6tbrc
SUBROUTINE C6TBRC( MXNSHL , MXJSHL , IZ1 , NBOT , & NTOP , TEV , QTHRC , FTHRCJ & ) C C C----------------------------------------------------------------------- C C ****************** FORTRAN77 SUBROUTINE: C6TBRC ********************* C C PURPOSE: SETS UP A TABLE OF RADIATIVE RECOMBINATION RATE C COEFFICIENTS FOR A BARE NUCLEUS, HELIUM-LIKE OR NEON-LIKE C ION TO EXCITED NLJ LEVELS. C C CALLING PROGRAM: ADAS306 C C INPUT : (I*4) MXNSHL = MAXIMUM VALUE OF PRINCIPAL QUANTUM NUMBER. C INPUT : (I*4) MXJSHL = MAXIMUM NUMBER OF J SUB-SHELLS. C INPUT : (I*4) IZ1 = ION CHARGE. C INPUT : (I*4) NBOT = MINIMUM PRINCIPAL QUANTUM NUMBER. C INPUT : (I*4) NTOP = MAXIMUM PRINCIPAL QUANTUM NUMBER. C INPUT : (R*8) TEV = ELECTRON TEMPERATURE. C UNITS: EV C C OUTPUT: (R*8) QTHRC() = RECOMBINATION RATE COEFFICIENT TO LEVEL N. C UNITS: CM3 SEC-1 C DIMENSION: N-SHELL C OUTPUT: (R*8) FTHRCJ(,)= FRACTION OF RECOMBINATION RATE OF LEVEL N C TO STATE NLJ. C 1ST DIMENSION: J-SHELL INDEX WHERE C 1 GIVES J=L+0.5 C 2 GIVES J=L-0.5 C 2ND DIMENSION: REFERENCED BY I4IDFL(). C C PARAM : (R*8) P1 = BOLTZMANN CONSTANT. C UNITS: EV K-1 C PARAM : (R*8) P2 = C C (I*4) N = PRINCIPAL QUANTUM NUMBER OF BOUND ELECTRON. C (I*4) L = ORBITAL QUANTUM NUMBER OF BOUND ELECTRON. C (I*4) L1 = ORBITAL QUANTUM NUMBER OF FREE ELECTRON. C (I*4) LP = TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM C NUMBER OF PARENT STATE. C (I*4) ISP = 2*SP+1 WHERE SP IS TOTAL SPIN OF PAREN T C STATE. C (I*4) LT = TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM C NUMBER OF BOUND SYSTEM. C (I*4) LT1 = TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM C NUMBER OF FREE SYSTEM. C (I*4) IS = 2*S+1 WHERE S IS TOTAL SPIN OF SYSTEM. C (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 (I*4) I = LOOP INDEX. C (I*4) J = LOOP INDEX. C (I*4) IDL = TABLE INDEX. C C (R*8) Z1 = REAL VALUE = IZ1. C (R*8) TE = ELECTRON TEMPERATURE. C UNITS: K C (R*8) V = EFFECTIVE PRINCIPAL QUANTUM NUMBER OF BOUND C ELECTRON. C (R*8) FACT = C (R*8) SUM = C (R*8) XL = REAL VALUE = L. C (R*8) WL = C (R*8) T = C (R*8) PREC1 = RADIATIVE RECOMBINATION INTEGRAL. C (R*8) PION1 = PHOTOIONISATION INTEGRAL. C (R*8) PSTIM1 = STIMULATED RECOMBINATION INTEGRAL. C (R*8) PREC2 = RADIATIVE RECOMBINATION INTEGRAL. C (R*8) PION2 = PHOTOIONISATION INTEGRAL. C (R*8) PSTIM2 = STIMULATED RECOMBINATION INTEGRAL. C C PARAM : (R*8) P1 = BOLTZMANN CONSTANT. C UNITS: EV K-1 C PARAM : (R*8) P2 = C C ROUTINES: C ROUTINE SOURCE BRIEF DESCRIPTION C ------------------------------------------------------------- C I4IDFL ADAS RETURNS UNIQUE INDEX GIVEN QUANTUM C NUMBERS N AND L. C CXPHOT ADAS CALCULATES PHOTO INTEGRALS USING GIIH C BOUND-FREE GAUNT-FACTORS. 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 VERSION: 1.2 DATE: 29-05-96 C MODIFIED: WILLIAM OSBORN C - REMOVED UNUSED VARIABLES C C----------------------------------------------------------------------- C C----------------------------------------------------------------------- INTEGER IZ1, MXJSHL, MXNSHL, NBOT INTEGER NTOP REAL*8 FTHRCJ(MXJSHL,(MXNSHL*(MXNSHL+1))/2) REAL*8 QTHRC(MXNSHL), TEV