Esope source files

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