xyznor_ply.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "scr14_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	xyznor_ply (iply, xnorm, nodglob, weight, empsizpl)

Function/Subroutine Documentation

xyznor_ply()

subroutine xyznor_ply	(	integer	iply,
			xnorm,
		integer, dimension(*)	nodglob,
		integer, dimension(numnod)	weight,
		integer	empsizpl )

Definition at line 33 of file xyznor_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      "com04_c.inc"
#include      "scr14_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NODGLOB(*),WEIGHT(NUMNOD),EMPSIZPL
      INTEGER  IPLY
C     REAL
      my_real
     .   xnorm(3,*) 
C     REAL
      my_real
     .   s3000,s
      INTEGER I,INORM(3),K,P,IXNN,N,ARSZ,EMPL
      INTEGER, DIMENSION(:,:), ALLOCATABLE :: INORMAR
C-----------------------------------------------
      s3000 = three1000
      ixnn = s3000
      IF(fmt_ani==4)ixnn=0
C
      IF(nspmd == 1 ) THEN
        arsz = plynod(iply)%PLYNUMNODS
        ALLOCATE(inormar(3,arsz))
 
        DO n=1,plynod(iply)%PLYNUMNODS
          i = plynod(iply)%NODES(n)
          empl = plynod(iply)%PLYNODID(n)-empsizpl
          s = sqrt(xnorm(1,i)**2 + xnorm(2,i)**2 + xnorm(3,i)**2)
          IF(s/=zero)THEN
           s = s3000 / s
           inormar(1,empl) = xnorm(1,i) * s
           inormar(2,empl) = xnorm(2,i) * s
           inormar(3,empl) = xnorm(3,i) * s
          ELSE
           inormar(1,empl) = 0
           inormar(2,empl) = 0
           inormar(3,empl) = ixnn
          END IF
        END DO
        CALL write_s_c(inormar,3*arsz)
        DEALLOCATE(inormar)
        empsizpl = empsizpl + plynod(iply)%PLYNUMNODS
      ELSE
        CALL spmd_anim_ply_xyznor( nodglob,iply,xnorm,empsizpl)
      END IF
C-----------------------------------------------
      RETURN