C CLEBL3    SOURCE    PV        17/12/08    21:15:48     9660           
      SUBROUTINE CLEBL3(iwrk52,iele,igau)
C
C
C
C      Modele de Leblond et Devaux
C        descente au niveau des points de Gauss
C          06/97
C
C     SOURCE BRIOTTET
C
C
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C

-INC PPARAM
-INC CCOPTIO
-INC SMCHAML
-INC CCHAMP
-INC DECHE
C
      REAL*8 VWRK1(7),VWRK2(7)

           wrk52 = iwrk52
* DTPS fourni dans caracteristiques
           DTPS = valma0(14)
 10   continue
C
           T0=ture0(1)
C
           if ((tempf-temp0).gt.0) then
             TPOINT0=(turef(1) - ture0(1))/(tempf-temp0)
           else
             TPOINT0 = 0.
           endif
C
           ZA0=valma0(9)
C
           ZF0=valma0(10)
C
           ZB0=valma0(11)
C
           ZM0=valma0(12)
C
           VMS0=valma0(13)
C
            VWRK1(1)=T0
            VWRK1(2)=TPOINT0
            VWRK1(3)=ZA0
            VWRK1(4)=ZF0
            VWRK1(5)=ZB0
            VWRK1(6)=ZM0
            VWRK1(7)=VMS0
C
C     CALCUL DU NOUVEAU CHAMP DE PHASES
C
           CALL CLEBL4(iwrk52,VWRK1,DTPS,VWRK2)
C
C sortie
           xmatf(9)=VWRK2(3)
C
           xmatf(10)=VWRK2(4)
C
           xmatf(11)=VWRK2(5)
C
           xmatf(12)=VWRK2(6)
C
           xmatf(13)=VWRK2(7)
C
       do ii = 1,8
         xmatf(ii) = xmat(ii)
       enddo
C
      RETURN
      END