Télécharger la plaquette Cast3M


fr en

Télécharger ricoq.eso

Retour à la liste

Numérotation des lignes : 

ricoq

  1. C RICOQ     SOURCE    CB215821  26/08/24    21:18:16     12622          
  2.       SUBROUTINE RICOQ(RE,XA,HM,E,RIGL,REP,REM,P,XE,YE,ZE,XX,YY,RAI)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8(A-H,O-Z)
  5.       DIMENSION            E(3,3),RIGL(18,18),REP(9,9),REM(6,6),
  6.      1P(3,3),XE(3),YE(3),ZE(3),XX(3,3),YY(3,3)
  7.       DIMENSION RE(18,18),XA(3,3)
  8.       DATA EPSI/1.E-9/
  9. *
  10.       DO102I=1,18
  11.       DO9J=1,18
  12.       RE(I,J)=0.D0
  13.       RIGL(I,J)=0.D0
  14.     9 CONTINUE
  15.   102 CONTINUE
  16.       DO 1 I=1,3
  17.       XE(I)=XA(1,I)
  18.       YE(I)=XA(2,I)
  19.       ZE(I)=XA(3,I)
  20.     1 CONTINUE
  21.       CALL       PASSA(XE,YE,ZE,    P,X13,X23,Y13,Y23)
  22.       S=0.5*Y13*(X13-X23)
  23. *     SPECIAL    RIGI      PSI Z
  24.       RKZ= RAI*HM*S*EPSI
  25.       CALL RICOQ0(RIGL,X13,X23,Y13,Y23,RKZ)
  26. *
  27.       CALL RIMT6(REM,E,RAI,X13,X23,Y13,HM,S)
  28.       DO103I=1,6
  29.       DO100J=I,6
  30.       REM(J,I)=REM(I,J)
  31.   100 CONTINUE
  32.   103 CONTINUE
  33.       CALL RIPP9(REP,E,RAI,X13,X23,Y13,HM,S)
  34.       DO104I=1,9
  35.       DO101J=I,9
  36.       REP(J,I)=REP(I,J)
  37.   101 CONTINUE
  38.   104 CONTINUE
  39.       DO105K=1,3
  40.       K1=K-1
  41.       DO2L=K,3
  42.       L1=L-1
  43.       DO 106 I=1,2
  44.       M=6*K1+I
  45.       ML=2*K1+I
  46.       DO 3 J=1,2
  47.       N=6*L1+J
  48.       NL=2*L1+J
  49.       RIGL(M,N)=RIGL(M,N)+REM(ML,NL)
  50.     3 CONTINUE
  51.   106 CONTINUE
  52.       DO 107 I=1,3
  53.       M=6*K1+2+I
  54.       ML=3*K1+I
  55.       DO 4 J=1,3
  56.       N=6*L1+J+2
  57.       NL=3*L1+J
  58.       RIGL(M,N)=RIGL(M,N)+REP(ML,NL)
  59.     4 CONTINUE
  60.   107 CONTINUE
  61.     2 CONTINUE
  62.   105 CONTINUE
  63.       DO 108 K=1,6
  64.       K1=3*(K-1)
  65.       DO6L=K,6
  66.       L1=3*(L-1)
  67.       DO 109 I=1,3
  68.       IM=I+K1
  69.       DO 7 J=1,3
  70.       JM=J+L1
  71.       XX(I,J)=RIGL(IM,JM)
  72.     7 CONTINUE
  73.   109 CONTINUE
  74.       CALL PRODT(YY,XX,P,3,3)
  75.       DO 110 I=1,3
  76.       IM=I+K1
  77.       DO 8 J=1,3
  78.       JM=J+L1
  79.       RE(IM,JM)=YY(I,J)
  80.     8 CONTINUE
  81.   110 CONTINUE
  82.     6 CONTINUE
  83.   108 CONTINUE
  84.       DO 111 I=1,18
  85.       DO 10 J=1,I
  86.       RE(I,J)=RE(J,I)
  87.    10 CONTINUE
  88.   111 CONTINUE
  89. *     DO 555 I=1,18
  90. *     WRITE(6,77882) I
  91. *77882 FORMAT('0  LIGNE ',I4//)
  92. *     WRITE(6,77881) ( RE(I,J),J=1,18)
  93. *77881 FORMAT(1X,6(1X,1PE12.5)/)
  94. *555   CONTINUE
  95.       RETURN
  96.       END
  97.  
  98.  
  99.  

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