dmasani0.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    dmasani0       ../starter/source/output/anim/dmasani0.F
!||--- called by ------------------------------------------------------
!||    genani1        ../starter/source/output/anim/genani1.F
!||--- uses       -----------------------------------------------------
!||====================================================================
      SUBROUTINE dmasani0(ELBUF_TAB   ,IPARG   ,       
     2              MAS     ,PM      ,EL2FA   ,IPART   ,IPARTSP )       
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                                                            
      END