ADAS Subroutine b8totl
C SUBROUTINE B8TOTL( NDLEV , NDMET , & NORD , NMET , & IORDR , IMETR , ISTRN , & DENSX , & STCKMX , STACKX , & PLA1 , PLBAX , & PLAX , PLX , & PSAX , PSX & ) C----------------------------------------------------------------------- C C ****************** FORTRAN77 SUBROUTINE: B8TOTL ********************* C C PURPOSE: TO CALCULATE TOTAL/SPECIFIC LINE POWERS FOR METASTABLES C AND TOTAL/SPECIFIC EQUILIBRIUM LINE POWERS. C C DEVELOPMENT OF B6TOTL. C C NOTE: A SPECIFIC LINE IS EVALUATED WHICH TERMINATES ON EACH C METASTABLE. EACH IS RESOLVED INTO THE PART DRIVEN BY EACH C METASTABLE. THE EQUILIBRIUM POWER IN EACH OF THESE LINES IS C ALSO EVALUATED USING THE EQUILIBRIUM METASTABLE FRACTIONS. C C CALLING PROGRAM: ADAS208 C C SUBROUTINE: C C INPUT : (I*4) NDLEV = PARAMETER = MAX. NO. OF LEVELS ALLOWED C INPUT : (I*4) NDMET = PARAMETER = MAX. NO. OF METASTABLES ALLOWED C C INPUT : (I*4) NORD = NUMBER OF ORD. LEVELS (1 <= NORD <= 'NDLEV') C INPUT : (I*4) NMET = NUMBER OF METASTABLES (1 <= NMET <= 'NDMET') C C INPUT : (I*4) IORDR() = INDEX OF ORDINARY LEVELS IN COMPLETE LEVEL C LIST (ARRAY SIZE = 'NDLEV' ) C INPUT : (I*4) IMETR() = INDEX OF METASTABLE IN COMPLETE LEVEL LIST C (ARRAY SIZE = 'NDMET' ) C INPUT : (I*4) ISTRN() = SPECIFIC LINE POWER: SELECTED ELECTRON C IMPACT TRANSITION INDEX. (FOR USE WITH C 'IE1A()' , 'IE2A()' AND 'AA()' ARRAYS) C WHICH GIVES LARGEST POWER TO METASTABLE C DIMENSION: METASTABLE LINE COUNT INDEX C INPUT : (R*8) DENSX = ELECTRON DENSITY (UNITS: CM-3) C C INPUT : (R*8) STCKMX() = METASTABLE POPULATIONS STACK C AT FIXED TEMPERATURE AND DENSITY. C DIMENSION: METASTABLE INDEX C INPUT : (R*4) STACKX(,)= ORDINARY EXCITED LEVEL POPULAT'N DEPENDENCE C ON METASTABLE LEVEL. AT FIXED TEMPERATURE C AND DENSITY. C 1st DIMENSION: ORDINARY LEVEL INDEX C 2nd DIMENSION: METASTABLE INDEX C C INPUT : (R*8) PLA1() = DIRECT LINE POWER LOSS FOR EACH LEVEL. C (UNITS: ERGS SEC-1) (DIMENSION: LEVEL INDEX) C INPUT : (R*8) PLBAX() = HIGH N PROJECTED POWER BASED ON EXCITATIONS C FROM A PARTICULAR METASTABLE TO LEVELS C 'IPROJ' UPWARDS. (UNITS: ERGS CM3 SEC-1) C AT FIXED TEMPERATURE. C DIMENSION: METASTABLE INDEX C C OUTPUT: (R*8) PLAX = TOTAL EQUILIBRIUM LINE POWER COEFFICIENTS. C AT FIXED TEMPERATURE AND DENSITY. C (UNITS: ERGS CM3 SEC-1) C OUTPUT: (R*8) PLX() = TOTAL LINE POWERS FOR METASTABLES. THIS IS C THE SUM OF ALL EMISSIONS ORGINATING IN THE C COLLISIONAL-RADIATIVE SENSE FROM THE C METASTABLE. AT FIXED TEMPERATURE AND DENSITY C (UNITS: ERGS SEC-1 ) C DIMENSION: METASTABLE INDEX C C OUTPUT: (R*8) PSAX() = TOTAL EQUILIBRIUM LINE POWER COEFFICIENTS. C AT FIXED TEMPERATURE AND DENSITY. C (UNITS: ERGS CM3 SEC-1) C DIMENSION: METASTABLE INDEX C OUTPUT: (R*8) PSX(,) = TOTAL LINE POWERS FOR METASTABLES. THIS IS C THE SUM OF ALL EMISSIONS ORGINATING IN THE C COLLISIONAL-RADIATIVE SENSE FROM THE C METASTABLE. AT FIXED TEMPERATURE AND DENSITY C (UNITS: ERGS SEC-1 ) C 1ST. DIMENSION: METASTABLE INDEX FOR LINE C 2ND. DIMENSION: METASTABLE INDEX OF DRIVER C C (I*4) IM = METASTABLE LEVEL ARRAY INDEX C (I*4) IS = ORDINARY LEVEL ARRAY INDEX C (I*4) ISL = SPECIFIC LINE POWER INDEX C C C C ROUTINES: NONE C C C AUTHOR: PAUL E. BRIDEN (TESSELLA SUPPORT SERVICES PLC) C K1/0/37 C JET EXT. 5023 C C DATE: 09/10/90 C C UPDATE: 20/05/93-P BRIDEN: STACKX ARRAY CHANGED FROM REAL*8 -> REAL*4 C C UPDATE: 24/05/96 HP SUMMERS - EXTENSION FOR SPECIFIC LINE POWER C C*********************************************************************** C PUT UNDER S.C.C.S CONTROL: C C VERSION: 1.1 DATE: 16/07/95 C MODIFIED: WILLIAM OSBORN (TESSELLA SUPPORT SERVICES PLC) C - FIRST PUT UNDER S.C.C.S C C VERSION: 1.2 DATE: 20/07/07 C MODIFIED: Allan Whiteford C - Small modification to comments to allow for automatic C documentation preparation. C C----------------------------------------------------------------------- C C----------------------------------------------------------------------- INTEGER IMETR(NMET), IORDR(NORD), ISTRN(NDMET) INTEGER NDLEV, NDMET, NMET, NORD REAL*8 DENSX, PLA1(NDLEV), PLAX REAL*8 PLBAX(NDMET), PLX(NDMET) REAL*8 PSAX(NDMET), PSX(NDMET,NDMET) REAL STACKX(NDLEV,NDMET) REAL*8 STCKMX(NDMET)