lgmini_fxv.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    lgmini_fxv             ../starter/source/tools/lagmul/lgmini_fxv.F
!||--- called by ------------------------------------------------------
!||    lagm_ini               ../starter/source/tools/lagmul/lagm_ini.F
!||--- calls      -----------------------------------------------------
!||    ancmsg                 ../starter/source/output/message/message.F
!||    fretitl2               ../starter/source/starter/freform.F
!||--- uses       -----------------------------------------------------
!||    message_mod            ../starter/share/message_module/message_mod.F
!||====================================================================
      SUBROUTINE lgmini_fxv(IADLL  ,JLL    ,LLL    ,IBFV   ,VEL    ,
     .                      NC, MASS, INER, ITAB   ,NOM_OPT)
      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 I, 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
            IF(j>3.AND.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
          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
      END