C KREM3F    SOURCE    CHAT      05/01/13    01:06:16     5004           
      SUBROUTINE KREM3F(K2,KG,NG,KC1,NA2,NF2,NF3,NF1,NA3,NA1                    
     -                 ,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)                         
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)                                                 
C----------------------------------------------------------------------
C Calcul des facteurs de forme en 3D
C Sp appele par krem3g
C     3 FACES DIFFERENTES NON PARALLELES                                        
C     ----------------------------------                                        
C                                                                               
C     IL EXISTE UNE FACE GENEREE  PAR L'ARETE NA2                               
C     -------------------------------------------                               
C                                                                               
C----------------------------------------------------------------------
      DIMENSION KG(2),KH(2),KC1(2)                                              
-INC TFFOR3D

      SEGMENT SKCEL                                                             
      INTEGER  KBCEL(NR,NR),IINT(2,NSTAC),IS(NSTAC),JS(NSTAC)                   
      REAL*8   RMAX                                                             
      ENDSEGMENT                                                                
      IF (IFA(3,NA2).EQ.NF1) THEN                                               
                                                                                
      CALL KINITB(KBCEL,NRES,NRES,KG,NG)                                        
      CALL KREMPA(K2,KG,NG,NF2,NA1,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KREMPA(K2,KG,NG,NF2,NA2,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KREMPI(NIN,K2,KG,NG,NF2,C,U2,SHC3D,SKCEL,SKBUFF)                     
                                                                                
      CALL KINITB(KBCEL,NRES,NRES,KG,NG)                                        
      CALL KREMPA(K2,KG,NG,NF3,NA2,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KREMPA(K2,KG,NG,NF3,NA3,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KREMPI(NIN,K2,KG,NG,NF3,C,U2,SHC3D,SKCEL,SKBUFF)                     
                                                                                
      CALL KINITB(KBCEL,NRES,NRES,KG,NG)                                        
      CALL KREMPA(K2,KG,NG,NF1,NA1,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KREMPA(K2,KG,NG,NF1,NA3,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      KH(1) = (KG(1)/2 + KC1(1))/2                                              
      KH(2) = (KG(2)/2 + KC1(2))/2                                              
      NG = 0                                                                    
      KG(1) = 0                                                                 
      KG(2) = 0                                                                 
      CALL KREMPA(K2,KG,NG,NFN,NA2,3,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      KH(1) = (KG(1) + KH(1))/2                                                 
      KH(2) = (KG(2) + KH(2))/2                                                 
      CALL KREMPI(NIN,K2,KH,1,NF1,C,U2,SHC3D,SKCEL,SKBUFF)                      
                                                                                
      ELSE                                                                      
      NFF = IFA(3,NA2)                                                          
                                                                                
C     WRITE(6,*) ' KREMPY NFF ',NFF                                             
      CALL KINITB(KBCEL,NRES,NRES,KG,NG)                                        
      CALL KREMPA(K2,KG,NG,NF2,NA1,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KREMPA(K2,KG,NG,NF2,NA2,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KFINSO(NF1,NF3,KA,IM,KH)                                             
      CALL KFINSO(NFF,NF3,KA,IM,KG)                                             
      KH(1) = (KG(1) + KH(1))/2                                                 
      KH(2) = (KG(2) + KH(2))/2                                                 
      CALL KREMPI(NIN,K2,KH,1,NF2,C,U2,SHC3D,SKCEL,SKBUFF)                      
                                                                                
      CALL KINITB(KBCEL,NRES,NRES,KG,NG)                                        
      CALL KREMPA(K2,KG,NG,NF3,NA2,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KREMPA(K2,KG,NG,NF3,NA3,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KFINSO(NF1,NF2,KA,IM,KH)                                             
      CALL KFINSO(NFF,NF2,KA,IM,KG)                                             
      KH(1) = (KG(1) + KH(1))/2                                                 
      KH(2) = (KG(2) + KH(2))/2                                                 
      CALL KREMPI(NIN,K2,KH,1,NF3,C,U2,SHC3D,SKCEL,SKBUFF)                      
                                                                                
      CALL KINITB(KBCEL,NRES,NRES,KG,NG)                                        
      CALL KREMPA(K2,KG,NG,NF1,NA1,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KREMPA(K2,KG,NG,NF1,NA3,1,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KFINSO(NF2,NF3,KA,IM,KH)                                             
      KH(1) = (KG(1)/2 + KH(1))/2                                               
      KH(2) = (KG(2)/2 + KH(2))/2                                               
      CALL KREMPI(NIN,K2,KH,1,NF1,C,U2,SHC3D,SKCEL,SKBUFF)                      
C                                                                               
      CALL KINITB(KBCEL,NRES,NRES,KG,NG)                                        
      CALL KREMPA(K2,KG,NG,NFF,NA2,3,C,U2,SHC3D,SKCEL,SPROJA,SKBUFF)            
      CALL KFINSO(NF2,NF3,KA,IM,KG)                                             
      CALL KREMPI(NIN,K2,KG,1,NFF,C,U2,SHC3D,SKCEL,SKBUFF)                      
      ENDIF                                                                     
                                                                                
      RETURN                                                                    
      END