fvinject.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	fvinject (ibuf, elem, ibagjet, njet, igrsurf, itagel, nn, nel, id, tagvent, titr, itab)

Function/Subroutine Documentation

fvinject()

subroutine fvinject	(	integer, dimension(*)	ibuf,
		integer, dimension(3,nel)	elem,
		integer, dimension(nibjet,*)	ibagjet,
		integer	njet,
		type (surf_), dimension(nsurf)	igrsurf,
		integer, dimension(nel), intent(inout)	itagel,
		integer	nn,
		integer	nel,
		integer	id,
		integer, dimension(numnod)	tagvent,
		character(len=nchartitle)	titr,
		integer, dimension(numnod)	itab )

Definition at line 32 of file fvinject.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE groupdef_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"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IBUF(*), ELEM(3,NEL), IBAGJET(NIBJET,*), 
     .        NJET, 
     .        NN, NEL, ID, TAGVENT(NUMNOD),ITAB(NUMNOD)
      INTEGER, DIMENSION(NEL), INTENT(INOUT) :: ITAGEL
      CHARACTER(len=nchartitle) :: TITR
      TYPE (SURF_)   , DIMENSION(NSURF)   :: IGRSURF
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, II, ITABINV(NUMNOD), ITAGINJ(NN), IINJ, JINJ,
     .        ISU, NELSU,  IEL, NG1, NG2, NG3, NG4, N1, N2, N3, N4,
     .        NALL, NN1, NN2, IERROR, INOD
C
      DO i=1,numnod
         itabinv(i)=0
      ENDDO
C
      DO i=1,nn
         ii=ibuf(i)
         IF(itabinv(ii)==0)itabinv(ii)=i
C         ITAGINJ(I)=0
      ENDDO
C
      DO iinj=1,njet
         isu=ibagjet(14,iinj)
         IF(isu==0) cycle
         nelsu=igrsurf(isu)%NSEG
         inod=0
         DO i = 1, nn
            itaginj(i) = 0
         ENDDO
         DO iel=1,nelsu
            ng1 = igrsurf(isu)%NODES(iel,1)
            ng2 = igrsurf(isu)%NODES(iel,2)
            ng3 = igrsurf(isu)%NODES(iel,3)
            ng4 = igrsurf(isu)%NODES(iel,4)
            n1=itabinv(ng1)
            n2=itabinv(ng2)
            n3=itabinv(ng3)
            n4=itabinv(ng4)
            ierror = 0
            IF(n1==0.AND.tagvent(ng1)==0) THEN
              ierror = 1
            ELSE
              jinj=itaginj(n1)
              IF(jinj > 0 .AND. jinj < iinj) THEN
                  inod=inod+1
                  CALL ancmsg(msgid=1046,msgtype=msgwarning,
     .                        anmode=aninfo_blind_2,
     .                        i1=id,i2=itab(ng1),i3=jinj)
              ENDIF
              itaginj(n1)=iinj
            END IF
            IF(n2==0.AND.tagvent(ng2)==0) THEN
              ierror = 1
            ELSE
              jinj=itaginj(n2)
              IF(jinj > 0 .AND. jinj < iinj) THEN
                  inod=inod+1
                  CALL ancmsg(msgid=1046,msgtype=msgwarning,
     .                        anmode=aninfo_blind_2,
     .                        i1=id,i2=itab(ng2),i3=jinj)
              ENDIF
              itaginj(n2)=iinj
            END IF
            IF(n3==0.AND.tagvent(ng3)==0) THEN
              ierror = 1
            ELSE
              jinj=itaginj(n3)
              IF(jinj > 0 .AND. jinj < iinj) THEN
                  inod=inod+1
                  CALL ancmsg(msgid=1046,msgtype=msgwarning,
     .                        anmode=aninfo_blind_2,
     .                        i1=id,i2=itab(ng3),i3=jinj)
              ENDIF
              itaginj(n3)=iinj
            END IF
            IF(n4==0.AND.tagvent(ng4)==0) THEN
              ierror = 1
            ELSE
              jinj=itaginj(n4)
              IF(jinj > 0 .AND. jinj < iinj) THEN
                  inod=inod+1
                  CALL ancmsg(msgid=1046,msgtype=msgwarning,
     .                        anmode=aninfo_blind_2,
     .                        i1=id,i2=itab(ng4),i3=jinj)
              ENDIF
              itaginj(n4)=iinj
            END IF
            IF(ierror==1)THEN
              CALL ancmsg(msgid=632,msgtype=msgerror,anmode=aninfo,
     .                    i1=id,i2=iinj,c1=titr,c2='INFLATOR')
            ENDIF                
         ENDDO
         IF(inod > 0) THEN
           CALL ancmsg(msgid=1047,msgtype=msgwarning,anmode=aninfo,
     .                 i1=id,i2=iinj,i3=inod)
         ENDIF
C
         DO iel=1,nel
            n1=elem(1,iel)
            n2=elem(2,iel)
            n3=elem(3,iel)
            nall=itaginj(n1)*itaginj(n2)*itaginj(n3)
            IF (nall/=0) THEN
               nn1=itaginj(n2)-itaginj(n1)
               nn2=itaginj(n3)-itaginj(n1)
               IF (nn1==0.AND.nn2==0) THEN
                  IF (itagel(iel) > 0) THEN
                     CALL ancmsg(msgid=1616,msgtype=msgwarning,
     .                    anmode=aninfo_blind_2,
     .                    i1=id,i2=itagel(iel),i3=itaginj(n1))
                  ENDIF
                  itagel(iel)=itaginj(n1)
               ENDIF
            ENDIF
         ENDDO
      ENDDO
 
C
      RETURN