Télécharger cmaxtb.eso

Retour à la liste

Numérotation des lignes :

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

© Cast3M 2003 - Tous droits réservés.
Mentions légales