Télécharger la plaquette Cast3M


fr en

Télécharger cmaxr2.eso

Retour à la liste

Numérotation des lignes : 

cmaxr2

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

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