accele.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	accele (a, ar, v, ms, in, size_nale, nale, ms_2d, size_npby, npby)

Function/Subroutine Documentation

accele()

subroutine accele	(	intent(inout)	a,
		intent(inout)	ar,
		intent(inout)	v,
		intent(inout)	ms,
		intent(inout)	in,
		integer, intent(in)	size_nale,
		integer, dimension(size_nale), intent(in)	nale,
		intent(inout)	ms_2d,
		integer, intent(in)	size_npby,
		integer, dimension(nnpby,size_npby), intent(in)	npby )

Parameters

[in]	size_nale	size of NALE array
[in]	size_npby	2nd dimension of NPBY array
[in]	nale	array of ALE node (1 if the node is related to ALE computation
[in]	npby	rigid body array

Definition at line 28 of file accele.F.

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      "inter22.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER, INTENT(IN) :: SIZE_NALE !< size of NALE array
        INTEGER, INTENT(IN) :: SIZE_NPBY !< 2nd dimension of NPBY array
        INTEGER, DIMENSION(SIZE_NALE), INTENT(IN) :: NALE !< array of ALE node (1 if the node is related to ALE computation
        INTEGER, DIMENSION(NNPBY,SIZE_NPBY), INTENT(IN) :: NPBY !< rigid body array
        my_real, DIMENSION(3,NUMNOD), INTENT(INOUT) :: a,v,ar !< acceleration, velocity and  angular acceleration
        my_real, DIMENSION(NUMNOD), INTENT(INOUT) :: ms,in,ms_2d !< mass and inertia
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: N,I,MSR
      my_real :: rtmp
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------     
      
      !------------------------------------------------!
      !       acceleration computation 
      !------------------------------------------------!       
        IF(n2d == 0) THEN
!$OMP DO SCHEDULE(guided)
            DO n=1,numnod
                IF(ms(n)>zero) THEN
                    rtmp = one / ms(n)
                    a(1,n) = a(1,n) * rtmp
                    a(2,n) = a(2,n) * rtmp
                    a(3,n) = a(3,n) * rtmp
                ELSE
                    a(1,n) = zero
                    a(2,n) = zero
                    a(3,n) = zero
                    IF(nmult/=0) THEN
                        v(1,n) = zero
                        v(2,n) = zero
                        v(3,n) = zero
                    ENDIF
                ENDIF
            ENDDO
!$OMP END DO
        ELSE
            DO i=1,nrbody        
                msr = npby(1,i)
                IF(msr>0) THEN
                    ms(msr) = ms_2d(msr)
                ENDIF
            ENDDO
!$OMP DO SCHEDULE(guided)
            DO n=1,numnod
                IF(ms(n)>zero) THEN
                    rtmp = one / ms(n)
                    a(1,n) = a(1,n) * rtmp
                    a(2,n) = a(2,n) * rtmp
                    a(3,n) = a(3,n) * rtmp
                ELSE
                    a(1,n) = zero
                    a(2,n) = zero
                    a(3,n) = zero
                    IF(nmult/=0) THEN
                        v(1,n) = zero
                        v(2,n) = zero
                        v(3,n) = zero
                    ENDIF
                ENDIF
            ENDDO
!$OMP END DO
        ENDIF
 
      !------------------------------------------------!
      !       angular acceleration computation 
      !------------------------------------------------!
        IF (iroddl/=0) THEN
!$OMP DO SCHEDULE(guided)
            DO n=1,numnod
                IF(in(n)>zero) THEN
                    rtmp = one / in(n)
                    ar(1,n) = ar(1,n) * rtmp
                    ar(2,n) = ar(2,n) * rtmp
                    ar(3,n) = ar(3,n) * rtmp
                ELSE
                    ar(1,n) = zero
                    ar(2,n) = zero
                    ar(3,n) = zero      
                ENDIF
            ENDDO
!$OMP END DO
        ENDIF
 
      !------------------------------------------------!
      !      FVM DOES NOT NEED ACCEL AND VEL           !
      !------------------------------------------------!      
        IF(int22>0)THEN
!$OMP DO SCHEDULE(guided)
            DO n=1,numnod
                IF(nale(n)/=0)THEN
                    a( 1:3,n) = zero
                    ar(1:3,n) = zero
                    ms(n)     = zero
                ENDIF
            ENDDO  
!$OMP END DO   
        ENDIF      
      !------------------------------------------------!    
 
        RETURN