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

	SUBROUTINE corr_bol_lejeune(t,fe,tfe,dtfedt,dtfedfe)

c subroutine PUBLIC du module mod_photo

c corrections bolometriques de Lejeune T., Cuisinier F., Buser R.,
c 1998, A&AS 130, 65

c adaptation P. Morel, Laboratoire Lagrange, O.C.A., CESAM2k
c interpolation en Teff,[Fe/H]

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

	USE mod_donnees, ONLY : nom_chemin
	USE mod_kind
	USE mod_numerique, ONLY : bsp2dn

	IMPLICIT NONE

	INTEGER, PARAMETER :: nt=37, nfe=10, nf=11, mt=2, mfe=2

	REAL (kind=dp), INTENT(in) :: t, fe
	REAL (kind=dp), INTENT(out), DIMENSION(nf) :: tfe, dtfedt,
	1 dtfedfe

 	REAL (kind=dp), SAVE, DIMENSION(nf,nt,nfe) :: data

	INTEGER :: i
	REAL (kind=dp), PARAMETER, DIMENSION(nt) :: t_eff=(/
	1 2.0d3, 2.2d3, 2.5d3, 2.8d3, 3.0d3, 3.2d3, 3.35d3,
	2 (3.5d3+2.5d2*i,i=0,nt-11), 10.5d3, 11.0d3, 11.5d3 /)
	REAL (kind=dp), PARAMETER, DIMENSION(nfe) :: fesh=(/
	1 (-3.5d0+0.5d0*i,i=0,nfe-1) /)

	REAL (kind=dp), SAVE, DIMENSION(nt+mt) :: t_efft
	REAL (kind=dp), SAVE, DIMENSION(nfe+mfe) :: fesht

	REAL(kind=dp) :: Teff, U_B, B_V, V_I, V_K, R_I, J_H, H_K, J_K,
	1 K_L, BcV, logg, Fe_H

	INTEGER, SAVE :: knott, knotfe, lt=mt, lfe=mfe
	INTEGER :: j

	LOGICAL :: init=.TRUE.

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

2000	FORMAT(8es10.3)

c	write(*,2000)t_eff ; write(*,2000)fesh ; pause

	IF(init)THEN
	 init=.FALSE.
	 OPEN(unit=1,file=TRIM(nom_chemin)//'cor_bol_lejeune.data',form='formatted')

	 DO j=1,nfe
	  DO i=1,nt
           READ(15,11,end=2)Teff,U_B,B_V,V_I,V_K,R_I,J_H,H_K,J_K,K_L,
	1  BcV,logg,Fe_H
11	   FORMAT(2x,F6.0,F7.3,1x,F7.3,1x,F7.3,1x,F7.3,1x,F7.3,1x,F7.3,
	1  1x,F7.3,1x,F7.3,1x,F7.3,1x,F7.3,1x,F7.3,1x,F4.1)
c          write(*,11)Teff,U_B,B_V,V_I,V_K,R_I,J_H,H_K,J_K,K_L,BcV,logg,Fe_H
           DATA(:,i,j)=(/ U_B,B_V,V_I,V_K,R_I,J_H,H_K,J_K,K_L,BcV,logg /)
          ENDDO
	 ENDDO
2	 CLOSE(unit=15)

	 CALL bsp2dn(nf,data,t_eff,t_efft,nt,knott,mt,lt,fesh,fesht,nfe,knotfe,mfe,lfe,
	1 .FALSE.,t_eff(1),fesh(1),tfe,dtfedt,dtfedfe)
c	 write(*,2000)tfe ; write(*,2000)dtfedt	; write(*,2000)dtfedfe

	 PRINT*
	 PRINT*,'corr.bol. BcV, log g, indices couleur, f(Teff,[Fe/H])'
	 PRINT*,'Lejeune T., Cuisinier F., Buser R., 1998, A&AS 130, 65'
	 PRINT*

	ENDIF

	CALL bsp2dn(nf,data,t_eff,t_efft,nt,knott,mt,lt,fesh,fesht,nfe,knotfe,mfe,lfe,
	1 .TRUE.,t,fe,tfe,dtfedt,dtfedfe)

	RETURN

	END SUBROUTINE corr_bol_lejeune