1 : C IFACE4 SOURCE PV 07/10/05 21:15:55 5890 2 : C---------------------------------------------------------------------| 3 : C | 4 : FUNCTION IFACE4(IP,JP,KP,LP) 5 : C | 6 : C CETTE FONCTION RECHERCHE LA FACETTE IFC QUI CONTIENT | 7 : C LES POINTS IP, JP, KP ET LP. | 8 : C SI JFC N'EXISTE PAS, ON RECHERCHE L'EXISTENCE EVENTUELLE | 9 : C D'UNE FACETTE TRIANGULAIRE (OU CARRE) CONTENANT 3 DES 4 POINTS | 10 : C ELLE EST RAMENEE EN NEGATIF | 11 : C - IRET=0 SI AUCUNE FACETTE N'EXISTE | 12 : C | 13 : C---------------------------------------------------------------------| 14 : C 15 : IMPLICIT INTEGER(I-N) 16 : IMPLICIT REAL*8(A-H,O-Z) 17 : -INC TDEMAIT 18 : -INC CCOPTIO 19 : C 20 : DO 100 I=1,40 21 : IF=NPF(I,IP) 22 : IF (IF.EQ.0) GOTO 140 23 : DO 110 J=1,40 24 : IF (NPF(J,JP).EQ.0) GOTO 100 25 : IF (IF.NE.NPF(J,JP)) GOTO 110 26 : DO 120 K=1,40 27 : IF (NPF(K,KP).EQ.0) GOTO 100 28 : IF (IF.NE.NPF(K,KP)) GOTO 120 29 : DO 130 L=1,40 30 : IF (NPF(L,LP).EQ.0) GOTO 200 31 : IF (IF.NE.NPF(L,LP)) GOTO 130 32 : IFACE4=IF 33 : RETURN 34 : 130 CONTINUE 35 : GOTO 200 36 : 120 CONTINUE 37 : 110 CONTINUE 38 : 100 CONTINUE 39 : C 40 : 140 IFC=IFACE3(IP,JP,KP) 41 : IF (IFC.EQ.0) GOTO 150 42 : IFACE4=-IFC 43 : RETURN 44 : C 45 : 150 IFC=IFACE3(IP,JP,LP) 46 : IF (IFC.EQ.0) GOTO 160 47 : IFACE4=-IFC 48 : RETURN 49 : C 50 : 160 IFC=IFACE3(IP,KP,LP) 51 : IF (IFC.EQ.0) GOTO 170 52 : IFACE4=-IFC 53 : RETURN 54 : C 55 : 170 IFC=IFACE3(JP,KP,LP) 56 : IF (IFC.EQ.0) GOTO 180 57 : IFACE4=-IFC 58 : RETURN 59 : C 60 : 180 IFACE4=0 61 : RETURN 62 : C 63 : 200 CONTINUE 64 : IFACE4=-IF 65 : RETURN 66 : C 67 : C FIN DE LA SUBROUTINE FACE4 68 : END 69 : 70 :
© Cast3M 2003 - All rights reserved.
Disclaimer