dmasani0.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	dmasani0 (elbuf_tab, iparg, mas, pm, el2fa, ipart, ipartsp)

Function/Subroutine Documentation

dmasani0()

subroutine dmasani0	(	type (elbuf_struct_), dimension(ngroup), target	elbuf_tab,
		integer, dimension(nparg,*)	iparg,
			mas,
			pm,
		integer, dimension(*)	el2fa,
		integer, dimension(lipart1,*)	ipart,
		integer, dimension(*)	ipartsp )

Definition at line 29 of file dmasani0.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      "com01_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,*),EL2FA(*),IPART(LIPART1,*),                      
     .        IPARTSP(*)                                                
      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                                                           
     .   off, VALUE                                                     
      INTEGER I, NG, NEL, N, MT,IALEL, IPRT,NN1                           
      REAL R4                                                           
      TYPE(G_BUFEL_) ,POINTER :: GBUF
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
      nn1 = 1                                                           
C-----------------------------------------------                        
C       face ext                                                        
C-----------------------------------------------                        
      DO 490 ng=1,ngroup                                                
        nel = iparg(2,ng)                                            
        nft = iparg(3,ng)
        ity = iparg(5,ng)
        lft=1                                                           
        llt=nel                                                         
C-----------------------------------------------                        
        IF(ity==51)THEN                                               
C          PARTICULES SPH.                                              
C-----------------------------------------------                        
           gbuf => elbuf_tab(ng)%GBUF                                                      
           ialel=iparg(7,ng)+iparg(11,ng)                               
           DO i=lft,llt                                                 
             n = i + nft                                                
             IF(el2fa(nn1+n)/=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(nn1+n)) = VALUE                                
             ENDIF                                                      
           ENDDO                                                        
        ENDIF                                                           
C-----------------------------------------------                        
C       FIN DE BOUCLE SUR LES GROUPES                                   
C-----------------------------------------------                        
 490  CONTINUE                                                          
C-----------------------------------------------                        
      RETURN