s10deri3.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	s10deri3 (vol, ngl, xx, yy, zz, px, py, pz, nx, rx, ry, rz, sx, sy, sz, tx, ty, tz, volu, voln, elbuf_str, volg)
subroutine	s10jacob (alph, beta, w, x1b, x2b, x3b, x4a, x5b, x6b, x7b, x8b, x9b, x10b, x8a, x9a, x10a, y1b, y2b, y3b, y4a, y5b, y6b, y7b, y8b, y9b, y10b, y8a, y9a, y10a, z1b, z2b, z3b, z4a, z5b, z6b, z7b, z8b, z9b, z10b, z8a, z9a, z10a, px1, px2, px3, px4, px5, px6, px7, px8, px9, px10, py1, py2, py3, py4, py5, py6, py7, py8, py9, py10, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pz9, pz10, nx1, nx2, nx3, nx4, nx5, nx6, nx7, nx8, nx9, nx10, vol, voldp)
subroutine	s10pijsav (px, py, pz, pij, nel)

Function/Subroutine Documentation

s10deri3()

subroutine s10deri3	(		vol,
		integer, dimension(*)	ngl,
		double precision, dimension(mvsiz,10)	xx,
		double precision, dimension(mvsiz,10)	yy,
		double precision, dimension(mvsiz,10)	zz,
			px,
			py,
			pz,
			nx,
			rx,
			ry,
			rz,
			sx,
			sy,
			sz,
			tx,
			ty,
			tz,
			volu,
			voln,
		type(elbuf_struct_), target	elbuf_str,
			volg )

Definition at line 35 of file s10deri3.F.

      USE message_mod
      USE elbufdef_mod            
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   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "vect01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
C     REAL 
      double precision
     .   xx(mvsiz,10), yy(mvsiz,10), zz(mvsiz,10)
      my_real
     .   vol(mvsiz,5),
     .   rx(*), ry(*), rz(*), sx(*), sy(*), sz(*), tx(*),ty(*),tz(*),
     .   nx(mvsiz,10,5),volu(*),voln(*),volg(mvsiz),
     .   px(mvsiz,10,5),py(mvsiz,10,5),pz(mvsiz,10,5)
      TYPE(ELBUF_STRUCT_), TARGET :: ELBUF_STR
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NGL(*), I,IP,N,K1,K2,K3,K4,K5,K6,K7,K8,K9,K10,
     .        M,IPERM(10,4),
     .        NBAD(MVSIZ)
      double precision
     .   xa(mvsiz,10),ya(mvsiz,10),za(mvsiz,10),
     .   xb(mvsiz,10),yb(mvsiz,10),zb(mvsiz,10),
     .   a4,b4,a4m1,b4m1,aa
      my_real
     .   alph,beta,w
      TYPE(L_BUFEL_) ,POINTER :: LBUF     
      DATA iperm/
     .            2, 4, 3, 1, 9,10, 6, 5, 8, 7,
     .            4, 1, 3, 2, 8, 7,10, 9, 5, 6,
     .            1, 4, 2, 3, 8, 9, 5, 7,10, 6,
     .            1, 2, 3, 4, 5, 6, 7, 8, 9,10/
C-----------------------------------------------
      IF(npt==4)THEN
        alph = zep5854102
        beta  = zep1381966
        w = fourth
      ELSEIF(npt==5)THEN
        alph = half
        beta  = one_over_6
        w = nine_over_20
      ELSE
        alph = -huge(alph)
        beta  = -huge(beta)
        w = -huge(w)
      ENDIF
      a4 = four * alph
      b4 = four * beta
      a4m1  = a4- one
      b4m1  = b4- one
C
      DO i=lft,llt
        rx(i) = xx(i,1) - xx(i,4)
        ry(i) = yy(i,1) - yy(i,4)
        rz(i) = zz(i,1) - zz(i,4)
        sx(i) = xx(i,2) - xx(i,4)
        sy(i) = yy(i,2) - yy(i,4)
        sz(i) = zz(i,2) - zz(i,4)
        tx(i) = xx(i,3) - xx(i,4)
        ty(i) = yy(i,3) - yy(i,4)
        tz(i) = zz(i,3) - zz(i,4)
        volg(i) =zero
      ENDDO
C
      DO n=1,4
        DO i=lft,llt
          xa(i,n) = a4m1*xx(i,n)
          ya(i,n) = a4m1*yy(i,n)
          za(i,n) = a4m1*zz(i,n)
C
          xb(i,n) = b4m1*xx(i,n)
          yb(i,n) = b4m1*yy(i,n)
          zb(i,n) = b4m1*zz(i,n)
        ENDDO
      ENDDO
C
      DO n=5,10
        DO i=lft,llt
          xa(i,n) = a4*xx(i,n)
          ya(i,n) = a4*yy(i,n)
          za(i,n) = a4*zz(i,n)
C
          xb(i,n) = b4*xx(i,n)
          yb(i,n) = b4*yy(i,n)
          zb(i,n) = b4*zz(i,n)
        ENDDO
      ENDDO
C
      DO ip=1,4
        k1 = iperm(1,ip)
        k2 = iperm(2,ip)
        k3 = iperm(3,ip)
        k4 = iperm(4,ip)
        k5 = iperm(5,ip)
        k6 = iperm(6,ip)
        k7 = iperm(7,ip)
        k8 = iperm(8,ip)
        k9 = iperm(9,ip)
        k10= iperm(10,ip)
        lbuf => elbuf_str%BUFLY(1)%LBUF(ip,1,1)
        CALL s10jacob(alph,beta,w,
     .   xb(1,k1),xb(1,k2),xb(1,k3),xa(1,k4),xb(1,k5),
     .   xb(1,k6),xb(1,k7),xb(1,k8),xb(1,k9),xb(1,k10),
     .   xa(1,k8),xa(1,k9),xa(1,k10),
     .   yb(1,k1),yb(1,k2),yb(1,k3),ya(1,k4),yb(1,k5),
     .   yb(1,k6),yb(1,k7),yb(1,k8),yb(1,k9),yb(1,k10),
     .   ya(1,k8),ya(1,k9),ya(1,k10),
     .   zb(1,k1),zb(1,k2),zb(1,k3),za(1,k4),zb(1,k5),
     .   zb(1,k6),zb(1,k7),zb(1,k8),zb(1,k9),zb(1,k10),
     .   za(1,k8),za(1,k9),za(1,k10),
     .   px(1,k1,ip) ,px(1,k2,ip),px(1,k3,ip),px(1,k4,ip),px(1,k5,ip),
     .   px(1,k6,ip) ,px(1,k7,ip),px(1,k8,ip),px(1,k9,ip),px(1,k10,ip),
     .   py(1,k1,ip) ,py(1,k2,ip),py(1,k3,ip),py(1,k4,ip),py(1,k5,ip),
     .   py(1,k6,ip) ,py(1,k7,ip),py(1,k8,ip),py(1,k9,ip),py(1,k10,ip),
     .   pz(1,k1,ip) ,pz(1,k2,ip),pz(1,k3,ip),pz(1,k4,ip),pz(1,k5,ip),
     .   pz(1,k6,ip) ,pz(1,k7,ip),pz(1,k8,ip),pz(1,k9,ip),pz(1,k10,ip),
     .   nx(1,k1,ip) ,nx(1,k2,ip),nx(1,k3,ip),nx(1,k4,ip),nx(1,k5,ip),
     .   nx(1,k6,ip) ,nx(1,k7,ip),nx(1,k8,ip),nx(1,k9,ip),nx(1,k10,ip),
     .   vol(1,ip)   ,lbuf%VOL0DP)
c
      ENDDO
C
C
      IF(npt==5)THEN
        alph = fourth
        beta  = fourth
        a4 = one
        b4 = one
        a4m1  = zero
        b4m1  = zero
        w = - four_over_5
        ip = 5
C
        DO n=1,4
         DO i=lft,llt
          xa(i,n) = zero
          ya(i,n) = zero
          za(i,n) = zero
         ENDDO
        ENDDO
CC
        lbuf => elbuf_str%BUFLY(1)%LBUF(npt,1,1)
        CALL s10jacob(alph,beta,w,
     .   xa(1,k1),xa(1,k2),xa(1,k3),xa(1,k4),xx(1,k5),
     .   xx(1,k6),xx(1,k7),xx(1,k8),xx(1,k9),xx(1,k10),
     .   xx(1,k8),xx(1,k9),xx(1,k10),
     .   ya(1,k1),ya(1,k2),ya(1,k3),ya(1,k4),yy(1,k5),
     .   yy(1,k6),yy(1,k7),yy(1,k8),yy(1,k9),yy(1,k10),
     .   yy(1,k8),yy(1,k9),yy(1,k10),
     .   za(1,k1),za(1,k2),za(1,k3),za(1,k4),zz(1,k5),
     .   zz(1,k6),zz(1,k7),zz(1,k8),zz(1,k9),zz(1,k10),
     .   zz(1,k8),zz(1,k9),zz(1,k10),
     .   px(1,k1,ip) ,px(1,k2,ip),px(1,k3,ip),px(1,k4,ip),px(1,k5,ip),
     .   px(1,k6,ip) ,px(1,k7,ip),px(1,k8,ip),px(1,k9,ip),px(1,k10,ip),
     .   py(1,k1,ip) ,py(1,k2,ip),py(1,k3,ip),py(1,k4,ip),py(1,k5,ip),
     .   py(1,k6,ip) ,py(1,k7,ip),py(1,k8,ip),py(1,k9,ip),py(1,k10,ip),
     .   pz(1,k1,ip) ,pz(1,k2,ip),pz(1,k3,ip),pz(1,k4,ip),pz(1,k5,ip),
     .   pz(1,k6,ip) ,pz(1,k7,ip),pz(1,k8,ip),pz(1,k9,ip),pz(1,k10,ip),
     .   nx(1,k1,ip) ,nx(1,k2,ip),nx(1,k3,ip),nx(1,k4,ip),nx(1,k5,ip),
     .   nx(1,k6,ip) ,nx(1,k7,ip),nx(1,k8,ip),nx(1,k9,ip),nx(1,k10,ip),
     .   vol(1,ip)   ,lbuf%VOL0DP)
      ENDIF
C      
      IF (ismstr==10.OR.ismstr==12) THEN
       DO ip=1,npt
        lbuf => elbuf_str%BUFLY(1)%LBUF(ip,1,1)
        CALL s10pijsav(px(1,1,ip),py(1,1,ip),pz(1,1,ip),lbuf%PIJ,llt)
       ENDDO
      ENDIF
C
      DO i=lft,llt
       nbad(i)=0
      ENDDO
 
      DO ip=1,npt
        DO i=lft,llt
 
          IF(ip<=4.AND.vol(i,ip)<=zero)THEN
              nbad(i)=nbad(i)+1
              vol(i,ip)=zero
          ENDIF
 
          volg(i) =volg(i) + vol(i,ip)
        ENDDO
      ENDDO
 
      DO i=lft,llt
       IF(nbad(i)>0)THEN
         CALL ancmsg(msgid=489,
     .               msgtype=msgerror,
     .               anmode=aninfo_blind_1,
     .               i1=ngl(i))
       ENDIF
      ENDDO
 
      DO i=lft,llt
        volu(i) =volg(i)/npt
        voln(i) =volg(i)
      END DO
C
      RETURN
C
C  1000 FORMAT(/' ZERO OR NEGATIVE VOLUME : 3D-ELEMENT NB',I10/)
C  2000 FORMAT(/' ZERO OR NEGATIVE VOLUME : DELETE 3D-ELEMENT NB',I10/)

s10jacob()

subroutine s10jacob	(		alph,
			beta,
			w,
		double precision, dimension(mvsiz)	x1b,
		double precision, dimension(mvsiz)	x2b,
		double precision, dimension(mvsiz)	x3b,
		double precision, dimension(mvsiz)	x4a,
		double precision, dimension(mvsiz)	x5b,
		double precision, dimension(mvsiz)	x6b,
		double precision, dimension(mvsiz)	x7b,
		double precision, dimension(mvsiz)	x8b,
		double precision, dimension(mvsiz)	x9b,
		double precision, dimension(mvsiz)	x10b,
		double precision, dimension(mvsiz)	x8a,
		double precision, dimension(mvsiz)	x9a,
		double precision, dimension(mvsiz)	x10a,
		double precision, dimension(mvsiz)	y1b,
		double precision, dimension(mvsiz)	y2b,
		double precision, dimension(mvsiz)	y3b,
		double precision, dimension(mvsiz)	y4a,
		double precision, dimension(mvsiz)	y5b,
		double precision, dimension(mvsiz)	y6b,
		double precision, dimension(mvsiz)	y7b,
		double precision, dimension(mvsiz)	y8b,
		double precision, dimension(mvsiz)	y9b,
		double precision, dimension(mvsiz)	y10b,
		double precision, dimension(mvsiz)	y8a,
		double precision, dimension(mvsiz)	y9a,
		double precision, dimension(mvsiz)	y10a,
		double precision, dimension(mvsiz)	z1b,
		double precision, dimension(mvsiz)	z2b,
		double precision, dimension(mvsiz)	z3b,
		double precision, dimension(mvsiz)	z4a,
		double precision, dimension(mvsiz)	z5b,
		double precision, dimension(mvsiz)	z6b,
		double precision, dimension(mvsiz)	z7b,
		double precision, dimension(mvsiz)	z8b,
		double precision, dimension(mvsiz)	z9b,
		double precision, dimension(mvsiz)	z10b,
		double precision, dimension(mvsiz)	z8a,
		double precision, dimension(mvsiz)	z9a,
		double precision, dimension(mvsiz)	z10a,
			px1,
			px2,
			px3,
			px4,
			px5,
			px6,
			px7,
			px8,
			px9,
			px10,
			py1,
			py2,
			py3,
			py4,
			py5,
			py6,
			py7,
			py8,
			py9,
			py10,
			pz1,
			pz2,
			pz3,
			pz4,
			pz5,
			pz6,
			pz7,
			pz8,
			pz9,
			pz10,
			nx1,
			nx2,
			nx3,
			nx4,
			nx5,
			nx6,
			nx7,
			nx8,
			nx9,
			nx10,
			vol,
		double precision, dimension(*)	voldp )

Definition at line 262 of file s10deri3.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-----------------------------------------------
C     REAL
      double precision
     .   x1b(mvsiz),x2b(mvsiz),x3b(mvsiz),x4a(mvsiz),x5b(mvsiz),
     .   x6b(mvsiz),x7b(mvsiz),x8b(mvsiz),x9b(mvsiz),x10b(mvsiz),
     .   x8a(mvsiz),x9a(mvsiz),x10a(mvsiz),
     .   y1b(mvsiz),y2b(mvsiz),y3b(mvsiz),y4a(mvsiz),y5b(mvsiz),
     .   y6b(mvsiz),y7b(mvsiz),y8b(mvsiz),y9b(mvsiz),y10b(mvsiz),
     .   y8a(mvsiz),y9a(mvsiz),y10a(mvsiz),
     .   z1b(mvsiz),z2b(mvsiz),z3b(mvsiz),z4a(mvsiz),z5b(mvsiz),
     .   z6b(mvsiz),z7b(mvsiz),z8b(mvsiz),z9b(mvsiz),z10b(mvsiz),
     .   z8a(mvsiz),z9a(mvsiz),z10a(mvsiz),voldp(*)
      my_real
     .   px1(mvsiz),px2(mvsiz),px3(mvsiz),px4(mvsiz),px5(mvsiz),
     .   px6(mvsiz),px7(mvsiz),px8(mvsiz),px9(mvsiz),px10(mvsiz),
     .   py1(mvsiz),py2(mvsiz),py3(mvsiz),py4(mvsiz),py5(mvsiz),
     .   py6(mvsiz),py7(mvsiz),py8(mvsiz),py9(mvsiz),py10(mvsiz),
     .   pz1(mvsiz),pz2(mvsiz),pz3(mvsiz),pz4(mvsiz),pz5(mvsiz),
     .   pz6(mvsiz),pz7(mvsiz),pz8(mvsiz),pz9(mvsiz),pz10(mvsiz),
     .   nx1(mvsiz),nx2(mvsiz),nx3(mvsiz),nx4(mvsiz),nx5(mvsiz),
     .   nx6(mvsiz),nx7(mvsiz),nx8(mvsiz),nx9(mvsiz),nx10(mvsiz),
     .   vol(mvsiz),alph,beta,w
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
C     REAL
      double precision
     .   a4m1,b4m1,aa
c      my_real
      double precision
     .   d,a4mb4,a4,b4,bb,ab,b2,
     .   dxdr,dxds,dxdt,dydr,dyds,dydt,dzdr,dzds,dzdt
c      my_real
      double precision
     .   drdx(mvsiz), dsdx(mvsiz), dtdx(mvsiz),
     .   drdy(mvsiz), dsdy(mvsiz), dtdy(mvsiz),
     .   drdz(mvsiz), dsdz(mvsiz), dtdz(mvsiz),
     .   det(mvsiz)
C-----------------------------------------------
      a4 = four * alph
      b4 = four * beta
      a4m1  = a4 - one
      b4m1  = b4 - one
      aa = (two*alph - one)*alph
      bb = (two*beta - one)*beta
      b2 = four*beta*beta
      ab = four*alph*beta      
C-----------------------------------------------
      DO i=lft,llt
        nx1(i) = bb
        nx2(i) = bb
        nx3(i) = bb
        nx4(i) = aa
        nx5(i) = b2
        nx6(i) = b2
        nx7(i) = b2
        nx8(i) = ab
        nx9(i) = ab
        nx10(i)= ab
      ENDDO
C
      DO i=lft,llt
       aa = x5b(i) + x6b(i) + x7b(i)
     .    - x4a(i) - x8b(i) - x9b(i) - x10b(i)
       dxdr = x1b(i) +x8a(i) - x6b(i) + aa
       dxds = x2b(i) +x9a(i) - x7b(i) + aa
       dxdt = x3b(i) +x10a(i)- x5b(i) + aa
C
       aa = y5b(i) + y6b(i) + y7b(i)
     .    - y4a(i) - y8b(i) - y9b(i) - y10b(i)
       dydr = y1b(i) +y8a(i) - y6b(i) + aa
       dyds = y2b(i) +y9a(i) - y7b(i) + aa
       dydt = y3b(i) +y10a(i)- y5b(i) + aa
C
       aa = z5b(i) + z6b(i) + z7b(i)
     .    - z4a(i) - z8b(i) - z9b(i) - z10b(i)
       dzdr = z1b(i) +z8a(i) - z6b(i) + aa
       dzds = z2b(i) +z9a(i) - z7b(i) + aa
       dzdt = z3b(i) +z10a(i)- z5b(i) + aa
C
       drdx(i)=dyds*dzdt-dzds*dydt
       dsdx(i)=dydt*dzdr-dzdt*dydr
       dtdx(i)=dydr*dzds-dzdr*dyds
C
       drdy(i)=dzds*dxdt-dxds*dzdt
       dsdy(i)=dzdt*dxdr-dxdt*dzdr
       dtdy(i)=dzdr*dxds-dxdr*dzds
C
       drdz(i)=dxds*dydt-dyds*dxdt
       dsdz(i)=dxdt*dydr-dydt*dxdr
       dtdz(i)=dxdr*dyds-dydr*dxds
C
       det(i) = dxdr * drdx(i)
     .        + dydr * drdy(i)
     .        + dzdr * drdz(i)
       voldp(i) = w * det(i) * one_over_6
       vol(i) = voldp(i)
C
      ENDDO
C      IF(IDTMIN(1)==1)THEN
C       DO I=LFT,LLT
C        IF(OFF(I) ==0.)THEN
C          DET(I)=1.
C        ELSEIF(DET(I)<=VOLMIN)THEN
C          DET(I)=1.
C          OFF(I)=0.
C          WRITE(ISTDO,2000) NGL(I)
C          WRITE(IOUT ,2000) NGL(I)
C        ELSEIF(DET(I)<=0.0)THEN
C          WRITE(ISTDO,1000) NGL(I)
C          WRITE(IOUT ,1000) NGL(I)
C        ENDIF
C       ENDDO
C      ELSEIF(IDTMIN(1)==2)THEN
C       DO I=LFT,LLT
C        IF(OFF(I) ==0.)THEN
C           DET(I)=1.
C        ELSEIF(DET(I)<=VOLMIN.OR.DET(I)<=0.0)THEN
C           DET(I)=1.
C           OFF(I)=0.
C           WRITE(ISTDO,2000) NGL(I)
C           WRITE(IOUT ,2000) NGL(I)
C        ENDIF
C       ENDDO
C      ELSE
C        DO I=LFT,LLT
C          IF(OFF(I) ==0.)DET(I)=1.
C        ENDDO
C      ENDIF
CC
      DO i=lft,llt
        a4 = four * alph
        b4 = four * beta
        a4m1  = a4 - one
        b4m1  = b4 - one
        d = one/max(em30,det(i))
C
        b4m1 = d * b4m1
        a4m1 = d * a4m1
        b4 = d * b4
        a4 = d * a4
        a4mb4  = a4 - b4
C
        px1(i) = b4m1 * drdx(i)
        py1(i) = b4m1 * drdy(i)
        pz1(i) = b4m1 * drdz(i)
C
        px2(i) = b4m1 * dsdx(i)
        py2(i) = b4m1 * dsdy(i)
        pz2(i) = b4m1 * dsdz(i)
C
        px3(i) = b4m1 * dtdx(i)
        py3(i) = b4m1 * dtdy(i)
        pz3(i) = b4m1 * dtdz(i)
C
        px4(i) =-a4m1 * (drdx(i)+dsdx(i)+dtdx(i))
        py4(i) =-a4m1 * (drdy(i)+dsdy(i)+dtdy(i))
        pz4(i) =-a4m1 * (drdz(i)+dsdz(i)+dtdz(i))
C
        px5(i) = b4 * (drdx(i)+dsdx(i))
        py5(i) = b4 * (drdy(i)+dsdy(i))
        pz5(i) = b4 * (drdz(i)+dsdz(i))
C
        px6(i) = b4 * (dsdx(i)+dtdx(i))
        py6(i) = b4 * (dsdy(i)+dtdy(i))
        pz6(i) = b4 * (dsdz(i)+dtdz(i))
C
        px7(i) = b4 * (dtdx(i)+drdx(i))
        py7(i) = b4 * (dtdy(i)+drdy(i))
        pz7(i) = b4 * (dtdz(i)+drdz(i))
C
        px8(i) = a4mb4  * drdx(i) - px6(i)
        py8(i) = a4mb4  * drdy(i) - py6(i)
        pz8(i) = a4mb4  * drdz(i) - pz6(i)
C
        px9(i) = a4mb4  * dsdx(i) - px7(i)
        py9(i) = a4mb4  * dsdy(i) - py7(i)
        pz9(i) = a4mb4  * dsdz(i) - pz7(i)
C
        px10(i)= a4mb4  * dtdx(i) - px5(i)
        py10(i)= a4mb4  * dtdy(i) - py5(i)
        pz10(i)= a4mb4  * dtdz(i) - pz5(i)
C
      ENDDO
C
      RETURN

s10pijsav()

subroutine s10pijsav	(		px,
			py,
			pz,
			pij,
		integer	nel )

Definition at line 479 of file s10deri3.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 NEL
C     REAL
      my_real
     .   px(mvsiz,10),py(mvsiz,10),pz(mvsiz,10),
     .   pij(nel,3,10)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,K
 
C-----------------------------------------------
      DO j=1,10
       DO i=1,nel
        pij(i,1,j) = px(i,j) 
        pij(i,2,j) = py(i,j) 
        pij(i,3,j) = pz(i,j)  
       END DO
      END DO
C-----------
      RETURN