resol_head.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/.
!||====================================================================
!||    resol_head             ../engine/source/engine/resol_head.F
!||--- called by ------------------------------------------------------
!||    radioss2               ../engine/source/engine/radioss2.F
!||--- calls      -----------------------------------------------------
!||    resol                  ../engine/source/engine/resol.F
!||    trace_in               ../engine/source/system/trace_back.F
!||    trace_out              ../engine/source/system/trace_back.F
!||--- uses       -----------------------------------------------------
!||    ale_connectivity_mod   ../common_source/modules/ale/ale_connectivity_mod.F
!||    cluster_mod            ../engine/share/modules/cluster_mod.F
!||    connectivity_mod       ../common_source/modules/connectivity.F90
!||    coupling_adapter_mod   ../engine/source/coupling/coupling_adapter.f90
!||    damp_mod               ../engine/share/modules/damp_mod.F
!||    diffusion_mod          ../engine/share/modules/diffusion_mod.F
!||    drape_mod              ../engine/share/modules/drape_mod.F
!||    dsgraph_mod            ../engine/share/modules/dsgraph_mod.F
!||    dt_mod                 ../engine/source/modules/dt_mod.F
!||    ebcs_mod               ../common_source/modules/boundary_conditions/ebcs_mod.F90
!||    eig_mod                ../engine/share/modules/eig_mod.F
!||    fxb_mod                ../engine/share/modules/fxb_mod.F
!||    glob_therm_mod         ../common_source/modules/mat_elem/glob_therm_mod.F90
!||    groupdef_mod           ../common_source/modules/groupdef_mod.F
!||    h3d_mod                ../engine/share/modules/h3d_mod.F
!||    intbuf_fric_mod        ../common_source/modules/interfaces/intbuf_fric_mod.F90
!||    interfaces_mod         ../common_source/modules/interfaces/interfaces_mod.f90
!||    loads_mod              ../common_source/modules/loads/loads_mod.F90
!||    mat_elem_mod           ../common_source/modules/mat_elem/mat_elem_mod.F90
!||    monvol_struct_mod      ../engine/share/modules/monvol_struct_mod.F
!||    multi_fvm_mod          ../common_source/modules/ale/multi_fvm_mod.F90
!||    names_and_titles_mod   ../common_source/modules/names_and_titles_mod.F
!||    nodal_arrays_mod       ../common_source/modules/nodal_arrays.f90
!||    output_mod             ../common_source/modules/output/output_mod.F90
!||    pblast_mod             ../common_source/modules/loads/pblast_mod.F90
!||    pinchtype_mod          ../common_source/modules/pinchtype_mod.F
!||    python_funct_mod       ../common_source/modules/python_mod.F90
!||    rbe3_mod               ../common_source/modules/constraints/rbe3_mod.F90
!||    restmod                ../engine/share/modules/restart_mod.F
!||    rwall_mod              ../common_source/modules/constraints/rwall_mod.f90
!||    segvar_mod             ../engine/share/modules/segvar_mod.F
!||    sensor_mod             ../common_source/modules/sensor_mod.F90
!||    skew_mod               ../common_source/modules/skew_mod.F90
!||    soltosph_mod           ../engine/share/modules/soltosph_mod.F
!||    sph_mod                ../engine/share/modules/sph_mod.F
!||    stack_mod              ../engine/share/modules/stack_mod.F
!||    state_mod              ../common_source/modules/state_mod.F
!||    table_glob_mod         ../engine/share/modules/table_glob_mod.F
!||    timer_mod              ../engine/source/system/timer_mod.f90
!||    unitab_mod             ../common_source/modules/unitab_mod.F
!||    user_windows_mod       ../common_source/modules/user_windows_mod.F
!||====================================================================
      SUBROUTINE resol_head( ELEMENT, NODES,
     .           ITASK       ,AF          ,IAF            ,IDATA      ,RDATA,
     .           GRAPHE      ,IFLOW       ,RFLOW          ,
     .           STACK       ,IRUNN_BIS   , TIMERS,
     .           MULTI_FVM   ,H3D_DATA    ,SUBSETS    ,IGRNOD          ,
     .           IGRBRIC     ,IGRQUAD     ,IGRSH4N        ,IGRSH3N    ,IGRTRUSS        ,
     .           IGRBEAM     ,IGRSPRING   ,IGRPART        ,IGRSURF    ,IGRSLIN         , 
     .           FORNEQS     ,PINCH_DATA  ,ALE_CONNECTIVITY,
     .           DRAPE_SH4N  ,DRAPE_SH3N  ,T_MONVOL       ,SENSORS    ,EBCS_TAB        ,
     .           DIFFUSION   ,SEGVAR      ,DYNAIN_DATA    ,DRAPEG     ,USER_WINDOWS    ,
     .           OUTPUT      ,INTERFACES  ,DT             ,LOADS      ,MAT_ELEM        ,
     .           PYTHON      ,NAMES_AND_TITLES, UNITAB    ,SKEWS      ,LIFLOW          ,
     .           LRFLOW      ,glob_therm ,PBLAST          ,RBE3       ,coupling        ,
     .           RWALL  )
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE connectivity_mod
      use nodal_arrays_mod
      USE restmod
      USE dsgraph_mod
      USE table_glob_mod
      USE sph_mod
      USE soltosph_mod
      USE cluster_mod
      USE stack_mod
      USE multi_fvm_mod
      USE h3d_mod
      USE intbuf_fric_mod
      USE groupdef_mod
      USE mat_elem_mod
      USE pinchtype_mod
      USE ale_connectivity_mod
      USE monvol_struct_mod
      USE drape_mod
      USE fxb_mod
      USE eig_mod
      USE sensor_mod
      USE ebcs_mod
      USE diffusion_mod
      USE segvar_mod
      USE loads_mod
      USE state_mod
      USE user_windows_mod
      USE output_mod
      USE interfaces_mod
      USE dt_mod
      USE python_funct_mod
      USE damp_mod
      USE names_and_titles_mod
      USE unitab_mod
      USE skew_mod
      use glob_therm_mod
      use pblast_mod
      use timer_mod
      use rbe3_mod
      use coupling_adapter_mod 
      use rwall_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C------------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "scr_fac_c.inc"
#include      "scr01_c.inc"
#include      "task_c.inc"
#include      "com_xfem1.inc"
C-----------------------------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      type(connectivity_), intent(inout) :: ELEMENT 
      type(nodal_arrays_), intent(inout) :: NODES
      type(timer_), intent(inout) :: TIMERS
      type(coupling_type), intent(inout) :: coupling
      INTEGER, INTENT(IN) :: LIFLOW  !< Size of IFLOW
      INTEGER, INTENT(IN) :: LRFLOW  !< Size of RFLOW
 
      INTEGER,INTENT(INOUT) :: IFLOW(LIFLOW)
      INTEGER,INTENT(INOUT) :: RFLOW(LRFLOW)
C
      INTEGER ITASK,IAF(*),IDATA(*),IRUNN_BIS
C     REAL
      my_real
     .   af(*),rdata(*),                          
     .   forneqs(*)
      TYPE(prgraph) :: GRAPHE(*)
      TYPE (STACK_PLY) :: STACK
      TYPE(multi_fvm_struct) :: MULTI_FVM
      TYPE(h3d_database) :: H3D_DATA
      TYPE(UNIT_TYPE_) :: UNITAB
!
      TYPE(subset_) ,DIMENSION(NSUBS)   :: SUBSETS
      TYPE(GROUP_)  ,DIMENSION(NGRNOD)  :: IGRNOD
      TYPE(group_)  ,DIMENSION(NGRBRIC) :: IGRBRIC
      TYPE(GROUP_)  ,DIMENSION(NGRQUAD) :: IGRQUAD
      TYPE(group_)  ,DIMENSION(NGRSHEL) :: IGRSH4N
      TYPE(GROUP_)  ,DIMENSION(NGRSH3N) :: IGRSH3N
      TYPE(GROUP_)  ,DIMENSION(NGRTRUS) :: IGRTRUSS
      TYPE(group_)  ,DIMENSION(NGRBEAM) :: IGRBEAM
      TYPE(group_)  ,DIMENSION(NGRSPRI) :: IGRSPRING
      TYPE(group_)  ,DIMENSION(NGRPART) :: IGRPART
      TYPE(surf_)   ,DIMENSION(NSURF)   :: IGRSURF
      TYPE(surf_)   ,DIMENSION(NSURF)   :: IGRSLIN
      TYPE(PINCH) :: PINCH_DATA
      TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
      TYPE(monvol_struct_), DIMENSION(NVOLU), INTENT(INOUT) :: T_MONVOL
      TYPE(drape_)                            :: DRAPE_SH3N(NUMELTG_DRAPE),DRAPE_SH4N(NUMELC_DRAPE)
      TYPE(drapeg_)                           :: DRAPEG
      TYPE(sensors_)  ,INTENT(INOUT)          :: SENSORS
      TYPE(t_ebcs_tab), INTENT(INOUT)         :: EBCS_TAB
      TYPE(t_diffusion)    ,INTENT(INOUT)     :: DIFFUSION
      TYPE(t_segvar)                          :: SEGVAR
      TYPE(dynain_database), INTENT(INOUT)    :: DYNAIN_DATA
      TYPE(user_windows_),INTENT(INOUT)       :: USER_WINDOWS
      TYPE(output_),INTENT(INOUT)             :: OUTPUT
      TYPE(interfaces_),INTENT(INOUT)         :: INTERFACES
      TYPE(dt_),INTENT(INOUT)                 :: DT
      TYPE(loads_)    ,INTENT(INOUT)          :: LOADS
      TYPE(mat_elem_) ,INTENT(INOUT)          :: MAT_ELEM
      TYPE(python_)   , INTENT(INOUT)         :: PYTHON
      TYPE(names_and_titles_),INTENT(INOUT)   :: NAMES_AND_TITLES !< NAMES_AND_TITLES host the input deck names and titles for outputs
      TYPE(skew_), INTENT(INOUT)              :: SKEWS
      type (glob_therm_) ,intent(inout)       :: glob_therm
      type (pblast_) ,intent(inout)           :: PBLAST
      type (rbe3_) ,intent(inout)             :: RBE3
      type (rwall_), intent(inout)            :: RWALL
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IUN,IUNUSED
      INTEGER ITRACE(10)
      DATA iun/1/,iunused/1/
      my_real
     .  fsavd(1)
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
C-----------------------------------------------
C     TRACE BACK  
C-----------------------------------------------
      itaskp1 = 1
      itaskp1_def = 1
      itrace(1)=0          
      itrace(2)=nthread
      CALL trace_in(2,itrace,zero)
C
      CALL resol(timers, element,nodes,coupling,
     .  af              ,iaf            ,
     1  iskwn          ,neth           ,ipart          ,nom_opt     ,kxx              ,ixx       ,
     2  ixtg            ,ixs            ,ixq            ,ixt            ,ixp         ,ixr       ,
     3  ifill          ,mat_elem       ,
     4  ims             ,npc            ,ibcl           ,
     5  ibfv            ,iecran         ,ilas           ,laccelm        ,nnlink         ,lnlink      ,
     6  iparg           ,dd_iad         ,igrv           ,iexlnk         ,
     7  ipari           ,iconx          ,npby           ,
     8  lpby            ,lrivet         ,nstrf          ,ljoint         ,nodpor         ,monvol      ,
     9  ilink     ,
     a  llink           ,linale         ,neflsw         ,nnflsw         ,icut           ,cluster     ,
     b  itask           ,iaf(if01)      ,
     c  thke        ,damp               ,
     d  pm             ,skews          ,geo            ,eani           ,bufmat      ,bufgeo           ,bufsf     ,
     e  w               ,veul           ,fill           ,dfill          ,alph           ,wb          ,dsave            ,asave     ,
     .  msnf            , 
     f  tf              ,forc           ,vel            ,fsav           ,fzero          ,xlas        ,accelm           ,
     g  grav            ,fr_wave        ,failwave       ,parts0         ,elbuf          ,sensors,
     h  rby             ,rivet          ,secbuf         ,volmon         ,lambda      ,
     i  wa              ,fv             ,partsav          ,
     j  uwa             ,val2           ,phi            ,segvar         ,r              ,crflsw      ,
     k  flsw            ,xcut           ,tani           ,secfcum     ,af(mf01),
     l  idata           ,rdata          ,
     m  iframe          ,kxsp           ,ixsp           ,nod2sp         ,ispsym         ,ispcond     ,
     n  xframe          ,spbuf          ,xspsym         ,vspsym         ,pv             ,
     o  fsavd           ,ibvel          ,lbvel          ,wasph          ,w16         ,
     p  isphio          ,lprtsph        ,lonfsph        ,vsphio         ,fbvel          ,lagbuf      ,ibcslag,
     q  iactiv          ,dampr          ,gjbufi         ,gjbufr         ,rbmpc          ,ibmpc       ,sphveln,
     t  nbrcvois        ,nbsdvois       ,lnrcvois       ,lnsdvois       ,nercvois       ,nesdvois    ,lercvois         ,
     u  lesdvois        ,npsegcom       ,lsegcom        ,nporgeo        ,
     v  ixtg1           ,npbyl          ,lpbyl          ,rbyl           ,igeo           ,ipm         ,
     w  madprt         ,madsh4         ,madsh3         ,madsol         ,madnod      ,madfail          ,
     x  iad_rby         ,fr_rby         ,iad_rby2       ,
     y  fr_rby2         ,iad_i2m        ,fr_i2m         ,addcni2        ,procni2        ,iadi2       ,fr_mv            ,
     z  iadmv2          ,fr_ll          ,fr_rl          ,iadcj          ,
     a  fr_cj           ,fr_sec         ,
     b  iad_sec         ,iad_cut        ,fr_cut         ,rg_cut         ,newfront       ,fr_mad      ,
     c  fxbipm          ,fxbrpm         ,fxbnod         ,fxbmod         ,fxbglm         ,fxbcpm      ,fxbcps           ,
     d  fxblm           ,fxbfls         ,fxbdls         ,fxbdep         ,fxbvit         ,fxbacc      ,fxbelm           ,
     e  fxbsig          ,fxbgrvi        ,fxbgrvr        ,eigipm         ,eigibuf        ,eigrpm      ,
     f  lnodpor         ,
     g  fr_i18          ,graphe         ,iflow          ,rflow,
     h  lgrav           ,dd_r2r         ,fasolfr        ,fr_lagf        ,
     j  llagf           ,icontact       ,rcontact       ,sh4tree        ,sh3tree     ,
     k  ipadmesh        ,padmesh        ,msc            ,mstg           ,inc            ,intg        ,ptg              ,
     l  iskwp           ,nskwp          ,isensp         ,nsensp         ,iaccp          ,naccp       ,ipart_state      ,
     m  acontact        ,pcontact       ,factiv         ,
     p  sh4trim         ,sh3trim        ,mscnd          ,incnd          ,
     q  ibfflux         ,fbfflux        ,rbym           ,irbym          ,lnrbym         ,icodrbym,
     r  ibcv            ,fconv          ,ibftemp        ,fbftemp        ,iad_rbym       ,fr_rbym,
     t  weight_rm       ,ms_ply         ,zi_ply         ,inod_pxfem     ,iel_pxfem      ,iadc_pxfem  ,
     u  adsky_pxfem     ,icode_ply      ,icodt_ply      ,iskew_ply      ,admsms      ,
     v  madclnod        ,nom_sect       ,mcpc           ,mcptg          ,dmelc          ,dmeltg      ,mssa             ,
     w  dmels           ,mstr           ,dmeltr         ,msp            ,dmelp          ,msrt        ,dmelrt           ,
     x  ibcr            ,fradia         ,res_sms        ,table          ,irbe2          ,lrbe2       ,iad_rbe2         ,
     y  fr_rbe2         ,phie           ,msf            ,
     z  procne_pxfem    ,iadsdp_pxfem   ,iadrcp_pxfem   ,icfield        ,lcfield        ,cfield      ,
     1  msz2            ,diag_sms       ,iloadp         ,lloadp      ,loadp            ,
     2  inod_crkxfem    ,iel_crkxfem    ,iadc_crkxfem   ,adsky_crkxfem  ,
     3  cne_crkxfem     ,procne_crkxfem ,iadsdp_crkxfem ,iadrcp_crkxfem ,ibufssg_io     ,
     4  ibc_ply         ,dmint2         ,ibordnode      ,
     5  mat_elem%ELBUF  ,por            ,nodedge        ,iad_edge       ,
     6  fr_edge         ,fr_nbedge      ,crknodiad      ,lgauge         ,gauge          ,
     7  igaup           ,ngaup          ,nodlevxf       ,dd_r2r_elem    ,
     8  nodglobxfe      ,sph2sol        ,sol2sph        ,irst           ,
     9  dmsph           ,wagap          ,mat_elem%XFEM_TAB,elcutc       ,nodenr         ,
     a  kxfenod2elc     ,enrtag         ,
     b  rthbuf          ,kxig3d         ,ixig3d         ,knot           ,wige           ,
     c  wsmcomp        ,stack          ,
     d  cputime_mp_glob ,cputime_mp     ,tab_ump        ,poin_ump       ,sol2sph_typ    ,
     e  irunn_bis       ,addcsrect      ,iad_frnor      ,fr_nor         ,procnor        ,
     f  interfaces%SPMD_ARRAYS%IAD_FREDG,interfaces%SPMD_ARRAYS%FR_EDG  ,drape_sh4n     ,drape_sh3n     ,tab_mat        , 
     g  nativ0_sms      ,multi_fvm      ,segquadfr      ,ms_2d          ,
     h  h3d_data        ,subsets        ,igrnod         ,igrbric        ,
     i  igrquad         ,igrsh4n        ,igrsh3n        ,igrtruss       ,igrbeam        ,
     j  igrspring       ,igrpart        ,igrsurf        ,forneqs        ,
     k  nloc_dmg        ,iskwp_l        ,knotlocpc      ,knotlocel      ,pinch_data     ,tag_skins6  ,ale_connectivity ,
     l  xcell           , xface         ,ne_nercvois    ,ne_nesdvois    ,ne_lercvois    ,ne_lesdvois ,ibcscyc          ,lbcscyc   ,
     m  t_monvol        ,id_global_vois ,face_vois      ,dynain_data    ,fcont_max   ,ebcs_tab         ,diffusion ,
     n  kloadpinter     ,loadpinter     ,dgaploadint    ,drapeg         ,user_windows ,output          ,interfaces,
     o  dt              ,loads          ,python         ,dpl0cld        ,vel0cld        ,
     p  ndamp_vrel      ,id_damp_vrel   ,fr_damp_vrel   ,ndamp_vrel_rbyg,names_and_titles,unitab,liflow,lrflow,
     r  glob_therm      ,pblast         ,rbe3           ,rwall          )
C
      CALL trace_out(2)
C
      RETURN
      END