errest.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	errest (nvar, naux, lcheck, mptr, nx, ny)

Function/Subroutine Documentation

errest()

subroutine errest	(		nvar,
			naux,
			lcheck,
			mptr,
			nx,
			ny
	)

for this grid at level lcheck: estimate the error by taking a large (2h,2k) step based on the values in the old storage loc., then take one regular (and for now wasted) step based on the new info. compare using an error relation for a pth order accurate integration formula. flag error plane as either bad (needs refinement), or good.

Definition at line 5 of file errest.f.

References amr_module::alloc, errf1(), amr_module::ibuff, amr_module::nghost, amr_module::node, prepbigstep(), prepregstep(), amr_module::store1, amr_module::store2, amr_module::storeaux, amr_module::storeflags, and amr_module::tempptr.

Referenced by flagger().

c
      use amr_module
      implicit double precision (a-h,o-z)
c
c   ### changed to stack based storage 2/23/13 
c   ### and broken into smaller routines to minimize 
c   ### stack space
     
      double precision valbgc(nvar,nx/2+2*nghost,ny/2+2*nghost)
      double precision auxbgc(naux,nx/2+2*nghost,ny/2+2*nghost)
     
 
c :::::::::::::::::::::::::: ERREST :::::::::::::::::::::::::::::::::::
c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
c
       mitot  = nx + 2*nghost
       mjtot  = ny + 2*nghost
       locnew = node(store1,mptr)
       locold = node(store2,mptr)
       locaux = node(storeaux,mptr)
       mi2tot = nx/2  + 2*nghost
       mj2tot = ny/2  + 2*nghost
c
c     prepare double the stencil size worth of boundary values,
c            then coarsen them for the giant step integration.
       midub = nx+4*nghost
       mjdub = ny+4*nghost
c
       call prepbigstep(nvar,naux,lcheck,mptr,nx,ny,midub,
     .                    mjdub,valbgc,auxbgc,mi2tot,mj2tot)

c
c  the one giant step based on old values is done. now take
c  one regular step based on new values. 
c  boundary values already in locbig, (see subr. flagger)
c
      locbig = node(tempptr,mptr)
      locaux = node(storeaux,mptr)
      call prepregstep(nvar,naux,lcheck,mptr,nx,ny,mitot,mjtot,
     .                 alloc(locbig),alloc(locaux))
c
c     ## locamrflags allocated in flagger. may previously have been used 
c     ## by flag2refine so make sure not to overwrite
      locamrflags = node(storeflags, mptr)    
      mbuff = max(nghost,ibuff+1)  
      mibuff = nx + 2*mbuff 
      mjbuff = ny + 2*mbuff 
      call errf1(alloc(locbig),nvar,valbgc,mptr,mi2tot,mj2tot,
     1           mitot,mjtot,alloc(locamrflags),mibuff,mjbuff)

c
      return

Here is the call graph for this function:

Here is the caller graph for this function: