C CALCEP33  SOURCE    BR232186  12/10/04    21:15:00     7520           
      SUBROUTINE CALCEP33 (EPSIM,XGAMMT,XGAMMC,XD1,XD2,XMUM,XLAMM,
     &                     XEPS33)
C
C====&===1=========2=========3=========4=========5=========6=========7==
C Commentaires :   Subroutine permettant de calculer la deformation ep33
C                  Appelé par la routine lcgldm.eso
C
C Auteur       : B. Richard (Dr - Ing.) - CEA/DEN/DANS/DM2S/SEMT/EMSI
C====&===1=========2=========3=========4=========5=========6=========7==
C
C----DECLARATION GENERALES----------------------------------------------
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      
      REAL*8 EPSIM(3,3),VAPM(3),VEPM(3,3)
C
C-----CALCUL DE LA DEFORMATION HORS PLAN
C
      XTRA = EPSIM(1,1) + EPSIM(2,2)

      CALL CALCGM(XTRA,XGAMMT,XGAMMC,XGM)

      CALL JACOB4(EPSIM,2,VAPM,VEPM)

      CALL CALCXIM (XTRA,XGAMMT,XGAMMC,XD1,XD2,XIM)

      XEPS33 = -1.0D0*(XLAMM*XIM)/(2.0D0*XMUM+XLAMM*XIM)*XTRA

      RETURN
      END