Télécharger dyne18.eso

Retour à la liste

Numérotation des lignes :

  1. C DYNE18 SOURCE CHAT 05/01/12 23:14:53 5004
  2. SUBROUTINE DYNE18(ICHP,KTQ,I1,I2,KCPR)
  3. IMPLICIT INTEGER(I-N)
  4. IMPLICIT REAL*8(A-H,O-Z)
  5. *--------------------------------------------------------------------*
  6. * *
  7. * Op{rateur DYNE : algorithme de Fu - de Vogelaere *
  8. * ________________________________________________ *
  9. * *
  10. * Remplissage du tableau des variables g{n{ralis{es. *
  11. * *
  12. * Param}tres: *
  13. * *
  14. * e ICHP : CHPOINT dans lequel on va chercher les valeurs *
  15. * es ITQ : Segment des variables de mouvement g{n{ralis{es *
  16. * e I1 : Indique si l'on travaille sur les d{placements, *
  17. * sur les vitesses, ou les acc{l{rations *
  18. * e I2 : Indice du pas *
  19. * e KCPR : Segment des points *
  20. * *
  21. * Auteur, date de cr{ation: *
  22. * *
  23. * Lionel VIVAN, le 26 Juillet 1989. *
  24. * *
  25. *--------------------------------------------------------------------*
  26. -INC CCOPTIO
  27. -INC SMCHPOI
  28. -INC SMCOORD
  29. -INC SMELEME
  30. *
  31. SEGMENT,ICPR(XCOOR(/1)/(IDIM+1))
  32. SEGMENT,ICOR(NP1)
  33. SEGMENT,MTQ
  34. REAL*8 Q1(NA1,4),Q2(NA1,4),Q3(NA1,4)
  35. REAL*8 WEXT(NA1,2),WINT(NA1,2)
  36. ENDSEGMENT
  37. *
  38. MTQ = KTQ
  39. ICPR = KCPR
  40. *
  41. MCHPOI = ICHP
  42. SEGACT,MCHPOI
  43. NSOUPO = IPCHP(/1)
  44. DO 10 II1=1,NSOUPO
  45. MSOUPO = IPCHP(II1)
  46. SEGACT,MSOUPO
  47. MELEME = IGEOC
  48. SEGACT,MELEME
  49. NP = NUM(/2)
  50. NP1 = NP
  51. SEGINI,ICOR
  52. DO 2 IP = 1,NP
  53. KNOE = ICPR(NUM(1,IP))
  54. IF (KNOE.NE.0) ICOR(IP) = KNOE
  55. 2 CONTINUE
  56. * end do
  57. SEGDES,MELEME
  58. NC = NOCOMP(/2)
  59. MPOVAL = IPOVAL
  60. SEGACT,MPOVAL
  61. N = VPOCHA(/1)
  62. IF (I1.EQ.1) THEN
  63. DO 20 J=1,N
  64. DO 20 K=1,NC
  65. IA = ICOR(J)
  66. IF (IA.NE.0) Q1(IA,I2) = VPOCHA(J,K)
  67. 20 CONTINUE
  68.  
  69. ELSE IF (I1.EQ.2) THEN
  70. DO 30 J=1,N
  71. DO 30 K=1,NC
  72. IA = ICOR(J)
  73. IF (IA.NE.0) Q2(IA,I2) = VPOCHA(J,K)
  74. 30 CONTINUE
  75.  
  76. ELSE IF (I1.EQ.3) THEN
  77. DO 40 J=1,N
  78. DO 40 K=1,NC
  79. IA = ICOR(J)
  80. IF (IA.NE.0) Q3(IA,I2) = VPOCHA(J,K)
  81. 40 CONTINUE
  82.  
  83. ELSE IF (I1.EQ.4) THEN
  84. DO 50 J=1,N
  85. DO 50 K=1,NC
  86. IA = ICOR(J)
  87. IF (IA.NE.0) WEXT(IA,I2) = VPOCHA(J,K)
  88. 50 CONTINUE
  89.  
  90. ELSE IF (I1.EQ.5) THEN
  91. DO 60 J=1,N
  92. DO 60 K=1,NC
  93. IA = ICOR(J)
  94. IF (IA.NE.0) WINT(IA,I2) = VPOCHA(J,K)
  95. 60 CONTINUE
  96.  
  97. ENDIF
  98. SEGDES,MPOVAL,MSOUPO
  99. SEGSUP,ICOR
  100. 10 CONTINUE
  101. * end do
  102. SEGDES,MCHPOI
  103. *
  104. END
  105.  
  106.  
  107.  

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