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