st_qaprint_driver.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/.
!||====================================================================
!||    st_qaprint_driver               ../starter/source/output/qaprint/st_qaprint_driver.F
!||--- called by ------------------------------------------------------
!||    lectur                          ../starter/source/starter/lectur.F
!||--- calls      -----------------------------------------------------
!||    st_qaprint_admesh               ../starter/source/output/qaprint/st_qaprint_admesh.F
!||    st_qaprint_ale_options_driver   ../starter/source/output/qaprint/st_qaprint_ale_options_driver.F
!||    st_qaprint_clusters             ../starter/source/output/qaprint/st_qaprint_clusters.F
!||    st_qaprint_composite_options    ../starter/source/output/qaprint/st_qaprint_composite_options.f
!||    st_qaprint_constraints          ../starter/source/output/qaprint/st_qaprint_constraints.F
!||    st_qaprint_dfs_detonators       ../starter/source/output/qaprint/st_qaprint_dfs_detonators.F
!||    st_qaprint_dfs_lasers           ../starter/source/output/qaprint/st_qaprint_dfs_lasers.F
!||    st_qaprint_ebcs                 ../starter/source/output/qaprint/st_qaprint_ebcs.F
!||    st_qaprint_element              ../starter/source/output/qaprint/st_qaprint_element.F
!||    st_qaprint_friction             ../starter/source/output/qaprint/st_qaprint_friction.F
!||    st_qaprint_general_controls     ../starter/source/output/qaprint/st_qaprint_general_controls.F
!||    st_qaprint_groups               ../starter/source/output/qaprint/st_qaprint_groups.F
!||    st_qaprint_initial_conditions   ../starter/source/output/qaprint/st_qaprint_initial_conditions.F
!||    st_qaprint_initial_state        ../starter/source/output/qaprint/st_qaprint_initial_state.F
!||    st_qaprint_inivol               ../starter/source/output/qaprint/st_qaprint_inivol.F
!||    st_qaprint_interfaces           ../starter/source/output/qaprint/st_qaprint_interfaces.F
!||    st_qaprint_internal_groups      ../starter/source/output/qaprint/st_qaprint_internal_groups.F
!||    st_qaprint_loads                ../starter/source/output/qaprint/st_qaprint_loads.F
!||    st_qaprint_madymo               ../starter/source/output/qaprint/st_qaprint_madymo.F
!||    st_qaprint_materials            ../starter/source/output/qaprint/st_qaprint_materials.F
!||    st_qaprint_model_tools          ../starter/source/output/qaprint/st_qaprint_model_tools.F
!||    st_qaprint_monvol               ../starter/source/output/qaprint/st_qaprint_monvol.F
!||    st_qaprint_multidomains         ../starter/source/output/qaprint/st_qaprint_multidomains.F
!||    st_qaprint_nodes                ../starter/source/output/qaprint/st_qaprint_nodes.F
!||    st_qaprint_output_databases     ../starter/source/output/qaprint/st_qaprint_output_databases.F
!||    st_qaprint_properties           ../starter/source/output/qaprint/st_qaprint_properties.F
!||    st_qaprint_reference_state      ../starter/source/output/qaprint/st_qaprint_reference_state.F
!||    st_qaprint_refsta               ../starter/source/output/qaprint/st_qaprint_refsta.F
!||    st_qaprint_seatbelts            ../starter/source/output/qaprint/st_qaprint_seatbelts.f
!||    st_qaprint_set                  ../starter/source/output/qaprint/st_qaprint_set.F
!||    st_qaprint_surf                 ../starter/source/output/qaprint/st_qaprint_surf.F
!||    st_qaprint_time_histories       ../starter/source/output/qaprint/st_qaprint_time_histories.f
!||    st_qaprint_transformations      ../starter/source/output/qaprint/st_qaprint_transformations.F
!||    st_qaprint_userwi               ../starter/source/output/qaprint/st_qaprint_userwi.F
!||--- uses       -----------------------------------------------------
!||    cluster_mod                     ../starter/share/modules1/cluster_mod.F
!||    defaults_mod                    ../starter/source/modules/defaults_mod.F90
!||    detonators_mod                  ../starter/share/modules1/detonators_mod.F
!||    drape_mod                       ../starter/share/modules1/drape_mod.F
!||    func2d_mod                      ../starter/share/modules1/func2d_mod.F
!||    inimap1d_mod                    ../starter/share/modules1/inimap1d_mod.F
!||    inimap2d_mod                    ../starter/share/modules1/inimap2d_mod.F
!||    monvol_struct_mod               ../starter/share/modules1/monvol_struct_mod.F
!||    submodel_mod                    ../starter/share/modules1/submodel_mod.F
!||    table_mod                       ../starter/share/modules1/table_mod.F
!||====================================================================
      SUBROUTINE st_qaprint_driver(IGEO        ,GEO        ,BUFGEO          ,IPM            ,PM               ,
     2                             BUFMAT      ,NOM_OPT    ,INOM_OPT        ,NUMLOADP       ,ILOADP           ,
     3                             LLOADP      ,LOADP      ,IBCL            ,FORC           ,IPRES            ,
     4                             PRES        ,NPBY       ,LPBY            ,RBY            ,IBCR             ,
     5                             FRADIA      ,IBCV       ,FCONV           ,IBFTEMP        ,FBFTEMP          ,
     6                             IGRV        ,LGRV       ,AGRV            ,IBFFLUX        ,FBFFLUX          ,
     7                             ITAB        ,V          ,VR              ,W              ,ICODE            ,
     8                             ISKEW       ,ICFIELD    ,LCFIELD         ,CFIELD         ,DAMPR            ,
     9                             TEMP        ,IBCSLAG    ,IPARI           ,INTBUF_TAB     ,CLUSTERS         ,
     A                             IBOX        ,IPMAS      ,IBFVEL          ,FBFVEL         ,NIMPACC          ,
     B                             LACCELM      ,ACCELM    ,NOM_SECT        ,NSTRF          ,SECBUF           ,
     C                             SKEW        ,ISKWN      ,XFRAME          ,T_MONVOL       ,T_MONVOL_METADATA,
     D                             I2RUPT      ,AREASL     ,INTBUF_FRIC_TAB ,NPFRICORTH     ,MAT_ELEM         ,
     E                             PFRICORTH   ,IREPFORTH  ,PHIFORTH        ,VFORTH         ,XREFC            ,
     F                             XREFTG      ,XREFS      ,TAGXREF         ,IXS            ,IXC              ,
     G                             IXTG        ,RWBUF      ,NPRW            ,LPRW           ,ITHVAR           ,
     H                             IPART       ,SUBSETS    ,IPARTTH         ,NTHGRPMX       ,NIMPDISP         ,
     M                             NIMPVEL     ,DETONATORS ,IBCSCYC         ,NPC            ,PLD              ,
     N                             TABLE       ,NPTS       ,IRBE3           ,LRBE3          ,FRBE3            ,
     O                             MGRBY       ,IXS10      ,ISOLNOD         ,IXR            ,R_SKEW           ,
     P                             IXP         ,IXT        ,X               ,THKE           ,SH4ANG           ,
     Q                             THKEC       ,SH3ANG     ,SET             ,LSUBMODEL      ,IGRNOD           ,
     R                             IGRPART     ,IGRBRIC    ,IGRSH4N         ,IGRSH3N        ,IGRQUAD          ,
     S                             IGRBEAM     ,IGRTRUSS   ,IGRSPRING       ,IGRSURF        ,IGRSLIN          ,
     T                             IXQ         ,ISPCOND    ,RTRANS          ,IRAND          ,ALEA             ,
     U                             XSEED       ,XLAS       ,LAS             ,IRBE2          ,LRBE2            ,
     V                             KXSP        ,IPARTSP    ,DRAPE           ,IXR_KJ         ,IACTIV           ,
     W                             FACTIV      ,UNITAB     ,NPBYL           ,LPBYL          ,RBYL             ,
     X                             XYZREF      ,SENSORS    ,FUNC2D          ,
     Y                             INICRACK    ,IPRELOAD   ,PRELOAD         ,IFLAG_BPRELOAD ,IBMPC,
     Z                             IBMPC2      ,IBMPC3     ,IBMPC4          ,RBMPC          ,LJOINT           ,
     A                             NNLINK      ,LNLINK     ,BUFSF           ,SBUFSF_        ,PM_STACK         ,
     B                             GEO_STACK   ,IGEO_STACK ,IPARG           ,IPADMESH       ,PADMESH          ,
     C                             LIFLOW      , LRFLOW    ,IFLOW           ,RFLOW          ,
     D                             SH4TREE     ,SH3TREE    ,SH4TRIM         ,SH3TRIM        ,QP_IPERTURB      ,
     E                             QP_RPERTURB ,LLINAL     ,LINALE          ,FVM_INIVEL     ,GJBUFI           ,
     F                             GJBUFR      ,MS         ,IN              ,LGAUGE         ,GAUGE            ,
     G                             KXX         ,IXX        ,IPARTX          ,IXRI           ,IXS16            ,
     H                             IEXMAD      ,FXBIPM     ,FXBFILE_TAB     ,EIGIPM         ,EIGRPM           ,
     I                             ISPHIO      ,VSPHIO     ,EBCS_TAB        ,INIMAP1D       ,INIMAP2D         ,
     J                             NSIGSH      ,SIGSH      ,NSIGI           ,SIGSP          ,NSIGS            ,
     K                             SIGI        ,NSIGBEAM   ,SIGBEAM         ,NSIGTRUSS      ,SIGTRUSS         ,
     L                             NSIGRS      ,SIGRS      ,MERGE_NODE_TAB  ,MERGE_NODE_TOL ,
     M                             IMERGE      ,NMERGE_TOT ,IEXLNK          ,DRAPEG         ,USER_WINDOWS     ,
     N                             OUTPUT      ,DEFAULTS   ,GLOB_THERM      ,PBLAST         ,IBEAM_VECTOR     ,
     O                             RBEAM_VECTOR,DAMP_RANGE_PART)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE qa_out_mod
      USE intbufdef_mod
      USE cluster_mod
      USE optiondef_mod
      USE monvol_struct_mod
      USE intbuf_fric_mod                     
      USE groupdef_mod      
      USE detonators_mod            
      USE table_mod
      USE setdef_mod
      USE submodel_mod
      USE drape_mod
      USE unitab_mod
      USE ale_ebcs_mod
      USE sensor_mod
      USE func2d_mod
      USE multi_fvm_mod
      USE ebcs_mod
      USE inimap1d_mod
      USE inimap2d_mod
      USE user_windows_mod
      USE output_mod
      USE mat_elem_mod
      USE defaults_mod
      USE names_and_titles_mod , ONLY : nchartitle
      use glob_therm_mod
      USE pblast_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      "com04_c.inc"
#include      "lagmult.inc"
#include      "param_c.inc"
#include      "scr03_c.inc"
#include      "scr17_c.inc"
#include      "scr23_c.inc"
#include      "tabsiz_c.inc"
#include      "sphcom.inc"
#include      "com_xfem1.inc"
#include      "fxbcom.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NOM_OPT(LNOPT1,SNOM_OPT1), INOM_OPT(SINOM_OPT)
      INTEGER, INTENT(IN) :: ITAB(NUMNOD),IXR_KJ(5,*)
      INTEGER, INTENT(IN) :: NNLINK(10,SNNLINK), LNLINK(SLNLINK)
      TYPE (CLUSTER_) ,DIMENSION(NCLUSTER) :: CLUSTERS
      INTEGER,INTENT(IN) :: NOM_SECT(SNOM_SECT),NSTRF(SNSTRF),IGEO_STACK(4* NPT_STACK+2,NS_STACK)
      my_real,INTENT(IN) :: secbuf(ssecbuf)    
      TYPE(monvol_struct_), DIMENSION(NVOLU), INTENT(IN) :: T_MONVOL
      TYPE(monvol_metadata_), INTENT(IN) :: T_MONVOL_METADATA
      TYPE(subset_), DIMENSION(NSUBS), INTENT(IN) :: SUBSETS
      TYPE(detonators_struct_) :: DETONATORS
      my_real, INTENT(IN) :: xlas(*)
      INTEGER, INTENT(IN) :: LAS(*)
      TYPE (DRAPE_)       :: DRAPE(NUMELC_DRAPE +  NUMELTG_DRAPE)
      TYPE (DRAPEG_)                         :: DRAPEG
      INTEGER, INTENT(IN) :: IACTIV(LACTIV,*)
      my_real, INTENT(IN) :: factiv(lractiv,*),geo_stack(6*npt_stack+1,ns_stack),pm_stack(20,ns_stack)
      TYPE (UNIT_TYPE_) ::UNITAB
      my_real, INTENT(IN) :: xyzref(sx)
      TYPE (SENSORS_) ,INTENT(IN) :: SENSORS
      TYPE(func2d_struct), DIMENSION(NFUNC2D), INTENT(IN) :: FUNC2D
      TYPE (INICRACK_) , DIMENSION(NINICRACK) :: INICRACK
      INTEGER, INTENT(IN) :: LLINAL
      INTEGER, DIMENSION(LLINAL), INTENT(IN) :: LINALE
      my_real, INTENT(IN) :: qp_rperturb(nperturb,4)
      INTEGER, INTENT(IN) :: QP_IPERTURB(NPERTURB,6)
      TYPE (FVM_INIVEL_STRUCT), DIMENSION(NINVEL), INTENT(IN) :: FVM_INIVEL
      INTEGER, INTENT(IN) :: LGAUGE(3,NBGAUGE)
      my_real, INTENT(IN) :: GAUGE(LLGAUGE,NBGAUGE)
      INTEGER, INTENT(IN) :: KXX(NIXX,*)
      INTEGER, INTENT(IN) :: IXX(*)
      INTEGER, INTENT(IN) :: IXRI(4,*)
      INTEGER, INTENT(IN) :: FXBIPM(NBIPM,NFXBODY),EIGIPM(*)
      my_real, INTENT(IN) :: EIGRPM(*)
      CHARACTER, DIMENSION(NFXBODY) :: FXBFILE_TAB*2148
      INTEGER ISPHIO(NISPHIO,NSPHIO)
      my_real
     .        VSPHIO(SVSPHIO)
      TYPE(t_ebcs_tab), INTENT(IN) :: EBCS_TAB
      TYPE(INIMAP1D_STRUCT), DIMENSION(NINIMAP1D), INTENT(IN) :: INIMAP1D
      TYPE(INIMAP2D_STRUCT), DIMENSION(NINIMAP2D), INTENT(IN) :: INIMAP2D
      INTEGER, INTENT(IN) :: NSIGSH,NSIGI,NSIGS,NSIGBEAM,NSIGTRUSS,NSIGRS
      my_real, INTENT(IN) :: SIGSH(MAX(1,NSIGSH),*),SIGSP(NSIGI,*),SIGI(NSIGS,*),
     .                       SIGBEAM(NSIGBEAM,*),SIGTRUSS(NSIGTRUSS,*),
     .                       SIGRS(NSIGRS,*)
      INTEGER, INTENT(IN) :: IMERGE(*),MERGE_NODE_TAB(*),NMERGE_TOT
      my_real, INTENT(IN) :: merge_node_tol(*)
      INTEGER,INTENT(IN) :: IPARG(NPARG,NGROUP)
      INTEGER, INTENT(IN) :: IBEAM_VECTOR(NUMELP)
      my_real, INTENT(IN) :: rbeam_vector(3,numelp)
      INTEGER ,INTENT(IN) :: DAMP_RANGE_PART(NPART)     !< flag to compute the damping range
      TYPE(user_windows_),INTENT(IN)      :: USER_WINDOWS
      TYPE(output_)      ,INTENT(IN)      :: OUTPUT
      TYPE(mat_elem_)    ,INTENT(IN)      :: MAT_ELEM
      TYPE(defaults_)    ,INTENT(IN)      :: DEFAULTS             !< Defaults mod
      TYPE(glob_therm_)  ,intent(in)      :: glob_therm
      TYPE(pblast_)      ,INTENT(IN)      :: PBLAST
C-----------------------------------------------
C      NOM_OPT(LNOPT1,SNOM_OPT1) 
C         * Possibly, NOM_OPT(1) = ID
C        NOM_OPT(LNOPT1-LTITL+1:LTITL) <=> TITLES of the OPTIONS
C--------------------------------------------------
C      SNOM_OPT1= NRBODY+NACCELM+NVOLU+NINTER+NINTSUB+
C     +           NRWALL+NJOINT+NSECT+NLINK+
C     +           NUMSKW+1+NUMFRAM+1+NFXBODY+NFLOW+NRBE2+
C     +           NRBE3+NSUBMOD+NFXVEL+NUMBCS+NUMMPC+
C     +           NGJOINT+NUNIT0+NFUNCT+NADMESH+
C     +           NSPHIO+NSPCOND+NRBYKIN+NEBCS+
C     +           NINICRACK+NODMAS+NBGAUGE+NCLUSTER+NINTERFRIC+
C     +           NRBMERGE
C-----------------------------------------------
C      INOM_OPT(SINOM_OPT)
C--------------------------------------------------
C      INOM_OPT(1) = NRBODY
C      INOM_OPT(2) = INOM_OPT(1) + NACCELM
C      INOM_OPT(3) = INOM_OPT(2) + NVOLU
C      INOM_OPT(4) = INOM_OPT(3) + NINTER
C      INOM_OPT(5) = INOM_OPT(4) + NINTSUB
C      INOM_OPT(6) = INOM_OPT(5) + NRWALL
C      INOM_OPT(7) = INOM_OPT(6) 
C      INOM_OPT(8) = INOM_OPT(7) + NJOINT
C      INOM_OPT(9) = INOM_OPT(8) + NSECT
C      INOM_OPT(10)= INOM_OPT(9) + NLINK
C      INOM_OPT(11)= INOM_OPT(10)+ NUMSKW+1+NUMFRAM+1+NSUBMOD
C      INOM_OPT(12)= INOM_OPT(11)+ NFXBODY
C      INOM_OPT(13)= INOM_OPT(12)+ NFLOW
C      INOM_OPT(14)= INOM_OPT(13)+ NRBE2
C      INOM_OPT(15)= INOM_OPT(14)+ NRBE3
C      INOM_OPT(16)= INOM_OPT(15)+ NFXVEL
C      INOM_OPT(17)= INOM_OPT(16)+ NUMBCS
C      INOM_OPT(18)= INOM_OPT(17)+ NUMMPC
C      INOM_OPT(19)= INOM_OPT(18)+ NGJOINT
C      INOM_OPT(20)= INOM_OPT(19)+ NUNIT0
C      INOM_OPT(21)= INOM_OPT(20)+ NFUNCT
C      INOM_OPT(22)= INOM_OPT(21)+ NADMESH
C      INOM_OPT(23)= INOM_OPT(22)+ NSPHIO
C      INOM_OPT(24)= INOM_OPT(23)+ NSPCOND
C      INOM_OPT(25)= INOM_OPT(24)+ NEBCS
C      INOM_OPT(26)= INOM_OPT(25)+ NINICRACK
C      INOM_OPT(27)= INOM_OPT(26)+ NODMAS
C      INOM_OPT(28)= INOM_OPT(27)+ NBGAUGE
C      INOM_OPT(29)= INOM_OPT(28)+ NCLUSTER
C      INOM_OPT(30)= INOM_OPT(29)+ NINTERFRIC
C      INOM_OPT(31)= INOM_OPT(30)+ NRBMERGE
C      .. TO BE MAINTAINED (cf doc/inom_opt.txt) ..
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NIMPDISP,NIMPVEL,NIMPACC
      INTEGER, INTENT(IN) :: IGEO(NPROPGI,NUMGEO),IPM(NPROPMI,NUMMAT)
      INTEGER, INTENT(IN) :: NUMLOADP, ILOADP(SIZLOADP,NLOADP), LLOADP(NUMLOADP)
      INTEGER, INTENT(IN) :: ICFIELD(SIZFIELD,NLOADC), LCFIELD(SLCFIELD)
      INTEGER, INTENT(IN) :: IBCL(NIBCLD,NCONLD-NPRELD), IPRES(NIBCLD,NPRELD)
      INTEGER, INTENT(IN) :: IGRV(NIGRV,NGRAV), LGRV(*)
      INTEGER, INTENT(IN) :: NPBY(NNPBY,NRBYKIN), NPBYL(NNPBY,NRBYLAG), LPBY(*), LPBYL(*)
      INTEGER, INTENT(IN) :: IBCR(GLOB_THERM%NIRADIA,GLOB_THERM%NUMRADIA)
      INTEGER, INTENT(IN) :: IBCV(GLOB_THERM%NICONV,GLOB_THERM%NUMCONV)
      INTEGER, INTENT(IN) :: IBFTEMP(GLOB_THERM%NIFT,GLOB_THERM%NFXTEMP)
      INTEGER, INTENT(IN) :: IBFFLUX(GLOB_THERM%NITFLUX,GLOB_THERM%NFXFLUX)
      INTEGER, INTENT(IN) :: IBFVEL(NIFV,NFXVEL)
      INTEGER, INTENT(IN) :: ICODE(NUMNOD), ISKEW(NUMNOD)
      INTEGER, INTENT(IN) :: IBCSLAG(5,NBCSLAG)
      INTEGER, INTENT(IN) :: IPARI(NPARI,NINTER)
      INTEGER, INTENT(IN) :: LACCELM(3, NACCELM)
      INTEGER, INTENT(IN) :: ISKWN(LISKN,*)
      INTEGER, INTENT(IN) :: NPFRICORTH , PFRICORTH(*) , IREPFORTH(*)
      INTEGER, INTENT(IN) :: TAGXREF(NUMNOD),IXC(NIXC,*),IXTG(NIXTG,*),IXS(NIXS,*), IXS10(6,*), 
     .                       IXS16(8,*)
      INTEGER, INTENT(IN) :: ISOLNOD(NUMELS),IXR(NIXR,*), R_SKEW(*),IXP(NIXP,*),IXT(NIXT,*),
     .                       IXQ(NIXQ,*),KXSP(NISP,*),IPARTSP(*),IPARTX(*)
      INTEGER, INTENT(IN) :: NPRW(NRWALL,NNPRW),LPRW(SLPRW)
      INTEGER, INTENT(IN) :: ITHVAR(SITHVAR),
     .                       IPART(LIPART1*(NPART+NTHPART)),
     .                       IPARTTH(18*(NPART+NTHPART)),NTHGRPMX,IBCSCYC(*)
      INTEGER, INTENT(IN) :: NPTS,NPC(*),IRBE3(*),LRBE3(*),IRBE2(*),LRBE2(*)  
      INTEGER, INTENT(IN) :: MGRBY(NMGRBY,SMGRBY)
      INTEGER, INTENT(IN) :: ISPCOND(NISPCOND,*),LJOINT(*),GJBUFI(LKJNI,*)
      INTEGER, INTENT(IN) :: IRAND(*)
      INTEGER, INTENT(IN) :: IPRELOAD(3,*), IFLAG_BPRELOAD(*)
      INTEGER, INTENT(IN) :: IBMPC(NUMMPC),IBMPC2(LMPC),IBMPC3(LMPC),IBMPC4(LMPC)
      INTEGER, INTENT(IN) :: IPADMESH(KIPADMESH,*)
      INTEGER, INTENT(IN) :: SH4TREE(KSH4TREE,*), SH3TREE(KSH3TREE,*),
     .                       SH4TRIM(*),SH3TRIM(*) 
      INTEGER, INTENT(IN) :: IEXMAD(*),IEXLNK(NR2R,*)
 
      my_real, INTENT(IN) ::
     .                        GEO(NPROPG,NUMGEO), BUFGEO(*), PM(NPROPM,NUMMAT), BUFMAT(*)
      my_real, INTENT(IN) ::
     .                        LOADP(LFACLOAD,NLOADP), CFIELD(LFACLOAD,NLOADC),
     .                        FORC(LFACCLD,NCONLD-NPRELD), PRES(LFACCLD,NPRELD),
     .                        AGRV(LFACGRV,NGRAV),PRELOAD(6,*)  
      my_real, INTENT(IN) :: RBY(NRBY,NRBYKIN),RBYL(NRBY,NRBYLAG),FRBE3(SFRBE3)
      my_real, INTENT(IN) :: GJBUFR(LKJNR,*),MS(*),IN(*)
      my_real, INTENT(IN) :: fradia(glob_therm%LFACTHER,glob_therm%NUMRADIA)
      my_real, INTENT(IN) :: fconv(glob_therm%LFACTHER,glob_therm%NUMCONV)
      my_real, INTENT(IN) :: fbftemp(glob_therm%LFACTHER,glob_therm%NFXTEMP)
      my_real, INTENT(IN) :: fbfflux(glob_therm%LFACTHER,glob_therm%NFXFLUX)
      my_real, INTENT(IN) :: fbfvel(lfxvelr,nfxvel)
      my_real, INTENT(IN) :: v(3,numnod), vr(svr) , ! SVR=3*NUMNOD*IRODDL
     .                        w(sw)                  ! SW=3*NUMNOD*IALE
      my_real, INTENT(IN) ::
     .                        dampr(nrdamp,*)
      my_real, INTENT(IN) ::
     .                        temp(numnod)
      my_real, INTENT(IN) ::
     .                        accelm(llaccelm, naccelm)
      my_real, INTENT(IN) ::
     .                        skew(lskew,*)
      my_real, INTENT(IN) ::
     .                        xframe(nxframe,*)
      my_real, INTENT(IN) ::
     .                        i2rupt(6,*),areasl(*)
      my_real, INTENT(IN) ::   phiforth(*), vforth(3,*)
      my_real, INTENT(IN) ::
     .                        xrefc(4,3,numelc),xreftg(3,3,numeltg),xrefs(8,3,numels8)
      my_real, INTENT(IN) ::
     .                        rwbuf(nrwlp,nrwall)
      my_real, INTENT(IN) ::
     .                        pld(*)
      my_real, INTENT(IN) ::
     .                        x(3,numnod),
     .                        thke(*),thkec(*),sh4ang(*),sh3ang(*)
      my_real, INTENT(IN) ::
     .                        rtrans(ntransf,*)
      my_real, INTENT(IN) ::
     .                        alea(*),xseed(*)
      my_real, INTENT(IN) ::
     .                        rbmpc(srbmpc)
 
      my_real, INTENT(IN)  ::  bufsf(sbufsf)
      INTEGER,INTENT(IN) :: SBUFSF_
 
      my_real, INTENT(IN) :: padmesh(kpadmesh,*)
     
      TYPE(intbuf_struct_) INTBUF_TAB(*)
      TYPE (BOX_)    ,DIMENSION(NBBOX)  ,INTENT(IN) :: IBOX
      TYPE (ADMAS_)  ,DIMENSION(NODMAS) ,INTENT(IN) :: IPMAS
      TYPE(intbuf_fric_struct_) INTBUF_FRIC_TAB(*)
      TYPE(ttable) TABLE(*)
!
      TYPE (SET_)    , DIMENSION(NSETS) :: SET
      TYPE(submodel_data),INTENT(IN)::LSUBMODEL(*)
      TYPE (GROUP_)  , INTENT(IN), DIMENSION(NGRNOD)  :: IGRNOD
      TYPE (GROUP_)  , INTENT(IN), DIMENSION(NGRPART) :: IGRPART
      TYPE (GROUP_)  , INTENT(IN), DIMENSION(NGRBRIC) :: IGRBRIC
      TYPE (GROUP_)  , INTENT(IN), DIMENSION(NGRSHEL) :: IGRSH4N
      TYPE (GROUP_)  , INTENT(IN), DIMENSION(NGRSH3N) :: IGRSH3N
      TYPE (GROUP_)  , INTENT(IN), DIMENSION(NGRQUAD) :: IGRQUAD
      TYPE (GROUP_)  , INTENT(IN), DIMENSION(NGRBEAM) :: IGRBEAM
      TYPE (GROUP_)  , INTENT(IN), DIMENSION(NGRTRUS) :: IGRTRUSS
      TYPE (GROUP_)  , INTENT(IN), DIMENSION(NGRSPRI) :: IGRSPRING
      TYPE (SURF_)   , INTENT(IN), DIMENSION(NSURF)   :: IGRSURF
      TYPE (SURF_)   , INTENT(IN), DIMENSION(NSLIN)   :: IGRSLIN
      INTEGER, INTENT(IN) :: LIFLOW, LRFLOW
      INTEGER, DIMENSION(LIFLOW), INTENT(IN) :: IFLOW
      my_real, DIMENSION(LRFLOW), INTENT(IN) :: rflow
C--------------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
 
      ! doqa is set in QA_OUT_MOD depending on DO_QA env variable setting
      ! we skip the specific QAPRINT treatment if doqa is != 1
      
      IF ( doqa /= 1  ) RETURN
 
      CALL st_qaprint_materials(mat_elem   ,ipm       ,pm        ,bufmat    )
 
      CALL st_qaprint_properties(igeo      ,geo       ,bufgeo    ,
     .                           pm_stack  ,geo_stack ,igeo_stack)
 
      CALL st_qaprint_initial_conditions(
     1     nom_opt   ,inom_opt  ,itab      ,v         ,vr        ,
     2     w         ,temp      ,inicrack  ,fvm_inivel, 
     3     inimap1d, inimap2d)
 
      CALL st_qaprint_constraints(nom_opt   ,inom_opt  ,npby      ,lpby      ,rby       ,
     2                            ibftemp   ,fbftemp   ,ibfflux   ,fbfflux   ,itab      ,
     3                            icode     ,iskew     ,ibcslag   ,ibfvel    ,fbfvel    ,
     4                            nimpdisp  ,nimpvel   ,nimpacc   ,rwbuf     ,nprw      ,
     5                            lprw      ,ibcscyc   ,irbe3     ,lrbe3     ,frbe3     ,
     6                            mgrby     ,ispcond   ,irbe2     ,lrbe2     ,npbyl     ,
     7                            lpbyl     ,rbyl      ,ibmpc     ,ibmpc2    ,ibmpc3    ,
     8                            ibmpc4    ,rbmpc     ,ljoint    ,nnlink    ,lnlink,
     9                            llinal    ,linale    ,gjbufi    ,gjbufr    ,ms        ,
     9                            in        ,fxbipm    ,fxbfile_tab,glob_therm)
 
      CALL st_qaprint_loads(nom_opt   ,inom_opt  ,numloadp  ,iloadp    ,lloadp    ,
     2                      loadp     ,ibcl      ,forc      ,ipres     ,pres      ,
     3                      ibcr      ,fradia    ,ibcv      ,fconv     ,igrv      ,
     4                      lgrv      ,agrv      ,icfield   ,lcfield   ,cfield    ,
     5                      ipreload  ,preload   ,iflag_bpreload,
     6                      liflow, lrflow, iflow,rflow     ,isphio    ,vsphio,
     7                      glob_therm  ,pblast  )
 
      CALL st_qaprint_general_controls(nom_opt   ,inom_opt  ,dampr, irand, alea, 
     1                                 xseed     ,unitab    ,qp_iperturb       ,
     2                                 qp_rperturb,eigipm   , eigrpm, defaults ,
     3                                 damp_range_part)
      CALL st_qaprint_interfaces(nom_opt   ,inom_opt  ,ipari  ,intbuf_tab,i2rupt ,  
     1                           areasl,glob_therm%INTHEAT)      
      CALL st_qaprint_clusters(nom_opt   ,inom_opt  ,clusters)
      CALL st_qaprint_output_databases(nom_opt, inom_opt, laccelm ,accelm,lgauge,gauge)
      CALL st_qaprint_model_tools(nom_opt, inom_opt, ibox, ipmas, nom_sect, nstrf, secbuf ,
     .                            skew   , iskwn   , xframe, npc,pld,table,npts, 
     .                            iactiv , factiv  ,sensors,func2d)
      CALL st_qaprint_monvol(t_monvol, t_monvol_metadata)
      CALL st_qaprint_friction(nom_opt,inom_opt,intbuf_fric_tab,npfricorth,pfricorth ,irepforth,phiforth,vforth) 
      CALL st_qaprint_reference_state(xrefc     ,xreftg    ,xrefs     ,tagxref   ,
     1                                ixs       ,ixc       ,ixtg      )
 
      CALL st_qaprint_time_histories(output%TH   ,ithvar   ,ipart    ,subsets  ,
     1                               ipartth  ,nthgrpmx )
 
      CALL st_qaprint_dfs_detonators(detonators)    
 
      CALL st_qaprint_dfs_lasers(xlas,las)   
 
      CALL st_qaprint_nodes(itab ,x,imerge,merge_node_tol, merge_node_tab,nmerge_tot,ms)
 
      CALL st_qaprint_element(ixs ,ixs10 ,ipm ,igeo ,itab ,
     1                        isolnod,ixr,r_skew,iskwn,ixp,
     2                        ixt    ,x  ,ixc   ,ixtg ,thke,
     3                        sh4ang ,thkec,sh3ang,kxsp,ipartsp,
     4                        ipart  ,ixr_kj,kxx, ixx ,ipartx,
     5                        ixri   ,ixs16 ,ixq, ibeam_vector,rbeam_vector)
      
!
      CALL st_qaprint_set(set    ,lsubmodel,itab     ,igrnod ,igrpart,
     .                    ipart  ,igrbric  ,igrsh4n  ,igrsh3n,igrquad,
     .                    igrbeam,igrtruss ,igrspring,igrsurf,igrslin,
     .                    ixc    ,ixtg     ,ixq      ,ixp    ,ixt    ,
     .                    ixr    ,ixs      )
      CALL st_qaprint_transformations(rtrans,lsubmodel,igrnod)
      CALL st_qaprint_composite_options(drape,drapeg)
      CALL st_qaprint_ale_options_driver
      CALL st_qaprint_inivol()  
      CALL st_qaprint_refsta(xrefc     ,xreftg  ,xrefs    ,tagxref   ,ixs   ,ixc   ,ixtg   ,itab, xyzref)   
      CALL st_qaprint_ebcs(ebcs_tab)  
      CALL st_qaprint_groups(igrnod ,igrpart ,igrbric  ,igrsh4n   ,igrsh3n   ,
     .                       igrquad,igrbeam ,igrtruss ,igrspring)  
      CALL st_qaprint_surf(igrsurf, igrslin, bufsf, sbufsf_)  
      CALL st_qaprint_admesh(ipart ,ipadmesh,padmesh,sh4tree,sh3tree,
     .                        ixc  ,ixtg    ,sh4trim,sh3trim    )
      CALL st_qaprint_madymo(iexmad ,itab   ,ipart  ,ixs    ,ixc    ,ixtg   )  
      CALL st_qaprint_seatbelts(itab)
      CALL st_qaprint_initial_state(
     .                    nsigsh   ,sigsh    ,nsigi   ,sigsp    ,nsigs   ,
     .                    sigi     ,nsigbeam ,sigbeam ,nsigtruss,sigtruss,
     .                    nsigrs   ,sigrs    )
      CALL st_qaprint_multidomains(ipart,iexlnk,igrnod)
      CALL st_qaprint_userwi(user_windows)
      CALL st_qaprint_internal_groups(iparg)
C-----------------------------------------------------------------------
 
      RETURN
      END