Esope source files

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