Search Site | Contact Details | FAQ

ADAS Subroutine c6tbqm

      SUBROUTINE C6TBQM( MXNSHL , MXJSHL , IZ0    , IZ1    ,
     &                   NBOT   , NTOP   , TEV    , DENS   ,
     &                   ZP     , TPV    , EMP    , TBLF   ,
     &                   TBQMEP , TBQMEM , TBQMIP , TBQMIM
     &                 )
C
C
C-----------------------------------------------------------------------
C
C  ****************** FORTRAN77 SUBROUTINE: C6TBQM *********************
C
C  PURPOSE:  SETS UP TABLES OF ELECTRON AND POSITIVE ION COLLISIONAL
C            RATE COEFFICIENTS BETWEEN NEARLY DEGENERATE LEVELS FOR
C            H-, LI-, AND NA-LIKE IONS.
C
C            THE RATES FOR THE SEPARATE TRANSITIONS NLJ->NL+1J' AND
C            NLJ->NL-1J' ARE OBTAINED.
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)  IZ0      = NUCLEAR CHARGE OF TARGET ION.
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  INPUT : (R*8)  DENS     = ELECTRON DENSITY.
C                            UNITS: CM-3
C  INPUT : (R*8)  ZP       = CHARGE OF COLLIDING POSITIVE ION.
C  INPUT : (R*8)  TPV      = TEMPERATURE (COLLIDING POSITIVE ION
C                            DISTRIBUTION).
C                            UNITS: EV
C  INPUT : (R*8)  EMP      = REDUCED MASS FOR COLLIDING POSITIVE ION.
C                            UNITS: ELECTRON MASSES
C  INPUT : (R*8)  ZP       = CHARGE OF COLLIDING POSITIVE ION.
C  INPUT : (R*8)  TPV      = POSITIVE ION TEMPERATURE.
C                            UNITS: EV
C  INPUT : (R*8)  EMP      = REDUCED MASS FOR COLLIDING POSITIVE ION.
C                            UNITS: ELECTRON MASSES
C  INPUT : (R*8)  TBLF()   = TABLE OF RADIATIVE LIFETIMES.
C                            UNITS: SECS
C                            DIMENSION: REFERENCED BY FUNC I4IDFL(N,L).
C
C  OUTPUT: (R*8)  TBQMEP(,)= ELECTRON RATE COEFFT. FOR NLJ->NL+1J'.
C                            1ST DIMENSION: J->J' TRANSITION INDEX.
C                            2ND DIMENSION: REFERENCED BY I4IDFL(N,L).
C  OUTPUT: (R*8)  TBQMEM(,)= ELECTRON RATE COEFFT. FOR NLJ->NL-1J'.
C                            1ST DIMENSION: J->J' TRANSITION INDEX.
C                            2ND DIMENSION: REFERENCED BY I4IDFL(N,L).
C  OUTPUT: (R*8)  TBQMIP(,)= POSITIVE ION RATE COEFFT. FOR NLJ->NL+1J'.
C                            1ST DIMENSION: J->J' TRANSITION INDEX.
C                            2ND DIMENSION: REFERENCED BY I4IDFL(N,L).
C  OUTPUT: (R*8)  TBQMIM(,)= POSITIVE ION RATE COEFFT. FOR NLJ->NL-1J'.
C                            1ST DIMENSION: J->J' TRANSITION INDEX.
C                            2ND DIMENSION: REFERENCED BY I4IDFL(N,L).
C
C  PARAM : (I*4)  MXJ     = 'MXJSHL'.
C
C          (I*4)  NI      = VALENCE ELECTRON PRINCIPAL QUANTUM NUMBER
C                           IN STATE I.
C          (I*4)  NJ      = VALENCE ELECTRON PRINCIPAL QUANTUM NUMBER
C                           IN STATE J.
C          (I*4)  LI      = VALENCE ELECTRON ORBITAL QUANTUM NUMBER IN
C                           STATE I.
C          (I*4)  LJ      = VALENCE ELECTRON ORBITAL QUANTUM NUMBER IN
C                           STATE J.
C          (I*4)  IDLI    = TABLE INDEX.
C          (I*4)  IDLJ    = TABLE INDEX.
C          (I*4)  I       = LOOP INDEX.
C          (I*4)  J       = LOOP INDEX.
C
C          (R*8)  GAE     = GAMA RATE PARAMETER FOR ELECTRON COLLISIONS.
C          (R*8)  GAP     = GAMA RATE PARAMETER FOR POSITIVE ION
C                           COLLISIONS.
C
C          (R*8)  QEP()   = ELECTRON RATE COEFFT. FOR NLJ->NL+1J'
C                           DIMENSION: J->J' TRANSITION INDEX.
C          (R*8)  QEM()   = ELECTRON RATE COEFFT. FOR NLJ->NL-1J'
C                           DIMENSION: J->J' TRANSITION INDEX.
C          (R*8)  QIP()   = POSITIVE ION RATE COEFFT. FOR NLJ->NL+1J'
C                           DIMENSION: J->J' TRANSITION INDEX.
C          (R*8)  QIM()   = POSITIVE ION RATE COEFFT. FOR NLJ->NL-1J'
C                           DIMENSION: J->J' TRANSITION INDEX.
C
C NOTES:
C       1) THE J->J' TRANSITION INDEX IS AS FOLLOWS:
C              1 : J=L+0.5 -> J'=L'+0.5
C              2 : J=L+0.5 -> J'=L'-0.5
C              3 : J=L-0.5 -> J'=L'+0.5
C              4 : J=L-0.5 -> J'=L'-0.5
C
C       2) BEFORE CALLING C6TBQM THE LIFETIME TABLE MUST BE FILLED
C          WITH A CALL TO C6TBLF.
C
C ROUTINES:
C          ROUTINE    SOURCE    BRIEF DESCRIPTION
C          -------------------------------------------------------------
C          I4UNIT     ADAS      RETURNS UNIT NO. FOR OUTPUT OF MESSAGES.
C          I4IDFL     ADAS      RETURNS UNIQUE INDEX GIVEN QUANTUM
C                               NUMBERS N AND L.
C          CXCRDG     ADAS      CALCULATES COLLISIONAL RATE COEFFICIENTS
C                               BETWEEN NEARLY DEGENERATE LEVELS OF
C                               H-, LI- OR NA-LIKE IONS.
C
C AUTHOR:  JONATHAN NASH (TESSELLA SUPPORT SERVICES PLC)
C          K1/0/81
C          JET EXT. 5183
C
C DATE:    02/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: 17-05-07
C MODIFIED: Allan Whiteford
C               - Updated comments as part of subroutine documentation
C                 procedure.
C
C-----------------------------------------------------------------------
C
C-----------------------------------------------------------------------
      INTEGER             IZ0,         IZ1,         MXJSHL,      MXNSHL
      INTEGER             NBOT,        NTOP
      REAL*8              DENS,        EMP
      REAL*8              TBLF((MXNSHL*(MXNSHL+1))/2)
      REAL*8              TBQMEM(2*MXJSHL,(MXNSHL*(MXNSHL+1))/2)
      REAL*8              TBQMEP(2*MXJSHL,(MXNSHL*(MXNSHL+1))/2)
      REAL*8              TBQMIM(2*MXJSHL,(MXNSHL*(MXNSHL+1))/2)
      REAL*8              TBQMIP(2*MXJSHL,(MXNSHL*(MXNSHL+1))/2)
      REAL*8              TEV,         TPV,         ZP
© Copyright 1995-2024 The ADAS Project
Comments and questions to: adas-at-adas.ac.uk