ADAS Subroutine cxcrip
SUBROUTINE CXCRIP( MXCOLL , MXTEMP , IZT , IZC , WI ,
& EI , WJ , EJ , EM , PHI ,
& NCOLL , EPS , OMEG , NTEMP , TVA ,
& RAT , QI , QJ , GA
& )
C
C
C-----------------------------------------------------------------------
C
C ****************** FORTRAN77 SUBROUTINE: CXCRIP *********************
C
C PURPOSE: CALCULATES ELECTRON AND POSITIVE ION COLLISION EXCITATION
C AND DEEXCITATION RATE COEFFICIENTS FOR DIPOLE TRANSITIONS
C IN THE IMPACT PARAMETER APPROXIMATION.
C
C (BURGESS AND SUMMERS,1976,MON.NOT.R.AST.SOC.,174,345)
C
C OPTIONALLY A SET OF INCIDENT PARTICLE ENERGIES AND
C COLLISION STRENGTHS MAY BE PROVIDED, IN WHICH CASE THE
C IMPACT PARAMETER THEORY IS USED TO CALCULATE THE COLLISION
C STRENGTHS AT HIGH ENERGY WITH VALUES SCALED TO THE HIGHEST
C ENERGY INPUT COLLISION STRENGTH.
C
C CALLING PROGRAM: CXCRDG
C
C INPUT : (I*4) MXCOLL = MAXIMUM NUMBER OF COLLISION STRENGTHS.
C INPUT : (I*4) MXTEMP = MAXIMUM NUMBER OF TEMPERATURES.
C INPUT : (I*4) IZT = TARGET ION CHARGE.
C INPUT : (I*4) IZC = CHARGE OF COLLIDING PARTICLE.
C INPUT : (R*8) WI = STATISTICAL WEIGHT OF STATE I.
C INPUT : (R*8) EI = BINDING ENERGY OF STATE I.
C UNITS: RYD
C INPUT : (R*8) WJ = STATISTICAL WEIGHT OF STATE J.
C INPUT : (R*8) EJ = BINDING ENERGY OF STATE J.
C UNITS: RYD
C INPUT : (R*8) EM = REDUCED MASS FOR COLLIDING PARTICLE.
C UNITS: ELECTRON MASSES
C INPUT : (R*8) PHI = FIJ/EIJ WHERE:
C FIJ = ABSORPTION OSCILLATOR STRENGTH;
C EIJ = EI-EJ = THE TRANSITION ENERGY (RYD).
C INPUT : (I*4) NCOLL = NUMBER OF TABULAR VALUES OF COLLISION
C STRENGTH.
C INPUT : (R*8) EPS() = INCIDENT ELECTRON ENERGIES.
C UNITS: RYD
C DIMENSION: COLLISION INDEX.
C INPUT : (R*8) OMEG() = COLLISION STRENGTHS.
C DIMENSION: COLLISION INDEX.
C INPUT : (I*4) NTEMP = NUMBER OF TEMPERATURES.
C INPUT : (R*8) TVA() = TEMPERATURES (INCIDENT PARTICLE
C DISTRIBUTION).
C UNITS: EV.
C DIMENSION: TEMPERATURE INDEX.
C
C OUTPUT: (R*8) RAT = RATIO OF OMEG(NCOLL) TO I.P. OMEGA.
C OUTPUT: (R*8) QI() = COLLISIONAL EXCITATION RATE COEFFICIENTS.
C UNITS: CM3 SEC-1
C DIMENSION: TEMPERATURE INDEX.
C OUTPUT: (R*8) QJ() = COLLISIONAL DEEXCITATION RATE COEFFICIENTS.
C UNITS: CM3 SEC-1
C DIMENSION: TEMPERATURE INDEX.
C OUTPUT: (R*8) GA() = GAMMA RATE PARAMETERS.
C UNITS:
C DIMENSION: TEMPERATURE INDEX.
C
C (I*4) I = LOOP INDEX.
C (I*4) K = LOOP INDEX.
C
C (R*8) Z1 = IZT+1.
C (R*8) Z2 = ABS( IZC ).
C (R*8) EIJ = TRANSITION ENERGY.
C UNITS: RYD
C (R*8) ELAM =
C (R*8) ZCOL =
C (R*8) Z2PHI =
C (R*8) SC =
C (R*8) T1 =
C (R*8) T2 =
C (R*8) EPSM =
C (R*8) R =
C (R*8) EIQ =
C (R*8) FLAG =
C (R*8) ATE =
C (R*8) S =
C (R*8) U1 =
C (R*8) G1 =
C (R*8) U2 =
C (R*8) G2 =
C (R*8) B =
C (R*8) S1 =
C (R*8) E =
C (R*8) P =
C
C ROUTINES:
C ROUTINE SOURCE BRIEF DESCRIPTION
C -------------------------------------------------------------
C CXEIQP ADAS
C
C AUTHOR: JONATHAN NASH (TESSELLA SUPPORT SERVICES PLC)
C K1/0/81
C JET EXT. 5183
C
C DATE: 08/10/93
C
C VERSION: 1.1 DATE: 20-06-95
C MODIFIED: TIM HAMMOND (Probably)
C - PUT UNDER S.C.C.S. CONTROL
C
C VERSION : 1.2 DATE: 17-05-07
C MODIFIED : Allan Whiteford
C - Updated comments as part of subroutine documentation
C procedure.
C-----------------------------------------------------------------------
C
C-----------------------------------------------------------------------
INTEGER IZC, IZT, MXCOLL, MXTEMP
INTEGER NCOLL, NTEMP
REAL*8 EI, EJ, EM
REAL*8 EPS(MXCOLL), GA(MXTEMP), OMEG(MXCOLL)
REAL*8 PHI, QI(MXTEMP), QJ(MXTEMP), RAT
REAL*8 TVA(MXTEMP), WI, WJ