parsor0.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/.
!||====================================================================
!||    parsor0     ../starter/source/output/anim/parsor0.F
!||--- called by ------------------------------------------------------
!||    genani1     ../starter/source/output/anim/genani1.F
!||--- calls      -----------------------------------------------------
!||====================================================================
      SUBROUTINE parsor0(IADD ,IPARG   ,MATER   ,EL2FA   ,              
     3                KXSP    ,IPARTSP )                                
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      "sphcom.inc"                                              
#include      "param_c.inc"                                             
C-----------------------------------------------                        
C   D u m m y   A r g u m e n t s                                       
C-----------------------------------------------                        
C     REAL                                                    
      INTEGER IADD(*),IPARG(NPARG,*),                                   
     .        mater(*),el2fa(*),                                        
     .        kxsp(nisp,*),ipartsp(*)                                   
C-----------------------------------------------                        
C     REAL                                                    
      my_real                                                           
     .   off                                                            
      INTEGER II, IE, NG, ITY, LFT, LLT, KPT, N, I,                     
     .        nel, iad, npar, nft, iprt,                                
     .        nn1,nn2,nn3,nn4,nn5,nn6,nn7,nn8,nn9,nn10,                 
     .        jj, inod,                                                 
     .        np(numsph),buf,bufp                                       
C-----------------------------------------------                        
C                                                                       
      nn1 = 1                                                           
      nn2 = nn1 + numsph+maxpjet                                        
      ie = 0                                                            
C-----------------------------------------------                        
C     MID                                                               
C-----------------------------------------------                        
      npar = 0                                                          
      jj = 0                                                            
C                                                                       
      DO 100 iprt = 1 , npart                                           
C                                                                       
       IF(mater(iprt)/=4) GOTO 100                                    
       npar = npar + 1                                                  
       DO 90 ng=1,ngroup                                                
        nel   =iparg(2,ng)                                              
        nft   =iparg(3,ng)                                              
        iad   =iparg(4,ng)                                              
        ity   =iparg(5,ng)                                              
        lft=1                                                           
        llt=nel                                                         
        IF(ity==51)THEN                                               
C-----------------------------------------------                        
C         PARTICULES SPH.                                               
C-----------------------------------------------                        
          DO 20 i=lft,llt                                               
           n = i + nft                                                  
           IF(ipartsp(n)/=iprt) GOTO 20                               
             inod=kxsp(3,n)                                             
             ii = inod-1                                                
             CALL write_i_c(ii,1)                                       
             jj = jj+1                                                  
             ie = ie + 1                                                
             el2fa(nn1+n) = ie                                          
 20      CONTINUE                                                       
        ELSE                                                            
        ENDIF                                                           
 90    CONTINUE                                                         
C-----------------------------------------------                        
C       PART ADRESS                                                     
C-----------------------------------------------                        
       iadd(npar) = ie                                                  
 100  CONTINUE                                                          
C                                                                       
C-----------------------------------------------                        
      RETURN                                                            
      END