Sources Esope | Cast3M

maxtrb
C MAXTRB    SOURCE    BP208322  17/03/01    21:17:54     9325           
      SUBROUTINE MAXTRB(WTRAV,WRK1,WRK5,WR12,NWA,NSTRS,NCHAIN,CMATE)
*
*      MODELE DE MAXWELL  : RECUPERATION DES VARIABLES INTERNES
*                           APRES ECOULEMENT
*
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
*
 
-INC PPARAM
-INC CCOPTIO
-INC SMEVOLL
-INC SMLREEL
-INC CCHAMP
*
      SEGMENT WRK1
        REAL*8 DDHOOK(LHOOK,LHOOK),SIG0(NSTRS),DEPST(NSTRS)
        REAL*8 SIGF(NSTRS),VAR0(NVARI),VARF(NVARI)
        REAL*8 DEFP(NSTRS),XCAR(ICARA)
      ENDSEGMENT
*
      SEGMENT WRK5
        REAL*8 EPIN0(NSTRS),EPINF(NSTRS),EPST0(NSTRS)
      ENDSEGMENT
*
      SEGMENT WR12
        REAL*8  EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
        REAL*8  EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
        REAL*8  EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
        REAL*8  SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
        REAL*8  SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
        REAL*8  SM8(NSTRS)
      ENDSEGMENT
*
 
 
 
      SEGMENT WTRAV
        REAL*8 DDAUX(LHOOK,LHOOK),VALMAT(NUMAT)
        REAL*8 VALCAR(NUCAR),DSIGT(NSTRS)
        REAL*8 TXR(IDIM,IDIM),DDHOMU(LHOOK,LHOOK)
        REAL*8 XLOC(3,3),XGLOB(3,3)
        REAL*8 D1HOOK(LHOOK,LHOOK),ROTHOO(LHOOK,LHOOK)
      ENDSEGMENT
*
 
      CHARACTER*8 CMATE
 
 
 
      DIMENSION NWA(9)
 
*
*    ON SE REMET DANS LES BONS AXES SI BESOIN
*
*
******* CAS D'UNE FORMULATION UNIDIRECTIONNELLE
*
       IF(CMATE.EQ.'UNIDIREC') THEN
         ICAS=2
         CALL MAXRO2(ICAS,WTRAV,WRK1,WRK5,WR12,NCHAIN)
         IF(IERR.NE.0) RETURN
       ENDIF
 
*
*      ON CREE LES PROGS ASSOCIES AUX VARIABLES INTERNES
*
             DO JC=1,NCHAIN
* zzzzzzzz  il faudra modifier le test sur var0
 
 
                 JG=NSTRS
                 SEGINI MLREEL
                 VARF(JC+1)=DBLE(MLREEL)
                 IF(JC.EQ.1) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM0(JD)
                    END DO
                 ELSE IF(JC.EQ.2) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM1(JD)
                    END DO
                 ELSE IF(JC.EQ.3) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM2(JD)
                    END DO
                 ELSE IF(JC.EQ.4) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM3(JD)
                    END DO
                 ELSE IF(JC.EQ.5) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM4(JD)
                    END DO
                 ELSE IF(JC.EQ.6) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM5(JD)
                    END DO
                 ELSE IF(JC.EQ.7) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM6(JD)
                    END DO
                 ELSE IF(JC.EQ.8) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM7(JD)
                    END DO
                 ELSE IF(JC.EQ.9) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM8(JD)
                    END DO
                 ENDIF
                 SEGDES MLREEL
 
             END DO
             DO JC=NCHAIN+1,9
                 JG=0
                 SEGINI MLREEL
                 VARF(JC+1)=DBLE(MLREEL)
                 SEGDES MLREEL
             END DO
*
      RETURN
      END