Sources Esope | Cast3M

kinpro
C KINPRO    SOURCE    CHAT      05/01/13    00:56:41     5004
      SUBROUTINE KINPRO (U1,SHC3D,SKBUFF)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C----------------------------------------------------------------------
C Calcul des facteurs de forme en 3D
C Sp appele pas FACGEN
C     DETERMINATION DES CELLULES VUES A PARTIR DE LA DIRECTION U1
C
C----------------------------------------------------------------------
      DIMENSION U1(1)
-INC TFFOR3D
C
      DO 1 NF = 1,KFC
      DO 2  I =  1,NRES
      DO 3  J =  1,NRES
          B = 0.
          DO 4    IES = 1,KES
          B = B + U1(IES)*KSI(NF,IES)*V(KRO(NF,IES),I,J)
 4        CONTINUE
          IF (B.GT.0.000001) THEN
             PSC(NF,I,J) = B
          ELSE
             PSC(NF,I,J) = -1.
          ENDIF
 3          CONTINUE
 2          CONTINUE
 1          CONTINUE
      RETURN
      END