disp_vel_saved_cload.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	disp_vel_saved_cload (v, d, vr, dr, ib, dpl0cld, vel0cld, nibcld, nconld, iroddl, numnod)

Function/Subroutine Documentation

disp_vel_saved_cload()

subroutine disp_vel_saved_cload	(	dimension(3,numnod), intent(in)	v,
		dimension(3,numnod), intent(in)	d,
		dimension(3,numnod), intent(in)	vr,
		dimension(3,numnod), intent(in)	dr,
		integer, dimension(nibcld,nconld), intent(in)	ib,
		dimension(6,nconld), intent(inout)	dpl0cld,
		dimension(6,nconld), intent(inout)	vel0cld,
		integer, intent(in)	nibcld,
		integer, intent(in)	nconld,
		integer, intent(in)	iroddl,
		integer, intent(in)	numnod )

Definition at line 28 of file disp_vel_saved_cload.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      !USE CONSTANT_MOD
      ! ,ONLY : THIRD,FOURTH
!      USE PINCHTYPE_MOD
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NIBCLD,NCONLD,IRODDL,NUMNOD
      INTEGER, INTENT(IN) :: IB(NIBCLD,NCONLD)
C     REAL
      my_real, INTENT(IN) :: 
     .   d(3,numnod) ,v(3,numnod) ,dr(3,numnod)  ,vr(3,numnod)
      my_real, INTENT(INOUT) :: 
     .   dpl0cld(6,nconld),vel0cld(6,nconld)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  NL,N1,N2,N3,N4,ISK,IDIR,IFUN
C-----------------------------------------------
C     IBCL(1:NIBCLD,1:NUMCLD)              IPRES = IBCL(1:NIBCLD,NUMCCLD+1,NUMCLD+NUMPRES)         
C     1: Node Number                       1st node number of the segment
C     2: NS = 10*Noskew+Idir               2nd node number of the segment
!! Save displacement and velocity in the previous cycle
      !DO NL=1,NCLOAD
      DO nl=1,nconld
        n1 = ib(1,nl)
        n2 = ib(2,nl)
        n3 = ib(3,nl)
        n4 = ib(4,nl)
        ifun = ib(9,nl)
        isk = ib(2,nl)/10
        idir  = ib(2,nl)-10*isk
        IF (n4 == -1) THEN
          IF (ifun > 1) THEN
            IF (idir<=3) THEN
              dpl0cld(idir,nl) = d(idir,n1)
              vel0cld(idir,nl) = v(idir,n1)
            ELSEIF (idir<=6 .AND. iroddl/=0) THEN
              dpl0cld(idir,nl) = dr(idir-3,n1)
              vel0cld(idir,nl) = vr(idir-3,n1)
            ENDIF
          ENDIF
        ENDIF  
      ENDDO
C
      RETURN