dmasanic.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dmasanic (elbuf_tab, x, d, geo, iparg, ixq, ixc, ixtg, mas, pm, el2fa, nbf)

Function/Subroutine Documentation

dmasanic()

subroutine dmasanic	(	type (elbuf_struct_), dimension(ngroup), target	elbuf_tab,
			x,
			d,
			geo,
		integer, dimension(nparg,*)	iparg,
		integer, dimension(nixq,*)	ixq,
		integer, dimension(nixc,*)	ixc,
		integer, dimension(nixtg,*)	ixtg,
			mas,
			pm,
		integer, dimension(*)	el2fa,
		integer	nbf )

Definition at line 29 of file dmasanic.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod            
C-----------------------------------------------                        
C   I m p l i c i t   T y p e s                                         
C-----------------------------------------------                        
#include      "implicit_f.inc"                                          
C-----------------------------------------------                        
C   C o m m o n   B l o c k s                                           
C-----------------------------------------------                        
#include      "com01_c.inc"                                             
#include      "com04_c.inc"                                             
#include      "param_c.inc"
C-----------------------------------------------                        
C   D u m m y   A r g u m e n t s                                       
C-----------------------------------------------                        
C     REAL                                                    
      my_real                                                           
     .   mas(*) ,pm(npropm,*),geo(npropg,*),x(3,*),            
     .   d(3,*)                                                         
      INTEGER IPARG(NPARG,*),IXC(NIXC,*),IXTG(NIXTG,*),   
     .        IXQ(NIXQ,*),EL2FA(*),NBF                                  
      TYPE (ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
C-----------------------------------------------                        
C   L o c a l   V a r i a b l e s                                       
C-----------------------------------------------                        
C     REAL                                                    
      my_real                                                           
     .   off,a0,thk0,rho0,xx1,xx2,xx3,yy1,yy2,yy3,zz1,zz2,zz3,VALUE                       
      INTEGER I, NG, NEL, NFT, ITY, LFT, IALEL,MT,LLT, 
     .        N,N1,N2,N3,N4,NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,NN9,NN10,                 
     .        NFT_FA,N_FA                        
      TYPE(G_BUFEL_)  ,POINTER :: GBUF     
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
      nn1 = 1                                                           
      nn2 = 1                                                           
      nn3 = 1                                                           
      nn4 = nn3 + numelq                                                
      nn5 = nn4 + numelc                                                
      nn6 = nn5 + numeltg                                               
      nn7 = nn6
      nn8 = nn7
      nn9 = nn8
C-----------------------------------------------                        
      DO 490 ng=1,ngroup                                                
        gbuf => elbuf_tab(ng)%GBUF
        nel   =iparg(2,ng)                                              
        ity   =iparg(5,ng)                                              
        nft   =iparg(3,ng)                                              
        lft=1                                                           
        llt=nel                                                         
        nft_fa = nft                                                    
C-----------------------------------------------                        
C       QUAD                                                            
C-----------------------------------------------                        
        IF(ity==2)THEN                                                
          ialel=(iparg(7,ng)+iparg(11,ng))                               
C                                                                       
          DO i=lft,llt                                                   
               n = i + nft                                              
               n_fa = i + nft_fa                                        
               IF(ialel==0)THEN                                       
                 mt=ixq(1,n)                                            
                 VALUE = pm(89,mt)* gbuf%VOL(i)                     
               ELSE                                                     
                 off = min(gbuf%OFF(i),one)                             
                 VALUE=  gbuf%RHO(i)*gbuf%VOL(i)*off                    
               ENDIF                                                    
               mas(el2fa(nn3+n_fa)) = VALUE                             
          ENDDO                                                          
C-----------------------------------------------                        
C       COQUES 4 N                                                      
C-----------------------------------------------                        
        ELSEIF(ity==3)THEN                                            
C                                                                       
          DO i=lft,llt                                                  
            n = i + nft                                                 
            n_fa = i + nft_fa                                           
            rho0 = pm(1,ixc(1,n))                                       
            thk0 = geo(1,ixc(6,n))                                      
            n1 = ixc(2,n)                                               
            n2 = ixc(3,n)                                               
            n3 = ixc(4,n)                                               
            n4 = ixc(5,n)                                               
            xx1 = x(1,n3)-d(1,n3)-x(1,n1)+d(1,n1)                       
            yy1 = x(2,n3)-d(2,n3)-x(2,n1)+d(2,n1)                       
            zz1 = x(3,n3)-d(3,n3)-x(3,n1)+d(3,n1)                       
            xx2 = x(1,n4)-d(1,n4)-x(1,n2)+d(1,n2)                       
            yy2 = x(2,n4)-d(2,n4)-x(2,n2)+d(2,n2)                       
            zz2 = x(3,n4)-d(3,n4)-x(3,n2)+d(3,n2)                       
            xx3 = yy1*zz2 - zz1*yy2                                     
            yy3 = zz1*xx2 - xx1*zz2                                     
            zz3 = xx1*yy2 - yy1*xx2                                     
            a0  = half*sqrt(xx3*xx3 + yy3*yy3 + zz3*zz3)              
            mas(el2fa(nn4+n_fa)) = rho0*thk0*a0                         
          ENDDO                                                         
C-----------------------------------------------                        
C       COQUES 3 N                                                      
C-----------------------------------------------                        
        ELSEIF(ity==7)THEN                                            
C                                                                       
          DO i=lft,llt                                                  
            n = i + nft                                                 
            n_fa = i + nft_fa                                           
            rho0 = pm(1,ixtg(1,n))                                      
            thk0 = geo(1,ixtg(5,n))                                     
            n1 = ixtg(2,n)                                              
            n2 = ixtg(3,n)                                              
            n3 = ixtg(4,n)                                              
            xx1 = x(1,n2)-d(1,n2)-x(1,n1)+d(1,n1)                       
            yy1 = x(2,n2)-d(2,n2)-x(2,n1)+d(2,n1)                       
            zz1 = x(3,n2)-d(3,n2)-x(3,n1)+d(3,n1)                       
            xx2 = x(1,n3)-d(1,n3)-x(1,n1)+d(1,n1)                       
            yy2 = x(2,n3)-d(2,n3)-x(2,n1)+d(2,n1)                       
            zz2 = x(3,n3)-d(3,n3)-x(3,n1)+d(3,n1)                       
            xx3 = yy1*zz2 - zz1*yy2                                     
            yy3 = zz1*xx2 - xx1*zz2                                     
            zz3 = xx1*yy2 - yy1*xx2                                     
            a0  = half*sqrt(xx3*xx3 + yy3*yy3 + zz3*zz3)              
            mas(el2fa(nn5+n_fa)) = rho0*thk0*a0                         
          ENDDO                                                         
C                                                                       
        ELSE                                                            
        ENDIF                                                           
C-----------------------------------------------                        
 490  CONTINUE                                                          
C-----------------------------------------------                        
C                                                                       
      RETURN