flow_depla.F File Reference

#include "implicit_f.inc"
#include "com08_c.inc"
#include "scr05_c.inc"
#include "vectorize.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	flow_depla (nale, v, d, nodft, nodlt)

Function/Subroutine Documentation

flow_depla()

subroutine flow_depla	(	integer, dimension(*)	nale,
			v,
			d,
		integer	nodft,
		integer	nodlt )

Definition at line 28 of file flow_depla.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      "com08_c.inc"
#include      "scr05_c.inc"
C-----------------------------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NODFT, NODLT, NALE(*)
      my_real d(3,*) ,v(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  :: N
      my_real :: vdt
      DOUBLE PRECISION VDT2
C-------------------------------------------------------------
C         CALCUL DES  DEPLACEMENTS TRANSLATIONS PUIS ROTATIONS
C-------------------------------------------------------------
 
      IF (iresp==1) THEN
#include      "vectorize.inc"
        DO n=nodft,nodlt
          IF(nale(n) == 3) THEN
             vdt2 = dt2*v(1,n)
             d(1,n)=d(1,n)+vdt2
cc             X(1,N)=X(1,N)+VDT2
cC             XDP(1,N)=XDP(1,N)+VDT2
C             X(1,N)=XDP(1,N)
             vdt2 = dt2*v(2,n)
             d(2,n)=d(2,n)+vdt2
cc             X(2,N)=X(2,N)+VDT2            
c             XDP(2,N)=XDP(2,N)+VDT2
c             X(2,N)=XDP(2,N)
             vdt2 = dt2*v(3,n)
             d(3,n)=d(3,n)+vdt2
cc             X(3,N)=X(3,N)+VDT2
C             XDP(3,N)=XDP(3,N)+VDT2
C             X(3,N)=XDP(3,N)
          ENDIF
        ENDDO
      ELSE
#include      "vectorize.inc"
        DO n=nodft,nodlt
          IF(nale(n) == 3) THEN
              vdt = dt2*v(1,n)
              d(1,n)=d(1,n)+vdt
cc              X(1,N)=X(1,N)+VDT
              vdt = dt2*v(2,n)
              d(2,n)=d(2,n)+vdt
cc              X(2,N)=X(2,N)+VDT
              vdt = dt2*v(3,n)
              d(3,n)=d(3,n)+vdt
cc              X(3,N)=X(3,N)+VDT
          ENDIF
        ENDDO
      ENDIF 
C
      RETURN