surftage.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	surftage (numel, ix, nix, nix1, nix2, ieltyp, tagbuf, isu, nseg, flag, indx, nindx, nseg_tot)

Function/Subroutine Documentation

surftage()

subroutine surftage	(	integer	numel,
		integer, dimension(nix,*)	ix,
		integer	nix,
		integer	nix1,
		integer	nix2,
		integer	ieltyp,
		integer, dimension(*)	tagbuf,
		type (surf_)	isu,
		integer	nseg,
		integer	flag,
		integer, dimension(*)	indx,
		integer	nindx,
		integer	nseg_tot )

Definition at line 30 of file surftage.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE groupdef_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NIX,IX(NIX,*),NIX1,NIX2,NUMEL,IELTYP,
     .      TAGBUF(*),NSEG,FLAG,INDX(*),NINDX,NSEG_TOT
!
      TYPE (SURF_) :: ISU
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER JJ,K,II,INV4(4),INV3(3)
      DATA inv4/4,3,2,1/
      DATA inv3/3,2,1/
C=======================================================================
       nseg=nseg+nindx
!
       IF (flag == 1) THEN
         DO ii = 1, nindx
            jj=indx(ii)
            nseg_tot = nseg_tot + 1
            IF(tagbuf(jj) == 1)THEN
              DO k=nix1,nix2
                isu%NODES(nseg_tot,k-1) = ix(k,jj)
              ENDDO
            ELSEIF(tagbuf(jj) == -1)THEN
              IF (ieltyp == 3) THEN
                DO k=nix2,nix1,-1
                  isu%NODES(nseg_tot,inv4(k-1)) = ix(k,jj)
                ENDDO
               ELSEIF (ieltyp == 7) THEN
                DO k=nix2,nix1,-1
                  isu%NODES(nseg_tot,inv3(k-1)) = ix(k,jj)
                ENDDO
               ENDIF ! IF (IELTYP == 3)
            ENDIF
            IF(nix2-nix1 == 2)THEN
              isu%NODES(nseg_tot,4) = isu%NODES(nseg_tot,3)
            ENDIF
            isu%ELTYP(nseg_tot) = ieltyp
            isu%ELEM(nseg_tot)  = jj
         ENDDO
       ENDIF
C-----------
      RETURN