Search Site | Contact Details | FAQ

ADAS Subroutine start7

       SUBROUTINE START7(IUTMP,IUPS1,IUPS2,STITLE,DSLPATH,
     &                   NBENG,NTEMP,NDENS, lbndl, lproj)
       IMPLICIT REAL*8(A-H,O-Z)
C-----------------------------------------------------------------------
C
C       *********  FORTRAN 77 ROUTINE : START7.F **********
C
C       PURPOSE : CALCULATION OF THE RESOLVED-NL POPULATION
C                 STRUCTURE.
C
C       NOTE    : THE RELEVANT QUANTUM NUMBERS ASSOCIATED WITH
C                 THE POPULATION STRUCTURE CALCULATION ARE
C                 STORED IN INTEGER ARRAYS WERE EACH ELEMENT
C                 IS 4 BYTES. THE FIRST 10 BITS OF THE 4 BYTE
C                 INTEGER ELEMENT OF THE ARRAY IS USED TO
C                 STORE THE TOTAL ANGULAR MOMENTUM QUANTUM
C                 NUMBER,L. THE NEXT TEN BITS IS USED TO STORE
C                 THE ORBITAL ANGULAR MOMENTUM QUANTUM
C                 NUMBER,l. THE LAST 12 BITS ARE USED TO STORE
C                 THE PRINCIPAL QUANTUM NUMBER,N.
C
C                 |<---------4 BYTE INTEGER--------->|
C                 |<---------32 BIT INTEGER--------->|
C                 |<----N----->|<----l--->|<---L---->|
C
C               [ |.....N......|.....l....|.....L....|  ]
C
C                 BIT OPERATORS ARE THEN EMPLOYED TO
C                 INTEROGATE ARRAYS,E.G IAND,ISHFR,
C                 USING HEXIDECIMAL MASKS.
C               
C
C       INPUT   :
C               
C       (CHR)   TITLE()    : NAME OF THE NEUTRAL BEAM SPECIES.
C       (I*4)   MN         :
C       (R*8)   RX3        :
C       (R*8)   DPT        :
C       (R*8)   EHCT       : CRITICAL ENERGY ?.
C       (I*4)   NHCT       : CRITICAL PRINCIPAL QUANTUM
C                            NUMBER ?.
C       (I*4)   LHCT       : CRITICAL ORBITAL QUANTUM
C                            NUMBER ?.
C       (I*4)   NIP        : RANGE OF DELTA N FOR IMPACT
C                            PARAMETER XSECTS. (LE.4)
C       (I*4)   NEX        :
C       (I*4)   IPRT       :
C       (I*4)   NDEL       :
C       (I*4)   INTD       : ORDER OF MAXWELL QUADRATURE
C                            FOR XSECTS. (LE.3)
C       (I*4)   IPRS       : 0 DEFAULT TO VAN REGEMORTER XSECTS.
C                              BEYOND NIP RANGE
C                            1 USE PERCIVAL-RICHARDS XSECTS.
C                              BEYOND NIP RANGE
C       (I*4)   ILOW       : 0 NO SPECIAL LOW LEVEL DATA ACCESSED.
C                            1 SPECIAL LOW LEVEL DATA ACCESSED.
C       (I*4)   IONIP      : 0 NO ION IMPACT COLLISIONS INCLUDED.
C                            1 ION IMPACT EXCITATION AND IONISATION
C                              INCLUDED.
C
C       (I*4)   NIONIP     : RANGE OF DELTA N FOR ION IMPACT
C                            EXCITATION XSECTS.
C       (I*4)   ILPRS      : 0 DEFAULT TO VAINSHTEIN ION IMPACT
C                              EXCITATION XSECTS.
C                            1 USE LODGE-PERCIVAL-RICHARDS ION
C                              IMPACT EXCITATION XSECTS.
C       (I*4)   IVDISP     : 0 ION IMPACT AT THERMAL MAXWELLIAN
C                              ENERGIES
C                            1 ION IMPACT AT DISPLACED THERMAL
C                              ENERGIES ACCORDING TO THE NEUTRAL
C                              BEAM ENERGY PARAMETER. IF(IVDISP=0 THEN
C                              SPECIAL LOW LEVEL DATA FOR ION IMPACT
C                              IS NOT SUBSTITUTED - ONLY VAINSHTEIN
C                              AND LODGE ET AL. OPTIONS ARE OPEN.
C                              ELECTRON IMPACT DATA SUBSTITUTION
C                              DOES OCCUR.
C       (R*4)   ZEFF       : EFFECTIVE CHARGE OF THE PLASMA.
C       (R*8)   DEDEG      : CRITICAL TRANSITION ENERGY (RYDBERGS)
C                            USED IN THE NEARLY DEGENERATE LEVEL
C                            TREATMENT.
C                            IF DE<=DEDEG THEN ASSUME ZERO A-VALUE
C                            AND NO SUPPLEMENTARY DATA.
C                            IF DE> DEDEG THEN ASSUME A-VALUE
C                            CALCULABLE AND SEARCH FOR SUPPLEMENTARY
C                            DATA.
C                            N.B. APPLIES TO DELTA N TRANSITIONS
C                            ONLY.
C
C       (I*4)   NL1        : PRINCIPAL QUANTUM NUMBER FROM WHICH
C                            THE RESOLVED-NL POPULATION STRUCTURE
C                            CALCULATION STARTS FROM
C                       
C       (I*4)   NL2        : PRINCIPAL QUANTUM NUMBER WHICH MARKS
C                            THE END OF THE RESOLVED-NL TREATMENT
C                            WITHIN THE POPULATION STRUCTURE
C                            CALCULATION AND INDICATES THE START
C                            OF THE BUNDLED-N APPROXIMATION.
C       (I*4)   NL3        : UPPER PRINCIPAL QUANTUM NUMBER OF
C                            THE BUNDLED-N APPROXIMATION.
C       (R*8)   Z0         : NUCLEAR CHARGE OF BEAM ATOM ?.
C       (R*8)   Z1         : ION CHARGE+1 OF BEAM ION ?.
C       (R*8)   ALF        : ADJUSTABLE PARAMETER ASSOCIATED
C                            WITH THE MODIFIED POTENTIAL USED
C                            WHEN SOLVING THE RADIAL WAVE
C                            EQUATION.
C       (R*8)   AMSZ0      :
C       (R*8)   AMSHYD     :
C       (I*4)   LP         :
C       (I*4)   ISP        :
C
C
C       OUTPUT  :
C       
C               (R*8)  ....ETC
C
C       GENERAL :
C
C       (I*4)   NLREP()    : ARRAY CONTAINING REPRESENTATIVE
C                            LEVELS.
C       (R*8)   ENL()      : EFFECTIVE PRINCIPAL QUANTUM
C                            NUMBER.
C       (R*8)   ENL2()     : RECIPROCAL OF THE EFFECTIVE PRINCIPAL
C                            QUANTUM NUMBER SQUARED.
C       (I*4)   KPF()      : ARRAY CONTAINING THE QUANTUM NUMBERS,
C                            N,l,L FROM NMIN TO NMAX, IN ORDER
C                            OF DECREASING BINDING ENERGY. SEE
C                            NOTE AT THE TOP OF PROGRAM.
C       (I*4)   KPB()      : ARRAY CONTAINING THE INDEX OF THE
C                            CORRESPONDING LEVEL IN KPF().
C       (R*8)   EGY        : IONISATION POTENTIAL ( RYDBERGS ).
C       (I*4)   IR         : COUNTER TO REFERENCE REPRESENTATIVE
C                            LEVELS.
C       (I*4)   I          : GENERAL COUNTER.
C       (R*8)   V          : EFFECTIVE PRINCIPAL QUANTUM NUMBER.
C       (R*8)   E          : RECIPROCAL OF THE EFFECTIVE PRINCIPAL
C                            QUANTUM NUMBER SQUARED.
C       (R*8)   EXE        : VARIABLE USED TO ASSIGN THE VALUE
C                            OF EXP(I/k*Te).
C       (R*8)   EXS        : VARIABLE USED TO ASSIGN THE VALUE
C                            OF EXP(I/k*TS)
C       (I*4)   K          : GENERAL COUNTER.
C       (R*8)   C1()       : COEFFICIENT OF THE QUANTUM DEFECT
C                            EXPANSION.
C       (R*8)   C2()       : COEFFICIENT OF THE QUANTUM DEFECT
C                            EXPANSION.
C       (R*8)   C3()       : COEFFICIENT OF THE QUANTUM DEFECT
C                            EXPANSION.
C
C
C
C
C       
C       ROUTINES:
C
C               ROUTINE        SOURCE     BRIEF DESCRIPTION
C               -----------------------------------------------
c               SUPPHE1         ADAS    OBTAINS FUNDAMENTAL DATA
C                                       FROM APPROPRIATE DATABASES.
C               OVLP            ADAS    ?????????????????????????
C               SETUP3          ADAS    ?????????????????????????
C               SPIJ            ADAS    ?????????????????????????
C               CCNST7          ADAS    ASSEMBLES ARRAYS USED TO
C                                       CONSTRUCT THE COLLISIONAL-
C                                       RADIATIVE MATRIX.
C               CCNSE4          ADAS    APPLIES MATRIX CONDENSATION
C                                       SCHEME TO ARRAYS USED TO
C                                       ASSEMBLE THE COLLISIONAL-
C                                       RADIATIVE MATRIX.
C               HYSCL           ADAS    ?????????????????????????
C
C
C       HISTORY : ORIGINALLY WRITTEN BY H.P.SUMMERS.
C
C       NOTE    : THE RESOLVED-NL CALCULATION WAS STRUCTURED
C                 IN SUCH A MANNER THAT THE CALCULATION WAS
C                 PERFORMED IN TWO STAGES. THIS TWO STAGE
C                 PROCCESS HAS BEEN REMOVED. CUBIC SPLINE
C                 INTERPLOATION IN L. THE DIMENSIONALITIES
C                 FOR EACH SPIN SYSTEM ARE AS FOLLOWS :
C
C                 NUMBER OF LEVELS <1000.
C                 NUMBER OF PRINCIPAL QUANTUM LEVELS<300.
C                 NUMBER OF RESOLVED PRINCIPAL QUANTUM LEVELS<40.
C                 NUMBER OF RESOLVED LEVELS <800
C                 NUMBER OF RESOLVED REPR.PRINC.QUANTUM LEVELS<11.
C                 NUMBER OF REPRESENTATIVE LEVELS<80.
C                 NUMBER OF PRINCIPAL QUANTUM REPRESENTATIVE LEVELS<30.
C
C
C       CONTACT : HARVEY ANDERSON
C                 UNIVERSITY OF STRATHCLYDE
C                 ANDERSON@PHYS.STRATH.AC.UK
C
C       DATE    : 2/2/98
C
C       NOTE :
C                 THE C-R MATRIX IS PASSED TO THE ROUTINE FINISH5.F
C                 VIA A SCRTACH FILE ON STREAM 12. IDEALLY THE VARIBLES
C                 SHOULD BE PASSED DIRECTLY TO THE ROUTINE.
C
C       
C       MODIFIED
C
C VERSION:      1.2                                     DATE:   21-10-99
C MODIFIED: RICHARD MARTIN
C               CHANGED HEXADECIMAL CONSTANTS TO Z'FFF00000' FORM.
C
C VERSION : 1.3
C MODIFIED: Martin O'Mullane      
C DATE    : 08-11-2004
C           Alter nmax in gamaf() from 200 to 500.        
C
C VERSION : 1.4
C DATE    : 18-11-2004
C MODIFIED: Martin O'Mullane      
C            - Align with Harvey Anderson's last version.
C            - Add lproj if projection output is requested.
C            - Add lbndl if adf36 output file is requested.
C            - The dsnps1 variable is replaced by iups2 in the 
C              parameter list.
C
C VERSION : 1.5
C DATE    : 16-05-07
C MODIFIED: Allan Whiteford    
C            - Moved parameter statement to below comment block
C              as part of subroutine documentation procedure.
C
C-----------------------------------------------------------------------
      CHARACTER*80        DSLPATH,     STITLE
      INTEGER             IUPS1,       IUPS2,       IUTMP,       NBENG
      INTEGER             NDENS,       NTEMP
      LOGICAL             LBNDL,       LPROJ
© Copyright 1995-2018 The ADAS Project
Comments and questions to: adas-at-adas.ac.uk