Esope source files

   1 : C KVERIF    SOURCE    CHAT      05/01/13    01:11:33     5004           
   2 :       SUBROUTINE KVERIF (SHC3D,SKBUFF)                                          
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
   5 : C----------------------------------------------------------------------
   6 : C Calcul des facteurs de forme en 3D
   7 : C Sp appele par FACGEN
   8 : C    VERIFICATION DU BILAN AUTOUR D'UNE DIRECTION DONNEE                        
   9 : C                                                                               
  10 : C----------------------------------------------------------------------
  11 : -INC TFFOR3D
  12 :                                                                                 
  13 :       F = 0                                                                     
  14 :       NC = 0                                                                    
  15 :       DO 500 NF = 1,KFC                                                         
  16 :       DO 501 I = 1,NRES                                                         
  17 :       DO 502 J = 1,NRES                                                         
  18 :          IF (PSC(NF,I,J).GT.-1.) THEN                                           
  19 :            NC = NC + 1                                                          
  20 :            F = F + PSC(NF,I,J)*G(I,J)                                           
  21 :          ENDIF                                                                  
  22 :  502  CONTINUE                                                                  
  23 :  501  CONTINUE                                                                  
  24 :  500  CONTINUE                                                                  
  25 :       ERR = DABS(1.-F)                                                          
  26 :       IF (ERR.GT.0.01.OR.KIMP.GE.4) THEN
  27 :            WRITE(6,1000) NC,F                                                   
  28 :  1000 FORMAT(1X,'VERIF: NCEL ',I6,' F ',F10.5)                                  
  29 :       ENDIF                                                                     
  30 :       RETURN                                                                    
  31 :       END                                                                       
  32 :  
  33 :