Esope source files

   1 : C BMLBST    SOURCE    CHAT      05/01/12    21:41:28     5004           
   2 : C BMLBST    SOURCE    BROC      85/07/01    23:05:03                            
   3 :       SUBROUTINE BMLBST(BLX,BLY,BLZ,BLT,DJAC,VML22,IFOU,LRE,NDDL,REL)           
   4 : C=======================================================================        
   5 : C                                                                               
   6 : C    INTEGRE LE PRODUIT BLT ML BL                                               
   7 : C       POUR LES ELEMENTS LIQUIDES                                              
   8 : C                                                                               
   9 : C          ROUTINE FORTRAN PUR                                                  
  10 : C=======================================================================        
  11 : C  INPUT                                                                        
  12 : C     BLX =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /X             
  13 : C     BLY =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /Y             
  14 : C     BLZ =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /Z             
  15 : C     BLT =MATRICE CONTENANT LES DERIVEES DES FONCTIONS DE FORME /TETA          
  16 : C     DJAC=JACOBIEN                                                             
  17 : C     VML22=-(COEFPI**2)/RHO                                                    
  18 : C     IFOU=IFOUR DE CCOPTIO                                                     
  19 : C     LRE =NOMBRE DE COLONNES DES MATRICES BL                                   
  20 : C     NDDL=NOMBRE DE DEGRES DE LIBERTE                                          
  21 : C  OUTPUT                                                                       
  22 : C     REL(LRE,LRE) =MATRICE MASSE                                               
  23 : C=======================================================================        
  24 :       IMPLICIT INTEGER(I-N)
  25 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  26 :       DIMENSION BLX(NDDL,1),BLY(NDDL,1),BLZ(NDDL,1),BLT(NDDL,1)                 
  27 :       DIMENSION REL(LRE,1)                                                      
  28 :       IFR=IFOU+4                                                                
  29 :       GOTO (666,10,10,10,30,40) ,IFR                                            
  30 : C                                                                               
  31 : C     ELEMENTS LIQUIDES 2DIM   CONT PLANES OU DEFO PLANES OU AXIS               
  32 : C                                                                               
  33 :   10  CONTINUE                                                                  
  34 :       DO 100 IX=1,LRE                                                           
  35 :       DO 100 IY=1,IX                                                            
  36 :       REL(IX,IY)=REL(IX,IY)+DJAC*VML22*(BLX(2,IY)*BLX(2,IX)+                    
  37 :      1                        BLY(2,IY)*BLY(2,IX))                              
  38 :       REL(IY,IX)=REL(IX,IY)                                                     
  39 : 100   CONTINUE                                                                  
  40 :       GOTO  666                                                                 
  41 : C                                                                               
  42 : C     ELEMENTS LIQUIDES BIDIM FOURIER                                           
  43 : C                                                                               
  44 :   30  CONTINUE                                                                  
  45 :       DO 300 IX=1,LRE                                                           
  46 :       DO 300 IY=1,IX                                                            
  47 :       REL(IX,IY)=REL(IX,IY)+DJAC*VML22*(BLX(2,IY)*BLX(2,IX)+                    
  48 :      1                        BLY(2,IY)*BLY(2,IX)+                              
  49 :      2                        BLT(2,IY)*BLT(2,IX))                              
  50 :       REL(IY,IX)=REL(IX,IY)                                                     
  51 : 300   CONTINUE                                                                  
  52 :       GOTO  666                                                                 
  53 : C                                                                               
  54 : C     ELEMENTS LIQUIDES TRIDIM                                                  
  55 : C                                                                               
  56 :   40  CONTINUE                                                                  
  57 :       DO 400 IX=1,LRE                                                           
  58 :       DO 400 IY=1,IX                                                            
  59 :       REL(IX,IY)=REL(IX,IY)+DJAC*VML22*(BLX(2,IY)*BLX(2,IX)+                    
  60 :      1                        BLY(2,IY)*BLY(2,IX)+                              
  61 :      2                        BLZ(2,IY)*BLZ(2,IX))                              
  62 :       REL(IY,IX)=REL(IX,IY)                                                     
  63 : 400   CONTINUE                                                                  
  64 :       GOTO  666                                                                 
  65 : C                                                                               
  66 :  666  CONTINUE                                                                  
  67 :       RETURN                                                                    
  68 :       END                                                                       
  69 :