cderii.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	cderii (px1, px2, py1, py2, x1g, x2g, x3g, x4g, y1g, y2g, y3g, y4g, z1g, z2g, z3g, z4g, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, x2l, x3l, x4l, y2l, y3l, y4l)

Function/Subroutine Documentation

cderii()

subroutine cderii	(	intent(out)	px1,
		intent(out)	px2,
		intent(out)	py1,
		intent(out)	py2,
		intent(in)	x1g,
		intent(in)	x2g,
		intent(in)	x3g,
		intent(in)	x4g,
		intent(in)	y1g,
		intent(in)	y2g,
		intent(in)	y3g,
		intent(in)	y4g,
		intent(in)	z1g,
		intent(in)	z2g,
		intent(in)	z3g,
		intent(in)	z4g,
		intent(in)	e1x,
		intent(in)	e2x,
		intent(in)	e3x,
		intent(in)	e1y,
		intent(in)	e2y,
		intent(in)	e3y,
		intent(in)	e1z,
		intent(in)	e2z,
		intent(in)	e3z,
		intent(out)	x2l,
		intent(out)	x3l,
		intent(out)	x4l,
		intent(out)	y2l,
		intent(out)	y3l,
		intent(out)	y4l )

Definition at line 30 of file cderii.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-----------------------------------------------
      my_real, DIMENSION(MVSIZ), INTENT(IN)  :: 
     .         e1x ,e2x ,e3x ,e1y ,e2y ,e3y ,e1z ,e2z ,e3z ,
     .         x1g,x2g,x3g,x4g,y1g,y2g,y3g,y4g,z1g,z2g,z3g,z4g
c
      my_real, DIMENSION(MVSIZ), INTENT(OUT) :: px1,px2,py1,py2,
     .         x2l,x3l,x4l,y2l,y3l,y4l
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "vect01_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real, DIMENSION(MVSIZ) :: x21g,y21g,z21g,x31g,y31g,z31g,x41g,
     .                             y41g,z41g,z2l,z3l,z4l
C=======================================================================
      DO i=lft,llt
        x21g(i) = x2g(i)-x1g(i)
        y21g(i) = y2g(i)-y1g(i)
        z21g(i) = z2g(i)-z1g(i)
        x31g(i) = x3g(i)-x1g(i)
        y31g(i) = y3g(i)-y1g(i)
        z31g(i) = z3g(i)-z1g(i)
        x41g(i) = x4g(i)-x1g(i)
        y41g(i) = y4g(i)-y1g(i)
        z41g(i) = z4g(i)-z1g(i)
      ENDDO
C
      DO i=lft,llt
        x2l(i)=e1x(i)*x21g(i)+e1y(i)*y21g(i)+e1z(i)*z21g(i)
        y2l(i)=e2x(i)*x21g(i)+e2y(i)*y21g(i)+e2z(i)*z21g(i)
c        Z2L(I)=E3X(I)*X21G(I)+E3Y(I)*Y21G(I)+E3Z(I)*Z21G(I)
        x3l(i)=e1x(i)*x31g(i)+e1y(i)*y31g(i)+e1z(i)*z31g(i)
        y3l(i)=e2x(i)*x31g(i)+e2y(i)*y31g(i)+e2z(i)*z31g(i)
c        Z3L(I)=E3X(I)*X31G(I)+E3Y(I)*Y31G(I)+E3Z(I)*Z31G(I)
        x4l(i)=e1x(i)*x41g(i)+e1y(i)*y41g(i)+e1z(i)*z41g(i)
        y4l(i)=e2x(i)*x41g(i)+e2y(i)*y41g(i)+e2z(i)*z41g(i)
c        Z4L(I)=E3X(I)*X41G(I)+E3Y(I)*Y41G(I)+E3Z(I)*Z41G(I)
      ENDDO
c
      DO i=lft,llt
        px1(i)= half *(y2l(i)-y4l(i))
        px2(i)= half * y3l(i)
        py1(i)=-half *(x2l(i)-x4l(i))
        py2(i)=-half * x3l(i)
      ENDDO
c-----------
      RETURN