hm_preread_impdisp.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "scr17_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	hm_preread_impdisp (nimpdisp, igrnod, ipart, ipartr, unitab, lsubmodel)

Function/Subroutine Documentation

hm_preread_impdisp()

subroutine hm_preread_impdisp	(	integer, intent(out)	nimpdisp,
		type (group_), dimension(ngrnod), intent(in)	igrnod,
		integer, dimension(lipart1,*)	ipart,
		integer, dimension(*)	ipartr,
		type(unit_type_), intent(in)	unitab,
		type(submodel_data), dimension(*), intent(in)	lsubmodel )

Definition at line 39 of file hm_preread_impdisp.F.

C============================================================================
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE groupdef_mod
      USE submodel_mod
      USE hm_option_read_mod
      USE unitab_mod
      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      "com04_c.inc"
#include      "scr17_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ,INTENT(OUT)   :: NIMPDISP
      INTEGER IPART(LIPART1,*), IPARTR(*)
      TYPE(UNIT_TYPE_) ,INTENT(IN) :: UNITAB
      TYPE (GROUP_)      , DIMENSION(NGRNOD)  ,INTENT(IN) :: IGRNOD
      TYPE(SUBMODEL_DATA), DIMENSION(*)       ,INTENT(IN) :: LSUBMODEL
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I,N,OPTID,NFDISP,NFGEO,IGS,GRNOD_ID,PART_ID,NNOD,JPART
      CHARACTER(LEN=NCHARKEY)   :: KEY
      CHARACTER(LEN=NCHARTITLE) :: TITR
      LOGICAL :: IS_AVAILABLE
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER  NGR2USRN
      EXTERNAL ngr2usrn
c--------------------------------------------------
c     Count number of nodes with imposed displacement => NIMPDISP
C======================================================================|
      is_available = .false.
c
      nimpdisp = 0
c--------------------------------------------------
c     /IMPDISP
c--------------------------------------------------
      CALL hm_option_count('/IMPDISP',nfdisp)
c
      CALL hm_option_start('/IMPDISP')
c
      DO i=1,nfdisp
c
        CALL hm_option_read_key(lsubmodel,
     .                          option_id   = optid,
     .                          option_titr = titr,
     .                          keyword2    = key)
c
        IF (key(1:4) /= 'FGEO') THEN
          CALL hm_get_intv('entityid' ,grnod_id ,is_available,lsubmodel)
          igs = ngr2usrn(grnod_id,igrnod,ngrnod,nnod)
          IF (igs > 0) nimpdisp = nimpdisp + nnod                                         
        ENDIF
c
      ENDDO ! DO I=1,NFDISP
c--------------------------      
c     /IMPDISP/FGEO    
c--------------------------         
      CALL hm_option_count('/IMPDISP/FGEO',nfgeo)
c
      CALL hm_option_start('/IMPDISP/FGEO')
c     
      DO i=1,nfgeo
c
        CALL hm_option_read_key(lsubmodel,
     .                          option_id   = optid,
     .                          option_titr = titr,
     .                          keyword2    = key)
        IF (key(1:4) == 'FGEO') THEN
          CALL hm_get_intv('rad_spring_part' ,part_id ,is_available,lsubmodel)
          IF (part_id > 0) THEN
            jpart = 0
            DO n=1,npart
              IF (ipart(4,n) == part_id) jpart = n
            ENDDO
            IF (jpart == 0) THEN
             CALL ancmsg(msgid=1077, msgtype=msgerror,
     .                   anmode=aninfo,
     .                   i1=optid,
     .                   c1=titr,
     .                   i2=part_id)
            ENDIF
            DO n=1,numelr
              IF (ipartr(n) == jpart) nimpdisp = nimpdisp + 1
            ENDDO
          ENDIF
c
          CALL hm_get_intv('distribution_table_count' ,nnod ,is_available,lsubmodel)
c
          nimpdisp = nimpdisp + nnod
        ENDIF
      ENDDO ! DO I=1,NFGEO
c-----------
      RETURN