delnumbc.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	delnumbc (iparg, ixq, ixc, ixtg, el2fa, nbf, inum, nelcut, nbpart, idcmax)

Function/Subroutine Documentation

delnumbc()

subroutine delnumbc	(	integer, dimension(nparg,*)	iparg,
		integer, dimension(nixq,*)	ixq,
		integer, dimension(nixc,*)	ixc,
		integer, dimension(nixtg,*)	ixtg,
		integer, dimension(*)	el2fa,
		integer	nbf,
		integer, dimension(*)	inum,
		integer	nelcut,
		integer	nbpart,
		integer	idcmax )

Definition at line 30 of file delnumbc.F.

      use element_mod , only : nixq,nixc,nixtg
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      "param_c.inc"                                             
C-----------------------------------------------                        
C   D u m m y   A r g u m e n t s                                       
C-----------------------------------------------                        
      INTEGER IXC(NIXC,*),IXTG(NIXTG,*),IPARG(NPARG,*),                 
     .        IXQ(NIXQ,*),EL2FA(*),INUM(*),
     .        NBF,NELCUT,NBPART,IDCMAX                              
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
C     REAL                                                    
      my_real                                                           
     .   off                                                            
      INTEGER IE, NG, ITY, LFT, LLT, N, I, J,                      
     .        NEL, IAD, NFT,                          
     .        NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,NN9
 
C-----------------------------------------------                        
      nn1 = 1                                                           
      nn2 = 1                                                           
      nn3 = 1                                                           
      nn4 = nn3 + numelq                                                
      nn5 = nn4 + numelc                                                
      nn6 = nn5 + numeltg                                               
      nn7 = nn6                                                         
      nn8 = nn7                                                         
      nn9 = nn8                                                         
C                                                                       
      DO ng=1,ngroup                                                
        nel   =iparg(2,ng)                                              
        nft   =iparg(3,ng)                                              
        iad   =iparg(4,ng)                                              
        ity   =iparg(5,ng)                                              
        lft=1                                                           
        llt=nel                                                         
C-----------------------------------------------                        
C       QUAD                                                            
C-----------------------------------------------                        
        IF(ity==2)THEN                                                
         DO i=lft,llt                                                   
          n = i + nft                                                   
          inum(el2fa(nn3+n)) = ixq(nixq,n)                              
         ENDDO                                                          
C-----------------------------------------------                        
C       COQUES                                                          
C-----------------------------------------------                        
        ELSEIF(ity==3)THEN                                            
         DO i=lft,llt                                                   
          n = i + nft                                                   
          inum(el2fa(nn4+n)) = ixc(nixc,n)                              
         ENDDO                                                          
C-----------------------------------------------                        
C       3-NODE SHELLS                                                 
C-----------------------------------------------                        
        ELSEIF(ity==7)THEN                                            
         DO i=lft,llt                                                   
          n = i + nft                                                   
          inum(el2fa(nn5+n)) = ixtg(nixtg,n)                            
         ENDDO     
        ELSE                                                            
        ENDIF                                                           
      ENDDO                                                          
                                                     
C-----------------------------------------------                        
      idcmax=0
      DO i=1,nbf
        idcmax=max(idcmax,inum(i))
      ENDDO
      CALL write_i_c(inum,nbf)                                          
      DO i=1,nelcut                                                     
        CALL write_i_c(0,1)                                             
      ENDDO                                                             
C                                                                       
      RETURN