Search Site | Contact Details | FAQ

ADAS Subroutine gpcalcx


       REAL*8 FUNCTION GPCALCX(XV0)                                     
       IMPLICIT REAL*8(A-H,O-Z)
C-----------------------------------------------------------------------
C  PURPOSE: ROUTINE TO PROVIDE BURGESS GENERAL PROGRAM RESULTS AT A
C  SERIES OF TEMPERATURES AND AT ZERO DENSITY.                           
C
C  EQUAL THE GENERAL FORMULA RESULTS AS FAR AS POSSIBLE BY MODIFICATION 
C  OF BETHE CORRECTIONS VIA A SINGLE SCALING PARAMETER CORFAC.          
C   THE CORRECTION FACTORS USED IN THE GENERAL PROGRAM                  
C  ARE OBTAINED BY ADJUSTMENT OF STANDARD SETS FOR SPECIFIC TYPES OF    
C  TRANSITION.  THE ADJUSTMENT IS                                       
C    (NEW COR(J))=EXP(-CORFAC/(L**DF+0.5))*(STANDARD COR(J)             
C  THE STANDARD COR'S ARE AS FOLLOWS:                                   
C   TYPE    TRANSITION                     COR'S                     DF 
C     1   NI=1,NJ>=2,LJ=LI+1:       0.05,0.30,0.50,0.90              2.0
C     2   NI=2,NJ=3,LJ=LI+1:        0.01,0.02,0.20,0.40,0.70,0.90    1.0
C     3   NI=2,NJ=3,LJ=LI-1:        0.01,0.01,0.01,0.08,0.30,0.70    1.0
C     4   NJ-NI=0, LJ=LI+1 :        0.30,0.35,0.40,0.45,0.70,0.90    0.5
C     5   NJ-NI=0, LJ=LI-1 :        0.30,0.35,0.40,0.45,0.70,0.90    0.5
C     6   NJ-NI>0, LJ=LI+1 :        0.01,0.02,0.20,0.40,0.70,0.90    1.0
C     7   NJ-NI>0, LJ=LI-1 :        0.01,0.01,0.01,0.08,0.30,0.70    1.0
C
C     (1)  INCLUDE NCUT AND EXTEND ARRAY SIZES
C     (2)  IMPLIMENTATION OF NCUT,LOW TEMPERATURE CHECK, CORRECTION
C          INVOLVING V1
C     (3)  NRAT INCREASED FROM 10 TO 15
C
C  ********** H.P. SUMMERS, JET             11 JUNE 1987  ************* 
C  ********** H.P. SUMMERS, JET     MOD.(1) 24 AUG  1989  ************* 
C  ********** W.J. DICKSON, JET     MOD.(2) 14 DEC  1989  *************
C  ********** P.E. BRIDEN ,TESSELLA MOD.(3) 23 AUG  1994  *************
C  INPUT                                                                
C      MAXT=NUMBER OF TEMPERATURES                                      
C      TEA(I)=ELECTRON TEMPERATURES (K)                                 
C      Z1=RECOMBINING ION CHARGE                                        
C      N0=LOWEST ACCESSIBLE N-SHELL BY RECOMBINATION                    
C      V0=EFFECTIVE PRINCIPAL QUANTUM NUMBER OF LOWEST ACCESSIBLE SHELL 
C      NI=LOWER PRINCIPAL QUANTUM NUMBER OF PARENT TRANSITION           
C      LI=LOWER ANGULAR QUANTUM NUMBER OF PARENT TRNASITION.            
C      WI=LOWER PARENT STATE STATISTICAL WEIGHT.                        
C      NJ=UPPER PRINCIPAL QUANTUM NUMBER OF PARENT TRANSITION           
C      LJ=UPPER ANGULAR QUANTUM NUMBER OF PARENT TRNASITION.            
C      WJ=UPPER PARENT STATE STATISTICAL WEIGHT.                        
C      EIJ=PARENT TRANSITION ENERGY (RYD)                               
C      FIJ=ABSORPTION OSCILLATOR STRENGTH OF PARENT TRANSITION          
C      EDISPG=UNIFORM ENERGY DISPLACEMENT FOR GENERAL FORMULA           
C      SCALEG=UNIFORM SCALING OF GENERAL FORMULA                        
C      PHFRAC=INITIAL ESTIMATE OF PHASE SPACE FACTOR                    
C      CORFAC=INITIAL ESTIMATE OF BETHE CORRECTION SCALER               
C      NCUT  =HIGH N CUT-OFF (APPLICABLE TO METASTABLE INITIAL STATES)  
C  OUTPUT                                                               
C      ALFO(I)=GENERAL PROGRAM DIELECTRONIC COEFFICIENTS (CM+3 SEC-1)   
C      PHFRAC=REVISED PHASE SPACE FACTOR                                
C      CORFAC=REVISED BETHE CORRECTION SCALER                           
C UNIX-IDL PORT:
C	   WILLIAM OSBORN, TESSELLA SUPPORT SERVICES PLC.
C
C DATE:    19TH APRIL 1996
C
C VERSION: 1.1				DATE: 19-04-96
C MODIFIED: WILLIAM OSBORN
C	    - NRAT INCREASED FROM 15 TO 100 IN LINE WITH GPCALL
C
C VERSION: 1.2				DATE: 25-04-96
C MODIFIED: WILLIAM OSBORN
C	    - CONVERTED TO A FUNCTION FOR FMIN TO USE
C
C VERSION: 1.3				DATE: 16-05-07
C MODIFIED: Allan Whiteford
C	    - Old IBM statement labels in columns 73-80 removed. 
C	    - Updated comments as part of subroutine documentation
C             procedure.
C
C-----------------------------------------------------------------------
      REAL*8              XV0
© Copyright 1995-2018 The ADAS Project
Comments and questions to: adas-at-adas.ac.uk