hm_read_rbody_lagmul.F File Reference

#include "implicit_f.inc"
#include "analyse_name.inc"
#include "units_c.inc"
#include "param_c.inc"
#include "lagmult.inc"
#include "scr17_c.inc"
#include "com04_c.inc"
#include "r2r_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	hm_read_rbody_lagmul (rbyl, npbyl, lpbyl, igrnod, lsubmodel, itab, itabm1, ikine, ikine1lag, nom_opt)

Function/Subroutine Documentation

hm_read_rbody_lagmul()

subroutine hm_read_rbody_lagmul	(		rbyl,
		integer, dimension(nnpby,*)	npbyl,
		integer, dimension(*)	lpbyl,
		type (group_), dimension(ngrnod)	igrnod,
		type(submodel_data), dimension(*), intent(in)	lsubmodel,
		integer, dimension(*)	itab,
		integer, dimension(*)	itabm1,
		integer, dimension(*)	ikine,
		integer, dimension(*)	ikine1lag,
		integer, dimension(lnopt1,*)	nom_opt )

Definition at line 44 of file hm_read_rbody_lagmul.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE groupdef_mod
      USE message_mod
      USE r2r_mod
      USE submodel_mod
      USE hm_option_read_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   A n a l y s e   M o d u l e
C-----------------------------------------------
#include      "analyse_name.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "units_c.inc"
#include      "param_c.inc"
#include      "lagmult.inc"
#include      "scr17_c.inc"
#include      "com04_c.inc"
#include      "r2r_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NPBYL(NNPBY,*),LPBYL(*), 
     .        ITAB(*), ITABM1(*),IKINE(*),IKINE1LAG(*)
      my_real rbyl(nrby,*)
      INTEGER NOM_OPT(LNOPT1,*)
C-----------------------------------------------
      TYPE (GROUP_)  , DIMENSION(NGRNOD)  :: IGRNOD
      TYPE(SUBMODEL_DATA),INTENT(IN)::LSUBMODEL(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER J,K,L,N,NR,MM,ID,IGU,IGS,
     .        NSKEW,NSL,MSL,SUB_INDEX,NRB,NRB_R2R
      CHARACTER MESS*40
      CHARACTER(LEN=NCHARTITLE) :: TITR
      CHARACTER(LEN=NCHARKEY) :: KEY
      LOGICAL IS_AVAILABLE
      DATA mess/'rigid body definitions'/
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER  USR2SYS,NODGRNR6
      EXTERNAL USR2SYS,NODGRNR6
C-----------------------------------
C   NPBYL(1,N) = MAIN NODE
C   NPBYL(2,N) = NUMBER OF SECONDARY NODES + MAIN
C   NPBYL(3,N) = 
C   NPBYL(4,N) = 
C   NPBYL(5,N) = 
C   NPBYL(6,N) = IDENTIFICATEUR
C   NPBYL(7,N) = 
C   NPBYL(8,N) = 
C======================================================================|
      WRITE(IOUT,1000)
C---
      IS_AVAILABLE = .FALSE.
      CALL HM_OPTION_START('/rbody')
C---
      K   = 0
      NRB = 0
      NRB_R2R = 0
C
      DO NR=1,NRBODY
C
C--------------------------------------------------
C EXTRACT DATAS OF /RBODY/... LINE
C--------------------------------------------------
C
       NRB_R2R = NRB_R2R + 1     
       IF (NSUBDOM > 0) THEN
         IF(TAGRBY(NRB_R2R) == 0) CALL HM_SZ_R2R(TAGRBY,NRB_R2R,LSUBMODEL)
       ENDIF
C
       KEY=''
       CALL HM_OPTION_READ_KEY(LSUBMODEL,
     .                       OPTION_ID = ID,
     .                       OPTION_TITR = TITR,
     .                       KEYWORD2 = KEY,
     .                       SUBMODEL_INDEX = SUB_INDEX)
       IF(KEY(1:6)=='lagmul')THEN     
        NRB = NRB + 1
        IF (NSUBDOM > 0) THEN ! TAGRBY is allocated only if NSUBDOM>0
          IF(TAGRBY(NRB) == 0) CALL HM_SZ_R2R(TAGRBY,NRB,LSUBMODEL)
        ENDIF
C---
        NOM_OPT(1,NRB)=ID
        CALL FRETITL(TITR,NOM_OPT(LNOPT1-LTITR+1,NRB),LTITR)
C---
        CALL HM_GET_INTV('node_id',MM,IS_AVAILABLE,LSUBMODEL)
        MM  = USR2SYS(MM,ITABM1,MESS,ID)
C---
        CALL HM_GET_INTV('grnd_id',IGU,IS_AVAILABLE,LSUBMODEL)
        NSL = NODGRNR6(MM,IGU,IGS,LPBYL(K+1),IGRNOD,ITABM1,MESS,ID)
        MSL = NSL+1
C---
        LPBYL(K+MSL) = MM
C---
        IF (NSL == 0) THEN
          CALL ANCMSG(MSGID=352,
     .                MSGTYPE=MSGWARNING,
     .                ANMODE=ANINFO_BLIND_2,
     .                I1=ID,
     .                C1=TITR)
        ENDIF
        CALL ANODSET(MM, CHECK_RB_M)
        DO J=1, NSL
          CALL ANODSET(LPBYL(J+K), CHECK_RB_S)
        ENDDO
C---
        DO J=1,MSL
          CALL KINSET(512,ITAB(LPBYL(J+K)),IKINE(LPBYL(J+K)),7,0,
     .                IKINE1LAG(LPBYL(J+K)))
        ENDDO
C---
        NPBYL(1,NRB) = MM
        NPBYL(2,NRB) = MSL
        NPBYL(6,NRB) = ID
        LAG_NCL = LAG_NCL + NSL*6
        LAG_NKL = LAG_NKL + NSL*21
C---
        WRITE(IOUT,1100)ID,TRIM(TITR),ITAB(MM),MSL
        WRITE(IOUT,1101)
        WRITE(IOUT,1102) (ITAB(LPBYL(J+K)),J=1,NSL)
        K = K + 3*MSL
       END IF ! IF(KEY(1:6)=='lagmul')THEN  
      ENDDO
C---
      RETURN
C------------------------------
1000  FORMAT(
     . /'      rigid body definitions(lagrange multipliers)'
     . /'      -------------------------------------------- '/)
1100  FORMAT( /5X,'rigid body id ',I10,1X,A,
     .       /10X,'primary node                            ',I10
     .       /10X,'number of nodes                         ',I10)
1101  FORMAT( 10X,'secondary nodes ')
1102  FORMAT(  9X,10I10)