jacob_f0.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	jacob_f0 (det0, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, div0, nel)

Function/Subroutine Documentation

jacob_f0()

subroutine jacob_f0	(		det0,
			px1,
			px2,
			px3,
			px4,
			py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			vx1,
			vx2,
			vx3,
			vx4,
			vx5,
			vx6,
			vx7,
			vx8,
			vy1,
			vy2,
			vy3,
			vy4,
			vy5,
			vy6,
			vy7,
			vy8,
			vz1,
			vz2,
			vz3,
			vz4,
			vz5,
			vz6,
			vz7,
			vz8,
		intent(inout)	div0,
		integer, intent(in)	nel )

Definition at line 32 of file jacob_f0.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(IN) :: NEL
C     REAL
      my_real
     .   vx1(*), vx2(*), vx3(*), vx4(*), vx5(*), vx6(*), vx7(*), vx8(*),
     .   vy1(*), vy2(*), vy3(*), vy4(*), vy5(*), vy6(*), vy7(*), vy8(*),
     .   vz1(*), vz2(*), vz3(*), vz4(*), vz5(*), vz6(*), vz7(*), vz8(*),
     .   px1(*), px2(*), px3(*), px4(*),
     .   py1(*), py2(*), py3(*), py4(*),
     .   pz1(*), pz2(*), pz3(*), pz4(*), det0(*)
      my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: div0
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
C                                                                     12
      my_real
     .   dxx(mvsiz), dxy(mvsiz), dxz(mvsiz),
     .   dyx(mvsiz), dyy(mvsiz), dyz(mvsiz),
     .   dzx(mvsiz), dzy(mvsiz), dzz(mvsiz)
C-----------------------------------------------
       CALL sdefot3(
     1   px1,     px2,     px3,     px4,
     2   py1,     py2,     py3,     py4,
     3   pz1,     pz2,     pz3,     pz4,
     4   vx1,     vx2,     vx3,     vx4,
     5   vx5,     vx6,     vx7,     vx8,
     6   vy1,     vy2,     vy3,     vy4,
     7   vy5,     vy6,     vy7,     vy8,
     8   vz1,     vz2,     vz3,     vz4,
     9   vz5,     vz6,     vz7,     vz8,
     a   dxx,     dxy,     dxz,     dyx,
     b   dyy,     dyz,     dzx,     dzy,
     c   dzz,     nel)
C-----------
       CALL jacob_j33(
     1   det0,    dxx,     dxy,     dxz,
     2   dyx,     dyy,     dyz,     dzx,
     3   dzy,     dzz,     nel)
C       
       div0(1:nel) = dxx(1:nel)+ dyy(1:nel)+ dzz(1:nel)
      RETURN