i9cor3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	i9cor3 (x, v, a, irect, msr, nsv, irtl, ix1, ix2, ix3, ix4, x1, x2, x3, x4, xi, y1, y2, y3, y4, yi, z1, z2, z3, z4, zi, lft, llt, nft)

Function/Subroutine Documentation

i9cor3()

subroutine i9cor3	(		x,
			v,
			a,
		integer, dimension(4,*)	irect,
		integer, dimension(*)	msr,
		integer, dimension(*)	nsv,
		integer, dimension(*)	irtl,
		integer, dimension(mvsiz), intent(out)	ix1,
		integer, dimension(mvsiz), intent(out)	ix2,
		integer, dimension(mvsiz), intent(out)	ix3,
		integer, dimension(mvsiz), intent(out)	ix4,
		intent(out)	x1,
		intent(out)	x2,
		intent(out)	x3,
		intent(out)	x4,
		intent(out)	xi,
		intent(out)	y1,
		intent(out)	y2,
		intent(out)	y3,
		intent(out)	y4,
		intent(out)	yi,
		intent(out)	z1,
		intent(out)	z2,
		intent(out)	z3,
		intent(out)	z4,
		intent(out)	zi,
		integer, intent(inout)	lft,
		integer, intent(inout)	llt,
		integer, intent(inout)	nft )

Definition at line 28 of file i9cor3.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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(INOUT) :: LFT
      INTEGER, INTENT(INOUT) :: LLT
      INTEGER, INTENT(INOUT) :: NFT
      INTEGER IRECT(4,*), MSR(*), NSV(*), IRTL(*)
      INTEGER, DIMENSION(MVSIZ), INTENT(OUT):: IX1,IX2,IX3,IX4
      my_real, DIMENSION(MVSIZ), INTENT(OUT) ::  x1,x2,x3,x4,xi
      my_real, DIMENSION(MVSIZ), INTENT(OUT) ::  y1,y2,y3,y4,yi
      my_real, DIMENSION(MVSIZ), INTENT(OUT) ::  z1,z2,z3,z4,zi
      my_real
     .   x(3,*),v(3,*),a(3,*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com08_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, IL, IG, L, NN
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      DO i=lft,llt
        il=i+nft
        ig=nsv(il)
        xi(i)=x(1,ig) + (v(1,ig) + a(1,ig)*dt12)*dt2
        yi(i)=x(2,ig) + (v(2,ig) + a(2,ig)*dt12)*dt2
        zi(i)=x(3,ig) + (v(3,ig) + a(3,ig)*dt12)*dt2
      ENDDO
      
      DO i=lft,llt
      il=i+nft
      ig=nsv(il)
      l=irtl(il)
      nn=msr(irect(1,l))
      ix1(i)=nn
      x1(i)=x(1,nn) + (v(1,nn) + a(1,nn)*dt12)*dt2
      y1(i)=x(2,nn) + (v(2,nn) + a(2,nn)*dt12)*dt2
      z1(i)=x(3,nn) + (v(3,nn) + a(3,nn)*dt12)*dt2
      nn=msr(irect(2,l))
      ix2(i)=nn
      x2(i)=x(1,nn) + (v(1,nn) + a(1,nn)*dt12)*dt2
      y2(i)=x(2,nn) + (v(2,nn) + a(2,nn)*dt12)*dt2
      z2(i)=x(3,nn) + (v(3,nn) + a(3,nn)*dt12)*dt2
      nn=msr(irect(3,l))
      ix3(i)=nn
      x3(i)=x(1,nn) + (v(1,nn) + a(1,nn)*dt12)*dt2
      y3(i)=x(2,nn) + (v(2,nn) + a(2,nn)*dt12)*dt2
      z3(i)=x(3,nn) + (v(3,nn) + a(3,nn)*dt12)*dt2
      nn=msr(irect(4,l))
      ix4(i)=nn
      x4(i)=x(1,nn) + (v(1,nn) + a(1,nn)*dt12)*dt2
      y4(i)=x(2,nn) + (v(2,nn) + a(2,nn)*dt12)*dt2
      z4(i)=x(3,nn) + (v(3,nn) + a(3,nn)*dt12)*dt2
      ENDDO
C-----------------------------------------------
      RETURN