inivel.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/.
C
!||====================================================================
!||    inivel                 ../starter/source/initial_conditions/general/inivel/inivel.F
!||--- called by ------------------------------------------------------
!||    lectur                 ../starter/source/starter/lectur.F
!||--- calls      -----------------------------------------------------
!||    usrtos                 ../starter/source/system/sysfus.F
!||--- uses       -----------------------------------------------------
!||    message_mod            ../starter/share/message_module/message_mod.F
!||    submodel_mod           ../starter/share/modules1/submodel_mod.F
!||====================================================================
      SUBROUTINE inivel(V       ,VR      ,SVR ,ITABM1  )
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE unitab_mod
      USE submodel_mod
      USE message_mod
      USE multi_fvm_mod
      USE groupdef_mod
      USE names_and_titles_mod , ONLY : ncharline
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      "scr16_c.inc"
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "units_c.inc"
#include      "scr03_c.inc"
#include      "param_c.inc"
#include      "sphcom.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER  ,INTENT(IN) :: ITABM1(2*NUMNOD)
      INTEGER  ,INTENT(IN) :: SVR
      my_real,INTENT(INOUT) :: v(3,numnod),vr(3,svr)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,N
      my_real
     .   v1, v2, v3
      CHARACTER(LEN=NCHARLINE) :: LLINE
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER  USRTOS
C=======================================================================
C--------------------------------------------------
C     V INI INI Ynn FILE
C--------------------------------------------------
      IF (isigi==3.OR.isigi==4.OR.isigi==5) THEN
  120   READ(iin4,fmt='(A)',END=199,ERR=199)lline
  122   IF(lline(1:31)/='/NODAL     /VECTOR    /VELOCITY')GOTO 120
        READ(iin4,fmt='(A)',END=199,ERR=199)lline
C
  125   READ(iin4,fmt='(A)',END=130,ERR=199)lline
        IF(lline(1:1)=='#')GOTO 125
        IF(lline(1:1)=='/')GOTO 130
          IF (ioutp_fmt==2) THEN
            READ(lline,'(I8,3F16.0)')n,v1,v2,v3
          ELSE
            READ(lline,'(I10,3F20.0)')n,v1,v2,v3
          ENDIF
          i = usrtos(n,itabm1)
          IF(i/=0)THEN
            v(1,i) = v1
            v(2,i) = v2
            v(3,i) = v3
          END IF
        GOTO 125
 130    CONTINUE
 199    CONTINUE
        rewind(iin4)
C-------------------
  220   READ(iin4,fmt='(A)',END=299,ERR=299)lline
  222   IF(lline(1:31)/='/NODAL     /VECTOR    /VELO_ROT')GOTO 220
        READ(iin4,fmt='(A)',END=299,ERR=299)lline
C
  225   READ(iin4,fmt='(A)',END=230,ERR=299)lline
        IF(lline(1:1)=='#')GOTO 225
        IF(lline(1:1)=='/')GOTO 230
          IF (ioutp_fmt==2) THEN
            READ(lline,'(I8,3F16.0)')n,v1,v2,v3
          ELSE
            READ(lline,'(I10,3F20.0)')n,v1,v2,v3
          ENDIF
          i = usrtos(n,itabm1)
          IF(i/=0)THEN
            vr(1,i) = v1
            vr(2,i) = v2
            vr(3,i) = v3
          END IF
        GOTO 225
 230    CONTINUE
 299    CONTINUE
        rewind(iin4)
      ENDIF
!
!-----------
      RETURN
      END