Esope source files

   1 : C MAXTRB    SOURCE    AM        08/02/14    21:26:42     6050           
   2 :       SUBROUTINE MAXTRB(WTRAV,WRK1,WRK5,WR12,NWA,NSTRS,NCHAIN,CMATE)
   3 : *
   4 : *      MODELE DE MAXWELL  : RECUPERATION DES VARIABLES INTERNES
   5 : *                           APRES ECOULEMENT
   6 : *
   7 :       IMPLICIT INTEGER(I-N)
   8 :       IMPLICIT REAL*8(A-H,O-Z)
   9 : *
  10 : -INC CCOPTIO
  11 : -INC SMEVOLL
  12 : -INC SMLREEL
  13 : -INC CCHAMP
  14 : *
  15 :       SEGMENT WRK1
  16 :         REAL*8 DDHOOK(LHOOK,LHOOK),SIG0(NSTRS),DEPST(NSTRS)
  17 :         REAL*8 SIGF(NSTRS),VAR0(NVARI),VARF(NVARI)
  18 :         REAL*8 DEFP(NSTRS),XCAR(ICARA)
  19 :       ENDSEGMENT
  20 : *
  21 :       SEGMENT WRK5
  22 :         REAL*8 EPIN0(NSTRS),EPINF(NSTRS),EPST0(NSTRS)
  23 :       ENDSEGMENT
  24 : *
  25 :       SEGMENT WR12
  26 :         REAL*8  EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
  27 :         REAL*8  EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
  28 :         REAL*8  EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
  29 :         REAL*8  SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
  30 :         REAL*8  SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
  31 :         REAL*8  SM8(NSTRS)
  32 :       ENDSEGMENT
  33 : *
  34 : 
  35 : 
  36 : 
  37 :       SEGMENT WTRAV
  38 :         REAL*8 DDAUX(LHOOK,LHOOK),VALMAT(NUMAT)
  39 :         REAL*8 VALCAR(NUCAR),DSIGT(NSTRS)
  40 :         REAL*8 TXR(IDIM,IDIM),DDHOMU(LHOOK,LHOOK)
  41 :         REAL*8 XLOC(IDIM,IDIM),XGLOB(IDIM,IDIM)
  42 :         REAL*8 D1HOOK(LHOOK,LHOOK),ROTHOO(LHOOK,LHOOK)
  43 :       ENDSEGMENT
  44 : *
  45 : 
  46 :       CHARACTER*8 CMATE
  47 : 
  48 : 
  49 : 
  50 :       DIMENSION NWA(9)
  51 : 
  52 : *
  53 : *    ON SE REMET DANS LES BONS AXES SI BESOIN
  54 : *
  55 : *
  56 : ******* CAS D'UNE FORMULATION UNIDIRECTIONNELLE
  57 : *
  58 :        IF(CMATE.EQ.'UNIDIREC') THEN
  59 :          ICAS=2
  59 :          CALL MAXRO2(ICAS,WTRAV,WRK1,WRK5,WR12,NCHAIN)
  60 :          IF(IERR.NE.0) RETURN
  61 :        ENDIF
  62 : 
  63 : *
  64 : *      ON CREE LES PROGS ASSOCIES AUX VARIABLES INTERNES
  65 : *
  66 :              DO JC=1,NCHAIN
  67 : * zzzzzzzz  il faudra modifier le test sur var0
  68 : 
  69 :               
  70 :                  JG=NSTRS
  71 :                  SEGINI MLREEL
  72 :                  VARF(JC+1)=DBLE(MLREEL)
  73 :                  IF(JC.EQ.1) THEN
  74 :                     DO JD=1,NSTRS
  75 :                        PROG(JD)=SM0(JD)
  76 :                     END DO
  77 :                  ELSE IF(JC.EQ.2) THEN
  78 :                     DO JD=1,NSTRS
  79 :                        PROG(JD)=SM1(JD)
  80 :                     END DO
  81 :                  ELSE IF(JC.EQ.3) THEN
  82 :                     DO JD=1,NSTRS
  83 :                        PROG(JD)=SM2(JD)
  84 :                     END DO
  85 :                  ELSE IF(JC.EQ.4) THEN
  86 :                     DO JD=1,NSTRS
  87 :                        PROG(JD)=SM3(JD)
  88 :                     END DO
  89 :                  ELSE IF(JC.EQ.5) THEN
  90 :                     DO JD=1,NSTRS
  91 :                        PROG(JD)=SM4(JD)
  92 :                     END DO
  93 :                  ELSE IF(JC.EQ.6) THEN
  94 :                     DO JD=1,NSTRS
  95 :                        PROG(JD)=SM5(JD)
  96 :                     END DO
  97 :                  ELSE IF(JC.EQ.7) THEN
  98 :                     DO JD=1,NSTRS
  99 :                        PROG(JD)=SM6(JD)
 100 :                     END DO
 101 :                  ELSE IF(JC.EQ.8) THEN
 102 :                     DO JD=1,NSTRS
 103 :                        PROG(JD)=SM7(JD)
 104 :                     END DO
 105 :                  ELSE IF(JC.EQ.9) THEN
 106 :                     DO JD=1,NSTRS
 107 :                        PROG(JD)=SM8(JD)
 108 :                     END DO
 109 :                  ENDIF
 110 :                  SEGDES MLREEL
 111 :             
 112 :              END DO
 113 :              DO JC=NCHAIN+1,9
 114 :                  JG=0
 115 :                  SEGINI MLREEL
 116 :                  VARF(JC+1)=DBLE(MLREEL)
 117 :                  SEGDES MLREEL
 118 :              END DO
 119 : *
 120 :       RETURN
 121 :       END
 122 : 
 123 :  
 124 :  
 125 :  
 126 :  
 127 :  
 128 :  
 129 :  
 130 :  
 131 :  
 132 :  
 133 :  
 134 :  
 135 :  
 136 :  
 137 :  
 138 :  
 139 :  
 140 :  
 141 :  
 142 :