OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i7buc1.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "units_c.inc"
#include "com04_c.inc"
#include "vect07_c.inc"
#include "scr06_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine i7buc1 (x, irect, nsv, bumult, nseg, nmn, nrtm, mwa, nsn, cand_e, cand_n, gap, xyzm, noint, i_stok, dist, tzinf, maxbox, minbox, msr, stf, stfn, multimp, istf, iddlevel, itab, gap_s, gap_m, igap, gapmin, gapmax, inacti, gap_s_l, gap_m_l, i_mem, id, titr, it19, prov_n, prov_e, nsvg, ix1, ix2, ix3, ix4, n11, n12, n13, pene, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, x0, y0, z0, nx1, ny1, nz1, nx2, ny2, nz2, nx3, ny3, nz3, nx4, ny4, nz4, p1, p2, p3, p4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, stif)

Function/Subroutine Documentation

◆ i7buc1()

subroutine i7buc1 ( x,
integer, dimension(4,*) irect,
integer, dimension(*) nsv,
bumult,
integer, dimension(*) nseg,
integer nmn,
integer nrtm,
integer, dimension(*) mwa,
integer nsn,
integer, dimension(*) cand_e,
integer, dimension(*) cand_n,
gap,
xyzm,
integer noint,
integer i_stok,
dist,
tzinf,
maxbox,
minbox,
integer, dimension(*) msr,
stf,
stfn,
integer multimp,
integer istf,
integer iddlevel,
integer, dimension(*) itab,
gap_s,
gap_m,
integer igap,
gapmin,
gapmax,
integer inacti,
gap_s_l,
gap_m_l,
integer i_mem,
integer id,
character(len=nchartitle) titr,
integer it19,
integer, dimension(mvsiz), intent(in) prov_n,
integer, dimension(mvsiz), intent(in) prov_e,
integer, dimension(mvsiz), intent(in) nsvg,
integer, dimension(mvsiz), intent(in) ix1,
integer, dimension(mvsiz), intent(in) ix2,
integer, dimension(mvsiz), intent(in) ix3,
integer, dimension(mvsiz), intent(in) ix4,
intent(in) n11,
intent(in) n12,
intent(in) n13,
intent(in) pene,
intent(in) x1,
intent(in) x2,
intent(in) x3,
intent(in) x4,
intent(in) y1,
intent(in) y2,
intent(in) y3,
intent(in) y4,
intent(in) z1,
intent(in) z2,
intent(in) z3,
intent(in) z4,
intent(in) xi,
intent(in) yi,
intent(in) zi,
intent(in) x0,
intent(in) y0,
intent(in) z0,
intent(in) nx1,
intent(in) ny1,
intent(in) nz1,
intent(in) nx2,
intent(in) ny2,
intent(in) nz2,
intent(in) nx3,
intent(in) ny3,
intent(in) nz3,
intent(in) nx4,
intent(in) ny4,
intent(in) nz4,
intent(in) p1,
intent(in) p2,
intent(in) p3,
intent(in) p4,
intent(in) lb1,
intent(in) lb2,
intent(in) lb3,
intent(in) lb4,
intent(in) lc1,
intent(in) lc2,
intent(in) lc3,
intent(in) lc4,
intent(in) stif )

Definition at line 38 of file i7buc1.F.

58C-----------------------------------------------
59C M o d u l e s
60C-----------------------------------------------
61 USE message_mod
63C-----------------------------------------------
64C I m p l i c i t T y p e s
65C-----------------------------------------------
66#include "implicit_f.inc"
67C-----------------------------------------------
68C G l o b a l P a r a m e t e r s
69C-----------------------------------------------
70#include "mvsiz_p.inc"
71C-----------------------------------------------
72C C o m m o n B l o c k s
73C-----------------------------------------------
74#include "units_c.inc"
75#include "com04_c.inc"
76#include "vect07_c.inc"
77#include "scr06_c.inc"
78C-----------------------------------------------
79C D u m m y A r g u m e n t s
80C-----------------------------------------------
81 INTEGER NMN, NRTM, NSN, NOINT,I_STOK,MULTIMP,ISTF,IGAP,
82 . INACTI
83 INTEGER IRECT(4,*),NSV(*),NSEG(*),MWA(*)
84 INTEGER CAND_E(*),CAND_N(*),MSR(*),MAXSIZ,IDDLEVEL
85 INTEGER ITAB(*),IT19
87 . stf(*),stfn(*),x(3,*),xyzm(6,*),gap_s(*),gap_m(*),
88 . dist,bumult,gap,tzinf,maxbox,minbox,gapmin,gapmax,
89 . gap_s_l(*),gap_m_l(*)
90 INTEGER ID
91 CHARACTER(LEN=NCHARTITLE) :: TITR
92 INTEGER, DIMENSION(MVSIZ), INTENT(IN) :: PROV_N,PROV_E,NSVG
93 INTEGER, DIMENSION(MVSIZ), INTENT(IN) :: IX1,IX2,IX3,IX4
94 my_real, DIMENSION(MVSIZ), INTENT(IN) :: n11,n12,n13,pene
95 my_real, DIMENSION(MVSIZ), INTENT(IN) :: x1,x2,x3,x4
96 my_real, DIMENSION(MVSIZ), INTENT(IN) :: y1,y2,y3,y4
97 my_real, DIMENSION(MVSIZ), INTENT(IN) :: z1,z2,z3,z4
98 my_real, DIMENSION(MVSIZ), INTENT(IN) :: xi,yi,zi
99 my_real, DIMENSION(MVSIZ), INTENT(IN) :: x0,y0,z0
100 my_real, DIMENSION(MVSIZ), INTENT(IN) :: nx1,ny1,nz1
101 my_real, DIMENSION(MVSIZ), INTENT(IN) :: nx2,ny2,nz2
102 my_real, DIMENSION(MVSIZ), INTENT(IN) :: nx3,ny3,nz3
103 my_real, DIMENSION(MVSIZ), INTENT(IN) :: nx4,ny4,nz4
104 my_real, DIMENSION(MVSIZ), INTENT(IN) :: p1,p2,p3,p4
105 my_real, DIMENSION(MVSIZ), INTENT(IN) :: lb1,lb2,lb3,lb4
106 my_real, DIMENSION(MVSIZ), INTENT(IN) :: lc1,lc2,lc3,lc4,stif
107C-----------------------------------------------
108C L o c a l V a r i a b l e s
109C-----------------------------------------------
110 INTEGER I, J, L, N1, N2, N3, N4, I_AMAX,I_MEM
111 INTEGER I_ADD, ADESTK, NB_NC, NB_EC, ADNSTK, IBID
112 INTEGER IP1, IP2, IP21, IP22, IP31,J_STOK,I_BID,NB_N_B
113 my_real
114 . dx1,dy1,dz1,
115 . dx3,dy3,dz3,
116 . dx4,dy4,dz4,
117 . dx6,dy6,dz6,
118 . dd1,dd2,dd3,dd4,dd,dd0,xmin,ymin,zmin,tb,
119 . xmax,ymax,zmax,minbox_st,maxbox_st,gapsmax,
120 . bid,tzinf_st,marge,marge_st,gapv(mvsiz),
121 . xmax_m,ymax_m,zmax_m,xmin_m,ymin_m,zmin_m,
122 . xmax_s,ymax_s,zmax_s,xmin_s,ymin_s,zmin_s
123 LOGICAL :: TYPE18
124C-----------------------------------------------
125C S o u r c e L i n e s
126C-----------------------------------------------
127
128 type18=.false.
129 IF(inacti==7)type18=.true.
130
131
132C 1-CALCUL TAILLE DES ZONES INFLUENCES
133C DD EST LA LONGEUR MOYENNE ELEMENT
134 dd=zero
135 DO 10 l=1,nrtm
136C CONNECIVITES ELEMENT
137 n1=irect(1,l)
138 n2=irect(2,l)
139 n3=irect(3,l)
140 n4=irect(4,l)
141C LONGUEUR COTE 1
142 dx1=(x(1,n1)-x(1,n2))
143 dy1=(x(2,n1)-x(2,n2))
144 dz1=(x(3,n1)-x(3,n2))
145 dd1=sqrt(dx1**2+dy1**2+dz1**2)
146C LONGUEUR COTE 2
147 dx3=(x(1,n1)-x(1,n4))
148 dy3=(x(2,n1)-x(2,n4))
149 dz3=(x(3,n1)-x(3,n4))
150 dd2=sqrt(dx3**2+dy3**2+dz3**2)
151C LONGUEUR COTE 3
152 dx4=(x(1,n3)-x(1,n2))
153 dy4=(x(2,n3)-x(2,n2))
154 dz4=(x(3,n3)-x(3,n2))
155 dd3=sqrt(dx4**2+dy4**2+dz4**2)
156C LONGUEUR COTE 4
157 dx6=(x(1,n4)-x(1,n3))
158 dy6=(x(2,n4)-x(2,n3))
159 dz6=(x(3,n4)-x(3,n3))
160 dd4=sqrt(dx6**2+dy6**2+dz6**2)
161 dd=dd+ (dd1+dd2+dd3+dd4)
162 10 CONTINUE
163C TAILLE ZINF = .1*TAILLE MOYENNE ELEMENT DE CHAQUE COTE
164C TAILLE BUCKET MIN = TZINF * BUMULT
165 dd0=dd/nrtm/four
166 dd = dd0
167C DD = MAX(DD0,ONEP251*GAP)
168C remove this additional test which breaks igap performance
169C in case this factor is needed, then needs to compute mean gap instead
170C TZINF = BUMULT*DD
171 marge = bumult*dd
172C calcul de TZINF en fct de la marge et non le contraire
173 tzinf = marge + gap
174C MARGE_ST : marge independante du BUMULT en input pour trouver les memes pene initiales (cas delete ou chgt coord.)
175 marge_st = bmul0*dd
176C Si IDDLEVEL = 0, alors 1er passage ici, avec la MARGE de l'ENGINE pour trouver les memes candidats
177 IF(iddlevel==0) marge_st = marge
178 tzinf_st = marge_st + gap
179C MAXBOX= ZEP9*(DD - TZINF)
180C DD + 2*TZINF : taille element augmentee de zone influence
181 maxbox= half*(dd + 2*tzinf)
182 minbox= half*maxbox
183 maxbox_st= half*(dd + 2*tzinf_st)
184 minbox_st= half*maxbox_st
185C MIS A ZERO POUR FAIRE SEARCH COMPLET CYCLE 0 ENGINE
186 dist = zero
187C--------------------------------
188C CALCUL DES BORNES DU DOMAINE
189C--------------------------------
190
191 xmin_m=ep30
192 xmax_m=-ep30
193 ymin_m=ep30
194 ymax_m=-ep30
195 zmin_m=ep30
196 zmax_m=-ep30
197
198 DO 20 i=1,nmn
199 j=msr(i)
200 xmin_m= min(xmin_m,x(1,j))
201 ymin_m= min(ymin_m,x(2,j))
202 zmin_m= min(zmin_m,x(3,j))
203 xmax_m= max(xmax_m,x(1,j))
204 ymax_m= max(ymax_m,x(2,j))
205 zmax_m= max(zmax_m,x(3,j))
206 20 CONTINUE
207
208 xmin_s=ep30
209 xmax_s=-ep30
210 ymin_s=ep30
211 ymax_s=-ep30
212 zmin_s=ep30
213 zmax_s=-ep30
214
215 DO 25 i=1,nsn
216 j=nsv(i)
217 xmin_s= min(xmin_s,x(1,j))
218 ymin_s= min(ymin_s,x(2,j))
219 zmin_s= min(zmin_s,x(3,j))
220 xmax_s= max(xmax_s,x(1,j))
221 ymax_s= max(ymax_s,x(2,j))
222 zmax_s= max(zmax_s,x(3,j))
223 25 CONTINUE
224C
225 xmin=min(xmin_m-tzinf_st,xmin_s)
226 ymin=min(ymin_m-tzinf_st,ymin_s)
227 zmin=min(zmin_m-tzinf_st,zmin_s)
228 xmax=max(xmax_m+tzinf_st,xmax_s)
229 ymax=max(ymax_m+tzinf_st,ymax_s)
230 zmax=max(zmax_m+tzinf_st,zmax_s)
231C
232C-----2- TRI PAR BOITES DES ELEMENTS ET DES NOEUDS
233C
234 nb_n_b = 1
235 i_mem = 0
236C SI ON A PAS ASSEZ DE MEMOIRE POUR LES PILES ON RECOMMENCE LE TRI
237C EN INCREMENTANT LE NB_N_B (NOMBRE DE NOEUDS PAR BOITE FINIE)
238 100 CONTINUE
239C POINTEUR NOM TAILLE
240C P1........Elt Bas Pile NRTM
241C P2........Elt PILE 2*NRTM
242C P21.......BPN NSN
243C P22.......PN NSN
244C P31.......ADDI 2000
245C
246C POUR ONE MAILLAGE DE TOPOLOGIE CARRE LA TAILLE DE P2 PEUT ETRE ESTIMEE A:
247C 4n: NUMELC + 6*SQRT(NUMELC) + 8*LOG2(NUMELC)
248C 3n: NUMELTG + 6*SQRT(2*NUMELTG) + 8*LOG2(NUMELTG)
249C
250C 4n: NUMELC + 6*SQRT(NUMELC) + 12*LOG(NUMELC) +
251C 3n: NUMELTG + 6*SQRT(2*NUMELTG) + 12*LOG(NUMELTG)
252C
253C NUMELC + NUMELTG + 6*SQRT(NUMELC+2*NUMELTG) + 12*LOG(NUMELC+NUMELTG)
254C
255C POUR ONE MAILLAGE DE TOPOLOGIE LINEAIRE LA TAILLE DE P2 PEUT ETRE ESTIMEE A:
256C NUMELC + NUMELTG + (NUMELC+NUMELTG)*(1 + 1/2 + 1/4 +...) + LOG2(NUMELC+NUMELTG)
257C 3*(NUMELC+NUMELTG) + LOG2(NUMELC+NUMELTG)
258C 3*(NUMELC+NUMELTG) + 300
259C
260 maxsiz = max(numnod,nrtm+100)
261 ip1 = 1
262 ip2 = ip1+maxsiz
263C IP21= IP2+2*MAXSIZ
264 ip21= ip2+3*maxsiz
265 ip22= ip21+numnod
266 ip31= ip22+numnod
267C IFIN= IP22+2000
268C
269C-----INITIALISATION DES ADRESSES ET X,Y,Z
270C
271C ADDE ADDN X Y Z
272C 1 1 XMIN YMIN ZMIN
273C 1 1 XMAX YMAX ZMAX
274C
275 mwa(ip31) = 0
276 mwa(ip31+1) = 0
277 mwa(ip31+2) = 0
278 mwa(ip31+3) = 0
279 i_add = 1
280 i_amax = 1
281 xyzm(1,i_add) = xmin
282 xyzm(2,i_add) = ymin
283 xyzm(3,i_add) = zmin
284 xyzm(4,i_add) = xmax
285 xyzm(5,i_add) = ymax
286 xyzm(6,i_add) = zmax
287 i_stok = 0
288 j_stok = 0
289 adnstk = 0
290 adestk = 0
291 nb_nc = nsn
292 nb_ec = nrtm
293C
294C-----COPIE DES NOS DE SEGMENTS ET DE NOEUDS DANS MWA(IP1) ET IP21
295C
296 DO 120 i=1,nb_ec
297 mwa(ip1+i-1) = i
298 120 CONTINUE
299 DO 140 i=1,nb_nc
300 mwa(ip21+i-1) = i
301 140 CONTINUE
302C
303C-----DEBUT DE LA PHASE DE TRI
304C
305C TANT QUE IL RESTE UNE ADRESSE A TRIER
306C------------------
307 200 CONTINUE
308C------------------
309C SEPARER B ET N EN TWO
310 CALL i7tri(
311 1 mwa(ip1) ,mwa(ip2) ,mwa(ip21),mwa(ip22),mwa(ip31+2*(i_add-2)),
312 2 irect ,x ,nb_nc ,nb_ec ,xyzm ,
313 3 i_add ,nsv ,i_amax ,xmax ,ymax ,
314 4 zmax ,3*maxsiz ,i_stok ,i_mem ,nb_n_b ,
315 5 cand_n ,cand_e ,nsn ,noint ,tzinf_st ,
316 6 maxbox_st,minbox_st,stf ,stfn ,j_stok ,
317 7 multimp ,istf , itab ,gap ,gap_s ,
318 8 gap_m ,igap ,gapmin ,gapmax ,marge_st ,
319 9 gap_s_l ,gap_m_l ,id ,titr ,
320 1 ix1 ,ix2 ,ix3,ix4 ,nsvg ,
321 2 prov_n ,prov_e ,n11,n12 ,n13 ,
322 3 pene ,x1 ,x2 ,x3 ,x4 ,
323 4 y1 ,y2 ,y3 ,y4 ,z1 ,
324 5 z2 ,z3 ,z4 ,xi ,yi ,
325 6 zi ,x0 ,y0 ,z0 ,nx1 ,
326 7 ny1 ,nz1 ,nx2,ny2 ,nz2 ,
327 8 nx3 ,ny3 ,nz3,nx4 ,ny4 ,
328 9 nz4 ,p1 ,p2 ,p3 ,p4 ,
329 1 lb1 ,lb2 ,lb3,lb4 ,lc1 ,
330 2 lc2 ,lc3 ,lc4,stif)
331C------------------
332 IF (i_mem == 2)THEN
333 RETURN
334 ENDIF
335C I_MEM = 1 ==> PAS ASSEZ DE MEMOIRE PILE
336C I_MEM = 2 ==> PAS ASSEZ DE MEMOIRE CANDIDATS
337 IF(i_mem==1)THEN
338 nb_n_b = nb_n_b + 1
339 i_mem = 0
340 GO TO 100
341 ELSE IF(i_mem==2) THEN
342 marge_st = three_over_4*marge_st
343 tzinf_st = marge_st + gap
344 i_mem = 0
345 IF(marge_st<em03) THEN
346C WRITE(ISTDO,*)' ***ERROR INFINITE LOOP DETECTED(I7BUC)'
347C WRITE(IOUT,*)' **ERROR INFINITE LOOP DETECTED(I7BUC)'
348C CALL ARRET(2)
349 CALL ancmsg(msgid=83,
350 . msgtype=msgerror,
351 . anmode=aninfo,
352 . i1=id,
353 . c1=titr)
354 ENDIF
355 GO TO 100
356 ENDIF
357 IF(i_add/=0) GO TO 200
358C FIN BOUCLE TANT QUE
359C---------------------------------
360 IF(j_stok/=0)THEN
361 lft = 1
362 llt = j_stok
363 CALL i7cor3(x ,irect,nsv ,prov_e ,prov_n,
364 . stf ,stfn ,gapv ,igap ,gap ,
365 . gap_s,gap_m,istf ,gapmin ,gapmax,
366 . gap_s_l,gap_m_l ,zero ,ix1 ,ix2 ,
367 5 ix3 ,ix4 ,nsvg,x1 ,x2 ,
368 6 x3 ,x4 ,y1 ,y2 ,y3 ,
369 7 y4 ,z1 ,z2 ,z3 ,z4 ,
370 8 xi ,yi ,zi ,stif ,zero ,
371 9 llt)
372 CALL i7dst3(ix3,ix4,x1 ,x2 ,x3 ,
373 1 x4 ,y1 ,y2 ,y3 ,y4 ,
374 2 z1 ,z2 ,z3 ,z4 ,xi ,
375 3 yi ,zi ,x0 ,y0 ,z0 ,
376 4 nx1,ny1,nz1,nx2,ny2,
377 5 nz2,nx3,ny3,nz3,nx4,
378 6 ny4,nz4,p1 ,p2 ,p3 ,
379 7 p4 ,lb1,lb2,lb3,lb4,
380 8 lc1,lc2,lc3,lc4,llt)
381 CALL i7pen3(marge_st,gapv,n11,n12,n13 ,
382 1 pene ,nx1 ,ny1,nz1,nx2,
383 2 ny2 ,nz2 ,nx3,ny3,nz3,
384 3 nx4 ,ny4 ,nz4,p1 ,p2 ,
385 4 p3 ,p4,llt)
386 IF(i_stok+j_stok<multimp*nsn) THEN
387 CALL i7cmp3(i_stok,cand_e ,cand_n,1,pene,
388 1 prov_n,prov_e)
389 ELSE
390 i_bid = 0
391 CALL i7cmp3(i_bid,cand_e,cand_n,0,pene,
392 1 prov_n,prov_e)
393 IF(i_stok+i_bid<multimp*nsn) THEN
394 CALL i7cmp3(i_stok,cand_e,cand_n,1,pene,
395 1 prov_n,prov_e)
396 ELSE
397 i_mem = 2
398 RETURN
399 marge_st = three_over_4*marge_st
400 tzinf_st = marge_st + gap
401 i_mem = 0
402C ne pas dimunuer la taille des boite
403C MINBOX= THREE_OVER_4*MINBOX
404C MAXBOX= THREE_OVER_4*MAXBOX
405 IF(marge_st<em03) THEN
406C WRITE(ISTDO,*)' ***ERROR INFINITE LOOP DETECTED(I7BUC)'
407C WRITE(IOUT,*)' **ERROR INFINITE LOOP DETECTED(I7BUC)'
408C CALL ARRET(2)
409 CALL ancmsg(msgid=83,
410 . msgtype=msgerror,
411 . anmode=aninfo,
412 . i1=id,
413 . c1=titr)
414 ENDIF
415 GO TO 100
416 ENDIF
417 ENDIF
418 ENDIF
419C
420 IF(.NOT.type18)THEN
421 IF ((nsn/=0).AND.(it19<=0)) THEN
422 WRITE(iout,*)' POSSIBLE IMPACT NUMBER:',i_stok,' (<=', 1+(i_stok-1)/nsn,'*NSN)'
423 ELSEIF(nsn==0) THEN
424 CALL ancmsg(msgid=552,
425 . msgtype=msgwarning,
426 . anmode=aninfo_blind_2,
427 . i1=id,
428 . c1=titr)
429 ENDIF
430 endif!(.NOT.TYPE18)
431
432C
433C MISE A ZERO DE TAG POUR I7PWR3
434C
435 DO i=1,numnod+numfakenodigeo
436 mwa(i)=0
437 ENDDO
438C
439 RETURN
#define my_real
Definition cppsort.cpp:32
subroutine i7cmp3(i_stok, cand_e, cand_n, iflag, pene, prov_n, prov_e)
Definition i7cmp3.F:82
subroutine ymax(idn, fac, npc, pld, stiffmin, stiffmax, stiffini, stiffavg)
Definition law100_upd.F:272
#define min(a, b)
Definition macros.h:20
#define max(a, b)
Definition macros.h:21
initmumps id
integer, parameter nchartitle
subroutine i7cor3(x, irect, nsv, cand_e, cand_n, stf, stfn, gapv, igap, gap, gap_s, gap_m, istf, gapmin, gapmax, gap_s_l, gap_m_l, drad, ix1, ix2, ix3, ix4, nsvg, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, stif, dgapload, last)
Definition i7cor3.F:43
subroutine i7dst3(ix3, ix4, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, x0, y0, z0, nx1, ny1, nz1, nx2, ny2, nz2, nx3, ny3, nz3, nx4, ny4, nz4, p1, p2, p3, p4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, last)
Definition i7dst3.F:46
subroutine i7pen3(marge, gapv, n1, n2, n3, pene, nx1, ny1, nz1, nx2, ny2, nz2, nx3, ny3, nz3, nx4, ny4, nz4, p1, p2, p3, p4, last)
Definition i7pen3.F:43
subroutine i7tri(bpe, pe, bpn, pn, add, irect, x, nb_nc, nb_ec, xyzm, i_add, nsv, i_amax, xmax, ymax, zmax, maxsiz, i_stok, i_mem, nb_n_b, cand_n, cand_e, nsn, noint, tzinf, maxbox, minbox, stf, stfn, j_stok, multimp, istf, itab, gap, gap_s, gap_m, igap, gapmin, gapmax, marge, gap_s_l, gap_m_l, id, titr, ix1, ix2, ix3, ix4, nsvg, prov_n, prov_e, n11, n12, n13, pene, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, x0, y0, z0, nx1, ny1, nz1, nx2, ny2, nz2, nx3, ny3, nz3, nx4, ny4, nz4, p1, p2, p3, p4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, stif)
Definition i7tri.F:58
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:889