auxcoarsen.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	auxcoarsen (auxdub, midub, mjdub, auxbgc, mi2tot, mj2tot, naux, auxtype)

Function/Subroutine Documentation

auxcoarsen()

subroutine auxcoarsen	(	dimension(naux,midub, mjdub)	auxdub,
			midub,
			mjdub,
		dimension(naux,mi2tot,mj2tot)	auxbgc,
			mi2tot,
			mj2tot,
			naux,
		character*10, dimension(naux)	auxtype
	)

Coarsen the fine grid auxiliary data (with double the usual number of ghost cells to prepare coarsened data for error estimation).

Definition at line 6 of file auxcoarsen.f.

Referenced by prepbigstep().

      
       implicit double precision (a-h, o-z)

       dimension     auxdub(naux,midub, mjdub)
       dimension     auxbgc(naux,mi2tot,mj2tot)
       character*10  auxtype(naux)

c :::::::::::::::::::::::: coarsen ::::::::::::::::::::::::::::::::
c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::

       do 50 iaux = 1, naux

       if (auxtype(iaux) .eq. "center" .or. 
     .     auxtype(iaux) .eq. "capacity") then
         do 20 j = 1, mj2tot
            jfine = 2*(j-1) + 1
            do 20 i = 1, mi2tot
               ifine = 2*(i-1) + 1
               auxbgc(iaux,i,j) = (auxdub(iaux,ifine,jfine) +
     &                             auxdub(iaux,ifine+1,jfine)+
     &                             auxdub(iaux,ifine,jfine+1) +
     &                             auxdub(iaux,ifine+1,jfine+1))/4.d0
20       continue

       elseif (auxtype(iaux) .eq. "xleft") then 
         do 10 j = 1, mj2tot
            jfine = 2*(j-1) + 1
            do 10 i = 1, mi2tot
               ifine = 2*(i-1) + 1
               auxbgc(iaux,i,j) = (auxdub(iaux,ifine,jfine) +
     &                             auxdub(iaux,ifine,jfine+1)) /2.d0 
10       continue

       elseif (auxtype(iaux) .eq. "yleft") then 
         do 15 j = 1, mj2tot
            jfine = 2*(j-1) + 1
            do 15 i = 1, mi2tot
               ifine = 2*(i-1) + 1
               auxbgc(iaux,i,j) = (auxdub(iaux,ifine,jfine) +
     &                             auxdub(iaux,ifine+1,jfine))/2.d0
15       continue

       endif

50     continue

       return

Here is the caller graph for this function: