OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s10deri3.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| s10deri3 ../starter/source/elements/solid/solide10/s10deri3.F
25!||--- called by ------------------------------------------------------
26!|| inirig_mat ../starter/source/elements/initia/inirig_mat.F
27!|| s10init3 ../starter/source/elements/solid/solide10/s10init3.F
28!||--- calls -----------------------------------------------------
29!|| ancmsg ../starter/source/output/message/message.F
30!|| s10jacob ../starter/source/elements/solid/solide10/s10deri3.F
31!|| s10pijsav ../starter/source/elements/solid/solide10/s10deri3.F
32!||--- uses -----------------------------------------------------
33!|| message_mod ../starter/share/message_module/message_mod.F
34!||====================================================================
35 SUBROUTINE s10deri3(VOL,NGL,
36 . XX, YY, ZZ, PX,PY,PZ, NX,
37 . RX, RY, RZ, SX, SY, SZ , TX, TY, TZ,VOLU,VOLN,
38 . ELBUF_STR,VOLG)
39 USE message_mod
40 USE elbufdef_mod
41C-----------------------------------------------
42C I m p l i c i t T y p e s
43C-----------------------------------------------
44#include "implicit_f.inc"
45C-----------------------------------------------
46C G l o b a l P a r a m e t e r s
47C-----------------------------------------------
48#include "mvsiz_p.inc"
49C-----------------------------------------------
50C C o m m o n B l o c k s
51C-----------------------------------------------
52#include "vect01_c.inc"
53C-----------------------------------------------
54C D u m m y A r g u m e n t s
55C-----------------------------------------------
56C REAL
57 double precision
58 . xx(mvsiz,10), yy(mvsiz,10), zz(mvsiz,10)
59 my_real
60 . vol(mvsiz,5),
61 . rx(*), ry(*), rz(*), sx(*), sy(*), sz(*), tx(*),ty(*),tz(*),
62 . nx(mvsiz,10,5),volu(*),voln(*),volg(mvsiz),
63 . px(mvsiz,10,5),py(mvsiz,10,5),pz(mvsiz,10,5)
64 TYPE(elbuf_struct_), TARGET :: ELBUF_STR
65C-----------------------------------------------
66C L o c a l V a r i a b l e s
67C-----------------------------------------------
68 INTEGER NGL(*), I,IP,N,K1,K2,K3,K4,K5,K6,K7,K8,K9,K10,
69 . IPERM(10,4),
70 . NBAD(MVSIZ)
71 double precision
72 . xa(mvsiz,10),ya(mvsiz,10),za(mvsiz,10),
73 . xb(mvsiz,10),yb(mvsiz,10),zb(mvsiz,10),
74 . a4,b4,a4m1,b4m1
75 my_real
76 . alph,beta,w
77 TYPE(l_bufel_) ,POINTER :: LBUF
78 DATA iperm/
79 . 2, 4, 3, 1, 9,10, 6, 5, 8, 7,
80 . 4, 1, 3, 2, 8, 7,10, 9, 5, 6,
81 . 1, 4, 2, 3, 8, 9, 5, 7,10, 6,
82 . 1, 2, 3, 4, 5, 6, 7, 8, 9,10/
83C-----------------------------------------------
84 IF(npt==4)THEN
85 alph = zep5854102
86 beta = zep1381966
87 w = fourth
88 ELSEIF(npt==5)THEN
89 alph = half
90 beta = one_over_6
91 w = nine_over_20
92 ELSE
93 alph = -huge(alph)
94 beta = -huge(beta)
95 w = -huge(w)
96 ENDIF
97 a4 = four * alph
98 b4 = four * beta
99 a4m1 = a4- one
100 b4m1 = b4- one
101C
102 DO i=lft,llt
103 rx(i) = xx(i,1) - xx(i,4)
104 ry(i) = yy(i,1) - yy(i,4)
105 rz(i) = zz(i,1) - zz(i,4)
106 sx(i) = xx(i,2) - xx(i,4)
107 sy(i) = yy(i,2) - yy(i,4)
108 sz(i) = zz(i,2) - zz(i,4)
109 tx(i) = xx(i,3) - xx(i,4)
110 ty(i) = yy(i,3) - yy(i,4)
111 tz(i) = zz(i,3) - zz(i,4)
112 volg(i) =zero
113 ENDDO
114C
115 DO n=1,4
116 DO i=lft,llt
117 xa(i,n) = a4m1*xx(i,n)
118 ya(i,n) = a4m1*yy(i,n)
119 za(i,n) = a4m1*zz(i,n)
120C
121 xb(i,n) = b4m1*xx(i,n)
122 yb(i,n) = b4m1*yy(i,n)
123 zb(i,n) = b4m1*zz(i,n)
124 ENDDO
125 ENDDO
126C
127 DO n=5,10
128 DO i=lft,llt
129 xa(i,n) = a4*xx(i,n)
130 ya(i,n) = a4*yy(i,n)
131 za(i,n) = a4*zz(i,n)
132C
133 xb(i,n) = b4*xx(i,n)
134 yb(i,n) = b4*yy(i,n)
135 zb(i,n) = b4*zz(i,n)
136 ENDDO
137 ENDDO
138C
139 DO ip=1,4
140 k1 = iperm(1,ip)
141 k2 = iperm(2,ip)
142 k3 = iperm(3,ip)
143 k4 = iperm(4,ip)
144 k5 = iperm(5,ip)
145 k6 = iperm(6,ip)
146 k7 = iperm(7,ip)
147 k8 = iperm(8,ip)
148 k9 = iperm(9,ip)
149 k10= iperm(10,ip)
150 lbuf => elbuf_str%BUFLY(1)%LBUF(ip,1,1)
151 CALL s10jacob(alph,beta,w,
152 . xb(1,k1),xb(1,k2),xb(1,k3),xa(1,k4),xb(1,k5),
153 . xb(1,k6),xb(1,k7),xb(1,k8),xb(1,k9),xb(1,k10),
154 . xa(1,k8),xa(1,k9),xa(1,k10),
155 . yb(1,k1),yb(1,k2),yb(1,k3),ya(1,k4),yb(1,k5),
156 . yb(1,k6),yb(1,k7),yb(1,k8),yb(1,k9),yb(1,k10),
157 . ya(1,k8),ya(1,k9),ya(1,k10),
158 . zb(1,k1),zb(1,k2),zb(1,k3),za(1,k4),zb(1,k5),
159 . zb(1,k6),zb(1,k7),zb(1,k8),zb(1,k9),zb(1,k10),
160 . za(1,k8),za(1,k9),za(1,k10),
161 . px(1,k1,ip) ,px(1,k2,ip),px(1,k3,ip),px(1,k4,ip),px(1,k5,ip),
162 . px(1,k6,ip) ,px(1,k7,ip),px(1,k8,ip),px(1,k9,ip),px(1,k10,ip),
163 . py(1,k1,ip) ,py(1,k2,ip),py(1,k3,ip),py(1,k4,ip),py(1,k5,ip),
164 . py(1,k6,ip) ,py(1,k7,ip),py(1,k8,ip),py(1,k9,ip),py(1,k10,ip),
165 . pz(1,k1,ip) ,pz(1,k2,ip),pz(1,k3,ip),pz(1,k4,ip),pz(1,k5,ip),
166 . pz(1,k6,ip) ,pz(1,k7,ip),pz(1,k8,ip),pz(1,k9,ip),pz(1,k10,ip),
167 . nx(1,k1,ip) ,nx(1,k2,ip),nx(1,k3,ip),nx(1,k4,ip),nx(1,k5,ip),
168 . nx(1,k6,ip) ,nx(1,k7,ip),nx(1,k8,ip),nx(1,k9,ip),nx(1,k10,ip),
169 . vol(1,ip) ,lbuf%VOL0DP)
170c
171 ENDDO
172C
173C
174 IF(npt==5)THEN
175 alph = fourth
176 beta = fourth
177 a4 = one
178 b4 = one
179 a4m1 = zero
180 b4m1 = zero
181 w = - four_over_5
182 ip = 5
183C
184 DO n=1,4
185 DO i=lft,llt
186 xa(i,n) = zero
187 ya(i,n) = zero
188 za(i,n) = zero
189 ENDDO
190 ENDDO
191CC
192 lbuf => elbuf_str%BUFLY(1)%LBUF(npt,1,1)
193 CALL s10jacob(alph,beta,w,
194 . xa(1,k1),xa(1,k2),xa(1,k3),xa(1,k4),xx(1,k5),
195 . xx(1,k6),xx(1,k7),xx(1,k8),xx(1,k9),xx(1,k10),
196 . xx(1,k8),xx(1,k9),xx(1,k10),
197 . ya(1,k1),ya(1,k2),ya(1,k3),ya(1,k4),yy(1,k5),
198 . yy(1,k6),yy(1,k7),yy(1,k8),yy(1,k9),yy(1,k10),
199 . yy(1,k8),yy(1,k9),yy(1,k10),
200 . za(1,k1),za(1,k2),za(1,k3),za(1,k4),zz(1,k5),
201 . zz(1,k6),zz(1,k7),zz(1,k8),zz(1,k9),zz(1,k10),
202 . zz(1,k8),zz(1,k9),zz(1,k10),
203 . px(1,k1,ip) ,px(1,k2,ip),px(1,k3,ip),px(1,k4,ip),px(1,k5,ip),
204 . px(1,k6,ip) ,px(1,k7,ip),px(1,k8,ip),px(1,k9,ip),px(1,k10,ip),
205 . py(1,k1,ip) ,py(1,k2,ip),py(1,k3,ip),py(1,k4,ip),py(1,k5,ip),
206 . py(1,k6,ip) ,py(1,k7,ip),py(1,k8,ip),py(1,k9,ip),py(1,k10,ip),
207 . pz(1,k1,ip) ,pz(1,k2,ip),pz(1,k3,ip),pz(1,k4,ip),pz(1,k5,ip),
208 . pz(1,k6,ip) ,pz(1,k7,ip),pz(1,k8,ip),pz(1,k9,ip),pz(1,k10,ip),
209 . nx(1,k1,ip) ,nx(1,k2,ip),nx(1,k3,ip),nx(1,k4,ip),nx(1,k5,ip),
210 . nx(1,k6,ip) ,nx(1,k7,ip),nx(1,k8,ip),nx(1,k9,ip),nx(1,k10,ip),
211 . vol(1,ip) ,lbuf%VOL0DP)
212 ENDIF
213C
214 IF (ismstr==10.OR.ismstr==12) THEN
215 DO ip=1,npt
216 lbuf => elbuf_str%BUFLY(1)%LBUF(ip,1,1)
217 CALL s10pijsav(px(1,1,ip),py(1,1,ip),pz(1,1,ip),lbuf%PIJ,llt)
218 ENDDO
219 ENDIF
220C
221 DO i=lft,llt
222 nbad(i)=0
223 ENDDO
224
225 DO ip=1,npt
226 DO i=lft,llt
227
228 IF(ip<=4.AND.vol(i,ip)<=zero)THEN
229 nbad(i)=nbad(i)+1
230 vol(i,ip)=zero
231 ENDIF
232
233 volg(i) =volg(i) + vol(i,ip)
234 ENDDO
235 ENDDO
236
237 DO i=lft,llt
238 IF(nbad(i)>0)THEN
239 CALL ancmsg(msgid=489,
240 . msgtype=msgerror,
241 . anmode=aninfo_blind_1,
242 . i1=ngl(i))
243 ENDIF
244 ENDDO
245
246 DO i=lft,llt
247 volu(i) =volg(i)/npt
248 voln(i) =volg(i)
249 END DO
250C
251 RETURN
252C
253C 1000 FORMAT(/' ZERO OR NEGATIVE VOLUME : 3D-ELEMENT NB',I10/)
254C 2000 FORMAT(/' ZERO OR NEGATIVE VOLUME : DELETE 3D-ELEMENT NB',I10/)
255 END
256!||====================================================================
257!|| s10jacob ../starter/source/elements/solid/solide10/s10deri3.F
258!||--- called by ------------------------------------------------------
259!|| s10deri3 ../starter/source/elements/solid/solide10/s10deri3.F
260!|| s10jaci3 ../starter/source/elements/solid/solide10/s10jaci3.F
261!||====================================================================
262 SUBROUTINE s10jacob(ALPH,BETA,W,
263 . X1B,X2B,X3B,X4A,X5B,X6B,X7B,X8B,X9B,X10B,X8A,X9A,X10A,
264 . Y1B,Y2B,Y3B,Y4A,Y5B,Y6B,Y7B,Y8B,Y9B,Y10B,Y8A,Y9A,Y10A,
265 . Z1B,Z2B,Z3B,Z4A,Z5B,Z6B,Z7B,Z8B,Z9B,Z10B,Z8A,Z9A,Z10A,
266 . PX1,PX2,PX3,PX4,PX5,PX6,PX7,PX8,PX9,PX10,
267 . PY1,PY2,PY3,PY4,PY5,PY6,PY7,PY8,PY9,PY10,
268 . PZ1,PZ2,PZ3,PZ4,PZ5,PZ6,PZ7,PZ8,PZ9,PZ10,
269 . NX1,NX2,NX3,NX4,NX5,NX6,NX7,NX8,NX9,NX10,
270 . VOL,VOLDP)
271C-----------------------------------------------
272C I m p l i c i t T y p e s
273C-----------------------------------------------
274#include "implicit_f.inc"
275C-----------------------------------------------
276C G l o b a l P a r a m e t e r s
277C-----------------------------------------------
278#include "mvsiz_p.inc"
279C-----------------------------------------------
280C D u m m y A r g u m e n t s
281C-----------------------------------------------
282C REAL
283 double precision
284 . x1b(mvsiz),x2b(mvsiz),x3b(mvsiz),x4a(mvsiz),x5b(mvsiz),
285 . x6b(mvsiz),x7b(mvsiz),x8b(mvsiz),x9b(mvsiz),x10b(mvsiz),
286 . x8a(mvsiz),x9a(mvsiz),x10a(mvsiz),
287 . y1b(mvsiz),y2b(mvsiz),y3b(mvsiz),y4a(mvsiz),y5b(mvsiz),
288 . y6b(mvsiz),y7b(mvsiz),y8b(mvsiz),y9b(mvsiz),y10b(mvsiz),
289 . y8a(mvsiz),y9a(mvsiz),y10a(mvsiz),
290 . z1b(mvsiz),z2b(mvsiz),z3b(mvsiz),z4a(mvsiz),z5b(mvsiz),
291 . z6b(mvsiz),z7b(mvsiz),z8b(mvsiz),z9b(mvsiz),z10b(mvsiz),
292 . z8a(mvsiz),z9a(mvsiz),z10a(mvsiz),voldp(*)
293 my_real
294 . px1(mvsiz),px2(mvsiz),px3(mvsiz),px4(mvsiz),px5(mvsiz),
295 . px6(mvsiz),px7(mvsiz),px8(mvsiz),px9(mvsiz),px10(mvsiz),
296 . py1(mvsiz),py2(mvsiz),py3(mvsiz),py4(mvsiz),py5(mvsiz),
297 . py6(mvsiz),py7(mvsiz),py8(mvsiz),py9(mvsiz),py10(mvsiz),
298 . pz1(mvsiz),pz2(mvsiz),pz3(mvsiz),pz4(mvsiz),pz5(mvsiz),
299 . pz6(mvsiz),pz7(mvsiz),pz8(mvsiz),pz9(mvsiz),pz10(mvsiz),
300 . nx1(mvsiz),nx2(mvsiz),nx3(mvsiz),nx4(mvsiz),nx5(mvsiz),
301 . nx6(mvsiz),nx7(mvsiz),nx8(mvsiz),nx9(mvsiz),nx10(mvsiz),
302 . vol(mvsiz),alph,beta,w
303C-----------------------------------------------
304C C o m m o n B l o c k s
305C-----------------------------------------------
306#include "vect01_c.inc"
307C-----------------------------------------------
308C L o c a l V a r i a b l e s
309C-----------------------------------------------
310 INTEGER I
311C REAL
312 DOUBLE PRECISION
313 . A4M1,B4M1,AA
314c my_real
315 DOUBLE PRECISION
316 . D,A4MB4,A4,B4,BB,AB,B2,
317 . DXDR,DXDS,DXDT,DYDR,DYDS,DYDT,DZDR,DZDS,DZDT
318c my_real
319 DOUBLE PRECISION
320 . DRDX(MVSIZ), DSDX(MVSIZ), DTDX(MVSIZ),
321 . drdy(mvsiz), dsdy(mvsiz), dtdy(mvsiz),
322 . drdz(mvsiz), dsdz(mvsiz), dtdz(mvsiz),
323 . det(mvsiz)
324C-----------------------------------------------
325 a4 = four * alph
326 b4 = four * beta
327 a4m1 = a4 - one
328 b4m1 = b4 - one
329 aa = (two*alph - one)*alph
330 bb = (two*beta - one)*beta
331 b2 = four*beta*beta
332 ab = four*alph*beta
333C-----------------------------------------------
334 DO i=lft,llt
335 nx1(i) = bb
336 nx2(i) = bb
337 nx3(i) = bb
338 nx4(i) = aa
339 nx5(i) = b2
340 nx6(i) = b2
341 nx7(i) = b2
342 nx8(i) = ab
343 nx9(i) = ab
344 nx10(i)= ab
345 ENDDO
346C
347 DO i=lft,llt
348 aa = x5b(i) + x6b(i) + x7b(i)
349 . - x4a(i) - x8b(i) - x9b(i) - x10b(i)
350 dxdr = x1b(i) +x8a(i) - x6b(i) + aa
351 dxds = x2b(i) +x9a(i) - x7b(i) + aa
352 dxdt = x3b(i) +x10a(i)- x5b(i) + aa
353C
354 aa = y5b(i) + y6b(i) + y7b(i)
355 . - y4a(i) - y8b(i) - y9b(i) - y10b(i)
356 dydr = y1b(i) +y8a(i) - y6b(i) + aa
357 dyds = y2b(i) +y9a(i) - y7b(i) + aa
358 dydt = y3b(i) +y10a(i)- y5b(i) + aa
359C
360 aa = z5b(i) + z6b(i) + z7b(i)
361 . - z4a(i) - z8b(i) - z9b(i) - z10b(i)
362 dzdr = z1b(i) +z8a(i) - z6b(i) + aa
363 dzds = z2b(i) +z9a(i) - z7b(i) + aa
364 dzdt = z3b(i) +z10a(i)- z5b(i) + aa
365C
366 drdx(i)=dyds*dzdt-dzds*dydt
367 dsdx(i)=dydt*dzdr-dzdt*dydr
368 dtdx(i)=dydr*dzds-dzdr*dyds
369C
370 drdy(i)=dzds*dxdt-dxds*dzdt
371 dsdy(i)=dzdt*dxdr-dxdt*dzdr
372 dtdy(i)=dzdr*dxds-dxdr*dzds
373C
374 drdz(i)=dxds*dydt-dyds*dxdt
375 dsdz(i)=dxdt*dydr-dydt*dxdr
376 dtdz(i)=dxdr*dyds-dydr*dxds
377C
378 det(i) = dxdr * drdx(i)
379 . + dydr * drdy(i)
380 . + dzdr * drdz(i)
381 voldp(i) = w * det(i) * one_over_6
382 vol(i) = voldp(i)
383C
384 ENDDO
385C IF(IDTMIN(1)==1)THEN
386C DO I=LFT,LLT
387C IF(OFF(I) ==0.)THEN
388C DET(I)=1.
389C ELSEIF(DET(I)<=VOLMIN)THEN
390C DET(I)=1.
391C OFF(I)=0.
392C WRITE(ISTDO,2000) NGL(I)
393C WRITE(IOUT ,2000) NGL(I)
394C ELSEIF(DET(I)<=0.0)THEN
395C WRITE(ISTDO,1000) NGL(I)
396C WRITE(IOUT ,1000) NGL(I)
397C ENDIF
398C ENDDO
399C ELSEIF(IDTMIN(1)==2)THEN
400C DO I=LFT,LLT
401C IF(OFF(I) ==0.)THEN
402C DET(I)=1.
403C ELSEIF(DET(I)<=VOLMIN.OR.DET(I)<=0.0)THEN
404C DET(I)=1.
405C OFF(I)=0.
406C WRITE(ISTDO,2000) NGL(I)
407C WRITE(IOUT ,2000) NGL(I)
408C ENDIF
409C ENDDO
410C ELSE
411C DO I=LFT,LLT
412C IF(OFF(I) ==0.)DET(I)=1.
413C ENDDO
414C ENDIF
415CC
416 DO i=lft,llt
417 a4 = four * alph
418 b4 = four * beta
419 a4m1 = a4 - one
420 b4m1 = b4 - one
421 d = one/max(em30,det(i))
422C
423 b4m1 = d * b4m1
424 a4m1 = d * a4m1
425 b4 = d * b4
426 a4 = d * a4
427 a4mb4 = a4 - b4
428C
429 px1(i) = b4m1 * drdx(i)
430 py1(i) = b4m1 * drdy(i)
431 pz1(i) = b4m1 * drdz(i)
432C
433 px2(i) = b4m1 * dsdx(i)
434 py2(i) = b4m1 * dsdy(i)
435 pz2(i) = b4m1 * dsdz(i)
436C
437 px3(i) = b4m1 * dtdx(i)
438 py3(i) = b4m1 * dtdy(i)
439 pz3(i) = b4m1 * dtdz(i)
440C
441 px4(i) =-a4m1 * (drdx(i)+dsdx(i)+dtdx(i))
442 py4(i) =-a4m1 * (drdy(i)+dsdy(i)+dtdy(i))
443 pz4(i) =-a4m1 * (drdz(i)+dsdz(i)+dtdz(i))
444C
445 px5(i) = b4 * (drdx(i)+dsdx(i))
446 py5(i) = b4 * (drdy(i)+dsdy(i))
447 pz5(i) = b4 * (drdz(i)+dsdz(i))
448C
449 px6(i) = b4 * (dsdx(i)+dtdx(i))
450 py6(i) = b4 * (dsdy(i)+dtdy(i))
451 pz6(i) = b4 * (dsdz(i)+dtdz(i))
452C
453 px7(i) = b4 * (dtdx(i)+drdx(i))
454 py7(i) = b4 * (dtdy(i)+drdy(i))
455 pz7(i) = b4 * (dtdz(i)+drdz(i))
456C
457 px8(i) = a4mb4 * drdx(i) - px6(i)
458 py8(i) = a4mb4 * drdy(i) - py6(i)
459 pz8(i) = a4mb4 * drdz(i) - pz6(i)
460C
461 px9(i) = a4mb4 * dsdx(i) - px7(i)
462 py9(i) = a4mb4 * dsdy(i) - py7(i)
463 pz9(i) = a4mb4 * dsdz(i) - pz7(i)
464C
465 px10(i)= a4mb4 * dtdx(i) - px5(i)
466 py10(i)= a4mb4 * dtdy(i) - py5(i)
467 pz10(i)= a4mb4 * dtdz(i) - pz5(i)
468C
469 ENDDO
470C
471 RETURN
472 END
473!||====================================================================
474!|| s10pijsav ../starter/source/elements/solid/solide10/s10deri3.F
475!||--- called by ------------------------------------------------------
476!|| s10deri3 ../starter/source/elements/solid/solide10/s10deri3.F
477!|| s10jaci3 ../starter/source/elements/solid/solide10/s10jaci3.F
478!||====================================================================
479 SUBROUTINE s10pijsav(PX,PY,PZ, PIJ,NEL)
480C-----------------------------------------------
481C I m p l i c i t T y p e s
482C-----------------------------------------------
483#include "implicit_f.inc"
484C-----------------------------------------------
485C G l o b a l P a r a m e t e r s
486C-----------------------------------------------
487#include "mvsiz_p.inc"
488C-----------------------------------------------
489C D u m m y A r g u m e n t s
490C-----------------------------------------------
491 INTEGER NEL
492C REAL
493 my_real
494 . PX(MVSIZ,10),PY(MVSIZ,10),PZ(MVSIZ,10),
495 . PIJ(NEL,3,10)
496C-----------------------------------------------
497C L o c a l V a r i a b l e s
498C-----------------------------------------------
499 INTEGER I, J
500
501C-----------------------------------------------
502 DO j=1,10
503 DO i=1,nel
504 pij(i,1,j) = px(i,j)
505 pij(i,2,j) = py(i,j)
506 pij(i,3,j) = pz(i,j)
507 END DO
508 END DO
509C-----------
510 RETURN
511 END
my_real
#define my_real
Definition cppsort.cpp:32
max
#define max(a, b)
Definition macros.h:21
message_mod
Definition message_mod.F:1257
s10jacob
subroutine s10jacob(alph, beta, w, x1b, x2b, x3b, x4a, x5b, x6b, x7b, x8b, x9b, x10b, x8a, x9a, x10a, y1b, y2b, y3b, y4a, y5b, y6b, y7b, y8b, y9b, y10b, y8a, y9a, y10a, z1b, z2b, z3b, z4a, z5b, z6b, z7b, z8b, z9b, z10b, z8a, z9a, z10a, px1, px2, px3, px4, px5, px6, px7, px8, px9, px10, py1, py2, py3, py4, py5, py6, py7, py8, py9, py10, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pz9, pz10, nx1, nx2, nx3, nx4, nx5, nx6, nx7, nx8, nx9, nx10, vol, voldp)
Definition s10deri3.F:271
s10pijsav
subroutine s10pijsav(px, py, pz, pij, nel)
Definition s10deri3.F:480
s10deri3
subroutine s10deri3(vol, ngl, xx, yy, zz, px, py, pz, nx, rx, ry, rz, sx, sy, sz, tx, ty, tz, volu, voln, elbuf_str, volg)
Definition s10deri3.F:39
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:895