1 : C TRJIEL SOURCE CHAT 05/01/13 03:49:52 5004 2 : SUBROUTINE TRJIEL(IEL1,IEL2,ICONT,NF,IFACEL,IZCENT,IELTFA) 3 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 4 : C 5 : C CHERCHE LE NUMERO DE L ELEMENT IEL2 AYANT LA FACE NF COMMUNE AVEC IEL1 6 : C A PARTIR DES ELEMENTS DE LA TABLE GENEREE PAR KDOM 7 : C 8 : C ENTRES 9 : C IEL1 : NUMERO DE L ELEMENT 10 : C ICONT : NUMERO LOCAL DE LA FACE DANS L ELEMENT IEL1 11 : C IFACEL : POINTEUR DU MAILLAGE "FACEL" (KDOM) 12 : C IZCENT : POINTEUR DU MAILLAGE "CENTRE" (KDOM) 13 : C IELTFA : POINTEUR DU MAILLAGE "ELTFA" (KDOM) 14 : C ON SUPPOSE CES SEGMENTS ACTIFS 15 : C SORTIES 16 : C IEL2 : NUMERO DE L ELEMENT AYANT LA FACE NF COMMUNE AVEC IEL1 17 : C NF : SERA LE NUMERO DE LA FACE ICONT DANS LE MAILLAGE GENERAL 18 : CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 19 : IMPLICIT INTEGER(I-N) 20 : IMPLICIT REAL*8 (A-H,O-Z) 21 : -INC SMELEME 22 : POINTEUR IELTFA.MELEME ,IFACEL.MELEME,IZCENT.MELEME 23 : C 24 : NEL0=0 25 : IEL2=0 25 : CALL MELNEL(IEL1,IELTFA,IPT1,NEL0,1) 26 : NF=IPT1.NUM(ICONT,IEL1-NEL0) 27 : C write(6,*)' ieltrj ', nf,iel1,icont,nel0 28 : NFC=IFACEL.NUM(/2) 29 : NEL=IZCENT.NUM(/2) 30 : DO 50 K=1,NFC 31 : IF(IFACEL.NUM(2,K).EQ.NF)THEN 32 : IF(IFACEL.NUM(1,K).EQ.IZCENT.NUM(1,IEL1))THEN 33 : J2=3 34 : ELSEIF(IFACEL.NUM(3,K).EQ.IZCENT.NUM(1,IEL1))THEN 35 : J2=1 36 : ELSE 37 : GO TO 35 38 : ENDIF 39 : DO 30 L=1,NEL 40 : IF(IFACEL.NUM(J2,K).EQ.IZCENT.NUM(1,L))THEN 41 : IF(L.EQ.IEL1)GO TO 25 42 : IEL2=L 43 : ENDIF 44 : 25 CONTINUE 45 : 30 CONTINUE 46 : 35 CONTINUE 47 : ENDIF 48 : C 49 : 50 CONTINUE 50 : RETURN 51 : END 52 :
© Cast3M 2003 - All rights reserved.
Disclaimer