Download kinpro.eso

Back to the list

   1 : C KINPRO    SOURCE    CHAT      05/01/13    00:56:41     5004           
   2 :       SUBROUTINE KINPRO (U1,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 pas FACGEN
   8 : C     DETERMINATION DES CELLULES VUES A PARTIR DE LA DIRECTION U1               
   9 : C                                                                               
  10 : C----------------------------------------------------------------------
  11 :       DIMENSION U1(1)                                                           
  12 : -INC TFFOR3D
  13 : C                                                                               
  14 :       DO 1 NF = 1,KFC                                                           
  15 :       DO 2  I =  1,NRES                                                         
  16 :       DO 3  J =  1,NRES                                                         
  17 :           B = 0.                                                                
  18 :           DO 4    IES = 1,KES                                                   
  19 :           B = B + U1(IES)*KSI(NF,IES)*V(KRO(NF,IES),I,J)                        
  20 :  4        CONTINUE                                                              
  21 :           IF (B.GT.0.000001) THEN                                               
  22 :              PSC(NF,I,J) = B                                                    
  23 :           ELSE                                                                  
  24 :              PSC(NF,I,J) = -1.                                                  
  25 :           ENDIF                                                                 
  26 :  3          CONTINUE                                                            
  27 :  2          CONTINUE                                                            
  28 :  1          CONTINUE                                                            
  29 :       RETURN                                                                    
  30 :       END                                                                       
  31 :  
  32 :  

© Cast3M 2003 - All rights reserved.
Disclaimer