flow_accele.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "vectorize.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	flow_accele (nale, ms, a, v, nodft, nodlt)

Function/Subroutine Documentation

flow_accele()

subroutine flow_accele	(	integer, dimension(*)	nale,
			ms,
			a,
			v,
		integer	nodft,
		integer	nodlt )

Definition at line 28 of file flow_accele.F.

C----6---------------------------------------------------------------7---------8
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"
C-----------------------------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NODFT, NODLT,NALE(*)
C     REAL
      my_real
     .   a(3,*) ,v(3,*) ,ms(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  J,K,N
      my_real
     .         rtmp
C
#include      "vectorize.inc"
      DO n=nodft,nodlt
        IF(nale(n) == 3) THEN
         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
        ELSE ! lagrangian or not Law77 node.
           a(1,n) = zero
           a(2,n) = zero
           a(3,n) = zero
        ENDIF  
      ENDDO
C--------------------------------------------
      RETURN