euldx.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "com08_c.inc"
#include "scr05_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	euldx (v, d, x, ddp, nale, nodft, nodlt)

Function/Subroutine Documentation

euldx()

subroutine euldx	(		v,
			d,
			x,
		real(kind=8), dimension(3,*), intent(inout)	ddp,
		integer, dimension(numnod)	nale,
		integer	nodft,
		integer	nodlt )

Definition at line 28 of file euldx.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      "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 NALE(NUMNOD),NODFT,NODLT
      my_real x(3,numnod), d(3,numnod), v(3,numnod)  
      REAL(kind=8), dimension(3,*), INTENT(INOUT) :: ddp
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  N
      my_real xlag
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
        IF(n2d == 1)THEN
          DO n=nodft,nodlt
            IF(x(2,n)+dt2*v(2,n) < zero)v(2,n)=zero
          ENDDO
        ENDIF
        DO n=nodft,nodlt
          xlag=min(1,iabs(iabs(nale(n))-2))
          d(1,n)=d(1,n)+dt2*v(1,n)*xlag
          d(2,n)=d(2,n)+dt2*v(2,n)*xlag
          d(3,n)=d(3,n)+dt2*v(3,n)*xlag
          x(1,n)=x(1,n)+dt2*v(1,n)*xlag
          x(2,n)=x(2,n)+dt2*v(2,n)*xlag
          x(3,n)=x(3,n)+dt2*v(3,n)*xlag
        ENDDO
!   -------------------------
!   single precision
        IF(iresp==1) THEN
            DO n=nodft,nodlt
                ddp(1,n)= d(1,n)
                ddp(2,n)= d(2,n)
                ddp(3,n)= d(3,n)
            ENDDO
        ENDIF    
C-----------------------------------------------       
      RETURN