ADAS Subroutine bxmcca
SUBROUTINE BXMCCA( NDLEV , IL ,
& LPSEL , LISEL ,
& DENE , DENP ,
& CRA ,
& CRCE , CRCP , CIE ,
& CC
& )
C-----------------------------------------------------------------------
C
C ****************** FORTRAN77 SUBROUTINE: BXMCCA *********************
C
C PURPOSE: TO CONSTRUCT WHOLE RATE MATRIX 'CC' FOR TRANSITIONS BETWEEN
C ALL ENERGY LEVELS AT A FIXED TEMPERATURE AND GIVEN DENSITY
C 'DENE/DENP'.
C
C CALLING PROGRAM: ADAS205/ADAS206
C
C SUBROUTINE:
C
C INPUT : (I*4) NDLEV = MAXIMUM NUMBER OF ENERGY LEVELS ALLOWED
C INPUT : (I*4) IL = NUMBER OF ENERGY LEVELS
C
C INPUT : (L*4) LPSEL = .TRUE. => INCLUDE PROTON COLLISIONS
C .FALSE. => DO NOT INCLUDE PROTON COLLISIONS
C INPUT : (L*4) LISEL = .TRUE. => INCLUDE IONISATION RATES
C .FALSE. => DO NOT INCLUDE IONISATION RATES
C
C INPUT : (R*8) DENE = ELECTRON DENSITY (UNITS: CM-3)
C INPUT : (R*8) DENP = PROTON DENSITY (UNITS: CM-3)
C
C INPUT : (R*8) CRA(,) = A-VALUE (sec-1) MATRIX COVERING ALL
C TRANSITIONS.
C 1st DIMENSION: ENERGY LEVEL INDEX
C 2nd DIMENSION: ENERGY LEVEL INDEX
C (NOTE: DIAGONAL ELEMENTS REPRESENT THE
C NEGATIVE SUM OF THEIR RESPECTIVE
C COLUMNS.)
C
C INPUT : (R*8) CRCE(,) = ELECTRON IMPACT TRANSITIONS:
C EXCIT'N/DE-EXCIT'N RATE COEFFT MATRIX
C COVERING ALL TRANSITIONS (cm**3/s).
C VALUES FOR GIVEN TEMPERATURE.
C 1st DIMENSION: ENERGY LEVEL INDEX
C 2nd DIMENSION: ENERGY LEVEL INDEX
C (NOTE: DIAGONAL ELEMENTS REPRESENT THE
C NEGATIVE SUM OF THEIR RESPECTIVE
C COLUMNS.)
C INPUT : (R*8) CRCP(,) = PROTON IMPACT TRANSITIONS:
C EXCIT'N/DE-EXCIT'N RATE COEFFT MATRIX
C COVERING ALL TRANSITIONS (cm**3/s).
C VALUES FOR GIVEN TEMPERATURE.
C 1st DIMENSION: ENERGY LEVEL INDEX
C 2nd DIMENSION: ENERGY LEVEL INDEX
C (NOTE: DIAGONAL ELEMENTS REPRESENT THE
C NEGATIVE SUM OF THEIR RESPECTIVE
C COLUMNS.)
C INPUT : (R*8) CIE() = IONISATION RATE COEFFICIENT VECTOR FOR
C FIXED TEMPERATURE.
C DIMENSION: ENERGY LEVEL INDEX
C
C OUTPUT: (R*8) CC(,) = RATE MATRIX COVERING ALL TRANSITIONS
C (UNITS: SEC-1)
C VALUES FOR GIVEN TEMPERATURE AND DENSITY.
C 1st DIMENSION: ENERGY LEVEL INDEX
C 2nd DIMENSION: ENERGY LEVEL INDEX
C
C (I*4) IS1 = ENERGY LEVEL ARRAY INDEX
C (I*4) IS2 = ENERGY LEVEL ARRAY INDEX
C
C
C ROUTINES: NONE
C
C
C AUTHOR: PAUL E. BRIDEN (TESSELLA SUPPORT SERVICES PLC)
C K1/0/81
C JET EXT. 4569
C
C DATE: 09/10/90
C
C-----------------------------------------------------------------------
C
C-----------------------------------------------------------------------
INTEGER IL, NDLEV
LOGICAL LISEL, LPSEL
REAL*8 CC(NDLEV,NDLEV), CIE(NDLEV)
REAL*8 CRA(NDLEV,NDLEV), CRCE(NDLEV,NDLEV)
REAL*8 CRCP(NDLEV,NDLEV), DENE, DENP