s10get_x3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s10get_x3 (x, xdp, dr, numnod, xx, yy, zz, nc, isrot, iresp, nel)

Function/Subroutine Documentation

s10get_x3()

subroutine s10get_x3	(	intent(in)	x,
		double precision, dimension(3,numnod), intent(in)	xdp,
		intent(in)	dr,
		integer, intent(in)	numnod,
		double precision, dimension(mvsiz,10), intent(inout)	xx,
		double precision, dimension(mvsiz,10), intent(inout)	yy,
		double precision, dimension(mvsiz,10), intent(inout)	zz,
		integer, dimension(mvsiz,10), intent(in)	nc,
		integer, intent(in)	isrot,
		integer, intent(in)	iresp,
		integer, intent(in)	nel )

Definition at line 28 of file s10get_x3.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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: ISROT,NUMNOD,IRESP,NEL
      INTEGER, DIMENSION(MVSIZ,10), INTENT(IN) :: NC
     
      my_real, DIMENSION(3,NUMNOD), INTENT(IN) :: x,dr
      DOUBLE PRECISION, DIMENSION(3,NUMNOD), INTENT(IN) ::XDP
      DOUBLE PRECISION, DIMENSION(MVSIZ,10), INTENT(INOUT) ::
     .                                      XX, YY, ZZ
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, IPERM1(10),IPERM2(10),N,N1,N2,NN,IUN,MXT_1
C     REAL
      my_real
     .   dvx,dvy,dvz,dx,dy,dz
      DATA iperm1/0,0,0,0,1,2,3,1,2,3/
      DATA iperm2/0,0,0,0,2,3,1,4,4,4/
C-----------------------------------------------
C----------------------------
C     NODAL COORDINATES     |
C----------------------------
      IF (iresp==1)THEN
        DO n=1,4
           DO i=1,nel
             nn = nc(i,n)
             xx(i,n)=xdp(1,nn)
             yy(i,n)=xdp(2,nn)
             zz(i,n)=xdp(3,nn)
           END DO
        END DO
C        
        IF(isrot==0.OR.isrot==2)THEN
          DO n=5,10
            DO i=1,nel
              nn = nc(i,n)
              IF(nn/=0)THEN
                 xx(i,n)=xdp(1,nn)
                 yy(i,n)=xdp(2,nn)
                 zz(i,n)=xdp(3,nn)      
              ELSE
                n1=iperm1(n)
                n2=iperm2(n)
                xx(i,n) = half*(xx(i,n1)+xx(i,n2))
                yy(i,n) = half*(yy(i,n1)+yy(i,n2))
                zz(i,n) = half*(zz(i,n1)+zz(i,n2))
              END IF
            END DO
          END DO
        ELSEIF(isrot==1)THEN
          DO n=5,10
            DO i=1,nel
              nn = nc(i,n)
              n1=iperm1(n)
              n2=iperm2(n)
              dx = (yy(i,n2)-yy(i,n1))*(dr(3,nc(i,n2))-dr(3,nc(i,n1)))
     .           - (zz(i,n2)-zz(i,n1))*(dr(2,nc(i,n2))-dr(2,nc(i,n1)))
              dy = (zz(i,n2)-zz(i,n1))*(dr(1,nc(i,n2))-dr(1,nc(i,n1)))
     .           - (xx(i,n2)-xx(i,n1))*(dr(3,nc(i,n2))-dr(3,nc(i,n1)))
              dz = (xx(i,n2)-xx(i,n1))*(dr(2,nc(i,n2))-dr(2,nc(i,n1)))
     .           - (yy(i,n2)-yy(i,n1))*(dr(1,nc(i,n2))-dr(1,nc(i,n1)))
 
              xx(i,n) = half*(xx(i,n1)+xx(i,n2)) + one_over_8 * dx
              yy(i,n) = half*(yy(i,n1)+yy(i,n2)) + one_over_8 * dy
              zz(i,n) = half*(zz(i,n1)+zz(i,n2)) + one_over_8 * dz
            END DO
          END DO
        END IF
      ELSE !DP
        DO n=1,4
           DO i=1,nel
             nn = nc(i,n)
             xx(i,n)=x(1,nn)
             yy(i,n)=x(2,nn)
             zz(i,n)=x(3,nn)
           END DO
        END DO
        IF(isrot==0.OR.isrot==2)THEN
          DO n=5,10
            DO i=1,nel
              nn = nc(i,n)
              IF(nn/=0)THEN
                 xx(i,n)=x(1,nn)
                 yy(i,n)=x(2,nn)
                 zz(i,n)=x(3,nn)      
              ELSE
                n1=iperm1(n)
                n2=iperm2(n)
                xx(i,n) = half*(xx(i,n1)+xx(i,n2))
                yy(i,n) = half*(yy(i,n1)+yy(i,n2))
                zz(i,n) = half*(zz(i,n1)+zz(i,n2))
              END IF
            END DO
          END DO
        ELSEIF(isrot==1)THEN
          DO n=5,10
            DO i=1,nel
              nn = nc(i,n)
              n1=iperm1(n)
              n2=iperm2(n)
              dx = (yy(i,n2)-yy(i,n1))*(dr(3,nc(i,n2))-dr(3,nc(i,n1)))
     .           - (zz(i,n2)-zz(i,n1))*(dr(2,nc(i,n2))-dr(2,nc(i,n1)))
              dy = (zz(i,n2)-zz(i,n1))*(dr(1,nc(i,n2))-dr(1,nc(i,n1)))
     .           - (xx(i,n2)-xx(i,n1))*(dr(3,nc(i,n2))-dr(3,nc(i,n1)))
              dz = (xx(i,n2)-xx(i,n1))*(dr(2,nc(i,n2))-dr(2,nc(i,n1)))
     .           - (yy(i,n2)-yy(i,n1))*(dr(1,nc(i,n2))-dr(1,nc(i,n1)))
 
              xx(i,n) = half*(xx(i,n1)+xx(i,n2)) + one_over_8 * dx
              yy(i,n) = half*(yy(i,n1)+yy(i,n2)) + one_over_8 * dy
              zz(i,n) = half*(zz(i,n1)+zz(i,n2)) + one_over_8 * dz
            END DO
          END DO
        END IF
      END IF ! (IRESP==1)THEN
C-----------
      RETURN