c*********************************************************************

	REAL (kind=dp) FUNCTION rate(cr,tr)

c routine PRIVATE du module mod_nuc
	
c calcul du tau de réaction suivant fkthr Theoretical statistical rate calculated
c with SMOKER code rate=EXP(a1+a2/t9+a3/t913+a4*t913+a5*t9+a6*t953+a7*lnt9)
c Stellar Reaction Rates by F. Thielemann, et al.
c Reaclib Data Tables of reaction rates based on F. Thielemann, et al., Adv. Nucl.
c Astro., 525 (1987). 1991 updated version (updated to Z=46 in 1995 by Ch.
c Freiburghaus). Downloaded from quasar.physik.unibas.ch on 1-18-1995 (RBF).

c entrées:
c	cr: table des coefficients
c	tr: table des t9
	
c Auteur: P.Morel, Département Cassiopée, O.C.A., CESAM2k

c-------------------------------------------------------------------

	USE mod_kind

	IMPLICIT NONE

	REAL (kind=dp), INTENT(in), DIMENSION(7) :: cr
	REAL (kind=dp), INTENT(in), DIMENSION(5) :: tr	

c---------------------------------------------------------------

2000	FORMAT(8es10.3)

c LOG(rate)=a1+a2/t9+a3/t913+a4*t913+a5*t9+a6*t953+a7*lnt9		
	rate=cr(1)+cr(2)/tr(1)+cr(3)/tr(2)+cr(4)*tr(2)+cr(5)*tr(3)
	1 +cr(6)*tr(4)+cr(7)*tr(5)
c	WRITE(*,2000)a ; WRITE(*,2000)t ; WRITE(*,2000)rate
 	
 	RETURN
 	
 	END FUNCTION rate