Télécharger insi.eso

Retour à la liste

Numérotation des lignes :

  1. C INSI SOURCE CHAT 05/01/13 00:39:12 5004
  2. C$$$$ INSI
  3. C INSI SOURCE ISPRA 90/06/12
  4. SUBROUTINE INSI
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. C=======================================================================
  8. C OPERATEUR INSI
  9. C
  10. C D*EVOLUTION V*EVOLUTION = INSI A*EVOLUTION (METH*MOT)
  11. C
  12. C=======================================================================
  13. C PROGRAMMEUR : P.P.
  14. C=======================================================================
  15. C
  16. CHARACTER *72 TI
  17. CHARACTER*12 MOTX,MOTY
  18. C
  19. PARAMETER (NMOCLE=2)
  20. CHARACTER*4 MOTCLE(NMOCLE)
  21. C
  22. -INC CCOPTIO
  23. -INC SMEVOLL
  24. -INC SMLREEL
  25. C
  26. POINTEUR IACCE.MLREEL,IVITE.MLREEL,IDEPL.MLREEL,ITEMP.MLREEL
  27. POINTEUR JACCE.MEVOLL,JVITE.MEVOLL,JDEPL.MEVOLL
  28. POINTEUR KACCE.KEVOLL,KVITE.KEVOLL,KDEPL.KEVOLL
  29. C
  30. C 1) LECTURE DES DONNEES GIBIANE
  31. C
  32. C 1.1) LISTE DES MOTS CLEF
  33. C
  34. DATA MOTCLE/'SIMP','LINE'/
  35. C
  36. C 1.2) DEFAUTS
  37. C
  38. IMETH=1
  39. C
  40. C 1.3) LECTURE DE L'OBJET EVOLUTIO CONTENANT L'ACCELERATION
  41. C
  42. CALL LIROBJ('EVOLUTIO',JACCE,1,IRET)
  43. IF(IRET.EQ.0) RETURN
  44. C
  45. C 1.4) LECTURE DU MOT-CLEF
  46. C (OPTIONEL)
  47. C
  48. CALL LIRMOT(MOTCLE,NMOCLE,IVAL,0)
  49. C
  50. IF(IVAL.NE.0)THEN
  51. IMETH=IVAL
  52. ENDIF
  53. C
  54. C
  55. C 2) INITIALISATION DIVERSES
  56. C
  57. SEGACT, JACCE
  58. N=JACCE.IEVOLL(/1)
  59. C
  60. TI=JACCE.IEVTEX
  61. SEGINI, JVITE
  62. SEGINI, JDEPL
  63. JVITE.IEVTEX='Integration ('//MOTCLE(IMETH)//') de:'//TI(1:50)
  64. JDEPL.IEVTEX='Integration ('//MOTCLE(IMETH)//') de:'//TI(1:50)
  65. C
  66. C 3) LOOP DE CALCUL
  67. C
  68. DO 100 IE1=1,N
  69. C
  70. C 3.1) INITIALISATION ET ACTIVATION
  71. C
  72. KACCE=JACCE.IEVOLL(IE1)
  73. SEGACT, KACCE
  74. C
  75. ITEMP=KACCE.IPROGX
  76. IACCE=KACCE.IPROGY
  77. ICOUL=KACCE.NUMEVX
  78. SEGDES, KACCE
  79. C
  80. SEGACT, IACCE
  81. SEGACT, ITEMP
  82. NPT=ITEMP.PROG(/1)
  83. JG=NPT
  84. C
  85. SEGINI, KVITE
  86. JVITE.IEVOLL(IE1)=KVITE
  87. C
  88. WRITE(TI,'(A17,1X,A6,1X,I2)')'Signal en vitesse',
  89. > 'numero',IE1
  90. WRITE(MOTY,'(9HVitesse ,1X,I2)')IE1
  91. SEGINI, IVITE
  92. C
  93. KVITE.KEVTEX=TI
  94. KVITE.NUMEVX=ICOUL
  95. KVITE.NUMEVY='REEL'
  96. KVITE.TYPX='LISTREEL'
  97. KVITE.IPROGX=ITEMP
  98. KVITE.TYPY='LISTREEL'
  99. KVITE.IPROGY=IVITE
  100. KVITE.NOMEVY=MOTY(1:12)
  101. SEGDES, KVITE
  102. C
  103. SEGINI, KDEPL
  104. JDEPL.IEVOLL(IE1)=KDEPL
  105. C
  106. WRITE(TI,'(A21,1X,A6,1X,I2)')'Signal en deplacement',
  107. > 'numero',IE1
  108. WRITE(MOTY,'(9HDeplacem.,1X,I2)')IE1
  109. SEGINI, IDEPL
  110. C
  111. KDEPL.KEVTEX=TI
  112. KDEPL.NUMEVX=ICOUL
  113. KDEPL.NUMEVY='REEL'
  114. KDEPL.TYPX='LISTREEL'
  115. KDEPL.IPROGX=ITEMP
  116. KDEPL.TYPY='LISTREEL'
  117. KDEPL.IPROGY=IDEPL
  118. KDEPL.NOMEVY=MOTY(1:12)
  119. SEGDES, KDEPL
  120. C
  121. C 3.2) CALCUL DE V(I) ET D(I)
  122. C
  123. C 3.2.1) METHODE SIMPLIFIEE
  124. C
  125. IF(IMETH.EQ.1)THEN
  126. IVITE.PROG(1)=0.D0
  127. IDEPL.PROG(1)=0.D0
  128. DO 20 IE2=2,NPT
  129. DT=ITEMP.PROG(IE2)-ITEMP.PROG(IE2-1)
  130. IVITE.PROG(IE2)=IVITE.PROG(IE2-1)
  131. > +DT*(IACCE.PROG(IE2-1)+IACCE.PROG(IE2))/2
  132. IDEPL.PROG(IE2)=IDEPL.PROG(IE2-1)
  133. > +DT*(IVITE.PROG(IE2-1)+IVITE.PROG(IE2))/2
  134. 20 CONTINUE
  135. ENDIF
  136. C
  137. C 3.2.2) METHODE LINEAIRE
  138. C
  139. IF(IMETH.EQ.2)THEN
  140. IVITE.PROG(1)=0.D0
  141. IDEPL.PROG(1)=0.D0
  142. DO 25 IE2=2,NPT
  143. DT=ITEMP.PROG(IE2)-ITEMP.PROG(IE2-1)
  144. IVITE.PROG(IE2)=IVITE.PROG(IE2-1)
  145. > +DT*(IACCE.PROG(IE2-1)+IACCE.PROG(IE2))/2
  146. IDEPL.PROG(IE2)=IDEPL.PROG(IE2-1)
  147. > +DT*(IVITE.PROG(IE2-1)
  148. > +DT*(IACCE.PROG(IE2)+2*IACCE.PROG(IE2-1))/6)
  149. 25 CONTINUE
  150. ENDIF
  151. C
  152. C 3.3) FIN ACTIVATION
  153. C
  154. SEGDES, IACCE
  155. SEGDES, ITEMP
  156. SEGDES, IVITE
  157. SEGDES, IDEPL
  158. C
  159. 100 CONTINUE
  160. C
  161. C
  162. C
  163. SEGDES, JACCE
  164. SEGDES, JVITE
  165. SEGDES, JDEPL
  166. C
  167. C 4) RETOUR A GIBIANE
  168. C
  169. CALL ECROBJ('EVOLUTIO',JVITE)
  170. CALL ECROBJ('EVOLUTIO',JDEPL)
  171. C
  172. C
  173. 666 CONTINUE
  174. RETURN
  175. END
  176.  
  177.  
  178.  
  179.  
  180.  

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