anioffs.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	anioffs (elbuf_tab, iparg, ioff, el2fa, nbf, nbpart, isph3d)

Function/Subroutine Documentation

anioffs()

subroutine anioffs	(	type (elbuf_struct_), dimension(ngroup), target	elbuf_tab,
		integer, dimension(nparg,*)	iparg,
		integer, dimension(*)	ioff,
		integer, dimension(*)	el2fa,
		integer	nbf,
		integer	nbpart,
		integer	isph3d )

Definition at line 30 of file anioffs.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      "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-----------------------------------------------                        
C     REAL                                                    
      INTEGER IPARG(NPARG,*),EL2FA(*),NBF,IOFF(*),                      
     .        NBPART,ISPH3D                                                    
      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-----------------------------------------------                        
C     REAL                                                    
      INTEGER I, NG, NEL, NFT, ITY, LFT, NPT,                      
     .        N, LLT, MLW,ISOLNOD,NN1,NN2,NN3 
      INTEGER IOFFBUF(NBF),RBUF                                         
      TYPE(G_BUFEL_)  ,POINTER :: GBUF     
C=======================================================================
      nn1 = 1                                                           
      nn2 = 1                                                           
      nn3 = nn2 + numels                                                
C-----------------------------------------------                        
C                                                                       
      DO 490 ng=1,ngroup                                                
        mlw   =iparg(1,ng)                                              
        nel   =iparg(2,ng)                                              
        ity   =iparg(5,ng)                                              
        nft   =iparg(3,ng)                                              
        isolnod = iparg(28,ng)                                          
        lft=1                                                           
        llt=nel                                                         
C-----------------------------------------------                        
C       SOLIDES 16N                                                     
C-----------------------------------------------                        
        IF (ity==1.AND.isolnod==16) THEN  
          gbuf =>  elbuf_tab(ng)%GBUF                           
          IF (mlw == 0 .OR. mlw == 13)THEN                                              
            DO i=lft,llt                                                
              n = nft + i                                               
              ioff(el2fa(nn2+n)) = 0                                    
              ioff(el2fa(nn2+n)+1) = 0                                  
              ioff(el2fa(nn2+n)+2) = 0                                  
              ioff(el2fa(nn2+n)+3) = 0                                  
            ENDDO                                                       
          ELSE                                                          
            DO i=lft,llt                                                
              n = nft + i                                               
              ioff(el2fa(nn2+n))   = nint(min(gbuf%OFF(i),one))
              ioff(el2fa(nn2+n)+1) = nint(min(gbuf%OFF(i),one))
              ioff(el2fa(nn2+n)+2) = nint(min(gbuf%OFF(i),one))
              ioff(el2fa(nn2+n)+3) = nint(min(gbuf%OFF(i),one))
            ENDDO                                                       
          ENDIF                                                         
C-----------------------------------------------                        
C       AUTRES SOLIDES                                                  
C-----------------------------------------------                        
        ELSEIF(ity==1)THEN                                            
          gbuf =>  elbuf_tab(ng)%GBUF                           
          IF(mlw==0)THEN                                              
            DO i=lft,llt                                                
              n = nft + i                                               
              ioff(el2fa(nn2+n)) = 0                                    
            ENDDO                                                       
          ELSE                                                          
            DO i=lft,llt                                                
              n = nft + i                                               
              ioff(el2fa(nn2+n)) = nint(min(gbuf%OFF(i),one))        
            ENDDO                                                       
          ENDIF                                                         
        ELSEIF(isph3d==1.AND.ity==51)THEN                           
          gbuf =>  elbuf_tab(ng)%GBUF                           
C-----------------------------------------------                        
C         TETRAS SPH.                                                  
C-----------------------------------------------                        
          IF(mlw==0)THEN                                              
            DO i=lft,llt                                                
              n = nft + i                                               
              ioff(el2fa(nn3+n)) = 0                                    
            ENDDO                                                       
          ELSE                                                          
            DO i=lft,llt                                                
              n = nft + i                                               
              ioff(el2fa(nn3+n)) = nint(min(gbuf%OFF(i),one))      
            ENDDO                                                       
          ENDIF                                                         
        ELSE                                                            
        ENDIF                                                           
C-----------------------------------------------                        
 490  CONTINUE                                                          
C-----------------------------------------------                        
      CALL write_c_c(ioff,nbf)                                          
C-----------                                                                       
      RETURN