Search Site | Contact Details | FAQ

ADAS Subroutine c6tbrc

      SUBROUTINE C6TBRC( MXNSHL , MXJSHL , IZ1    , NBOT   ,
     &                   NTOP   , TEV    , QTHRC  , FTHRCJ
     &                 )
C
C
C-----------------------------------------------------------------------
C
C  ****************** FORTRAN77 SUBROUTINE: C6TBRC *********************
C
C  PURPOSE:  SETS UP A TABLE OF RADIATIVE RECOMBINATION RATE
C            COEFFICIENTS FOR A BARE NUCLEUS, HELIUM-LIKE OR NEON-LIKE
C            ION TO EXCITED NLJ LEVELS.
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)  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
C  OUTPUT: (R*8)  QTHRC()  = RECOMBINATION RATE COEFFICIENT TO LEVEL N.
C                            UNITS: CM3 SEC-1
C                            DIMENSION: N-SHELL
C  OUTPUT: (R*8)  FTHRCJ(,)= FRACTION OF RECOMBINATION RATE OF LEVEL N
C                            TO STATE NLJ.
C                            1ST DIMENSION: J-SHELL INDEX WHERE
C                                           1 GIVES J=L+0.5
C                                           2 GIVES J=L-0.5
C                            2ND DIMENSION: REFERENCED BY I4IDFL().
C
C  PARAM : (R*8)  P1       = BOLTZMANN CONSTANT.
C                            UNITS: EV K-1
C  PARAM : (R*8)  P2       =
C
C          (I*4)  N        = PRINCIPAL QUANTUM NUMBER OF BOUND ELECTRON.
C          (I*4)  L        = ORBITAL QUANTUM NUMBER OF BOUND ELECTRON.
C          (I*4)  L1       = ORBITAL QUANTUM NUMBER OF FREE ELECTRON.
C          (I*4)  LP       = TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM
C                           NUMBER OF PARENT STATE.
C          (I*4)  ISP      = 2*SP+1 WHERE SP IS TOTAL SPIN OF PAREN     T
C                            STATE.
C          (I*4)  LT       = TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM
C                            NUMBER OF BOUND SYSTEM.
C          (I*4)  LT1      = TOTAL ORBITAL ANGULAR MOMENTUM QUANTUM
C                            NUMBER OF FREE SYSTEM.
C          (I*4)  IS       = 2*S+1 WHERE S IS TOTAL SPIN OF SYSTEM.
C          (I*4)  IRES     = LEVEL OF RESOLUTION.
C                          = 1 :
C                          = 2 : ABOVE LT1 SUM.
C                          = 3 : ABOVE LT SUM.
C                          = 4 : ABOVE S SUM.
C                          = 5 : UNRESOLVED GBF.
C          (I*4)  I        = LOOP INDEX.
C          (I*4)  J        = LOOP INDEX.
C          (I*4)  IDL      = TABLE INDEX.
C
C          (R*8)  Z1       = REAL VALUE = IZ1.
C          (R*8)  TE       = ELECTRON TEMPERATURE.
C                            UNITS: K
C          (R*8)  V        = EFFECTIVE PRINCIPAL QUANTUM NUMBER OF BOUND
C                            ELECTRON.
C          (R*8)  FACT     =
C          (R*8)  SUM      =
C          (R*8)  XL       = REAL VALUE = L.
C          (R*8)  WL       =
C          (R*8)  T        =
C          (R*8)  PREC1    = RADIATIVE RECOMBINATION INTEGRAL.
C          (R*8)  PION1    = PHOTOIONISATION INTEGRAL.
C          (R*8)  PSTIM1   = STIMULATED RECOMBINATION INTEGRAL.
C          (R*8)  PREC2    = RADIATIVE RECOMBINATION INTEGRAL.
C          (R*8)  PION2    = PHOTOIONISATION INTEGRAL.
C          (R*8)  PSTIM2   = STIMULATED RECOMBINATION INTEGRAL.
C
C  PARAM : (R*8)  P1       = BOLTZMANN CONSTANT.
C                            UNITS: EV K-1
C  PARAM : (R*8)  P2       =
C
C ROUTINES:
C          ROUTINE    SOURCE    BRIEF DESCRIPTION
C          -------------------------------------------------------------
C          I4IDFL     ADAS      RETURNS UNIQUE INDEX GIVEN QUANTUM
C                               NUMBERS N AND L.
C          CXPHOT     ADAS      CALCULATES PHOTO INTEGRALS USING GIIH
C                               BOUND-FREE GAUNT-FACTORS.
C
C AUTHOR:  JONATHAN NASH (TESSELLA SUPPORT SERVICES PLC)
C          K1/0/81
C          JET EXT. 5183
C
C DATE:    05/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: 29-05-96
C MODIFIED: WILLIAM OSBORN
C               - REMOVED UNUSED VARIABLES
C
C-----------------------------------------------------------------------
C
C-----------------------------------------------------------------------
      INTEGER             IZ1,         MXJSHL,      MXNSHL,      NBOT
      INTEGER             NTOP
      REAL*8              FTHRCJ(MXJSHL,(MXNSHL*(MXNSHL+1))/2)
      REAL*8              QTHRC(MXNSHL),            TEV
© Copyright 1995-2018 The ADAS Project
Comments and questions to: adas-at-adas.ac.uk