Télécharger cmaxr2.eso

Retour à la liste

Numérotation des lignes :

  1. C CMAXR2 SOURCE BP208322 17/03/01 21:15:44 9325
  2. SUBROUTINE CMAXR2(wrk52,wrk53,wrk54,ICAS,WR12,NCHAIN)
  3. C MAXRO2 SOURCE AM 00/12/13 21:38:11 4045
  4. * SUBROUTINE MAXRO2(ICAS,WTRAV,WRK1,WRK5,WR12,NCHAIN)
  5. *
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. *
  9. -INC CCOPTIO
  10. -INC SMEVOLL
  11. -INC SMLREEL
  12. -INC CCHAMP
  13. -INC DECHE
  14. *
  15. SEGMENT WR12
  16. REAL*8 EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
  17. REAL*8 EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
  18. REAL*8 EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
  19. REAL*8 SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
  20. REAL*8 SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
  21. REAL*8 SM8(NSTRS)
  22. ENDSEGMENT
  23. *
  24. SEGMENT WTRA2
  25. REAL*8 A(NDIM,NDIM)
  26. REAL*8 RR(NDIM,NDIM),RT(NDIM,NDIM),RTRA(NDIM,NDIM)
  27. ENDSEGMENT
  28. *
  29. IDIM=XGLOB(/1)
  30. NDIM=IDIM
  31. IF(IFOUR.EQ.1) NDIM=IDIM+1
  32.  
  33. NSTRS=SIG0(/1)
  34.  
  35. SEGINI WTRA2
  36.  
  37. * ICAS = 1 : RT A R
  38. * ICAS = 2 : R A RT
  39.  
  40. DO I=1,IDIM
  41. DO J=1,IDIM
  42. RR(I,J)=XGLOB(I,J)
  43. ENDDO
  44. ENDDO
  45. IF(IDIM.EQ.2.AND.IFOUR.EQ.1) RR(3,3)=1.D0
  46. CALL TRSPOD(RR,NDIM,NDIM,RT)
  47. *
  48. * ROTATION DES TENSEURS
  49. *
  50. IDEFO=0
  51. CALL MAXRO3(ICAS,IDEFO,SM0,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  52. CALL MAXRO3(ICAS,IDEFO,SM1,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  53. CALL MAXRO3(ICAS,IDEFO,SM2,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  54. CALL MAXRO3(ICAS,IDEFO,SM3,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  55. CALL MAXRO3(ICAS,IDEFO,SM4,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  56. IF(NCHAIN.GE.6) THEN
  57. CALL MAXRO3(ICAS,IDEFO,SM5,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  58. ENDIF
  59. IF(NCHAIN.GE.7) THEN
  60. CALL MAXRO3(ICAS,IDEFO,SM6,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  61. ENDIF
  62. IF(NCHAIN.GE.8) THEN
  63. CALL MAXRO3(ICAS,IDEFO,SM7,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  64. ENDIF
  65. IF(NCHAIN.GE.9) THEN
  66. CALL MAXRO3(ICAS,IDEFO,SM8,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  67. ENDIF
  68.  
  69.  
  70. IF(ICAS.EQ.1) THEN
  71.  
  72. CALL MAXRO3(ICAS,IDEFO,SIG0,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  73. IDEFO=1
  74. CALL MAXRO3(ICAS,IDEFO,DEPST,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  75. CALL MAXRO3(ICAS,IDEFO,EPIN0,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  76. *
  77.  
  78. ELSE IF(ICAS.EQ.2) THEN
  79.  
  80. CALL MAXRO3(ICAS,IDEFO,SIGF,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  81. IDEFO=1
  82. CALL MAXRO3(ICAS,IDEFO,DEFP ,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  83. CALL MAXRO3(ICAS,IDEFO,EPINF,A,RR,RT,RTRA,NDIM,IDIM,IFOUR)
  84.  
  85. ENDIF
  86.  
  87. SEGSUP WTRA2
  88.  
  89. RETURN
  90. END
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  

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