33 1 TAGBUF,NPART,FLAG,IDOM)
38#include "implicit_f.inc"
42 INTEGER NIX,IX(NIX,*),NIX1,NIX2,NUMEL,IPARTE(*),
43 1 tagbuf(*),npart,flag,idom,tag_el
47 INTEGER J,L,NUM_KJ,RES
55 IF (tagbuf(iparte(j))==0)
THEN
57 IF (tagbuf(ix(l,j)+npart)==1)
THEN
58 tagbuf(ix(l,j)+npart)=0
69 IF (tagbuf(iparte(j))==1)
THEN
71 IF (tagbuf(ix(l,j)+npart)<2)
THEN
72 tagbuf(ix(l,j)+npart)=1
83 IF (tagbuf(iparte(j))==0)
THEN
85 IF (tagbuf(ix(l,j)+npart) == 1)
THEN
86 tagbuf(ix(l,j)+npart)=1+idom
87 ELSEIF (tagbuf(ix(l,j)+npart)>1)
THEN
88 IF (tagbuf(ix(l,j)+npart)/=(1+idom))
THEN
93 . i2=tagbuf(ix(l,j)+npart)-1,
99 tagbuf(ix(l,j)+npart)=-idom
105 ELSEIF (flag==2)
THEN
111 tag_el=iparte(j+npart)
114 IF (tagbuf(ix(l,j)+npart)<=2)
THEN
115 tagbuf(ix(l,j)+npart)=2*iparte(j+npart)
118 ELSEIF (tag_el==-1)
THEN
120 IF (tagbuf(ix(l,j)+npart)<0)
THEN
121 tagbuf(ix(l,j)+npart)=0
127 ELSEIF (flag==3)
THEN
132 IF (tagbuf(iparte(j))==1)
THEN
135 IF (tagbuf(ix(l,j)+npart)<1)
THEN
136 tagbuf(ix(l,j)+npart)=1
144 IF (tagbuf(iparte(j))==0)
THEN
147 IF (tagbuf(ix(l,j)+npart)==0)
THEN
148 tagbuf(ix(l,j)+npart)=-1
155 ELSEIF ((flag==4).AND.(numel>0))
THEN
159 num_kj = ix(1,numel+1)
163 IF (tagbuf(iparte(res))==1)
THEN
166 IF (tagbuf(ix(l,j)+npart)<1)
THEN
167 tagbuf(ix(l,j)+npart)=1
176 IF (tagbuf(iparte(res))==0)
THEN
179 IF (tagbuf(ix(l,j)+npart)==0)
THEN
180 tagbuf(ix(l,j)+npart)= -1
266 use element_mod ,
only : nixs
270#include "implicit_f.inc"
274#include "com04_c.inc"
278 INTEGER IXS(NIXS,*),IPARTS(*),IXS10(6,*),IXS20(12,*),IXS16(8,*),
279 1 TAGBUF(*),FLAG,IDOM
295 IF (tagbuf(iparts(j)) == 0)
THEN
297 IF (tagbuf(ixs(l,j)+np)<2) tagbuf(ixs(l,j)+np)=0
304 IF (tagbuf(iparts(j)) == 0)
THEN
306 IF (tagbuf(ixs(l,j)+np)<2) tagbuf(ixs(l,j)+np)=0
309 IF (ixs10(l,i) /= 0)
THEN
310 IF (tagbuf(ixs10(l,i)+np)<2) tagbuf(ixs10(l,i)+np)=0
317 j = i + numels8 + numels10
318 IF (tagbuf(iparts(j)) == 0)
THEN
320 IF (tagbuf(ixs(l,j)+np)<2) tagbuf(ixs(l,j)+np)=0
323 IF (ixs20(l,i) /= 0)
THEN
324 IF (tagbuf(ixs20(l,i)+np)<2) tagbuf(ixs20(l,i)+np)=0
331 j = i + numels8 + numels10 + numels20
332 IF (tagbuf(iparts(j)) == 0)
THEN
334 IF (tagbuf(ixs(l,j)+np)<2) tagbuf(ixs(l,j)+np)=0
337 IF (ixs16(l,i) /= 0)
THEN
338 IF (tagbuf(ixs16(l,i)+np)<2) tagbuf(ixs16(l,i)+np)=0
348 ELSEIF (flag==0)
THEN
351 IF (tagbuf(iparts(j)) == 1)
THEN
353 IF (tagbuf(ixs(l,j)+np)<2) tagbuf(ixs(l,j)+np)=1
360 IF (tagbuf(iparts(j)) == 1)
THEN
362 IF (tagbuf(ixs(l,j)+np)<2) tagbuf(ixs(l,j)+np)=1
365 IF (ixs10(l,i) /= 0)
THEN
366 IF (tagbuf(ixs10(l,i)+np)<2) tagbuf(ixs10(l,i)+np)=1
373 j = i + numels8 + numels10
374 IF (tagbuf(iparts(j)) == 1)
THEN
376 IF (tagbuf(ixs(l,j)+np)<2) tagbuf(ixs(l,j)+np)=1
379 IF (ixs20(l,i) /= 0)
THEN
380 IF (tagbuf(ixs20(l,i)+np)<2) tagbuf(ixs20(l,i)+np)=1
387 j = i + numels8 + numels10 + numels20
388 IF (tagbuf(iparts(j)) == 1)
THEN
390 IF (tagbuf(ixs(l,j)+np)<2) tagbuf(ixs(l,j)+np)=1
393 IF (ixs16(l,i) /= 0)
THEN
394 IF (tagbuf(ixs16(l,i)+np)<2) tagbuf(ixs16(l,i)+np)=1
404 ELSEIF (flag==1)
THEN
407 IF (tagbuf(iparts(j)) == 0)
THEN
409 IF (tagbuf(ixs(l,j)+np)==1) tagbuf(ixs(l,j)+np)=1+idom
410 IF (tagbuf(ixs(l,j)+np)>1)
THEN
411 IF (tagbuf(ixs(l,j)+np)/=(1+idom))
THEN
416 . i2=tagbuf(ixs(l,j)+np)-1,
421 IF (tagbuf(ixs(l,j)+np)<1) tagbuf(ixs(l,j)+np)=-idom
429 IF (tagbuf(iparts(j)) == 0)
THEN
431 IF (tagbuf(ixs(l,j)+np)==1) tagbuf(ixs(l,j)+np)=1+idom
432 IF (tagbuf(ixs(l,j)+np)>1)
THEN
433 IF (tagbuf(ixs(l,j)+np)/=(1+idom))
THEN
438 . i2=tagbuf(ixs(l,j)+np)-1,
443 IF (TAGBUF(IXS(L,J)+NP)<1) TAGBUF(IXS(L,J)+NP)=-IDOM
446 IF (IXS10(L,I) /= 0) THEN
447 IF (TAGBUF(IXS10(L,I)+NP)==1) TAGBUF(IXS10(L,I)+NP)=1+IDOM
448 IF (TAGBUF(IXS10(L,I)+NP)>1) THEN
449 IF (TAGBUF(IXS10(L,I)+NP)/=(1+IDOM)) THEN
451 CALL ANCMSG(MSGID=838,
454 . I2=TAGBUF(IXS10(L,I)+NP)-1,
459 IF (TAGBUF(IXS10(L,I)+NP)<1) TAGBUF(IXS10(L,I)+NP)=-IDOM
467 J = I + NUMELS8 + NUMELS10
468 IF (TAGBUF(IPARTS(J)) == 0)THEN
470 IF (TAGBUF(IXS(L,J)+NP)==1) TAGBUF(IXS(L,J)+NP)=1+IDOM
471 IF (TAGBUF(IXS(L,J)+NP)>1) THEN
472 IF (TAGBUF(IXS(L,J)+NP)/=(1+IDOM)) THEN
474 CALL ANCMSG(MSGID=838,
477 . I2=TAGBUF(IXS(L,J)+NP)-1,
482 IF (TAGBUF(IXS(L,J)+NP)<1) TAGBUF(IXS(L,J)+NP)=-IDOM
485 IF (IXS20(L,I) /= 0) THEN
486 IF (TAGBUF(IXS20(L,I)+NP)==1) TAGBUF(IXS20(L,I)+NP)=1+IDOM
487 IF (TAGBUF(IXS20(L,I)+NP)>1) THEN
488 IF (TAGBUF(IXS20(L,I)+NP)/=(1+IDOM)) THEN
490 CALL ANCMSG(MSGID=838,
493 . I2=TAGBUF(IXS20(L,I)+NP)-1,
498 IF (TAGBUF(IXS20(L,I)+NP)<1) TAGBUF(IXS20(L,I)+NP)=-IDOM
506 J = I + NUMELS8 + NUMELS10 + NUMELS20
507 IF (TAGBUF(IPARTS(J)) == 0) THEN
509 IF (TAGBUF(IXS(L,J)+NP)==1) TAGBUF(IXS(L,J)+NP)=1+IDOM
510 IF (TAGBUF(IXS(L,J)+NP)>1) THEN
511 IF (TAGBUF(IXS(L,J)+NP)/=(1+IDOM)) THEN
513 CALL ANCMSG(MSGID=838,
516 . I2=TAGBUF(IXS(L,J)+NP)-1,
521 IF (TAGBUF(IXS(L,J)+NP)<1) TAGBUF(IXS(L,J)+NP)=-IDOM
524 IF (IXS16(L,I) /= 0) THEN
525 IF (TAGBUF(IXS16(L,I)+NP)==1) TAGBUF(IXS16(L,I)+NP)=1+IDOM
526 IF (TAGBUF(IXS16(L,I)+NP)>1) THEN
527 IF (TAGBUF(IXS16(L,I)+NP)/=(1+IDOM)) THEN
529 CALL ANCMSG(MSGID=838,
532 . I2=TAGBUF(IXS16(L,I)+NP)-1,
537 IF (TAGBUF(IXS16(L,I)+NP)<1) TAGBUF(IXS16(L,I)+NP)=-IDOM
548 ELSEIF (FLAG==2) THEN
551 IF (IPARTS(J+NP)/=0)THEN
553 IF (TAGBUF(IXS(L,J)+NP)<3) TAGBUF(IXS(L,J)+NP)=2*IPARTS(J+NP)
560 IF (IPARTS(J+NP)/=0)THEN
562 IF (TAGBUF(IXS(L,J)+NP)<3) TAGBUF(IXS(L,J)+NP)=2*IPARTS(J+NP)
565 IF (IXS10(L,I) /= 0) THEN
566 IF (TAGBUF(IXS10(L,I)+NP)<3) TAGBUF(IXS10(L,I)+NP)=2*IPARTS(J+NP)
573 J = I + NUMELS8 + NUMELS10
574 IF (IPARTS(J+NP)/=0)THEN
576 IF (TAGBUF(IXS(L,J)+NP)<3) TAGBUF(IXS(L,J)+NP)=2*IPARTS(J+NP)
579 IF (IXS20(L,I) /= 0) THEN
580 IF (TAGBUF(IXS20(L,I)+NP)<3) TAGBUF(IXS20(L,I)+NP)=2*IPARTS(J+NP)
587 J = I + NUMELS8 + NUMELS10 + NUMELS20
588 IF (IPARTS(J+NP)/=0)THEN
590 IF (TAGBUF(IXS(L,J)+NP)<3) TAGBUF(IXS(L,J)+NP)=2*IPARTS(J+NP)
593 IF (IXS16(L,I) /= 0) THEN
594 IF (TAGBUF(IXS16(L,I)+NP)<3) TAGBUF(IXS16(L,I)+NP)=2*IPARTS(J+NP)
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)
subroutine tagnod_r2r(ix, nix, nix1, nix2, numel, iparte, tagbuf, npart, flag, idom)