Télécharger joiloc.eso

Retour à la liste

Numérotation des lignes :

joiloc
  1. C JOILOC SOURCE AM 11/05/03 21:17:18 6955
  2. SUBROUTINE JOILOC(REL,BPSS,RELO,P,LRE,IDIM)
  3. C-----------------------------------------------------------------------
  4. C CALCUL DU VECTEUR DEPLACEMENTS DANS LE REPERE GLOBAL LOCAL
  5. C
  6. C ENTREE
  7. C REL VECTEUR EXPRIME DANS LE REPERE GLOBAL
  8. C BPSS MATRICE DE PASSAGE (IDIM x IDIM)
  9. C SORTIE
  10. C REL VECTEUR EXPRIME DANS LE REPERE LOCAL
  11. C-----------------------------------------------------------------------
  12. IMPLICIT INTEGER(I-N)
  13. IMPLICIT REAL*8(A-H,O-Z)
  14. DIMENSION REL(LRE),BPSS(3,3),RELO(LRE),P(LRE,LRE)
  15. *
  16. CALL ZERO(P,LRE,LRE)
  17. CALL ZERO(RELO,LRE,1)
  18. *
  19. * ASSEMBLAGE DE PASSAGE P
  20. *
  21. IF(IDIM.EQ.3) THEN
  22. DO K=0,3
  23. DO I=1,3
  24. DO J=1,3
  25. P(I+3*K,J+3*K)=BPSS(I,J)
  26. ENDDO
  27. ENDDO
  28. ENDDO
  29. ELSE IF(IDIM.EQ.2) THEN
  30. DO K=0,1
  31. DO I=1,3
  32. DO J=1,3
  33. P(I+3*K,J+3*K)=BPSS(I,J)
  34. ENDDO
  35. ENDDO
  36. ENDDO
  37. ENDIF
  38. *
  39. C
  40. C REG = P*REL
  41. C
  42. DO I=1,LRE
  43. DO J=1,LRE
  44. RELO(I) = RELO(I) + P(I,J)*REL(J)
  45. ENDDO
  46. ENDDO
  47. *
  48. CALL ZERO(REL,LRE,1)
  49. C
  50. C REL = RELO
  51. C
  52. DO I=1,LRE
  53. REL(I) = RELO(I)
  54. ENDDO
  55. *
  56. RETURN
  57. END
  58.  
  59.  
  60.  

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