fixcapaq.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	fixcapaq (val, aux, mitot, mjtot, valc, auxc, mic, mjc, nvar, naux, levc, setflags)

Function/Subroutine Documentation

fixcapaq()

subroutine fixcapaq	(	dimension(nvar,mitot,mjtot)	val,
		dimension(naux,mitot,mjtot)	aux,
			mitot,
			mjtot,
		dimension(nvar,mic,mjc)	valc,
		dimension(naux,mic,mjc)	auxc,
			mic,
			mjc,
			nvar,
			naux,
			levc,
		dimension(mitot,mjtot)	setflags
	)

Definition at line 6 of file fixcapaq.f.

References amr_module::intratx, amr_module::intraty, amr_module::mcapa, amr_module::needs_to_be_set, amr_module::nghost, and setflags().

Referenced by filval().

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

c
c :::::::::::::::::::::::  FIXCAPAQ ::::::::::::::::::::::::::::::
c  new fine grid solution q was linearly interpolated. but want
c  to conserve kappa*q, not q. calculate the discrepancy
c  in kappa*q using this q, and modify q to account for it.
c ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::

      dimension   val(nvar,mitot,mjtot), valc(nvar,mic,mjc)
      dimension   aux(naux,mitot,mjtot), auxc(naux,mic,mjc)
      dimension   setflags(mitot,mjtot)

      dcapamax = 0.d0
      lratiox  = intratx(levc)
      lratioy  = intraty(levc)

      do 10 ic = 2, mic-1
      do 10 jc = 2, mjc-1


       do 15 ivar = 1, nvar

       capaqfine = 0.d0

       do 20 ico = 1, lratiox
       ifine = (ic-2)*lratiox + nghost + ico
       do 20 jco = 1, lratioy
         jfine = (jc-2)*lratioy + nghost + jco
         capaqfine = capaqfine + aux(mcapa,ifine,jfine)*
     &                           val(ivar,ifine,jfine)
20     continue

       dcapaq = auxc(mcapa,ic,jc)*valc(ivar,ic,jc)-
     &          capaqfine/(lratiox*lratioy)
       dcapamax = dmax1(dcapamax,dabs(dcapaq))
      
       do 30 ico = 1, lratiox
       ifine = (ic-2)*lratiox + nghost + ico
       do 30 jco = 1, lratioy
         jfine = (jc-2)*lratioy + nghost + jco

         if (setflags(ifine,jfine) .eq. needs_to_be_set) then 
         ! was set by coarsegrid, need to check for adjustment
           val(ivar,ifine,jfine) = val(ivar,ifine,jfine) +
     &                             dcapaq/aux(mcapa,ifine,jfine)
         endif
30     continue

15     continue

10     continue

c      write(6,*)" max discrepancy ", dcapamax

       return

Here is the call graph for this function:

Here is the caller graph for this function: