i17xsave.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "com04_c.inc"
#include "com08_c.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i17xsave (x, v, a, ixs, ixs16, nmes, nme_t, esh_t, nelem, neles, eminxm, itask, xsav, xmin, ymin, zmin, xmax, ymax, zmax)

Function/Subroutine Documentation

i17xsave()

subroutine i17xsave	(		x,
			v,
			a,
		integer, dimension(nixs,*)	ixs,
		integer, dimension(8,*)	ixs16,
		integer	nmes,
		integer	nme_t,
		integer	esh_t,
		integer, dimension(*)	nelem,
		integer, dimension(*)	neles,
			eminxm,
		integer	itask,
			xsav,
			xmin,
			ymin,
			zmin,
			xmax,
			ymax,
			zmax )

Definition at line 31 of file i17xsave.F.

      use element_mod , only : nixs
C saving of xsav and calculation of domain bounds
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "com08_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NMES,NME_T,ESH_T,ITASK,
     .        NELEM(*),NELES(*),IXS(NIXS,*),IXS16(8,*)
      my_real
     .        xmax, ymax, zmax, xmin, ymin, zmin,
     .        x(3,*), v(3,*), a(3,*), eminxm(6,*), xsav(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NMES_F, NMES_L, I, J, K
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
C
      nmes_f = 1 + itask*nmes / nthread
      nmes_l = (itask+1)*nmes / nthread
C
      xmin=ep30
      xmax=-ep30
      ymin=ep30
      ymax=-ep30
      zmin=ep30
      zmax=-ep30
C
      DO k=1,8
        DO i=1+esh_t,nme_t+esh_t                      
          j=ixs(k+1,nelem(i))
          xsav(1,j) = x(1,j)+dt2*(v(1,j)+dt12*a(1,j))
          xsav(2,j) = x(2,j)+dt2*(v(2,j)+dt12*a(2,j))
          xsav(3,j) = x(3,j)+dt2*(v(3,j)+dt12*a(3,j))
          j=ixs16(k,nelem(i)-numels8-numels10-numels20)
          xsav(1,j) = x(1,j)+dt2*(v(1,j)+dt12*a(1,j))
          xsav(2,j) = x(2,j)+dt2*(v(2,j)+dt12*a(2,j))
          xsav(3,j) = x(3,j)+dt2*(v(3,j)+dt12*a(3,j))
C
          xmin = min( xmin , eminxm(1,i) )
          ymin = min( ymin , eminxm(2,i) )
          zmin = min( zmin , eminxm(3,i) )
          xmax = max( xmax , eminxm(4,i) )
          ymax = max( ymax , eminxm(5,i) )
          zmax = max( zmax , eminxm(6,i) )
        ENDDO
        DO i=nmes_f,nmes_l
          j=ixs(k+1,neles(i))
          xsav(1,j) = x(1,j)+dt2*(v(1,j)+dt12*a(1,j))
          xsav(2,j) = x(2,j)+dt2*(v(2,j)+dt12*a(2,j))
          xsav(3,j) = x(3,j)+dt2*(v(3,j)+dt12*a(3,j))
          j=ixs16(k,neles(i)-numels8-numels10-numels20)
          xsav(1,j) = x(1,j)+dt2*(v(1,j)+dt12*a(1,j))
          xsav(2,j) = x(2,j)+dt2*(v(2,j)+dt12*a(2,j))
          xsav(3,j) = x(3,j)+dt2*(v(3,j)+dt12*a(3,j))
        ENDDO
      ENDDO
C
      RETURN