Esope source files

   1 : C NMLNST    SOURCE    CHAT      05/01/13    01:57:46     5004           
   2 :       SUBROUTINE NMLNST(XGENE,DJAC,VML12,LRE,NST,REL)                           
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C      INTEGRE LE PRODUIT NT ML 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     VML12=-(COEFPR*COEFPI)/(RHO*C**2)                                         
  14 : C     LRE=TAILLE DE LA MATRICE DE MASSE                                         
  15 : C     NST=NOMBRE DE LIGNES DE LA MATRICE NGENE                                  
  16 : C     REL=MATRICE MASSE                                                         
  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 NTMLN POUR LES DDL IX ET IY                                     
  25 : C                                                                               
  26 :       REL(IX,IY)=REL(IX,IY)+DJAC*VML12*(XGENE(1,IY)*XGENE(2,IX)+                
  27 :      1                          XGENE(2,IY)*XGENE(1,IX))                        
  28 :       REL(IY,IX)=REL(IX,IY)                                                     
  29 :   300 CONTINUE                                                                  
  30 :       RETURN                                                                    
  31 :       END                                                                       
  32 :