ADAS Subroutine xxceia
SUBROUTINE XXCEIA( EIA ) C----------------------------------------------------------------------- C C ****************** FORTRAN77 SUBROUTINE: XXCEIA ********************* C C PURPOSE: CONVERTS IONISATION POTENTIALS FROM WAVE NUMBERS TO C RYDBERGS AND FILLS IN ANY MISSING VALUES UP TO AN ION CHARGE C OF 50. C C (IF 'EIA()' IS ALL ZERO - RETURN) C (PRIOR TO 16/08/90 WAS KNOWN AS 'E3EIA' - SLIGHTLY AMENDED) C C CALLING PROGRAMS: GENERAL USE C C SUBROUTINE: C C I/O : (R*8) EIA() = IONISATION POTENTIALS: ()=ION CHARGE C INPUT : UNITS - WAVE NUMBERS (CM-1) C OUTPUT: UNITS - RYDBERGS C C (R*8) WN2RYD = WAVE NUMBER (CM-1) TO RYDBERG CONVERSION C C (I*4) I4UNIT = FUNCTION (SEE ROUTINE SECTION BELOW) C (I*4) IPOT = NUMBER OF IONISATION POTENTIAL VALUES C PRESENT IN THE INPUT FILE. C (I*4) I = GENERAL USE C C (R*8) A1 = EXTRAPOLATION EQUATION COEFFICIENT C (R*8) A2 = EXTRAPOLATION EQUATION COEFFICIENT C (R*8) A3 = EXTRAPOLATION EQUATION COEFFICIENT C (R*8) XI = VALUE AT WHICH EXTRAPOLATION IS REQUIRED C C ROUTINES: C ROUTINE SOURCE BRIEF DESCRIPTION C ------------------------------------------------------------- C I4UNIT ADAS INTEGER*4 FUNCTION - C FETCH UNIT NUMBER FOR OUTPUT OF MESSAGES C C NOTE: C THE NUMBER OF IONISATION COEFFICIENTS PRESENT WILL BE EITHER C 30 OR 50. IF 30 THEN THE VALUES FROM 31 TO 50 NEED TO BE C EXTRAPOLATED. THE EXTRAPOLATION EQUATION IS BASED ON THE C VALUES OF EIA(20), EIA(25) AND EIA(50). C C THE EXTRAPOLATION EQUATION BEING: C C EIA(XI) = A1 + A2*XI + A2*XI*XI C C ROUTINES: NONE C C AUTHOR: PAUL E. BRIDEN (TESSELLA SUPPORT SERVICES PLC) C K1/0/81 C JET EXT. 4569 C C DATE: 16/08/90 C C UPDATE: 11/08/93 HP SUMMERS - DO NOT STOP IF IPOT.LE.30, BUT DO NOT C EXTRAPOLATE C C UNIX-IDL PORT: C C VERSION: 1.1 DATE: 06-09-95 C MODIFIED: TIM HAMMOND (TESSELLA SUPPORT SERVICES PLC) C - FIRST RELEASE (NO CHANGES) C C VERSION: 1.2 DATE: 20-07-07 C MODIFIED: Allan Whiteford C - Small modifications to comments to allow for C automatic documentation preparation. C C----------------------------------------------------------------------- C----------------------------------------------------------------------- REAL*8 EIA(50)