ADAS Subroutine cxsgei
SUBROUTINE CXSGEI( MXNSHL , IZ1 , ZT , ZT1 , THETA ,
& VEL , NA , LA , NB ,
& XSECNA , FRACLA
& )
C
C
C-----------------------------------------------------------------------
C
C ****************** FORTRAN77 SUBROUTINE: CXSGEI *********************
C
C PURPOSE: USES THE EIKONAL APPROXIMATION TO CALCULATE CHARGE EXCHANGE
C CROSS-SECTIONS INTO N'L' RESOLVED EXCITED STATES OF
C HYDROGENIC IONS IN CAPTURE FROM HYDROGEN NL STATES.
C
C THE SUBROUTINE CAN BE USED FOR CAPTURE FROM HELIUM BY
C APPROPRIATE CHOICE OF ZT AND ZT1.
C
C FOR HYDROGEN, ZT = ZT1 = 1.0
C FOR HELIUM , ZT = ZT1 = 1.6875
C
C EICHLER (1981) PHYS.REV.A,23,498.
C
C CALLING PROGRAM: ADAS308.
C
C INPUT : (I*4) MXNSHL = MAXIMUM VALUE OF N QUANTUM NUMBER.
C INPUT : (I*4) IZ1 = HYDROGENIC ION CHARGE (PROJECTILE).
C INPUT : (R*8) ZT = SCREENING CHARGE FOR THE 1S ELECTRON OF THE
C TARGET ATOM IN THE INITIAL STATE.
C INPUT : (R*8) ZT1 = EFFECTIVE CHARGE FOR THE 1S ELECTRON OF THE
C TARGET ATOM IN THE FINAL STATE.
C INPUT : (R*8) THETA = PARAMETER TO GIVE CORRECT BINDING ENERGY
C FOR INITIAL TARGET STATE.
C INPUT : (R*8) VEL = COLLISION SPEED.
C UNITS: AT. UNITS
C INPUT : (I*4) NA = PRINCIPAL QUANTUM NUMBER OF INTIAL STATE OF
C TARGET FROM WHICH CAPTURE IS MADE.
C INPUT : (I*4) LA = ORBITAL QUANTUM NUMBER OF INTIAL STATE OF
C TARGET FROM WHICH CAPTURE IS MADE.
C INPUT : (I*4) NB = PRINCIPAL QUANTUM NUMBER OF FINAL STATES OF
C PROJECTILE TO WHICH CAPTURE IS MADE.
C
C OUTPUT: (R*8) XSECNA = N-RESOLVED CROSS-SECTION FOR CAPTURE.
C UNITS: AT. UNITS
C OUTPUT: (R*8) FRACLA() = L-RESOLVED CROSS-SECTION AS A FRACTION OF
C CORRESPONDING N-RESOLVED CROSS-SECTION.
C DIMENSION: REFERENCED BY L QUANTUM NUMBER.
C
C PARAM : (I*4) MXN = 'MXNSHL'.
C PARAM : (R*8) PI = PI.
C
C (I*4) LB = ORBITAL QUANTUM NUMBER OF FINAL STATE OF
C PROJECTILE.
C (I*4) I = LOOP INDEX.
C (I*4) K = LOOP INDEX.
C (I*4) IU = LOOP INDEX.
C (I*4) MB = LOOP INDEX.
C (I*4) MB1 = LOOP INDEX.
C (I*4) IND =
C (I*4) IL = LOOP INDEX.
C (I*4) IG = LOOP INDEX.
C (I*4) IWM = LOOP INDEX.
C (I*4) IW1 = LOOP INDEX.
C (I*4) IW2 = LOOP INDEX.
C (I*4) IW3 = LOOP INDEX.
C (I*4) IW4 = LOOP INDEX.
C (I*4) IBA = LOOP INDEX.
C (I*4) IVA = LOOP INDEX.
C (I*4) ITA = LOOP INDEX.
C (I*4) IAM =
C (I*4) IA1 = LOOP INDEX.
C (I*4) IA2 = LOOP INDEX.
C (I*4) IA3 =
C (I*4) ISA =
C (I*4) IBB = LOOP INDEX.
C (I*4) IVB = LOOP INDEX.
C (I*4) ITB = LOOP INDEX.
C (I*4) IAM =
C (I*4) IB1 = LOOP INDEX.
C (I*4) IB2 = LOOP INDEX.
C (I*4) IB3 =
C (I*4) ISB =
C (I*4) IP = LOOP INDEX.
C
C (R*8) Z1 = REAL VALUE = IZ1.
C (R*8) XNA = REAL VALUE = N.
C (R*8) XNB = REAL VALUE = N.
C (R*8) ETA =
C (R*8) EPS =
C (R*8) PM =
C (R*8) PP =
C (R*8) QNA =
C (R*8) QNB =
C (R*8) QP =
C (R*8) ZT12 = ZT12**2.
C (R*8) T1 =
C (R*8) T2 =
C (R*8) T3 =
C (R*8) T4 =
C (R*8) T5 =
C (R*8) T6 =
C (R*8) T7 =
C (R*8) T8 =
C (R*8) SUM =
C (R*8) SUM1 =
C (R*8) BB =
C (R*8) BB1 =
C
C (Z*16) CTA =
C (Z*16) CTB =
C (Z*16) CT1 =
C (Z*16) CSUM =
C
C (R*8) B() =
C (R*8) F() =
C
C ROUTINES:
C ROUTINE SOURCE BRIEF DESCRIPTION
C -------------------------------------------------------------
C I4UNIT ADAS RETURNS UNIT NO. FOR OUTPUT OF MESSAGES.
C I4JGAM ADAS RETURNS VALUE FROM 'JGAM' TABLE.
C R8GAM ADAS RETURNS VALUE FROM 'GAM' TABLE.
C Z16CD ADAS
C
C AUTHOR: JONATHAN NASH (TESSELLA SUPPORT SERVICES PLC)
C K1/0/81
C JET EXT. 5183
C
C DATE: 06/10/93
C
C-----------------------------------------------------------------------
C
C-----------------------------------------------------------------------
INTEGER IZ1, LA, MXNSHL, NA
INTEGER NB
REAL*8 FRACLA(MXNSHL), THETA, VEL
REAL*8 XSECNA, ZT, ZT1