Esope source files

   1 : C NKSNST    SOURCE    CHAT      05/01/13    01:57:17     5004           
   2 :       SUBROUTINE NKSNST(XGENE,DJAC,VKS,LRE,NST,REL)                             
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C      INTEGRE LE PRODUIT NT KS N                                               
   6 : C          POUR LES ELEMENTS DE SURFACE LIBRE                                   
   7 : C                                                                               
   8 : C      ROUTINE FORTRAN PUR                                                      
   9 : C                                                                               
  10 : C     JACQUELINE BROCHARD  DECEMBRE 85                                          
  11 : C=======================================================================        
  12 : C     XGENE(NST,LRE)=MATRICE N                                                  
  13 : C     DJAC =JACOBIEN                                                            
  14 : C     VKS=RHO*G                                                                 
  15 : C     LRE=TAILLE DE LA MATRICE DE RIGIDITE                                      
  16 : C     NST=NOMBRE DE LIGNES DE LA MATRICE XGENE                                  
  17 : C     REL=MATRICE DE RIGIDITE                                                   
  18 : C=======================================================================        
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  21 :       DIMENSION REL(LRE,1),XGENE(NST,1)                                         
  22 :       DO   300 IX=1,LRE                                                         
  23 :       DO   300 IY=1,IX                                                          
  24 : C                                                                               
  25 : C     CALCUL DE NTKSN POUR LES DDL IX ET IY                                     
  26 : C                                                                               
  27 :       REL(IX,IY)=REL(IX,IY)+DJAC*VKS*XGENE(3,IY)*XGENE(3,IX)                    
  28 :       REL(IY,IX)=REL(IX,IY)                                                     
  29 :   300 CONTINUE                                                                  
  30 :       RETURN                                                                    
  31 :       END                                                                       
  32 :