nodalzvol_8F_source.html

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/.

!||====================================================================

!||    nodalzvol        ../engine/source/output/anim/generate/nodalzvol.F

!||--- called by ------------------------------------------------------

!||    genani           ../engine/source/output/anim/generate/genani.F

!||--- uses       -----------------------------------------------------

!||    elbufdef_mod     ../common_source/modules/mat_elem/elbufdef_mod.F90

!||    i22bufbric_mod   ../common_source/modules/interfaces/cut-cell-search_mod.F

!||    i22edge_mod      ../common_source/modules/interfaces/cut-cell-buffer_mod.F

!||    i22tri_mod       ../common_source/modules/interfaces/cut-cell-search_mod.F

!||    initbuf_mod      ../engine/share/resol/initbuf.F

!||====================================================================


      SUBROUTINE nodalzvol (IFUNC    , WA4, IFLOW, RFLOW, IPARG,

     .                       ELBUF_TAB, IX , NIX  , ITAB , NV46 )

C-----------------------------------------------

C   D e s c r i p t i o n

C-----------------------------------------------

C     This suroutine display polyhedra volume (old & new)

C   using free node as marker (inter22 only, see input card for grnod_id)

C-----------------------------------------------

C   P r e - C o n d i t i o n s

C-----------------------------------------------

C     Tested below during NG LOOP : IALEL > 0

C        where IALEL =IPARG(7,NG)+IPARG(11,NG)

C-----------------------------------------------

C   M o d u l e s

C-----------------------------------------------

      USE initbuf_mod

      USE elbufdef_mod

      USE i22bufbric_mod

      USE i22edge_mod

      USE i22tri_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      "vect01_c.inc"

#include      "param_c.inc"

#include      "inter22.inc"

C-----------------------------------------------

C   D u m m y   A r g u m e n t s

C-----------------------------------------------

      INTEGER,INTENT(IN) :: IFUNC, IFLOW(*),IPARG(NPARG,*),IX(NIX,*),ITAB(*),NIX,NV46

      my_real,INTENT(IN) :: rflow(*)

      REAL,INTENT(INOUT) :: WA4(*)


      TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB

C-----------------------------------------------

C   L o c a l   V a r i a b l e s

C-----------------------------------------------

      INTEGER                            :: IADI, IADR, I, ITYP, NINOUT, NNO, NEL, II1, II2,K1,K,

     .                                      ir1, ir2, j, jj, nno_l, nni_l, ii3, ii4, jjj, nni,

     .                                      ialel,nnod,ipos,iv,ngv,j1,j2,ibv, mlw,numnodcell,

     .                                      ng, kcvt, ii, nbf, nbl, ib, icell, nin, mcell,ncell,

     .                                      iphase,node_id,inew

      TYPE(g_bufel_)  ,POINTER           :: GBUF,GBUFv

      my_real, ALLOCATABLE, DIMENSION(:) :: count_vol

      my_real                            :: p,vf,d,v

      INTEGER,DIMENSION(:,:), POINTER    :: pAdjBRICK

C-----------------------------------------------

C   D e s c r i p t i o n

C-----------------------------------------------

C   This subroutine writes nodal VFRAC

C   /INTER/TYPE22  (only).

C-----------------------------------------------

C   P r e c o n d i t i o n s

C-----------------------------------------------

         IF(int22==0)RETURN

C-----------------------------------------------

C   S o u r c e   L i n e s

C-----------------------------------------------


         nnod   = nix-3   !8-node brick or 4-node quad

         inew = 29-ifunc  ! 1: new cell volume, 0:old


         !---------------------------------------------------------!

         !         /INTER/TYPE22                                   !

         !---------------------------------------------------------!

         !   CELL CENTROIDS ARE MARKED IN CUT CELL BUFFER          !

         !   PURPOSE IS TO DISPLAY CELL SCALR USING FREE NODE      !

         !---------------------------------------------------------!

         !---1. TAG FOR INTERSECTED BRICKS---!

         nbf = 1

         nbl = nb

         nin = 1

         !---1. COMPUTE NODAL PRESSURE---!

         DO ng = 1, ngroup

           mlw   = iparg(1,ng)

           nel   = iparg(2,ng)

           nft   = iparg(3,ng)

           ityp  = iparg(5,ng)

           ialel = iparg(7,ng)+iparg(11,ng)

           IF(ityp/=1 .AND. ityp/=2)cycle

           IF(ialel==0)cycle

           IF(mlw/=6 .AND. mlw/=11 .AND. mlw/=37.AND.mlw/=51)cycle

           gbuf => elbuf_tab(ng)%GBUF

           nin  = 1

           DO i=1,nel

             ib = nint(gbuf%TAG22(i))

             !---------------------------!

             ! NOT A CUT CELL            !

             !---------------------------!

             IF(ib==0)cycle

             !---------------------------!

             !        CUT CELL           !

             !---------------------------!

             icell   =  0

             ncell   = brick_list(nin,ib)%NBCUT

             DO WHILE (icell<=ncell) ! loop on polyhedron {1:NCELL} U {9}

               icell = icell +1

               IF (icell>ncell .AND. ncell/=0)icell=9

               node_id      = brick_list(nin,ib)%POLY(icell)%ID_FREE_NODE

               IF(node_id<=0)EXIT ! no enough nodes in the group or debug option not compatible with smp

               IF(inew == 1)THEN

                 wa4(node_id) = brick_list(nin,ib)%POLY(icell)%VNEW

               ELSE

                 wa4(node_id) = brick_list(nin,ib)%POLY(icell)%VOLD

               ENDIF

             ENDDO

           enddo!next I

         enddo!next NG


C-----------------------------------------------

      RETURN


      END

my_real
#define my_real
Definition cppsort.cpp:32

i22bufbric_mod
Definition cut-cell-search_mod.F:95

i22bufbric_mod::brick_list
type(brick_entity), dimension(:,:), allocatable, target brick_list
Definition cut-cell-search_mod.F:105

i22edge_mod
Definition cut-cell-buffer_mod.F:48

i22tri_mod
Definition cut-cell-search_mod.F:212

i22tri_mod::nb
integer nb
Definition cut-cell-search_mod.F:487

initbuf_mod
Definition initbuf.F:160

nodalzvol
subroutine nodalzvol(ifunc, wa4, iflow, rflow, iparg, elbuf_tab, ix, nix, itab, nv46)
Definition nodalzvol.F:36