xyznor.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	xyznor (xnorm, nodglob, weight)

Function/Subroutine Documentation

xyznor()

subroutine xyznor	(		xnorm,
		integer, dimension(*)	nodglob,
		integer, dimension(numnod)	weight )

Definition at line 31 of file xyznor.F.

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"
#include      "task_c.inc"
#include      "spmd_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NODGLOB(*),WEIGHT(NUMNOD)
C     REAL
      my_real
     .   xnorm(3,*) 
C     REAL
      my_real
     .   s3000,s
      INTEGER I,INORM(3),K,P,IXNN
C-----------------------------------------------
      s3000 = three1000
      ixnn = s3000
      IF(fmt_ani==4)ixnn=0
C
 
      IF(nspmd == 1) THEN
        DO i=1,numnod
         s = sqrt(xnorm(1,i)**2 + xnorm(2,i)**2 + xnorm(3,i)**2)
         IF(s/=zero)THEN
           s = s3000 / s
           inorm(1) = xnorm(1,i) * s
           inorm(2) = xnorm(2,i) * s 
           inorm(3) = xnorm(3,i) * s
         ELSE
           inorm(1) = 0
           inorm(2) = 0
           inorm(3) = ixnn
         END IF
           CALL write_s_c(inorm,3)
        END DO
 
      ELSE
         IF(ispmd==0) THEN
           CALL spmd_wrt_xyznor(xnorm,weight,nodglob,numnodg)
         ELSE
           CALL spmd_wrt_xyznor(xnorm,weight,nodglob,1)
         ENDIF
      END IF
 
      IF(nspmd == 1) THEN
        DO i=1,64*numelig3d
           inorm(1) = 0
           inorm(2) = 0
           inorm(3) = ixnn
           CALL write_s_c(inorm,3)
        END DO
 
      ENDIF
C-----------------------------------------------
      RETURN