putsp.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	putsp (lbase, level, nvar, naux)
	Reclaim list space in nodes cfluxptr and ffluxptr for all grids at level level More...

Function/Subroutine Documentation

putsp()

subroutine putsp	(		lbase,
			level,
			nvar,
			naux
	)

Reclaim list space in nodes cfluxptr and ffluxptr for all grids at level level

Definition at line 7 of file putsp.f.

References amr_module::cfluxptr, amr_module::ffluxptr, amr_module::intratx, amr_module::intraty, amr_module::levelptr, amr_module::lfine, amr_module::listsp, amr_module::lstart, amr_module::ndihi, amr_module::ndilo, amr_module::ndjhi, amr_module::ndjlo, amr_module::node, and reclam().

Referenced by cleanup(), and gfixup().

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

c
c ::::::::::::::::::::::::::::::: PUTSP :::::::::::::::::::::::::
c
c reclaim list space in nodes cfluxptr and ffluxptr for grids at level
c
c first compute max. space allocated in node cfluxptr.
c
c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
c
      if (level .eq. lfine) go to 30
c
      mptr  = lstart(level)
 20      call reclam(node(cfluxptr,mptr), 5*listsp(level))
         node(cfluxptr,mptr) = 0
      mptr  = node(levelptr,mptr)
      if (mptr .ne. 0) go to 20
c
 30    if (level .eq. lbase) go to 99
      mptr = lstart(level)
 40       nx    = node(ndihi,mptr) - node(ndilo,mptr) + 1
          ny    = node(ndjhi,mptr) - node(ndjlo,mptr) + 1
          ikeep = nx/intratx(level-1)
          jkeep = ny/intraty(level-1)
          lenbc = 2*(ikeep+jkeep)
c         twice perimeter since saving plus or minus fluxes 
c         plus coarse solution storage
          call reclam(node(ffluxptr,mptr), 2*nvar*lenbc+naux*lenbc)
          mptr  = node(levelptr,mptr)
          if (mptr .ne. 0) go to 40
c
 99    return

Here is the call graph for this function:

Here is the caller graph for this function: