Esope source files

   1 : C CMAXTB    SOURCE    KICH      09/09/11    21:16:24     6505           
   2 :       SUBROUTINE CMAXTB(WRK52,WRK53,WRK54,WR12,NWA,NCHAIN)
   3 : C MAXTRB    SOURCE    AM        00/12/13    21:38:55     4045           
   4 : *      SUBROUTINE MAXTRB(WTRAV,WRK1,WRK5,WR12,NWA,NSTRS,NCHAIN,CMATE)
   5 : *
   6 : *      MODELE DE MAXWELL  : RECUPERATION DES VARIABLES INTERNES
   7 : *                           APRES ECOULEMENT
   8 : *
   9 :       IMPLICIT INTEGER(I-N)
  10 :       IMPLICIT REAL*8(A-H,O-Z)
  11 : *
  12 : -INC CCOPTIO
  13 : -INC SMEVOLL
  14 : -INC SMLREEL
  15 : -INC CCHAMP
  16 : -INC DECHE
  17 : *
  18 : *
  19 :       SEGMENT WR12
  20 :         REAL*8  EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
  21 :         REAL*8  EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
  22 :         REAL*8  EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
  23 :         REAL*8  SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
  24 :         REAL*8  SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
  25 :         REAL*8  SM8(NSTRS)
  26 :       ENDSEGMENT
  27 : *
  28 : 
  29 :       DIMENSION NWA(9)
  30 : 
  31 : *
  32 : *    ON SE REMET DANS LES BONS AXES SI BESOIN
  33 : *
  34 : *
  35 : ******* CAS D'UNE FORMULATION UNIDIRECTIONNELLE
  36 : *
  37 :        IF(CMATE.EQ.'UNIDIREC') THEN
  38 :          ICAS=2
  38 :              CALL CMAXR2(wrk52,wrk53,wrk54,ICAS,WR12,NCHAIN)
  39 :          IF(IERR.NE.0) RETURN
  40 :        ENDIF
  41 : 
  42 : *
  43 : *      ON CREE LES PROGS ASSOCIES AUX VARIABLES INTERNES
  44 : *
  45 :              DO JC=1,NCHAIN
  46 : * zzzzzzzz  il faudra modifier le test sur var0
  47 : 
  48 :               
  49 :                  JG=NSTRS
  50 :                  SEGINI MLREEL
  51 :                  VARF(JC+1)=DBLE(MLREEL)
  52 :                  IF(JC.EQ.1) THEN
  53 :                     DO JD=1,NSTRS
  54 :                        PROG(JD)=SM0(JD)
  55 :                     END DO
  56 :                  ELSE IF(JC.EQ.2) THEN
  57 :                     DO JD=1,NSTRS
  58 :                        PROG(JD)=SM1(JD)
  59 :                     END DO
  60 :                  ELSE IF(JC.EQ.3) THEN
  61 :                     DO JD=1,NSTRS
  62 :                        PROG(JD)=SM2(JD)
  63 :                     END DO
  64 :                  ELSE IF(JC.EQ.4) THEN
  65 :                     DO JD=1,NSTRS
  66 :                        PROG(JD)=SM3(JD)
  67 :                     END DO
  68 :                  ELSE IF(JC.EQ.5) THEN
  69 :                     DO JD=1,NSTRS
  70 :                        PROG(JD)=SM4(JD)
  71 :                     END DO
  72 :                  ELSE IF(JC.EQ.6) THEN
  73 :                     DO JD=1,NSTRS
  74 :                        PROG(JD)=SM5(JD)
  75 :                     END DO
  76 :                  ELSE IF(JC.EQ.7) THEN
  77 :                     DO JD=1,NSTRS
  78 :                        PROG(JD)=SM6(JD)
  79 :                     END DO
  80 :                  ELSE IF(JC.EQ.8) THEN
  81 :                     DO JD=1,NSTRS
  82 :                        PROG(JD)=SM7(JD)
  83 :                     END DO
  84 :                  ELSE IF(JC.EQ.9) THEN
  85 :                     DO JD=1,NSTRS
  86 :                        PROG(JD)=SM8(JD)
  87 :                     END DO
  88 :                  ENDIF
  89 :                  SEGDES MLREEL
  90 :             
  91 :              END DO
  92 :              DO JC=NCHAIN+1,9
  93 :                  JG=0
  94 :                  SEGINI MLREEL
  95 :                  VARF(JC+1)=DBLE(MLREEL)
  96 :                  SEGDES MLREEL
  97 :              END DO
  98 : *
  99 :       RETURN
 100 :       END
 101 : 
 102 :  
 103 :  
 104 :  
 105 :  
 106 :  
 107 :  
 108 :  
 109 :  
 110 :  
 111 :  
 112 :  
 113 :  
 114 :  
 115 :  
 116 :  
 117 :  
 118 :  
 119 :  
 120 :  
 121 :  
 122 :  
 123 :  
 124 :  
 125 :  
 126 :