c3curv3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	c3curv3 (jft, jlt, vrl1, vrl2, vrl3, ixtg, wxy, ismstr, kxx, kyy, kxy, px1, py1, py2, eyz, ezx, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z)

Function/Subroutine Documentation

c3curv3()

subroutine c3curv3	(	integer	jft,
		integer	jlt,
			vrl1,
			vrl2,
			vrl3,
		integer, dimension(nixtg,*)	ixtg,
			wxy,
		integer	ismstr,
			kxx,
			kyy,
			kxy,
			px1,
			py1,
			py2,
			eyz,
			ezx,
			e1x,
			e1y,
			e1z,
			e2x,
			e2y,
			e2z,
			e3x,
			e3y,
			e3z )

Definition at line 31 of file c3curv3.F.

      use element_mod , only : nixtg 
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 JFT, JLT,IXTG(NIXTG,*),ISMSTR
      my_real vrl1(mvsiz,3), vrl2(mvsiz,3), vrl3(mvsiz,3),wxy(*),
     .  e1x(mvsiz),  e1y(mvsiz),  e1z(mvsiz),  e2x(mvsiz), 
     .  e2y(mvsiz),  e2z(mvsiz),  e3x(mvsiz),  e3y(mvsiz), 
     .  e3z(mvsiz),  ezx(mvsiz),  eyz(mvsiz),               
     .  px1(mvsiz),  py1(mvsiz),  py2(mvsiz),
     .  kxx(mvsiz),  kyy(mvsiz),  kxy(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real rx1(mvsiz),  rx2(mvsiz),  rx3(mvsiz),
     .        ry1(mvsiz),  ry2(mvsiz),  ry3(mvsiz)
      my_real rx12t,rx13t,rx23t,ry12t,ry13t,ry23t,rxavt,ryavt
C=======================================================================
      DO i=jft,jlt
        rx1(i)=e1x(i)*vrl1(i,1)+e1y(i)*vrl1(i,2)+e1z(i)*vrl1(i,3)
        ry1(i)=e2x(i)*vrl1(i,1)+e2y(i)*vrl1(i,2)+e2z(i)*vrl1(i,3)
        ry2(i)=e2x(i)*vrl2(i,1)+e2y(i)*vrl2(i,2)+e2z(i)*vrl2(i,3)
        rx2(i)=e1x(i)*vrl2(i,1)+e1y(i)*vrl2(i,2)+e1z(i)*vrl2(i,3)
        rx3(i)=e1x(i)*vrl3(i,1)+e1y(i)*vrl3(i,2)+e1z(i)*vrl3(i,3)
        ry3(i)=e2x(i)*vrl3(i,1)+e2y(i)*vrl3(i,2)+e2z(i)*vrl3(i,3)
      ENDDO
C
      DO i=jft,jlt
        rx12t  = rx1(i) - rx2(i)
        rx13t  = rx1(i) - rx3(i)
        rx23t  = rx2(i) - rx3(i)
C
        kyy(i)=-py1(i) * rx13t - py2(i) * rx23t
        kxy(i)= px1(i) * rx12t
C
        ry12t = ry1(i) - ry2(i)
        ry13t = ry1(i) - ry3(i)
        ry23t = ry2(i) - ry3(i)
C
        kxx(i) = px1(i) * ry12t
        kxy(i) = py1(i) * ry13t + py2(i) * ry23t - kxy(i)
C
        ryavt  = px1(i) * (
     .              px1(i) * (-rx1(i)+rx2(i)) 
     .            + (two*py1(i)+three*py2(i)) * ry1(i)
     .            + (three*py1(i)+two*py2(i)) * ry2(i)
 
     .            + ( py1(i) + py2(i) )   * ry3(i) )
        rxavt  = -px1(i) * (   
     .           + (two*py1(i)+py2(i)) * rx1(i)
     .           + (py1(i)+two*py2(i)) * rx2(i)
     .           + three*(py1(i)+py2(i)) * rx3(i) )
     .           + py1(i) * (py1(i)+two*py2(i)) * ry1(i)
     .           - py2(i) * (two*py1(i)+py2(i)) * ry2(i)
     .           + ( py2(i)**2 - py1(i)**2 )   * ry3(i)
C
        ezx(i) = ezx(i) + ryavt * third
        eyz(i) = eyz(i) + rxavt * third   
      ENDDO
c
      IF (ismstr==10) THEN
        DO i=jft,jlt
         rx12t = rx1(i) - rx2(i)
         ry13t = ry1(i) - ry3(i)
         ry23t = ry2(i) - ry3(i)
C
         wxy(i) = py1(i) * ry13t + py2(i) * ry23t + px1(i) * rx12t
        ENDDO 
      ENDIF
c-----------
      RETURN