ADAS Subroutine b8corp
C
SUBROUTINE B8CORP(NLEV , MAXT , MAXD , NMET ,
& TEVA , COEFF )
C-----------------------------------------------------------------------
C
C ******************* FORTRAN77 SUBROUTINE: B8CORP ********************
C
C PURPOSE: Corrects unphysical low temperature recombination
C contributions to a PEC.
C
C There is a low temperature problem in the production
C calculation which gives unphysical recombination
C contributions to the PECs. Generally the first 3-4
C temperatures in the ADAS 96 standard are affected. This
C routine replaces the first 4 temperatures from the
C recombination contribution with extrapolated values
C from the remaining data.
C
C
C CALLING PROGRAM: ADAS208 (B8WR11)
C
C
C INPUT : (I*4) NLEV = NUMBER OF LEVELS
C INPUT : (I*4) MAXT = NUMBER OF TEMPERATURES
C INPUT : (I*4) MAXD = NUMBER OF DENSITIES
C INPUT : (I*4) NMET = NUMBER OF METASTABLES
C
C I/O : (R*4) COEFF() = RECOMBINATION CONTRIBUTION TO THE PEC
C (STVR IN ADAS208 CALL)
C 1ST DIMENSION : LEVELS
C 2ND DIMENSION : TEMPERATURES
C 3RD DIMENSION : DENSITIES
C 4TH DIMENSION : METASTABLES
C
C
C ROUTINES:
C ROUTINE SOURCE BRIEF DESCRIPTION
C ------------------------------------------------------------
C XXSPLN ADAS SPLINE SUBROUTINE
C R8FUN1 ADAS REAL*8 FUNCTION: ( X -> X )
C
C AUTHOR: Martin O'Mullane
C
C
C DATE: 14-09-99
C
C
C VERSION: 1.1 DATE: 14-09-99
C MODIFIED: Martin O'Mullane
C - First version
C
C VERSION : 1.2
C DATE : 26-10-99
C MODIFIED: Martin O'Mullane
C - Consider each level separately for extrapolation.
C - Increase NDLEV in line with adas208.
C
C VERSION : 1.3
C DATE : 20-07-07
C MODIFIED: Allan Whiteford
C - Small modification to comments to allow for automatic
C documentation preparation.
C
C VERSION : 1.4
C DATE : 02-09-2007
C MODIFIED: Martin O'Mullane
C - ITAG was defined as real*8 rather than integer.
C
C-----------------------------------------------------------------------
C
C-----------------------------------------------------------------------
INTEGER MAXD, MAXT, NLEV, NMET
REAL COEFF(NDLEV,NDTEM,NDDEN,NDMET)
REAL*8 TEVA(NDTEM)