C IFACE4    SOURCE    PV        07/10/05    21:15:55     5890           
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 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