shell_surface_buffer_remesh.F File Reference

#include "implicit_f.inc"
#include "scr17_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	shell_surface_buffer_remesh (ix, nix, nix1, nix2, ieltyp, nseg, iad_surf, buftmpsurf, clause, iparte, ipart, kshtree, shtree)

Function/Subroutine Documentation

shell_surface_buffer_remesh()

subroutine shell_surface_buffer_remesh	(	integer, dimension(nix,*)	ix,
		integer	nix,
		integer	nix1,
		integer	nix2,
		integer	ieltyp,
		integer	nseg,
		integer	iad_surf,
		integer, dimension(*)	buftmpsurf,
		type (set_)	clause,
		integer, dimension(*)	iparte,
		integer, dimension(lipart1,*)	ipart,
		integer	kshtree,
		integer, dimension(kshtree,*)	shtree )

Definition at line 29 of file shell_surface_buffer_remesh.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE setdef_mod
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      "scr17_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NIX,NIX1,NIX2,IELTYP,NSEG,IAD_SURF
      INTEGER IX(NIX,*),BUFTMPSURF(*),IPARTE(*),
     .        IPART(LIPART1,*),KSHTREE,SHTREE(KSHTREE,*)
!
      TYPE (SET_) :: CLAUSE
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER J,JJ,K,IP,NLEV,MY_LEV,INV4(4),INV3(3),NELEM
      DATA inv4/4,3,2,1/
      DATA inv3/3,2,1/
C-----------------------------------------------
!
      nelem = 0
      IF (ieltyp == 3) nelem = clause%NB_SH4N
      IF (ieltyp == 7) nelem = clause%NB_SH3N
!
      DO j=1,nelem
        jj= -huge(jj)
        IF(ieltyp==3) jj = clause%SH4N(j)
        IF(ieltyp==7) jj = clause%SH3N(j)
!
        ip=iparte(jj)
        nlev  =ipart(10,ip)
        my_lev=shtree(3,jj)
        IF (my_lev < 0) my_lev=-(my_lev+1)
        IF (my_lev==nlev) THEN
          nseg=nseg+1
          DO k=nix1,nix2
            buftmpsurf(iad_surf) = ix(k,jj)
            iad_surf=iad_surf+1
          ENDDO
!
          IF (nix2-nix1 == 2) THEN
            buftmpsurf(iad_surf) = buftmpsurf(iad_surf-1)  ! SH3N (NOD4 == NOD3)
            iad_surf=iad_surf+1
          ENDIF
!
          buftmpsurf(iad_surf) = ieltyp
          iad_surf=iad_surf+1
          buftmpsurf(iad_surf) = jj
          iad_surf=iad_surf+1
        ENDIF ! IF (MY_LEV==NLEV)
      ENDDO ! DO J=1,NELEM
C------------
      RETURN