ADAS Subroutine c9cxee
SUBROUTINE C9CXEE( MXNENG , MXNSHL , NGRND , NTOT ,
& NBOT , NTOP , IRZ0 , IRZ1 ,
& RAMSNO , TEV , TIEV , DENS ,
& DENSZ , ZEFF , BMAG , BMENG ,
& ITHEOR , IBSTAT , IEMMS , NTU ,
& NTL , NMINF , NMAXF , NENRGY ,
& ENRGYA , ALPHAA , XSECNA , FRACLA ,
& ERATE
& )
C
C
C-----------------------------------------------------------------------
C
C ****************** FORTRAN77 SUBROUTINE: C9CXEE *********************
C
C PURPOSE: CALCULATES THE L-RESOLVED EFFECTIVE EMISSIVITY RATE
C COEFFICIENT FOR THE GIVEN TRANSITION.
C
C IT IS APPLICABLE TO IMPURITIES IN PLASMA TRAVERSED BY
C NEUTRAL BEAMS OF H OR HE.
C
C THE RECOMBINED TARGET ION IS TREATED AS H-LIKE.
C
C THE MODEL INCLUDES CAPTURE, N-N' LEVEL CASCADE, AND MIXING
C AMONG L LEVELS OF SAME N BY COLLISIONS.
C
C AN INTERNAL EIKONAL APPROXIMATION IS USED FOR CAPTURE FROM
C EXCITED H OR HE STATES, ALTHOUGH NORMALLY THE EXTERNAL DATA
C SET SHOULD BE USED.
C
C CALLING PROGRAM: ADAS309
C
C INPUT : (I*4) MXNENG = MAXIMUM NO. OF ENERGIES IN DATA SET.
C INPUT : (I*4) MXNSHL = MAXIMUM NUMBER OF N SHELLS.
C INPUT : (I*4) NGRND = PRINCIPAL QUANTUM NUMBER OF GROUND STATE.
C INPUT : (I*4) NTOT = PRINCIPAL QUANTUM NUMBER OF HIGHEST BOUND
C STATE.
C INPUT : (I*4) NBOT = MINIMUM PRINCIPAL QUANTUM NUMBER FOR
C RATE TABLES.
C INPUT : (I*4) NTOP = MAXIMUM PRINCIPAL QUANTUM NUMBER FOR
C RATE TABLES.
C INPUT : (I*4) IRZ0 = RECEIVER NUCLEAR CHARGE.
C INPUT : (I*4) IRZ1 = RECEIVER ION INITIAL CHARGE.
C INPUT : (R*8) RAMSNO = RECEIVER ATOMIC MASS.
C INPUT : (R*8) TEV = ELECTRON TEMPERATURE.
C UNITS: EV
C INPUT : (R*8) TIEV = ION TEMPERATURE.
C UNITS: EV
C INPUT : (R*8) DENS = ELECTRON DENSITY.
C UNITS: CM-3
C INPUT : (R*8) DENSZ = PLASMA ION DENSITY.
C UNITS: CM-3
C INPUT : (R*8) ZEFF = EFFECTIVE ION CHARGE.
C INPUT : (R*8) BMAG = PLASMA MAGNETIC INDUCTION.
C UNITS: TESLA
C INPUT : (R*8) BMENG = BEAM ENERGY.
C UNITS: EV/AMU
C INPUT : (I*4) ITHEOR = CHARGE EXCHANGE MODEL OPTION.
C 1 => USE INPUT DATA SET.
C 2 => USE EIKONAL MODEL.
C INPUT : (I*4) IBSTAT = DONOR STATE FOR EIKONAL MODEL.
C 1 => H(1S)
C 2 => H(2S)
C 3 => H(2P)
C 4 => HE(1S2)
C 5 => HE(1S2S)
C INPUT : (I*4) IEMMS = EMISSION MEASURE MODEL OPTION.
C 1 => CHARGE EXCHANGE.
C 2 => ELECTRON IMPACT EXCITATION.
C INPUT : (I*4) NTL = LOWER PRINCIPAL QUANTUM NUMBER OF
C TRANSITION.
C INPUT : (I*4) NTU = UPPER PRINCIPAL QUANTUM NUMBER OF
C TRANSITION.
C INPUT : (I*4) NMINF = LOWEST N-SHELL FOR WHICH DATA READ.
C INPUT : (I*4) NMAXF = HIGHEST N-SHELL FOR WHICH DATA READ.
C INPUT : (I*4) NENRGY = NUMBER OF ENERGIES READ FROM DATA SET.
C INPUT : (R*8) ENRGYA() = COLLISION ENERGIES READ FROM INPUT DATA
C SET.
C UNITS: EV/AMU
C DIMENSION: ENERGY INDEX
C INPUT : (R*8) ALPHAA() = EXTRAPOLATION PARAMETER ALPHA READ FROM
C INPUT DATA SET.
C DIMENSION: ENERGY INDEX
C INPUT : (R*8) XSECNA(,) = N-RESOLVED CHARGE EXCHANGE CROSS-SECTIONS
C READ FROM INPUT DATA SET.
C UNITS: CM2
C 1ST DIMENSION: ENERGY INDEX
C 2ND DIMENSION: N-SHELL
C INPUT : (R*8) FRACLA(,) = L-RESOLVED CHARGE EXCHANGE CROSS-SECTIONS.
C AFTER CXDATA: ABSOLUTE VALUES (CM2).
C AFTER CXFRAC: FRACTION OF N-RESOLVED
C DATA.
C 1ST DIMENSION: ENERGY INDEX
C 2ND DIMENSION: INDEXED BY I4IDFL(N,L)
C
C OUTPUT: (R*8) ERATE = EFFECTIVE EMISSIVITY RATE COEFFICIENT FOR
C REQUESTED TRANSITION
C SPECTRUM LINE.
C UNITS: CM3 SEC-1
C
C PARAM : (I*4) MXN = MXNSHL.
C PARAM : (I*4) MXJSHL = MAXIMUM NUMBER OF J SUB-SHELLS.
C PARAM : (I*4) MXBEAM = MAXIMUM NUMBER OF BEAM COMPONENTS.
C PARAM : (I*4) MXOBSL = MAXIMUM NUMBER OF OBSERVED SPECTRUM
C LINES.
C
C PARAM : (I*4) MXPRSL = MAXIMUM NUMBER OF SPECTRUM LINES TO
C PREDICT.
C
C PARAM : (R*8) EMP = REDUCED MASS FOR POSITIVE ION.
C UNITS: ELECTRON MASSES
C
C (I*4) NBEAM = NUMBER OF BEAM ENERGIES.
C (I*4) NOLINE = NUMBER OF OBSERVED SPECTRUM LINES.
C (I*4) NPLINE = NUMBER OF SPECTRUM LINES TO PREDICT.
C (I*4) NUMIN = MINIMUM UPPER PRINCIPAL QUANTUM NUMBER
C FOR OBSERVED SPECTRUM LINES.
C (I*4) NUMAX = MAXIMUM UPPER PRINCIPAL QUANTUM NUMBER
C FOR OBSERVED SPECTRUM LINES.
C
C (R*8) EM = EMMISSION MEASURE.
C UNITS: CM-5
C
C (I*4) NL() = LIST OF LOWER PRINCIPAL QUANTUM NUMBERS
C OF OBSERVED SPECTRUM LINES.
C DIMENSION: SPECTRUM LINE INDEX.
C (I*4) NU() = LIST OF UPPER PRINCIPAL QUANTUM NUMBERS
C OF OBSERVED SPECTRUM LINES.
C DIMENSION: SPECTRUM LINE INDEX.
C (I*4) NPL() = LIST OF LOWER PRINCIPAL QUANTUM NUMBERS
C OF SPECTRUM LINES TO PREDICT.
C DIMENSION: SPECTRUM LINE INDEX.
C (I*4) NPU() = LIST OF UPPER PRINCIPAL QUANTUM NUMBERS
C OF SPECTRUM LINES TO PREDICT.
C DIMENSION: SPECTRUM LINE INDEX.
C
C (R*8) BMFRA() = BEAM COMPONENT FRACTIONS.
C DIMENSION: COMPONENT INDEX.
C (R*8) BMENA() = BEAM ENERGY COMPONENTS.
C UNITS: EV/AMU
C (R*8) EMISA() = LIST OF EMISSIVITIES OF OBSERVED SPECTRUM
C LINES.
C UNITS: PH CM-2 SEC-1
C DIMENSION: SPECTRUM LINE INDEX.
C (R*8) TBLF() = TABLE OF RADIATIVE LIFETIMES.
C UNITS: SECS
C DIMENSION: REFERENCED BY I4IDFL(N,L).
C (R*8) TBQEX() = MEAN EXCITATION RATE COEFFICIENTS FOR
C NL-LEVELS AVERAGED OVER BEAM FRACTIONS.
C UNITS: CM3 SEC-1
C DIMENSION: REFERENCED BY I4IDFL(N,L).
C (R*8) QTHEX() = MEAN EXCITATION RATE COEFFICIENTS FOR
C N-LEVELS AVERAGED OVER BEAM FRACTIONS.
C UNITS: CM3 SEC-1
C DIMENSION: N SHELL INDEX.
C (R*8) FTHEX() = FRACTION OF N-LEVEL MEAN EXCITATION RATE
C COEFFICIENTS IN NL-LEVEL.
C DIMENSION: REFERENCED BY I4IDFL(N,L).
C (R*8) QTHCH() = MEAN CHARGE EXCHANGE COEFFICIENTS FOR
C N-LEVELS AVERAGED OVER BEAM FRACTIONS.
C UNITS: CM3 SEC-1
C DIMENSION: N SHELL INDEX.
C (R*8) FTHCH() = FRACTION OF N-LEVEL MEAN CHARGE EXCHANGE
C COEFFICIENTS IN NL-LEVEL.
C DIMENSION: REFERENCED BY I4IDFL(N,L).
C (R*8) TBQMEP() = ELECTRON COLLISIONAL RATE COEFFT. FOR
C NL->NL+1.
C DIMENSION: REFERENCED BY I4IDFL(N,L).
C (R*8) TBQMEM() = ELECTRON COLLISIONAL RATE COEFFT. FOR
C NL->NL-1.
C DIMENSION: REFERENCED BY I4IDFL(N,L).
C (R*8) TBQMIP() = POSITIVE ION COLLISIONAL RATE COEFFT. FOR
C NL->NL+1.
C DIMENSION: REFERENCED BY I4IDFL(N,L).
C (R*8) TBQMIM() = POSITIVE ION COLLISIONAL RATE COEFFT. FOR
C NL->NL-1.
C DIMENSION: REFERENCED BY I4IDFL(N,L).
C (R*8) QEX() =
C DIMENSION: N SHELL INDEX.
C (R*8) TOTPOP() = TOTAL COLLISION POP. FOR PREDICTED
C SPECTRUM LINE.
C UNITS: CM-2
C DIMENSION: PREDICTED LINE INDEX.
C (R*8) TOTEMI() = TOTAL COLLISION EMISSIVITIES FOR PREDICTED
C SPECTRUM LINE.
C UNITS: PH CM-2 SEC-1
C DIMENSION: PREDICTED LINE INDEX.
C (R*8) AVRGWL() = AVERAGE AIR WAVELENGTH FOR PREDICTED
C SPECTRUM LINE.
C UNITS: A
C DIMENSION: PREDICTED LINE INDEX.
C (R*8) QEFF() = EFF. RATE COEFFICIENT FOR PREDICTED
C SPECTRUM LINE.
C UNITS: CM3 SEC-1
C DIMENSION: PREDICTED LINE INDEX.
C
C (R*8) TBLPOP(,) = TABLE OF COLLISION POP. FOR PREDICTED
C SPECTRUM LINE.
C UNITS: CM-2
C 1ST DIMENSION: PREDICTED LINE INDEX.
C 2ND DIMENSION: REFERENCED BY I4IDLI().
C (R*8) TBLEMI(,) = TABLE OF COLLISION EMISSIVITIES FOR
C PREDICTED SPECTRUM LINE.
C UNITS: PH CM-2 SEC-1
C 1ST DIMENSION: PREDICTED LINE INDEX.
C 2ND DIMENSION: REFERENCED BY I4IDLI().
C (R*8) TBLWLN(,) = TABLE OF WAVELENGTHS FOR PREDICTED
C SPECTRUM LINE.
C UNITS: A
C 1ST DIMENSION: PREDICTED LINE INDEX.
C 2ND DIMENSION: REFERENCED BY I4IDLI().
C
C ROUTINES:
C ROUTINE SOURCE BRIEF DESCRIPTION
C -------------------------------------------------------------
C I4UNIT ADAS RETURNS UNIT NO. FOR OUTPUT OF MESSAGES.
C CXTBLF ADAS FILLS L-RESOLVED RADIATIVE LIFETIME
C TABLE.
C CXTBEX ADAS FILLS N AND L-RESOLVED ELECTRON IMPACT
C EXCITATION RATE TABLES.
C CXQEIK ADAS FILLS N AND L-RESOLVED CHARGE EXCHANGE
C RATE TABLES USING EIKONAL APPROXIMATION.
C CXQXCH ADAS FILLS N AND L-RESOLVED CHARGE EXCHANGE
C RATE TABLES USING INPUT DATA SET.
C C8TBQM ADAS FILLS N AND L-RESOLVED COLLISIONAL RATE
C TABLES.
C C9EMIS ADAS PREDICTS THE L-RESOLVED EMISSIVITY FOR
C REQUESTED TRANSITIONS.
C
C AUTHOR: JONATHAN NASH (TESSELLA SUPPORT SERVICES PLC)
C K1/0/87
C JET EXT. 5183
C
C DATE: 03/12/93
C
C UNIX PORT: TIM HAMMOND (TESSELLA SUPPORT SERVICES PLC)
C DATE: 10/07/95 VERSION 1.1
C
C-----------------------------------------------------------------------
C
C-----------------------------------------------------------------------
INTEGER IBSTAT, IEMMS, IRZ0, IRZ1
INTEGER ITHEOR, MXNENG, MXNSHL, NBOT
INTEGER NENRGY, NGRND, NMAXF, NMINF
INTEGER NTL, NTOP, NTOT, NTU
REAL*8 ALPHAA(MXNENG), BMAG, BMENG
REAL*8 DENS, DENSZ, ENRGYA(MXNENG)
REAL*8 ERATE
REAL*8 FRACLA(MXNENG,(MXNSHL*(MXNSHL+1))/2), RAMSNO
REAL*8 TEV, TIEV
REAL*8 XSECNA(MXNENG,MXNSHL), ZEFF