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