h3d_write_tensor_corner.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	h3d_write_tensor_corner (iok_part, is_written, tensor_corner, nel, offset, nft, value, nnod, is_written_tensor)

Function/Subroutine Documentation

h3d_write_tensor_corner()

subroutine h3d_write_tensor_corner	(	integer, dimension(*), intent(in)	iok_part,
		integer, dimension(*), intent(out)	is_written,
		intent(out)	tensor_corner,
		integer, intent(in)	nel,
		integer, intent(in)	offset,
		integer, intent(in)	nft,
		intent(inout)	value,
		integer, intent(in)	nnod,
		integer, dimension(*), intent(in)	is_written_tensor )

Definition at line 28 of file h3d_write_tensor_corner.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        my_real, DIMENSION(6,*), INTENT(INOUT) :: VALUE
        my_real, DIMENSION(6,*), INTENT(OUT) :: tensor_corner
        INTEGER, INTENT(IN) :: NEL,OFFSET,NFT,NNOD
        INTEGER, DIMENSION(*), INTENT(IN) :: IOK_PART
        INTEGER, DIMENSION(*), INTENT(OUT) :: IS_WRITTEN
        INTEGER, DIMENSION(*), INTENT(IN) :: IS_WRITTEN_TENSOR
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER :: I,J
C-----------------------------------------------
C
        DO i=1,nel
          DO j=1,nnod
            tensor_corner(1,nnod*(offset+nft+i-1)+j) = value(1,20*(i-1)+j)
            tensor_corner(2,nnod*(offset+nft+i-1)+j) = value(2,20*(i-1)+j)
            tensor_corner(3,nnod*(offset+nft+i-1)+j) = value(3,20*(i-1)+j)
            tensor_corner(4,nnod*(offset+nft+i-1)+j) = value(4,20*(i-1)+j)
            tensor_corner(5,nnod*(offset+nft+i-1)+j) = value(5,20*(i-1)+j)
            tensor_corner(6,nnod*(offset+nft+i-1)+j) = value(6,20*(i-1)+j)
            IF(iok_part(i) == 1  .AND. is_written_tensor(i) == 1) is_written(offset+nft+i)= 1
            value(1:6,20*(i-1)+j) = zero
          ENDDO
        ENDDO
 
        RETURN