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

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

!||    flow0              ../engine/source/fluid/flow0.F

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

!||    resol              ../engine/source/engine/resol.F

!||--- calls      -----------------------------------------------------

!||    daasolv            ../engine/source/fluid/daasolv.F

!||    daasolvp           ../engine/source/fluid/daasolvp.F

!||    incpflow           ../engine/source/fluid/incpflow.F

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

!||    output_mod         ../common_source/modules/output/output_mod.F90

!||    python_funct_mod   ../common_source/modules/python_mod.F90

!||    sensor_mod         ../common_source/modules/sensor_mod.F90

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


      SUBROUTINE flow0(OUTPUT,IFLOW  , RFLOW  , WIFLOW , WRFLOW, X      ,

     .                 V      , A      , NPC    , TF    , SENSOR_TAB,

     .                 NBGAUGE, LGAUGE , GAUGE  , NSENSOR,

     .                 IGRV   , AGRV   , NFUNCT , PYTHON, WFEXT)

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

C   M o d u l e s

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

      USE sensor_mod

      USE output_mod

      USE python_funct_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      "com08_c.inc"

#include      "param_c.inc"

#include      "scr07_c.inc"

#include      "flowcom.inc"

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

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

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

      TYPE(output_), INTENT(INOUT) :: OUTPUT

      INTEGER ,INTENT(IN) :: NSENSOR,NBGAUGE,NFUNCT

      INTEGER IFLOW(*), WIFLOW(*), NPC(*),LGAUGE(3,*)

      INTEGER IGRV(NIGRV,*)

      my_real rflow(*), wrflow(*), x(3,*), v(3,*), a(3,*), tf(*), gauge(llgauge,*), agrv(lfacgrv,*)

      TYPE (SENSOR_STR_) ,DIMENSION(NSENSOR) :: SENSOR_TAB

      TYPE (PYTHON_), INTENT(INOUT) :: PYTHON

      DOUBLE PRECISION,INTENT(INOUT) :: WFEXT

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, NNN, NEL_LOC,

     .        II1, II2, II3, II4, II5, II6, II7, IR1, IR2, IR3, IR4,

     .        IR5, IR6, IR7, IR8, IR9, IR10, IR11, II8, II9, II10, II11, IADIP, IADH, IADG,

     .        ii12, ii13, ii14, ipint, ndim, jform, freesurf

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

C Description of arrays ITYP=1

C     II1 -> IFLOW: integer parameters

C     II2  -> IBUF    : local-global surface node correspondences

C     II3 -> Elem: Local Connectivitis of the Triangles

C     II4  -> IINOUT  : integer variables for inlet-outlet surfaces

C     II5  -> IBUFI   : local-global internal node correspondences

C     II6  -> ITAGEL  : tag of elements on inlet-outlet surfaces

C---- SPMD tables

C     II7  -> IBUFR   : local node row index (0 if not on this proc)

C     II8  -> IBUFC   : local node column index (0 if not on this proc)

C     II9 -> IBUFL: List Noes surface on current processor for x, v, a

C     II10 -> CNP: Procedure processor for nodes surface for x, v, a

C     II11 -> Ibufil: Internal noded list on current processor for x, v, a

C     II12 -> CNPI: Award processor for internal nodes for x, v, a

C     II13 -> IBUFELR : local elem row index (0 if not on this proc)

C     II14 -> IBUFELC : local elem column index (0 if not on this proc)

C

C     IR1  -> RFLOW  : parametres reels

C     IR2  -> PHI    : potentiel

C     IR3  -> PRES   : pression

C     IR4 -> U: velocity

C     IR5  -> RINOUT : real variables for inlet-outlet surfaces

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

C Description of arrays ITYP=3

C     II1 -> IFLOW: integer parameters

C     II2  -> IBUF    : local-global surface node correspondences

C     II3 -> Elem: Local connectivitis of triangles/quadrangles+flag

C     II4  -> IBUF_L  : local-global surface node correspondences (to be removed in SPMD)

C     II5  -> SHELL_GA: local shell-gauge correspondences

C     II6  -> CNP     : SPMD number of processors for each node

C engine

C     II7  -> IPIV    : lapack resolution nel > nelmax

C     II8  -> IBUFELR : SPMD

C     II9  -> IBUFELC : SPMD

C

C     IR1  -> RFLOW  : real parameters

C     IR2  -> NORMAL : normal

C     IR3  -> TA     : arrival time

C     IR4  -> AREA   : element area

C     IR5  -> COSG   : direction cosine

C     IR6  -> DIST   : distance charge

C     IR7  -> MFLE   : fluide mass matrix (inverse) nel < nelmax

C     IR7  -> MFLE   : C**t B + B**t C  nel >= nelmax

C     IR8  -> ACCF   : acceleration point fluide

C     IR9  -> PS     : scattered pressure

C     IR10 -> PTI    : incident pressure time integral

C     IR11 -> CMAT   : Matrice C (nel >= nelmax)

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

      ii1 = 1

      ii2 = 1

      ii3 = 1

      ii4 = 1

      ii5 = 1

      ii6 = 1

      ii7 = 1

      ii8 = 1

      ii9 = 1

      ii10= 1

      ipint=0

      IF ((tt>=output%TANIM .AND. tt<=output%TANIM_STOP).OR.

     .   manim==4 .OR.manim==5 .OR.manim==6 .OR.manim==7.OR.

     .   manim==12.OR.manim==13.OR.manim==14.OR.manim==15)

     .   ipint=1

      iadi=0

      iadr=0

      DO i=1,nflow

         ityp=iflow(iadi+2)

         nno=iflow(iadi+5)

         nel=iflow(iadi+6)

         ii1=iadi+1

         ii2=ii1+niflow

         ii3=ii2+nno

         ir1=iadr+1

         ir2=ir1+nrflow

         IF (ityp==1) THEN

            ninout=iflow(iadi+4)

            nnn=iflow(iadi+7)

            ii4=ii3+3*nel

            ii5=ii4+ninout*niioflow

            ii6=ii5+nnn

            IF (nspmd > 1) THEN

               ii7=ii6+nel

               ii8=ii7+nno

               ii9=ii8+nno

               ii10=ii9+nno

               ii11=ii10+nno

               ii12=ii11+nnn

               ii13=ii12+nnn

               ii14=ii13+nel

            ELSE

               ii7=ii6+nel

               ii8=ii7

               ii9=ii8

               ii10=ii9+nno

               ii11=ii10

               ii12=ii11+nnn

               ii13=ii12

               ii14=ii13

            ENDIF

            ir3=ir2+nno+nnn

            ir4=ir3+nno+nnn

            ir5=ir4+3*(nno+nnn)


            iadip=iflow(iadi+10)

            iadh=iflow(iadi+11)

            iadg=iflow(iadi+20)


            CALL incpflow(

     .      nno,         nel,         ninout,      nnn,         iflow(ii1),

     .      iflow(ii2),  iflow(ii3),  iflow(ii4),  iflow(ii5),  iflow(ii7),

     .      iflow(ii8),  iflow(ii9),  rflow(ir1),  rflow(ir2),  rflow(ir3),

     .      rflow(ir4),  rflow(ir5),  x,           v,           a,

     .      npc,         tf        ,  nsensor    , sensor_tab ,

     .      iflow(ii10), iflow(ii6),  iflow(ii13), iflow(ii14), wiflow(iadip),

     .      wrflow(iadh),wrflow(iadg),iflow(ii11), iflow(ii12), ipint,

     .      python      ,wfext)


         ELSEIF(ityp == 3) THEN

            jform    = iflow(iadi+4)

            freesurf = iflow(iadi+25)

            IF(jform==1) THEN

               ii4=ii3+3*nel

               ii5=ii4+nno

               ii6=ii5

               ndim=3

            ELSEIF(jform==2) THEN

               ii4=ii3+5*nel

               ii5=ii4+nno

               ii6=ii5+nbgauge

               ndim=5

            ENDIF


            ir3 = ir2+nel*3

            ir4 = ir3+nel*freesurf

            ir5 = ir4+nel

            ir6 = ir5+nel*freesurf

            ir7 = ir6+nel*freesurf

            ir8 = ir7+nel*nel

            ir9 = ir8+nel

            ir10= ir9+nel

            ir11= ir10+nel


            IF(nspmd == 1) THEN

              ii7=ii6

              CALL daasolv(ndim, nno, nel,

     .                     iflow(ii1), iflow(ii2), iflow(ii3), iflow(ii5),

     .                     rflow(ir1), rflow(ir2), rflow(ir3), rflow(ir4), rflow(ir5), rflow(ir6),

     .                     rflow(ir7), rflow(ir8), rflow(ir9), rflow(ir10),x, v, a, npc, tf,

     .                     nbgauge, lgauge, gauge, nsensor, sensor_tab, igrv, agrv,

     .                     rflow(ir11),iflow(ii7), nfunct, python, wfext)

            ELSE

              ii7=ii6+nno

              ii8=ii7+nel

              ii9=ii8+nel

              nel_loc =iflow(iadi+20)

              iadip =iflow(iadi+10)

              iadh  =iflow(iadi+11)

              CALL daasolvp(ndim, nno, nel, nel_loc,

     .                      iflow(ii1), iflow(ii2), iflow(ii3), iflow(ii4), iflow(ii5), iflow(ii6),

     .                      rflow(ir1), rflow(ir2), rflow(ir3), rflow(ir4), rflow(ir5), rflow(ir6),

     .                      rflow(ir7), rflow(ir8), rflow(ir9), rflow(ir10),x, v, a, npc, tf,

     .                      nbgauge, lgauge, gauge, nsensor, sensor_tab, igrv, agrv,

     .                      rflow(ir11),wiflow(iadip), wrflow(iadh), iflow(ii8), iflow(ii9),

     .                      nfunct, python, wfext)

            ENDIF

         ENDIF

         iadr=iadr+iflow(iadi+15)

         iadi=iadi+iflow(iadi+14)

      ENDDO


      RETURN


      END


my_real
#define my_real
Definition cppsort.cpp:32

daasolv
subroutine daasolv(ndim, nno, nel, iflow, ibuf, elem, shell_ga, rflow, normal, ta, areaf, cosg, dist, mfle, accf, pm, pti, x, v, a, npc, tf, nbgauge, lgauge, gauge, nsensor, sensor_tab, igrv, agrv, cbem, ipiv, nfunct, python, wfext)
Definition daasolv.F:41

daasolvp
subroutine daasolvp(ndim, nno, nel, nel_loc, iflow, ibuf, elem, ibufl, shell_ga, cnp, rflow, normal, ta, areaf, cosg, dist, mfle, accf, pm, pti, x, v, a, npc, tf, nbgauge, lgauge, gauge, nsensor, sensor_tab, igrv, agrv, cbem, ipiv_l, mfle_l, ibufelr, ibufelc, nfunct, python, wfext)
Definition daasolvp.F:42

flow0
subroutine flow0(output, iflow, rflow, wiflow, wrflow, x, v, a, npc, tf, sensor_tab, nbgauge, lgauge, gauge, nsensor, igrv, agrv, nfunct, python, wfext)
Definition flow0.F:40

incpflow
subroutine incpflow(nno, nel, ninout, nni, iflow, ibuf, elem, iinout, ibufi, ibufr, ibufc, ibufl, rflow, phi, pres, u, rinout, x, v, a, npc, tf, nsensor, sensor_tab, cnp, itagel, ibufelr, ibufelc, ipiv, hbem, gbem, ibufil, cnpi, ipint, python, wfext)
Definition incpflow.F:46