inivel.F File Reference

#include "implicit_f.inc"
#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"

Go to the source code of this file.

Functions/Subroutines
subroutine	inivel (v, vr, svr, itabm1)

Function/Subroutine Documentation

inivel()

subroutine inivel	(	dimension(3,numnod), intent(inout)	v,
		dimension(3,svr), intent(inout)	vr,
		integer, intent(in)	svr,
		integer, dimension(2*numnod), intent(in)	itabm1 )

Definition at line 34 of file inivel.F.

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