hm_preread_impvel.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_impvel (nimpvel, igrnod, ipart, ipartr, nfvlag, unitab, lsubmodel)

Function/Subroutine Documentation

hm_preread_impvel()

subroutine hm_preread_impvel	(	integer, intent(out)	nimpvel,
		type (group_), dimension(ngrnod), intent(in)	igrnod,
		integer, dimension(lipart1,*)	ipart,
		integer, dimension(*)	ipartr,
		integer, intent(inout)	nfvlag,
		type(unit_type_), intent(in)	unitab,
		type(submodel_data), dimension(*), intent(in)	lsubmodel )

Definition at line 39 of file hm_preread_impvel.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)   :: NIMPVEL
      INTEGER ,INTENT(INOUT) :: NFVLAG
      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,NFVEL,NFGEO,IGS,GRNOD_ID,PART_ID,NNOD,JPART,SYS_TYPE
      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
c--------------------------------------------------
c     Count number of nodes with imposed velocity => NUMVEL
C======================================================================|
      is_available = .false.
c
      nimpvel = 0   ! initialize counter of nodes
c--------------------------------------------------
c     /IMPVEL and /IMPVEL/LAGMUL
c--------------------------------------------------
      CALL hm_option_count('/impvel',NFVEL)
c
      CALL HM_OPTION_START('/impvel')
c
      DO I=1,NFVEL
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('rad_system_input_type'   ,SYS_TYPE ,IS_AVAILABLE,LSUBMODEL)
          CALL HM_GET_INTV('entityid'                ,GRNOD_ID ,IS_AVAILABLE,LSUBMODEL)
          IGS = NGR2USRN(GRNOD_ID,IGRNOD,NGRNOD,NNOD)
          IF (IGS > 0) THEN
            NIMPVEL = NIMPVEL + NNOD                                         
            IF (KEY(1:6) == 'lagmul') NFVLAG = NFVLAG + NNOD
          END IF
        ENDIF
c
      ENDDO ! DO I=1,NFVEL
c--------------------------      
c     /IMPVEL/FGEO    
c--------------------------         
      CALL HM_OPTION_COUNT('/impvel/fgeo',NFGEO)
c
      CALL HM_OPTION_START('/impvel/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) NIMPVEL = NIMPVEL + 1
            ENDDO
          ENDIF
c
          CALL HM_GET_INTV('distribution_table_count' ,NNOD ,IS_AVAILABLE,LSUBMODEL)
          NIMPVEL = NIMPVEL + NNOD
        ENDIF
      ENDDO ! DO I=1,NFGEO
c-----------
      RETURN