delnumbs.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	delnumbs (iparg, ixs, el2fa, nbf, inum, kxsp, isph3d)

Function/Subroutine Documentation

delnumbs()

subroutine delnumbs	(	integer, dimension(nparg,*)	iparg,
		integer, dimension(nixs,*)	ixs,
		integer, dimension(*)	el2fa,
		integer	nbf,
		integer, dimension(*)	inum,
		integer, dimension(nisp,*)	kxsp,
		integer	isph3d )

Definition at line 30 of file delnumbs.F.

      use element_mod , only : nixs
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-----------------------------------------------                        
      INTEGER IPARG(NPARG,*),                                           
     .        IXS(NIXS,*),EL2FA(*),INUM(*),                             
     .        NBF,KXSP(NISP,*)                                    
      INTEGER ISPH3D ,NG, ITY, LFT, LLT, N, I, NEL, IAD, NFT,                
     .        NN1,NN2,NN3,ISOLNOD         
C=======================================================================                       
      nn1 = 1                                                           
      nn2 = 1                                                           
      nn3 = nn2 + numels                                                
C                                                                       
      DO 490 ng=1,ngroup                                                
        nel   =iparg(2,ng)                                              
        nft   =iparg(3,ng)                                              
        iad   =iparg(4,ng)                                              
        ity   =iparg(5,ng)                                              
        isolnod = iparg(28,ng)                                          
        lft=1                                                           
        llt=nel                                                         
C-----------------------------------------------                        
C       16-NODE SOLID ELEMENTS                                                     
C-----------------------------------------------                        
        IF(ity==1.AND.isolnod==16)THEN                              
          DO i=lft,llt                                                  
           n = i + nft                                                  
           IF(el2fa(nn2+n)/=0)THEN                                    
               inum(el2fa(nn2+n)) = ixs(nixs,n)                         
               inum(el2fa(nn2+n)+1) = ixs(nixs,n)                       
               inum(el2fa(nn2+n)+2) = ixs(nixs,n)                       
               inum(el2fa(nn2+n)+3) = ixs(nixs,n)                       
           ENDIF                                                        
          ENDDO                                                         
C-----------------------------------------------                        
C       OTHER SOLID ELEMENTS                                                  
C-----------------------------------------------                        
        ELSEIF(ity==1)THEN                                            
          DO i=lft,llt                                                  
           n = i + nft                                                  
           IF(el2fa(nn2+n)/=0)THEN                                    
               inum(el2fa(nn2+n)) = ixs(nixs,n)                         
           ENDIF                                                        
          ENDDO                                                         
        ELSEIF(isph3d==1.AND.ity==51)THEN                           
C-----------------------------------------------                        
C         TETRAS SPH.                                                  
C-----------------------------------------------                        
          DO i=lft,llt                                                  
           n = i + nft                                                  
           IF(el2fa(nn3+n)/=0)THEN                                    
               inum(el2fa(nn3+n)) = kxsp(nisp,n)                        
           ENDIF                                                        
          ENDDO                                                         
        ELSE                                                            
        ENDIF                                                           
 490  CONTINUE                                                          
C-----------------------------------------------                        
      CALL write_i_c(inum,nbf)                                          
C-----------------------------------------------                        
      RETURN