inisrf.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	inisrf (x, v, vr, npby, rby, igrsurf, bufsf)

Function/Subroutine Documentation

inisrf()

subroutine inisrf	(		x,
			v,
			vr,
		integer, dimension(nnpby,*)	npby,
			rby,
		type (surf_), dimension(nsurf)	igrsurf,
			bufsf )

Definition at line 34 of file inisrf.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE groupdef_mod
      USE message_mod
C=======================================================================
C    * Verifies that there are not two RIGID BODIES (or two kinematic
C      conditions) which define the movement of the same surface.
C    * Initializes the position and velocities of the point where forces are calculated
C      in the interfaces.
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      "com04_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NPBY(NNPBY,*)
C     REAL
      my_real
     .   x(3,*),v(3,*),vr(3,*),rby(nrby,*),bufsf(*)
      TYPE (SURF_)   , DIMENSION(NSURF)   :: IGRSURF
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER N, ISU, ADRSRF, IM
C-----------------------------------------------
      DO 100 n=1,nrbykin
        isu=npby(8,n)
        IF (isu==0) GOTO 100
        IF (igrsurf(isu)%NB_MADYMO/=0) GOTO 999
C
        igrsurf(isu)%NB_MADYMO=n
        igrsurf(isu)%TYPE_MADYMO=1
C
        im=npby(1,n)
        adrsrf=igrsurf(isu)%IAD_BUFR
        bufsf(adrsrf+16)=x(1,im)
        bufsf(adrsrf+17)=x(2,im)
        bufsf(adrsrf+18)=x(3,im)
        bufsf(adrsrf+19)=v(1,im)
        bufsf(adrsrf+20)=v(2,im)
        bufsf(adrsrf+21)=v(3,im)
        bufsf(adrsrf+22)  =vr(1,im)
        bufsf(adrsrf+23)=vr(2,im)
        bufsf(adrsrf+24)=vr(3,im)
 100  CONTINUE
      GOTO 9999
C-----------------------------------------------
 999  CONTINUE
C      WRITE(ISTDO,*)' ** ERROR : MOTION OF SURFACE MULTIPLY DEFINED'
C      WRITE(IOUT,*) ' ** ERROR : SURFACES MOTION DEFINITION :'
C      WRITE(IOUT,*) ' ** ERROR : MOTION OF SURFACE MULTIPLY DEFINED :'
C      WRITE(IOUT,*) '            MOTION OF SURFACE ',ISURF(1,ISU),
C     .                 ' IS DEFINED BY MORE THAN ONE RIGID BODY.'
C      IERR=IERR+1
      CALL ancmsg(msgid=146,
     .            msgtype=msgerror,
     .            anmode=aninfo_blind_1,
     .            i1=igrsurf(isu)%ID)
C-----------------------------------------------
 9999 CONTINUE
      RETURN