Télécharger recvit.eso

Retour à la liste

Numérotation des lignes :

  1. C RECVIT SOURCE CHAT 05/01/13 02:46:57 5004
  2. SUBROUTINE RECVIT(ITRES,IPOS,IPTR,XVIT)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. *--------------------------------------------------------------------*
  6. * *
  7. * Recherche le déplacements généralisé de vitesse de rotation *
  8. * *
  9. * Paramètres: *
  10. * *
  11. * e ITRES table résultat issue de l'opérateur Dyne *
  12. * e IPOS position de XTEMP dans le listreel des temps *
  13. * e IPTR point de repère du mode de rotation,dans la base modale*
  14. * es XVIT valeur de la vitesse *
  15. * Auteur, date de création: *
  16. * *
  17. * Samuel DURAND le 15 Octobre 1996. *
  18. * *
  19. *--------------------------------------------------------------------*
  20. -INC CCOPTIO
  21. -INC SMCHPOI
  22. -INC SMELEME
  23. LOGICAL L0,L1
  24. * Recherche du champoint des déplacements généralisés
  25. CALL ACCTAB(ITRES,'ENTIER',IPOS,X0,' ',L0,IP0,
  26. & 'TABLE',I1,X1,' ',L1,ITDEP)
  27. CALL ACCTAB(ITDEP,'MOT',I0,X0,'VITESSE',L0,IP0,
  28. & 'CHPOINT',I1,X1,' ',L1,ICHPT)
  29. MCHPOI = ICHPT
  30. IF (MCHPOI.EQ.0) THEN
  31. * le CHPOINT des vitesses modales est nul
  32. MOTERR(1:8) = 'RCDEPL'
  33. CALL ERREUR(170)
  34. RETURN
  35. ENDIF
  36. SEGACT MCHPOI
  37. NSOU = IPCHP(/1)
  38. DO 10 ISOU = 1,NSOU
  39. MSOUPO = IPCHP(ISOU)
  40. SEGACT MSOUPO
  41. * on cherche la valeur du champ correspondant au point de
  42. * rotation
  43. MELEME = IGEOC
  44. SEGACT MELEME
  45. MPOVAL = IPOVAL
  46. SEGACT MPOVAL
  47. N2 = NUM(/2)
  48. DO 12 I = 1,N2
  49. IF (NUM(1,I).EQ.IPTR) THEN
  50. XVIT=VPOCHA(I,1)
  51. ENDIF
  52. 12 CONTINUE
  53. SEGDES MPOVAL,MELEME,MSOUPO
  54. 10 CONTINUE
  55. SEGDES MCHPOI
  56. *
  57. END
  58.  
  59.  
  60.  
  61.  
  62.  

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