outval.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	outval (val, nvar, mitot, mjtot, mptr, outgrd, naux, aux)
	print solution and aux. More...

Function/Subroutine Documentation

outval()

subroutine outval	(	dimension(nvar,mitot,mjtot)	val,
			nvar,
			mitot,
			mjtot,
			mptr,
		logical	outgrd,
			naux,
		dimension(naux,mitot,mjtot)	aux
	)

print solution and aux.

variables to output.

Definition at line 5 of file outval.f.

References amr_module::cornxlo, amr_module::cornylo, amr_module::hxposs, amr_module::hyposs, amr_module::nestlevel, amr_module::nghost, amr_module::node, amr_module::outunit, and amr_module::rnode.

Referenced by outmsh().

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

      dimension  val(nvar,mitot,mjtot)
      dimension  aux(naux,mitot,mjtot)
      logical    outgrd


c ::::::::::::::::::::::OUTVAL :::::::::::::::::::::::::::::::
c print solution and aux. variables to output. 
c if cell outside domain, don't print soln. value - nothing
c currently in ghost cells.
c ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
c
      if (.not. outgrd) go to 99
      level = node(nestlevel,mptr)
      hx    = hxposs(level)
      hy    = hyposs(level)
      cornx = rnode(cornxlo,mptr) -  nghost*hx
      corny = rnode(cornylo,mptr) -  nghost*hy
c
      do 25 j=nghost+1,mjtot-nghost
      do 20 i=nghost+1,mitot-nghost

          x  = cornx + hx*(dble(i)-.5d0)
          y  = corny + hy*(dble(j)-.5d0)
          write(outunit,107) x,y,i,j,(val(ivar,i,j),ivar=1,nvar)
 107      format(2hx=,f6.3,2hy=,f6.3,3h,i=,i3,3h,j=,i3,' a=',
     *           e25.15)
c    *           5(e9.3,1x))
          if (naux.gt.0) write(outunit,108) (aux(iaux,i,j),iaux=1,naux)
 108      format(1x,'aux = ',7(e10.3,1x))

 20   continue
 25   continue

 99   return

Here is the caller graph for this function: