hm_read_inter_lagmul.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/.
!||====================================================================
!||    hm_read_inter_lagmul          ../starter/source/interfaces/reader/hm_read_inter_lagmul.F
!||--- called by ------------------------------------------------------
!||    hm_read_interfaces            ../starter/source/interfaces/reader/hm_read_interfaces.F
!||--- calls      -----------------------------------------------------
!||    ancmsg                        ../starter/source/output/message/message.F
!||    fretitl                       ../starter/source/starter/freform.F
!||    hm_read_inter_lagmul_type02   ../starter/source/interfaces/int02/hm_read_inter_lagmul_type02.F
!||    hm_read_inter_lagmul_type07   ../starter/source/interfaces/int07/hm_read_inter_lagmul_type07.F
!||    hm_read_inter_lagmul_type16   ../starter/source/interfaces/int16/hm_read_inter_lagmul_type16.F
!||    hm_read_inter_lagmul_type17   ../starter/source/interfaces/int17/hm_read_inter_lagmul_type17.f
!||--- uses       -----------------------------------------------------
!||    message_mod                   ../starter/share/message_module/message_mod.F
!||    submodel_mod                  ../starter/share/modules1/submodel_mod.F
!||====================================================================
      SUBROUTINE hm_read_inter_lagmul(
     1        IPARI    ,STFAC      ,FRIGAP    ,IGRNOD      ,IGRSURF    ,
     2        ILAGM    ,NI         ,NOM_OPT   ,TITR        ,DEF_INTER  ,
     3        IGRBRIC  ,UNITAB     ,LSUBMODEL ,NOINT       ,KEY        ,
     4        NPARI    ,NPARIR)
C============================================================================
C
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE unitab_mod
      USE message_mod
      USE groupdef_mod
      USE submodel_mod , ONLY : submodel_data, nsubmod
      USE names_and_titles_mod , ONLY : nchartitle, ncharkey
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      "scr17_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER,INTENT(IN) :: NPARI, NPARIR !< array sizes IPARI and FRIGAP
      INTEGER NOM_OPT(LNOPT1,*)
      INTEGER ISU1,ISU2,ILAGM,NI,NOINT
      INTEGER IPARI(NPARI),DEF_INTER(100)
      my_real stfac
      my_real frigap(nparir)
      CHARACTER(LEN=NCHARTITLE)::TITR
      CHARACTER(LEN=NCHARKEY)::KEY
C-----------------------------------------------
      TYPE (GROUP_)  ,TARGET, DIMENSION(NGRNOD)  :: IGRNOD
      TYPE (GROUP_)  ,TARGET, DIMENSION(NGRBRIC) :: IGRBRIC
      TYPE (SURF_)   ,TARGET , DIMENSION(NSURF)   :: IGRSURF
      TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB 
      TYPE(submodel_data) LSUBMODEL(NSUBMOD)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "units_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER UID,NTYP,USR_TYP
C=======================================================================
C     READING LAGRANGE MULTIPLIER INTERFACES (/INTER/LAGMUL/TYPE2, 
C         /INTER/LAGMUL/TYPE7, /INTER/LAGMUL/TYPE16,/INTER/LAGMUL/TYPE17)
C========================================================================
C
      ilagm = 1
      ipari(33) = ilagm
C
      nom_opt(1,ni)=noint
C
      WRITE(iout,1000) noint,trim(titr)
C---------------------------------------------------------------
      CALL fretitl(titr,nom_opt(lnopt1-ltitr+1,ni),ltitr)
C---------------------------------------------------------------
C LAGMUL INTERFACES READING ROUTINES
C-----------------------------
      ntyp = 0
      SELECT CASE(key(1:len_trim(key)))
c------------------------------
         CASE ('TYPE2')
c------------------------------
           ntyp = 2
           CALL hm_read_inter_lagmul_type02(
     1        ipari      ,frigap     ,noint      ,
     2        igrnod     ,igrsurf    ,def_inter  ,titr      ,unitab     ,
     3        lsubmodel  ,npari      ,nparir)
c------------------------------
         CASE ('TYPE7')
c------------------------------
           ntyp = 7
           CALL hm_read_inter_lagmul_type07(
     1        ipari      ,stfac      ,frigap     ,noint     ,
     2        igrnod     ,igrsurf    , titr       ,unitab     ,lsubmodel  )
c------------------------------
         CASE ('TYPE16')
c------------------------------
           ntyp = 16 
           CALL hm_read_inter_lagmul_type16(
     .          ipari      ,stfac      ,frigap     ,noint     ,
     .          igrnod     ,igrbric    ,lsubmodel )
c------------------------------
         CASE ('TYPE17')
c------------------------------
           ntyp = 17
           CALL hm_read_inter_lagmul_type17(
     .          ipari      ,stfac      ,frigap    ,noint     ,
     .          igrbric    ,lsubmodel )
C---------------------------------------------------------------
      END SELECT
C---------------------------------------------------------------
      usr_typ = ntyp
      IF (usr_typ/=2.AND.usr_typ/=7.AND.usr_typ/=16.AND.usr_typ/=17) THEN
         CALL ancmsg(msgid=486,msgtype=msgerror,anmode=aninfo,
     .               i1=noint,
     .               c1=titr,
     .               i2=ntyp)
      ENDIF
C
      IF (stfac == zero ) stfac = one_fifth
C  
C--------------------------------------------------------------
 1000 FORMAT(/1x,'  INTERFACE NUMBER :',i10,1x,a)
C------------
      RETURN
      END