hm_read_inter_type15.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	hm_read_inter_type15 (ipari, stfac, frigap, noint, igrsurf, titr, lsubmodel, unitab)

Function/Subroutine Documentation

hm_read_inter_type15()

subroutine hm_read_inter_type15	(	integer, dimension(*)	ipari,
			stfac,
			frigap,
		integer	noint,
		type (surf_), dimension(nsurf), target	igrsurf,
		character(len=nchartitle), intent(in)	titr,
		type(submodel_data), dimension(*), intent(in)	lsubmodel,
		type (unit_type_), intent(in)	unitab )

Definition at line 36 of file hm_read_inter_type15.F.

C============================================================================
C
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE groupdef_mod
      USE submodel_mod
      USE unitab_mod
      USE names_and_titles_mod , ONLY : nchartitle
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-----------------------------------------------
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 stfac
      my_real frigap(*)
      CHARACTER(LEN=NCHARTITLE),INTENT(IN) :: TITR
      TYPE(SUBMODEL_DATA),INTENT(IN)::LSUBMODEL(*)
      TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
C-----------------------------------------------
      TYPE (SURF_)   ,TARGET , DIMENSION(NSURF)   :: IGRSURF
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-----------------------------------------------
      CHARACTER(LEN=NCHARTITLE) :: TITR1
      INTEGER I,J,L, NTYP,IS1, IS2,ISU20,INTKG
      my_real fric,gap,startt,stopt,visc
      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=======================================================================
C     READING PENALTY INTERFACE /INTER/TYPE15
C=======================================================================
 
C Initializations
      is1=0
      is2=0
      intkg = 0
C
      fric = zero
      gap = zero
      startt = zero
      stopt=ep30
      visc = zero
C
      ntyp = 15
      ipari(15)=noint
      ipari(7)=ntyp
 
      is_available = .false.
C--------------------------------------------------
C EXTRACT DATAS (INTEGER VALUES)
C--------------------------------------------------
      CALL hm_get_intv('secondaryentityids',isu1,is_available,lsubmodel)
      CALL hm_get_intv('mainentityids',isu2,is_available,lsubmodel)
C--------------------------------------------------
C EXTRACT DATAS (REAL VALUES)
C--------------------------------------------------
      CALL hm_get_floatv('STIFF1',stfac,is_available,lsubmodel,unitab)
      CALL hm_get_floatv('FRIC',fric,is_available,lsubmodel,unitab)
 
C....* CHECKS *.............
 
      is1=1
      is2=4
      ingr2usr => igrsurf(1:nsurf)%ID
      IF(isu1/=0)isu1=ngr2usr(isu1,ingr2usr,nsurf)
        isu20=isu2
        isu2=ngr2usr(isu2,ingr2usr,nsurf)
        IF ( igrsurf(isu2)%TYPE/=100
     .    .AND.igrsurf(isu2)%TYPE/=101) THEN
           titr1 = igrsurf(isu20)%TITLE
             CALL ancmsg(msgid=111,
     .                   msgtype=msgerror,
     .                   anmode=aninfo,
     .                   i1=noint,
     .                   c1=titr,
     .                   i2=isu20,
     .                   c2=titr1)
      END IF
 
C.......* Storage IPARI FRIGAP *........
      ipari(45)=isu1
      ipari(46)=isu2
      ipari(13)=is1*10+is2
C
C.....* CHECKS *.....
C
      visc  =zero
      startt=zero
      stopt =ep30
 
C.....* Storage IPARI FRIGAP *.......
      frigap(1)=fric
      frigap(3)=startt
      frigap(11)=stopt
      frigap(14)=visc
C
      ipari(65) = intkg
      frigap(2)=gap
 
C
C------------------------------------------------------------
C     PRINTOUT
C------------------------------------------------------------
C
      WRITE(iout,1615)stfac,fric,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
C
C--------------------------------------------------------------
 1000 FORMAT(/1X,'  INTERFACE number :',I10,1X,A)
C------------
      RETURN
 
 
 1615 FORMAT(//
     .    '    type==15   elements to hyper-ellipsoid    ' //,
     .    '    INTERFACE stiffness factor. . . . . . . . . ',1PG20.13/,
     .    '    friction coefficient  . . . . . . . . . . . ',1PG20.13/,
     .    '    start time. . . . . . . . . . . . . . . . . ',1PG20.13/,
     .    '    stop time . . . . . . . . . . . . . . . . . ',1PG20.13/)