sgcoor10.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "scr18_c.inc"
#include "scr05_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sgcoor10 (xx, yy, zz, x, xdp, x0, y0, z0, vx0, vy0, vz0, sav, nc, d, off, off0, nel, mtn, ismstr)

Function/Subroutine Documentation

sgcoor10()

subroutine sgcoor10	(	double precision, dimension(mvsiz,10)	xx,
		double precision, dimension(mvsiz,10)	yy,
		double precision, dimension(mvsiz,10)	zz,
			x,
		double precision, dimension(3,*)	xdp,
		double precision, dimension(mvsiz,10)	x0,
		double precision, dimension(mvsiz,10)	y0,
		double precision, dimension(mvsiz,10)	z0,
			vx0,
			vy0,
			vz0,
		double precision, dimension(nel,30)	sav,
		integer, dimension(mvsiz,10)	nc,
			d,
			off,
			off0,
		integer	nel,
		integer, intent(in)	mtn,
		integer, intent(in)	ismstr )

Definition at line 28 of file sgcoor10.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr18_c.inc"
#include      "scr05_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: MTN
      INTEGER, INTENT(IN) :: ISMSTR
      INTEGER NC(MVSIZ,10),NEL
C     REAL
      my_real
     .  vx0(mvsiz,10),vy0(mvsiz,10),vz0(mvsiz,10),
     .  d(3,*) ,off(*),off0(*),x(3,*)
      DOUBLE PRECISION 
     .   XX(MVSIZ,10), YY(MVSIZ,10), ZZ(MVSIZ,10),XDP(3,*),
     .   X0(MVSIZ,10), Y0(MVSIZ,10), Z0(MVSIZ,10),SAV(NEL,30)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,NPE,N,N2,N3
C     REAL
C-----------------------------------------------
C----------------------------
C     NODAL COORDINATES INITIALES    |
C----------------------------
      npe = 10
      DO n=1,npe
       DO i=1,nel
        x0(i,n)=sav(i,n)
        y0(i,n)=sav(i,n+npe)
        z0(i,n)=sav(i,n+2*npe)
       ENDDO
      ENDDO 
C----------------------------
C     DISPLACEMENT     |
C----------------------------
        n = npe
        DO i=1,nel
        vx0(i,n)=xx(i,n)-x0(i,n)
        vy0(i,n)=yy(i,n)-y0(i,n)
        vz0(i,n)=zz(i,n)-z0(i,n)
        ENDDO
       DO n=1,npe
        DO i=1,nel
        vx0(i,n)=xx(i,n)-x0(i,n)-vx0(i,npe)
        vy0(i,n)=yy(i,n)-y0(i,n)-vy0(i,npe)
        vz0(i,n)=zz(i,n)-z0(i,n)-vz0(i,npe)
        ENDDO
       ENDDO 
C---   w/ /DT/BRICK/CST    
      IF (ismstr == 12.AND.idtmin(1)==3) THEN
        DO i=1,nel
         off0(i)=off(i)
        ENDDO
       IF (mtn==1) THEN
         IF(iresp==1)THEN
          DO i=1,nel
            IF (off(i) <= one ) cycle 
            n3 = nc(i,npe)
            vx0(i,npe)=xdp(1,n3)-x0(i,npe)
            vy0(i,npe)=xdp(2,n3)-y0(i,npe)
            vz0(i,npe)=xdp(3,n3)-z0(i,npe)
          ENDDO 
          DO n=1,npe
           DO i=1,nel
            IF (off(i) <= one ) cycle 
            n2 = nc(i,n)
            vx0(i,n) =xdp(1,n2)-x0(i,n)-vx0(i,npe)
            vy0(i,n) =xdp(2,n2)-y0(i,n)-vy0(i,npe)
            vz0(i,n) =xdp(3,n2)-z0(i,n)-vz0(i,npe)
           ENDDO
          ENDDO 
         ELSE      
          DO i=1,nel
            IF (off(i) <= one ) cycle 
            n3 = nc(i,npe)
            vx0(i,npe)=x(1,n3)-x0(i,npe)
            vy0(i,npe)=x(2,n3)-y0(i,npe)
            vz0(i,npe)=x(3,n3)-z0(i,npe)
          ENDDO 
          DO n=1,npe
           DO i=1,nel
            IF (off(i) <= one ) cycle 
            n2 = nc(i,n)
            vx0(i,n) =x(1,n2)-x0(i,n)-vx0(i,npe)
            vy0(i,n) =x(2,n2)-y0(i,n)-vy0(i,npe)
            vz0(i,n) =x(3,n2)-z0(i,n)-vz0(i,npe)
           ENDDO
          ENDDO 
         END IF !(IRESP==1)THEN
       ELSE
         DO n=1,npe
          DO i=1,nel
           IF (off(i) <= one ) cycle 
            vx0(i,n)=d(1,nc(i,n))-d(1,nc(i,npe))  
            vy0(i,n)=d(2,nc(i,n))-d(2,nc(i,npe))
            vz0(i,n)=d(3,nc(i,n))-d(3,nc(i,npe))
          ENDDO
         ENDDO
       END IF !(MTN==1) THEN
      END IF !(ISMSTR == 10.AND.IDTMIN(1)==3) THEN
C
      RETURN