h3d_write_tensor_corner.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    h3d_write_tensor_corner   ../engine/source/output/h3d/h3d_results/h3d_write_tensor_corner.F
!||--- called by ------------------------------------------------------
!||    h3d_solid_tensor_1        ../engine/source/output/h3d/h3d_results/h3d_solid_tensor_1.F
!||====================================================================
        SUBROUTINE h3d_write_tensor_corner(
     .                          IOK_PART  ,IS_WRITTEN ,TENSOR_CORNER ,NEL ,OFFSET ,NFT,
     .                          VALUE     ,NNOD  ,IS_WRITTEN_TENSOR)
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
        END SUBROUTINE h3d_write_tensor_corner