ADAS Subroutine cdintp

	SUBROUTINE CDINTP(CHOICE    ,OPTION   ,IECOUNT  ,INCOUNT  ,
     &			  MAXNE	    ,MAXTE    ,MAXEB    ,ZDATA    ,
     &			  GCRC1     ,GCRC2    ,GCRC3    ,GCRC4    ,
     &			  GCRC5     ,GCRC6    ,GCRC7    ,GCRC8    ,
     &			  GCRC9     ,IEREF    ,INREF    ,ITREF    ,
     &			  TVAL      ,RELMET   ,FI	,INA      ,
     &			  IEA       ,FII      ,FIII     ,NSPIN    ,
     &			  NNBNL     ,AVALUE   ,NERAY    ,ITA      ,
     &			  EBRAY     ,XDATA    ,YDATA    ,MULTI  )
C
C-----------------------------------------------------------------------
C
C  		****** FORTRAN77 ROUTINE: CDINTP 	********
C
C  	PURPOSE: INTERPOLATE BETWEEN THE EFFECTIVE CLOSE COUPLING
C		 COEFFICIENTS OF EFFECTIVE EMISSION COEFFICIENTS.
C
C
C	CALLING PROGRAM : ADAS313
C			
C	INPUT  :
C
C		(R*8)	FI()		: THE CONTRIBUTION FROM
C					  EXCITATION TO POPULATE A
C					  LEVEL RELATIVE TO THE
C					  GROUND STATE METASTABLE
C		(R*8)	FII()		: THE CONTRIBUTION FROM 
C					  EXCITATION TO POPULATE A
C					  LEVEL RELATIVE TO THE
C					  2(1)S METASTABLE.
C		(R*8)	FIII()		: THE CONTRIBUTION FROM
C					  EXCITATION TO POPULATE A
C					  LEVEL RELATIVE TO THE
C					  2(3)S METASTABLE.
C		(R*8)   GCRC1		: C-R COUPLING COEFFICIENT.
C		(R*8)   GCRC2		: C-R COUPLING COEFFICIENT.
C		(R*8)   GCRC3		: C-R COUPLING COEFFICIENT.
C		(R*8)   GCRC4		: C-R COUPLING COEFFICIENT.
C		(R*8)   GCRC5		: C-R COUPLING COEFFICIENT.
C		(R*8)   GCRC6		: C-R COUPLING COEFFICIENT.
C		(R*8)   GCRC7		: C-R COUPLING COEFFICIENT.
C		(R*8)   GCRC8		: C-R COUPLING COEFFICIENT.
C		(R*8)   GCRC9		: C-R COUPLING COEFFICIENT.
C		(R*8)	NNBNL()		: THE PRODUCT OF THE RELATIVE 
C			                  POPULATION OF A PARTICULAR
C			          	  LEVEL TO THE FIRST IONISATION 
C				  	  STAGE AND THE RECIPROCAL OF
C			          	  THE SAHA-BOLTZMANN B-FACTOR.
C		(R*8)   AVALUE		: TRANSITION PROBABILTY.
C		(R*8)	NERAY()		: ARRAY CONTAINING THE
C					  ELECTRON DENSITY (CM-3).
C		(R*8)	EBRAY()		: ARRAY CONTAINING THE BEAM 
C				          ENERGIES (eV/amu)
C		(I*4)	CHOICE		: SWITCH USED TO INDICATE IF
C					  COUPLING COEFFICIENCTS OR
C					  EMISSION DATA IS BEING
C					  HANDLED. CHOICE EQ 1 THEN
C					  COUPLING DATA. CHOICE EQ 2
C	    				  THEN EMISSION DATA.
C		(I*4)	OPTION		: SIWTCH USED TO SELECT THE
C					  ACTUAL CROSS COUPLING
C					  COEFFICIENT.
C		(I*4)	IECOUNT		: NUMBER OF BEAM ENERGIES.
C		(I*4)	INCOUNT		: NUMBER OF TARGET DENSITIES.
C		(I*4)	MAXNE		: THE MAXIMUM NUMBER OF
C					  TARGET DENSITIES
C		(I*4)	MAXTE		: MAXIMUM NUMBER OF
C					  TEMPERATURES.
C		(I*4)	MAXEB		: THE MAXIMUM NUMBER OF
C					  BEAM ENERGIES
C		(I*4)	IEREF		: INDEX TO THE REFERENCE
C					  BEAM ENERGY.
C		(I*4)	INREF		: INDEX TO THE REFERENCE
C					  TARGET DENSITY
C		(I*4)	ITREF		: INDEX TO THE REFERENCE
C					  TEMPERATURE.
C		(I*4)	TVAL		: ARRAY INDEX OF TEMPERATURE
C					  AT WHICH DATA HAS TO BE
C					  GENERATED.
C		(I*4)   RELMET		: SPECIFIES THE RELATIVE METASTABLE. 
C				    	   RELMET.EQ.1 1S2(1)S METASTABLE
C				    	   RELMET.EQ.2 1S2S(1)S METASTABLE
C				    	   RELMET.EQ.3 1S2S(3)S METASTABLE
C		(I*4)	INA()		: ARRAY CONTAINING THE INDEXES
C					  FOR THE DENSITY ARRAY.
C		(I*4)	IEA()		: ARRAY CONTAINING THE INDEXES
C					  FOR THE BEAM ENERGY ARRAY.
C		(I*4)	ITA()		: ARRAY CONTAINING THE INDEXES
C					  FOR THE TEMPERATURE ARRAY.
C		(I*4)   NSPIN		: NUMBER OF SPIN SYSTEMS.
C
C	OUTPUT :
C
C		(R*8)	ZDATA()		: ARRAY CONTAINING THE DERIVED
C					  DATA AS A FUNCTION OF BEAM
C					  ENERGY AND TARGET DENSITY
C					  FOR A SELECTED TEMPERATURE.
C		(R*8)	YDATA()		: ARRAY CONTAINING THE
C					  ELECTRON DENSITY (CM-3).
C		(R*8)	XDATA()		: ARRAY CONTAINING THE BEAM 
C				          ENERGIES (eV/amu) 
C
C	CONTACT : HARVEY ANDERSON
C		  UNIVERSITY OF STRATHCLYDE
C		  ANDERSON@PHYS.STRATH.AC.UK
C
C	DATE    : 07/5/98
C		 (FIRST VERSION)
C
C
C VERSION: 1.1						DATE: 16-03-99
C MODIFIED: RICHARD MARTIN
C		- PUT UNDER SCCS CONTROL.
C
C-----------------------------------------------------------------------
C
		
      INTEGER             CHOICE,      IEA(MAXEB),  IECOUNT,     IEREF
      INTEGER             INA(MAXNE),  INCOUNT,     INREF
      INTEGER             ITA(MAXTE),  ITREF,       MAXEB,       MAXNE
      INTEGER             MAXTE,       MULTI,       NSPIN,       OPTION
      INTEGER             RELMET,      TVAL
      REAL*8              AVALUE,      EBRAY(MAXEB)
      REAL*8              FI(MAXEB,MAXNE,MAXTE,NSPIN+1)
      REAL*8              FII(MAXEB,MAXNE,MAXTE,NSPIN+1)
      REAL*8              FIII(MAXEB,MAXNE,MAXTE,NSPIN+1)
      REAL*8              GCRC1(MAXEB,MAXNE,MAXTE)
      REAL*8              GCRC2(MAXEB,MAXNE,MAXTE)
      REAL*8              GCRC3(MAXEB,MAXNE,MAXTE)
      REAL*8              GCRC4(MAXEB,MAXNE,MAXTE)
      REAL*8              GCRC5(MAXEB,MAXNE,MAXTE)
      REAL*8              GCRC6(MAXEB,MAXNE,MAXTE)
      REAL*8              GCRC7(MAXEB,MAXNE,MAXTE)
      REAL*8              GCRC8(MAXEB,MAXNE,MAXTE)
      REAL*8              GCRC9(MAXEB,MAXNE,MAXTE), NERAY(MAXNE)
      REAL*8              NNBNL(MAXEB,MAXNE,MAXTE,NSPIN+1)
      REAL*8              XDATA(MAXEB),             YDATA(MAXNE)
      REAL*8              ZDATA(MAXEB,MAXNE)