xyznod_ply.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	xyznod_ply (iply, idply, nod_pxfem, x, zi_ply, nodglob, empsizpl)

Function/Subroutine Documentation

xyznod_ply()

subroutine xyznod_ply	(	integer	iply,
		integer	idply,
		integer, dimension(*)	nod_pxfem,
			x,
			zi_ply,
		integer, dimension(*)	nodglob,
		integer	empsizpl )

Definition at line 33 of file xyznod_ply.F.

C 
      USE plyxfem_mod        
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IPLY,IDPLY,NOD_PXFEM(*),NODGLOB(*),EMPSIZPL
      my_real
     .   x(3,*),zi_ply(nplyxfe,*),norm
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
      INTEGER I,N,II,ND,ARSZ,EMPL
      REAL, DIMENSION(:,:), ALLOCATABLE :: XNODARRAY
      my_real
     .   vn
C
       ii = idply
      IF (nspmd == 1) THEN
        arsz = plynod(iply)%PLYNUMNODS
        ALLOCATE (xnodarray(3,arsz))
 
        DO nd=1,plynod(iply)%PLYNUMNODS
             i = plynod(iply)%NODES(nd)
             empl = plynod(iply)%PLYNODID(nd)-empsizpl
             n = nod_pxfem(i)
             norm = sqrt(vn_nod(1,n)**2 + vn_nod(2,n)**2 + vn_nod(3,n)**2)
             norm = one/max(em20, norm)
             vn = vn_nod(1,n)*norm
             xnodarray(1,empl)=x(1,i)+zi_ply(n,iply)*vn+ply(iply)%U(1,n)
             vn = vn_nod(2,n)*norm
             xnodarray(2,empl)=x(2,i)+zi_ply(n,iply)*vn+ply(iply)%U(2,n)
             vn = vn_nod(3,n)*norm
             xnodarray(3,empl)=x(3,i)+zi_ply(n,iply)*vn+ply(iply)%U(3,n)
        END DO
        CALL write_r_c(xnodarray,3*arsz)
        DEALLOCATE(xnodarray)
        empsizpl = empsizpl + plynod(iply)%PLYNUMNODS
      ELSE
        CALL spmd_anim_ply_xyznod( nodglob,iply, idply,
     *                             nod_pxfem, x, zi_ply,empsizpl )
      ENDIF
C-----------------------------------------------
      RETURN