shiftset2.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	shiftset2 (rectflags, ilo, ihi, jlo, jhi, mbuff)
	For an input grid, flag cells near the previously flagged cells for creating buffer zone. More...

Function/Subroutine Documentation

shiftset2()

subroutine shiftset2	(	dimension(ilo-mbuff:ihi+mbuff,jlo-mbuff:jhi+mbuff)	rectflags,
			ilo,
			ihi,
			jlo,
			jhi,
			mbuff
	)

For an input grid, flag cells near the previously flagged cells for creating buffer zone.

Shift by + or - 1 in either direction to do bitwise calculus for proper nesting, buffering, etc. similar to cshift on CM machine. includes periodic buffering as well.

NEWER VERSION: DOES ALL DIRS AT SAME TIME rectflags array has been augmented by enough border cells to do buffering in place in the grid later will look to see if flagged pts are properly nested

Parameters

rectflags	array to be flagged
ilo	global i index of the left border of the grid being projected to (being flagged)
ihi	global i index of the right border of the grid being projected to (being flagged)
jlo	global j index of the lower border of the grid being projected to (being flagged)
jhi	global i index of the upper border of the grid being projected to (being flagged)
mbuff	width of the buffer zone

Definition at line 26 of file shiftset2.f.

References amr_module::ibuff.

Referenced by bufnst2().

       use amr_module
       implicit double precision (a-h, o-z)
       dimension rectflags(ilo-mbuff:ihi+mbuff,jlo-mbuff:jhi+mbuff)
       dimension copyflags(ilo-mbuff:ihi+mbuff,jlo-mbuff:jhi+mbuff)

       do j = jlo-mbuff, jhi+mbuff
       do i = ilo-mbuff, ihi+mbuff
          copyflags(i,j) = 0
       end do
       end do

c  note: not looking at ghost cells, only real cells should be flagged
c  but can buffer the flags into the ghost zone
       do j = jlo, jhi
       do i = ilo, ihi
          rflag = rectflags(i,j)
          if (rflag .gt. 0)  then
c            cell is flagged, buffer in all dirs by ibuff
c(note this is not nec same as mbuff)
c            use second array to avoid propagation
             mlo = i - ibuff
             mhi = i + ibuff
             klo = j - ibuff
             khi = j + ibuff
             do k = klo, khi 
             do m = mlo, mhi
                copyflags(m,k) = rflag  ! copy the flag (doesnt distinguish buffer flag from orig flag)
             end do
             end do
          endif

       end do
       end do


c   copy back. need flags in original array, not temp scratch array

       do 60 j = jlo-mbuff, jhi+mbuff
       do 60 i = ilo-mbuff, ihi+mbuff
         rectflags(i,j) = copyflags(i,j)
 60    continue


       return

Here is the caller graph for this function: