Télécharger recvit.eso

Retour à la liste

Numérotation des lignes :

recvit
  1. C RECVIT SOURCE CB215821 20/11/25 13:38:48 10792
  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.  
  21. -INC PPARAM
  22. -INC CCOPTIO
  23. -INC SMCHPOI
  24. -INC SMELEME
  25. LOGICAL L0,L1
  26. * Recherche du champoint des déplacements généralisés
  27. CALL ACCTAB(ITRES,'ENTIER',IPOS,X0,' ',L0,IP0,
  28. & 'TABLE',I1,X1,' ',L1,ITDEP)
  29. CALL ACCTAB(ITDEP,'MOT',I0,X0,'VITESSE',L0,IP0,
  30. & 'CHPOINT',I1,X1,' ',L1,ICHPT)
  31. MCHPOI = ICHPT
  32. IF (MCHPOI.EQ.0) THEN
  33. * le CHPOINT des vitesses modales est nul
  34. MOTERR(1:8) = 'RCDEPL'
  35. CALL ERREUR(170)
  36. RETURN
  37. ENDIF
  38. SEGACT MCHPOI
  39. NSOU = IPCHP(/1)
  40. DO 10 ISOU = 1,NSOU
  41. MSOUPO = IPCHP(ISOU)
  42. SEGACT MSOUPO
  43. * on cherche la valeur du champ correspondant au point de
  44. * rotation
  45. MELEME = IGEOC
  46. SEGACT MELEME
  47. MPOVAL = IPOVAL
  48. SEGACT MPOVAL
  49. N2 = NUM(/2)
  50. DO 12 I = 1,N2
  51. IF (NUM(1,I).EQ.IPTR) THEN
  52. XVIT=VPOCHA(I,1)
  53. ENDIF
  54. 12 CONTINUE
  55. SEGDES MPOVAL,MELEME,MSOUPO
  56. 10 CONTINUE
  57. SEGDES MCHPOI
  58. *
  59. END
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  

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