C IFACE4    SOURCE    JC220346  16/11/29    21:15:18     9221           
C---------------------------------------------------------------------|
C                                                                     |
       FUNCTION IFACE4(IP,JP,KP,LP)
C                                                                     |
C      CETTE FONCTION   RECHERCHE LA FACETTE IFC QUI CONTIENT         |
C      LES POINTS IP, JP, KP ET LP.                                   |
C      SI JFC N'EXISTE PAS, ON RECHERCHE L'EXISTENCE EVENTUELLE       |
C      D'UNE FACETTE TRIANGULAIRE (OU CARRE) CONTENANT 3 DES 4 POINTS |
C      ELLE EST RAMENEE EN NEGATIF                                    |
C      - IRET=0 SI AUCUNE FACETTE N'EXISTE                            |
C                                                                     |
C---------------------------------------------------------------------|
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
-INC TDEMAIT

-INC PPARAM
-INC CCOPTIO
C
       DO 100 I=1,40
           IF=NPF(I,IP)
              IF (IF.EQ.0) GOTO 140
              DO 110 J=1,40
                     IF (NPF(J,JP).EQ.0) GOTO 100
                     IF (IF.NE.NPF(J,JP)) GOTO 110
                     DO 120 K=1,40
                            IF (NPF(K,KP).EQ.0) GOTO 100
                            IF (IF.NE.NPF(K,KP)) GOTO 120
                            DO 130 L=1,40
                                   IF (NPF(L,LP).EQ.0) GOTO 200
                                   IF (IF.NE.NPF(L,LP)) GOTO 130
                                   IFACE4=IF
                                   RETURN
130                         CONTINUE
                         GOTO 200
120                  CONTINUE
110           CONTINUE
100    CONTINUE
C
140    IFC=IFACE3(IP,JP,KP)
       IF (IFC.EQ.0) GOTO 150
         IFACE4=-IFC
       RETURN
C
150    IFC=IFACE3(IP,JP,LP)
       IF (IFC.EQ.0) GOTO 160
         IFACE4=-IFC
       RETURN
C
160    IFC=IFACE3(IP,KP,LP)
       IF (IFC.EQ.0) GOTO 170
         IFACE4=-IFC
       RETURN
C
170    IFC=IFACE3(JP,KP,LP)
       IF (IFC.EQ.0) GOTO 180
         IFACE4=-IFC
       RETURN
C
 180   IFACE4=0
       RETURN
C
200    CONTINUE
         IFACE4=-IF
       RETURN
C
C      FIN DE LA SUBROUTINE FACE4
       END


 
