Télécharger pouvec.eso

Retour à la liste

Numérotation des lignes :

  1. C POUVEC SOURCE CHAT 05/01/13 02:20:50 5004
  2. SUBROUTINE POUVEC(U,V,P,ICOD)
  3. C-----------------------------------------------------------------------
  4. C ROUTINE DE PASSAGE DU VECTEUR DEPLACEMENT GLOBAL AU LOCAL
  5. C OU DU VECTEUR FORCE LOCAL AU GLOBAL
  6. C
  7. C ENTREE
  8. C ------ U DEPLACEMENTS GLOBAUX AUX 2 NOEUDS
  9. C OU FORCES LOCALES AUX 2 NOEUDS
  10. C P MATRICE DE PASSAGE
  11. C ICOD 1 POUR Q-LOC =P*Q-GLOB
  12. C 2 POUR F-GLOB=PT*F-LOC
  13. C
  14. C SORTIE
  15. C ------ V DEPLACEMENTS LOCAUX AUX 2 NOEUDS
  16. C OU FORCES GLOBALES AUX 2 NOEUDS
  17. C-----------------------------------------------------------------------
  18. IMPLICIT INTEGER(I-N)
  19. IMPLICIT REAL*8(A-H,O-Z)
  20. DIMENSION U(1),V(1),P(3,3)
  21. C
  22. GO TO (100,200),ICOD
  23. 100 CONTINUE
  24. DO 6 K=1,4
  25. IDEC=(K-1)*3
  26. DO 6 J=1,3
  27. JJ=J+IDEC
  28. V( JJ)=0.D00
  29. DO 6 I=1,3
  30. II=I+IDEC
  31. V(JJ)=V(JJ)+P(J,I)*U(II)
  32. 6 CONTINUE
  33. RETURN
  34. 200 CONTINUE
  35. DO 7 K=1,4
  36. IDEC=(K-1)*3
  37. DO 7 J=1,3
  38. JJ=J+IDEC
  39. V(JJ)=0.D00
  40. DO 7 I=1,3
  41. II=I+IDEC
  42. V(JJ)=V(JJ)+P(I,J)*U(II)
  43. 7 CONTINUE
  44. RETURN
  45. END
  46.  
  47.  

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