hm_read_inter_lagmul_type16.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "scr06_c.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	hm_read_inter_lagmul_type16 (ipari, stfac, frigap, noint, igrnod, igrbric, lsubmodel)

Function/Subroutine Documentation

hm_read_inter_lagmul_type16()

subroutine hm_read_inter_lagmul_type16	(	integer, dimension(*)	ipari,
			stfac,
			frigap,
		integer	noint,
		type (group_), dimension(ngrnod), target	igrnod,
		type (group_), dimension(ngrbric), target	igrbric,
		type(submodel_data), dimension(*)	lsubmodel )

Definition at line 35 of file hm_read_inter_lagmul_type16.F.

C============================================================================
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE groupdef_mod
      USE submodel_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ISU1,ISU2,NOINT
      INTEGER IPARI(*)
      my_real frigap(*) , stfac
      TYPE (GROUP_)  ,TARGET ,DIMENSION(NGRNOD)  :: IGRNOD
      TYPE (GROUP_)  ,TARGET, DIMENSION(NGRBRIC) :: IGRBRIC
      TYPE(SUBMODEL_DATA) LSUBMODEL(*)
C----------------------s-------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "scr06_c.inc"
#include      "units_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: ISU10,ISU20,IDUM,MULTIMP,ITIED,NTYP,IS1,IS2
      my_real :: startt,bumult,stopt,gap
      INTEGER, DIMENSION(:), POINTER :: INGR2USR
      LOGICAL IS_AVAILABLE
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER NGR2USR
C=======================================================================
C     reading /INTER/LAGMUL/TYPE16
C=======================================================================
 
C Initializations
      ntyp = 16
      startt = zero
      stopt  = ep20
      stfac  = one_fifth
      gap    = zero
c------------------------------------------------------------
c  Card1
      CALL hm_get_intv('secondaryentityids ',isu10,is_available,lsubmodel)
      CALL hm_get_intv('mainentityids',isu20,is_available,lsubmodel)
c
c  Card2
      CALL hm_get_intv('Itied' ,itied,is_available,lsubmodel)
c------------------------------------------------------------
c     CHECKS
c------------------------------------------------------------
      is1 = 0
      is2 = 5
      ingr2usr => igrnod(1:ngrnod)%ID
      isu1 = ngr2usr(isu10,ingr2usr,ngrnod)
      ingr2usr => igrbric(1:ngrbric)%ID
      isu2 = ngr2usr(isu20,ingr2usr,ngrbric)
c
      bumult = bmul0
      multimp=12
c
      IF (nspmd > 1) THEN
          CALL ancmsg(msgid=755,msgtype=msgerror,anmode=aninfo,c1='TYPE 16 LAGRANGE INTERFACE')
      END IF
c------------------------------------------------------------
c     STORAGE
c------------------------------------------------------------
      ipari(7)  = ntyp
      ipari(13) = is1*10+is2
      ipari(15) = noint
      ipari(23) = multimp
      ipari(30) = itied
      ipari(45) = isu1
      ipari(46) = isu2
C
      frigap(2) = gap
      frigap(3) = startt
      frigap(4) = bumult
      frigap(5) = one
      frigap(11)= stopt
C------------------------------------------------------------
C     PRINTOUT
C------------------------------------------------------------
      WRITE(iout,1000)isu10,isu20,itied,startt,stopt 
C 
      IF(is1 == 0)THEN
        WRITE(iout,'(6X,A)')'NO SECONDARY SURFACE INPUT'
      ELSEIF(is1 == 1)THEN
        WRITE(iout,'(6X,A)')'SECONDARY SURFACE INPUT BY SEGMENTS'
      ELSEIF(is1 == 2)THEN
        WRITE(iout,'(6X,A)')'SECONDARY SURFACE INPUT BY NODES'
      ELSEIF(is1 == 3)THEN
        WRITE(iout,'(6X,A)')'SECONDARY SURFACE INPUT BY SEGMENTS'
      ELSEIF(is1 == 4 )THEN
        WRITE(iout,'(6X,A)')'SECONDARY SIDE INPUT BY BRICKS'
      ELSEIF(is1 == 5 )THEN
        WRITE(iout,'(6X,A)')'SECONDARY SIDE INPUT BY SOLID ELEMENTS'        
      ENDIF
      IF(is2 == 0)THEN
        WRITE(iout,'(6X,A)')'NO MAIN SURFACE INPUT'
      ELSEIF(is2 == 1)THEN
        WRITE(iout,'(6X,A)')'MAIN SURFACE INPUT BY SEGMENTS'
      ELSEIF(is2 == 2)THEN
        WRITE(iout,'(6X,A)')'MAIN SURFACE INPUT BY NODES'
      ELSEIF(is2 == 3)THEN
        WRITE(iout,'(6X,A)')'MAIN SURFACE INPUT BY SEGMENTS'
      ELSEIF(is2 == 4)THEN
        WRITE(iout,'(6X,A)')'MAIN SURFACE REFERS ',
     .                      'TO HYPER-ELLIPSOIDAL SURFACE'
      ENDIF
      WRITE(iout,'(A)')' '
c-----------
      RETURN
c--------------------------------------------------------------
 1000 FORMAT(//
     .    '    TYPE 16   NODES TO 16 NODES THICK SHELL  ' //,
     .    '    NODE GROUP . . . . . . .  . . . . . . . . . ',i10/,
     .    '    SOLID ELEMENT GROUP. . . . . . .  . . . . . ',i10/,
     .    '    ITIED . . . . . . . . . . . . . . . . . . . ',i10/,
     .    '       0: SLIDING '/,
     .    '       1: TIED (during contact)'/,
     .    '       2: TIED (no rebound)'/,
     .    '    START TIME. . . . . . . . . . . . . . . . . ',1pg20.13/,
     .    '    STOP TIME . . . . . . . . . . . . . . . . . ',1pg20.13/)
c--------------------------------------------------------------