lgmini_fxv.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	lgmini_fxv (iadll, jll, lll, ibfv, vel, nc, mass, iner, itab, nom_opt)

Function/Subroutine Documentation

lgmini_fxv()

subroutine lgmini_fxv	(	integer, dimension(*)	iadll,
		integer, dimension(*)	jll,
		integer, dimension(*)	lll,
		integer, dimension(nifv,*)	ibfv,
			vel,
		integer	nc,
			mass,
			iner,
		integer, dimension(*)	itab,
		integer, dimension(lnopt1,*)	nom_opt )

Definition at line 33 of file lgmini_fxv.F.

      USE message_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-----------------------------------------------
#include      "param_c.inc"
#include      "com04_c.inc"
#include      "scr17_c.inc"
#include      "lagmult.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NC, IADLL(*), JLL(*), LLL(*), IBFV(NIFV,*),
     .        ITAB(*)
      my_real
     .        vel(lfxvelr,*), mass(*), iner(*)
      INTEGER NOM_OPT(LNOPT1,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER J, N, IK, ISK, NNO
      INTEGER ID
      CHARACTER(LEN=NCHARTITLE) :: TITR
C======================================================================| 
      DO n=1,nfxvel
        id=nom_opt(1,ibfv(12,n))
        CALL fretitl2(titr,nom_opt(lnopt1-ltitr+1,ibfv(12,n)),ltitr)
        IF (ibfv(8,n)/=0) THEN
          nno = iabs(ibfv(1,n))
          isk = ibfv(2,n)/10
          j=ibfv(2,n)-10*isk
          nc = nc + 1
          IF(j<=3.AND.mass(nno)==0)THEN
            CALL ancmsg(msgid=535,
     .                  msgtype=msgerror,
     .                  anmode=aninfo_blind_1,
     .                  c1='FIXED TRANSLATIONAL VELOCITY',
     .                  i1=id,
     .                  c2='FIXED TRANSLATIONAL VELOCITY',
     .                  c3=titr,c4=' ',
     .                  i2=itab(iabs(ibfv(1,n))))
            ELSEIF(j>3.AND.mass(nno)==0)THEN
            CALL ancmsg(msgid=535,
     .                  msgtype=msgerror,
     .                  anmode=aninfo_blind_1,
     .                  c1='FIXED ROTATIONAL VELOCITY',
     .                  i1=id,
     .                  c2='FIXED ROTATIONAL VELOCITY',
     .                  c3=titr,c4=' ',
     .                  i2=itab(iabs(ibfv(1,n))))
            ENDIF
Cm41u07---
            IF(nc>lag_ncf)THEN
            CALL ancmsg(msgid=468,
     .                  msgtype=msgerror,
     .                  anmode=aninfo,
     .                  i1=id,
     .                  c1='IMPOSED VELOCITY',
     .                  c2='IMPOSED VELOCITY',
     .                  c3=titr)
          ENDIF
          IF (isk<=1) THEN
            iadll(nc+1)=iadll(nc) + 1
            IF(iadll(nc+1)-1>lag_nkf)THEN
              CALL ancmsg(msgid=469,
     .                    msgtype=msgerror,
     .                    anmode=aninfo,
     .                    c1='imposed velocity',
     .                    I1=ID,
     .                    C2='imposed velocity',
     .                    C3=TITR)
            ENDIF
            IK = IADLL(NC)
            LLL(IK) = NNO
            JLL(IK) = J
.AND.            IF(J>3INER(NNO)==0)THEN
               CALL ANCMSG(MSGID=536,
     .                     MSGTYPE=MSGERROR,
     .                     ANMODE=ANINFO_BLIND_1,
     .                     C1='fixed rotational velocity',
     .                     I1=ID,
     .                     C2='fixed rotational velocity',
     .                     C3=TITR,
     .                     I2=ITAB(IABS(IBFV(1,N))))
              ENDIF
          ELSE
            IADLL(NC+1)=IADLL(NC) + 3
            IF(IADLL(NC+1)-1>LAG_NKF)THEN
              CALL ANCMSG(MSGID=469,
     .                    MSGTYPE=MSGERROR,
     .                    ANMODE=ANINFO,
     .                    I1=ID,
     .                    C1='imposed velocity',
     .                    C2='imposed velocity',
     .                    C3=TITR)
            ENDIF
            IF(J<=3)THEN
              IK = IADLL(NC)
              LLL(IK) = NNO
              JLL(IK) = 1
              IK = IK + 1
              LLL(IK) = NNO
              JLL(IK) = 2
              IK = IK + 1
              LLL(IK) = NNO
              JLL(IK) = 3
            ELSE
              IF(INER(NNO)==0)THEN
               CALL ANCMSG(MSGID=536,
     .                     MSGTYPE=MSGERROR,
     .                     ANMODE=ANINFO_BLIND_1,
     .                     C1='fixed rotational velocity',
     .                     I1=ID,
     .                     C2='fixed rotational velocity',
     .                     C3=TITR,
     .                     I2=ITAB(IABS(IBFV(1,N))))
                ENDIF
              IK = IADLL(NC)
              LLL(IK) = NNO
              JLL(IK) = 4
              IK = IK + 1
              LLL(IK) = NNO
              JLL(IK) = 5
              IK = IK + 1
              LLL(IK) = NNO
              JLL(IK) = 6
            ENDIF
          ENDIF
        ENDIF
      ENDDO
C---
      RETURN