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 :
© Cast3M 2003 - All rights reserved.
Disclaimer