Search Site | Contact Details | FAQ

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
© Copyright 1995-2018 The ADAS Project
Comments and questions to: adas-at-adas.ac.uk