Télécharger ricoq0.eso

Retour à la liste

Numérotation des lignes :

  1. C RICOQ0 SOURCE CHAT 05/01/13 02:56:14 5004
  2. SUBROUTINE RICOQ0(RIG,X13,X23,Y13,Y23,FAC)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. DIMENSION RIG(18,18),VEC(18)
  6. *
  7. S=(X13*Y23-X23*Y13)*2.D0
  8. FAC4=FAC*4.D0
  9. DO 100 I=1,18
  10. 100 VEC(I)=0.
  11. VEC(1)=-X23/S
  12. VEC(2)=-Y23/S
  13. VEC(7)=X13/S
  14. VEC(8)=Y13/S
  15. VEC(13)=(X23-X13)/S
  16. VEC(14)=(Y23-Y13)/S
  17. DO 102 I=1,3
  18. I6=6*I
  19. VEC(6) =0.D0
  20. VEC(12)=0.D0
  21. VEC(18)=0.D0
  22. VEC(I6)=1.D0
  23. DO 103 J=1,18
  24. DO 103 K=1,18
  25. RIG(J,K)=RIG(J,K)+VEC(J)*VEC(K)*FAC4
  26. 103 CONTINUE
  27. 102 CONTINUE
  28. *
  29. RETURN
  30. END
  31.  
  32.  

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