dmasanis.F File Reference

#include "implicit_f.inc"
#include "vect01_c.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "scr17_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dmasanis (elbuf_tab, iparg, ixs, mas, pm, el2fa, nbf, ipart, ipartsp, isph3d)

Function/Subroutine Documentation

dmasanis()

subroutine dmasanis	(	type (elbuf_struct_), dimension(ngroup), target	elbuf_tab,
		integer, dimension(nparg,*)	iparg,
		integer, dimension(nixs,*)	ixs,
			mas,
			pm,
		integer, dimension(*)	el2fa,
		integer	nbf,
		integer, dimension(lipart1,*)	ipart,
		integer, dimension(*)	ipartsp,
		integer	isph3d )

Definition at line 31 of file dmasanis.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      "vect01_c.inc"                                            
#include      "mvsiz_p.inc"                                            
#include      "com01_c.inc"                                             
#include      "com04_c.inc"                                             
#include      "scr17_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,*)                                              
      INTEGER IPARG(NPARG,*),IXS(NIXS,*),EL2FA(*),NBF,IPART(LIPART1,*),     
     .        IPARTSP(*),ISPH3D                                         
      TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
C-----------------------------------------------                        
C   L o c a l   V a r i a b l e s                                       
C-----------------------------------------------                        
      my_real                                                           
     .   EVAR(MVSIZ),                                                   
     .   OFF,VALUE
      INTEGER I, NG, NEL, ISS, ISC,N, J, MLW,                                          
     .        NN, K1, K2,JTURB,MT, IALEL,IPID,                     
     .        N1,N2,N3,N4,NN1,NN2,NN3,       
     .        OFFSET,NFT_FA,N_FA,                       
     .        INOD, ISOLNOD, IPRT,                                      
     .        JHBE, JIVF, JCLOSE, JPLASOL, IREP, IGTYP,          
     .        ICSEN, ISORTHG, IFAILURE, IINT                            
      TYPE(G_BUFEL_)  ,POINTER :: GBUF     
      REAL R4                                                           
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
      NN1 = 1                                                           
      NN2 = 1                                                           
      NN3 = NN2 + NUMELS                                                
C-----------------------------------------------                        
C       face ext                                                        
C-----------------------------------------------                        
      DO 490 NG=1,NGROUP                                                
        CALL INITBUF (IPARG    ,NG      ,                               
     2          MLW     ,NEL     ,NFT     ,IAD     ,ITY     ,           
     3          NPT     ,JALE    ,ISMSTR  ,JEUL    ,JTUR    ,           
     4          JTHE    ,JLAG    ,JMULT   ,JHBE    ,JIVF    ,           
     5          JPOR    ,JCVT    ,JCLOSE  ,JPLASOL ,           
     6          IREP    ,IINT    ,IGTYP   ,ISRAT   ,ISROT   ,           
     7          ICSEN   ,ISORTH  ,ISORTHG ,IFAILURE)                    
        ISOLNOD = IPARG(28,NG)                                          
        DO OFFSET = 0,NEL-1,NVSIZ                                       
          NFT   =IPARG(3,NG) + OFFSET                                   
          IAD   =IPARG(4,NG)                                            
          LFT=1                                                         
          LLT=MIN(NVSIZ,NEL-OFFSET)                                     
          NFT_FA = NFT                                                  
C-----------------------------------------------                        
C       SOLIDE                                                          
C-----------------------------------------------                        
.and.          IF (ITY==1  MLW > 0) THEN                                              
           IALEL=IPARG(7,NG)+IPARG(11,NG)                               
           GBUF => ELBUF_TAB(NG)%GBUF
           DO 130 I=LFT,LLT                                             
             N = I + NFT                                                
             N_FA = I + NFT_FA                                          
             IF(EL2FA(NN2+N_FA)/=0)THEN                               
               IF(IALEL==0)THEN                                       
                 MT=IXS(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                                                    
               IF(ISOLNOD==16)THEN                                    
                 VALUE = FOURTH*VALUE                                    
                 MAS(EL2FA(NN2+N_FA)) = VALUE                           
                 MAS(EL2FA(NN2+N_FA)+1) = VALUE                         
                 MAS(EL2FA(NN2+N_FA)+2) = VALUE                         
                 MAS(EL2FA(NN2+N_FA)+3) = VALUE                         
               ELSE                                                     
                 MAS(EL2FA(NN2+N_FA)) = VALUE                           
               ENDIF                                                    
             ENDIF                                                      
 130       CONTINUE                                                     
C                                                                       
.AND..and.        ELSEIF(ISPH3D==1ITY==51MLW > 0)THEN                           
C-----------------------------------------------                        
C        TETRAS SPH.                                                    
C-----------------------------------------------                        
           GBUF => ELBUF_TAB(NG)%GBUF
           IALEL=IPARG(7,NG)+IPARG(11,NG)                               
           DO 140 I=LFT,LLT                                             
             N = I + NFT                                                
             N_FA = I + NFT_FA                                          
             IF(EL2FA(NN3+N_FA)/=0)THEN                               
               IF(IALEL==0)THEN                                       
                 IPRT=IPARTSP(N)                                        
                 MT  =IPART(1,IPRT)                                     
                 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                             
             ENDIF                                                      
 140       CONTINUE                                                     
C                                                                       
        ELSE                                                            
        ENDIF                                                           
C-----------------------------------------------                        
C       FIN DE BOUCLE SUR LES OFFSET                                    
C-----------------------------------------------                        
       END DO                                                           
 490  CONTINUE                                                          
C-----------------------------------------------                        
C                                                                       
      RETURN