fluxsv.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	fluxsv (mptr, xfluxm, xfluxp, yfluxm, yfluxp, listbc,
	When a coarse grid cell is advanced, if it borders a coarse-fine interface, the flux or wave that emits from the interface and goes into the coarse cell is added to the corresponding location in node(ffluxptr, mkid) for conservative fix later, where mkid is grid number of the nested fine grid on the other side of that interface. More...

Function/Subroutine Documentation

fluxsv()

subroutine fluxsv	(		mptr,
		dimension(nvar,ndimx,ndimy)	xfluxm,
		dimension(nvar,ndimx,ndimy)	xfluxp,
		dimension(nvar,ndimx,ndimy)	yfluxm,
		dimension(nvar,ndimx,ndimy)	yfluxp,
		dimension(5,maxsp)	listbc
	)

When a coarse grid cell is advanced, if it borders a coarse-fine interface, the flux or wave that emits from the interface and goes into the coarse cell is added to the corresponding location in node(ffluxptr, mkid) for conservative fix later, where mkid is grid number of the nested fine grid on the other side of that interface.

Definition at line 11 of file fluxsv.f.

References amr_module::alloc, amr_module::dbugunit, amr_module::ffluxptr, amr_module::ndihi, amr_module::ndilo, amr_module::ndjhi, amr_module::ndjlo, amr_module::nestlevel, and amr_module::node.

Referenced by par_advanc().

     1                  ndimx,ndimy,nvar,maxsp,dtc,hx,hy)
c
      use amr_module
      implicit double precision (a-h,o-z)


      dimension xfluxp(nvar,ndimx,ndimy), yfluxp(nvar,ndimx,ndimy)
      dimension xfluxm(nvar,ndimx,ndimy), yfluxm(nvar,ndimx,ndimy)
      dimension listbc(5,maxsp)
c
c :::::::::::::::::::: fluxsv :::::::::::::::::::::::::
c
c  coarse grids should save their fluxes in cells adjacent to
c  their nested fine grids, for later conservation fixing.
c  listbc holds info for where to save which fluxes.
c  xflux holds 'f' fluxes, yflux holds 'g' fluxes.
c
c :::::::::::::::::::::::::::::;:::::::::::::::::::::::
 
 
      ispot   = 1
      level   = node(nestlevel,mptr)

 10      if (listbc(1,ispot).eq.0) go to 99          
c
         mkid     = listbc(4,ispot)
         intopl   = listbc(5,ispot)
         nx       = node(ndihi,mkid) - node(ndilo,mkid) + 1
         ny       = node(ndjhi,mkid) - node(ndjlo,mkid) + 1
         kidlst   = node(ffluxptr,mkid)
         i        = listbc(1,ispot)
         j        = listbc(2,ispot)
         inlist   = kidlst + nvar*(intopl-1) - 1
c
c side k(listbc 3) has which side of coarse cell has interface
c so can save appropriate fluxes.  (dont know why we didnt have
c which flux to save directly(i.e. put i+1,j to save that flux
c rather than putting in cell center coords).

         if (listbc(3,ispot) .eq. 1) then
c           ::::: cell i,j is on right side of a fine grid
            do 100 ivar = 1, nvar
               alloc(inlist + ivar) = -xfluxp(ivar,i,j)*dtc*hy

Here is the caller graph for this function: