domgrid.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	domgrid (lbase, lcheck)

Function/Subroutine Documentation

domgrid()

subroutine domgrid	(		lbase,
			lcheck
	)

loop over base grids to get proper nesting domain for grids at level lcheck but only upsize to the lcheck grids dimensions

Definition at line 5 of file domgrid.f.

References amr_module::alloc, amr_module::domflags2, amr_module::domflags_base, amr_module::ibuff, igetsp(), amr_module::levelptr, amr_module::lstart, amr_module::ndihi, amr_module::ndilo, amr_module::ndjhi, amr_module::ndjlo, amr_module::nghost, amr_module::node, and setdomflags().

c
      use amr_module
      implicit double precision (a-h,o-z)

      mbuff = max(nghost,ibuff+1)
c
c

      mptr = lstart(lcheck)
 10   continue
      ilo    = node(ndilo,mptr)
      ihi    = node(ndihi,mptr)
      jlo    = node(ndjlo,mptr)
      jhi    = node(ndjhi,mptr)
      nx     = node(ndihi,mptr) - node(ndilo,mptr) + 1
      ny     = node(ndjhi,mptr) - node(ndjlo,mptr) + 1
c     up to mbuff cells on each side might be flagged and
c     buffered, so to allow shrinkage, need yet one additional
c     cell on each side to be set from base grids
      mibuff = nx + 2*mbuff
      mjbuff = ny + 2*mbuff
      ibytesperdp = 8

c   bad names, for historical reasons. they are both smae size now
      locdomflags = igetsp( (mibuff*mjbuff)/ibytesperdp+1)
      locdom2 = igetsp( (mibuff*mjbuff)/ibytesperdp+1)


      node(domflags_base,mptr) = locdomflags
      node(domflags2,mptr) = locdom2
      call setdomflags(mptr,alloc(locdomflags),ilo,ihi,jlo,jhi,
     .                 mbuff,lbase,lcheck,mibuff,mjbuff)

      mptr = node(levelptr, mptr)
      if (mptr .ne. 0) go to 10

      return

Here is the call graph for this function: