anchor.diff:           if (maux.gt.0) then
               do 324 ma=1,maux
                if (auxtype(ma).eq."xleft") then
  -c                # Assuming velocity at left-face, this fix
  +c                # Assuming velocity at bottom-face, this fix
   c                # preserves conservation in incompressible flow:
  +c                # -----------------------------------------------
  +c                # See note for "Side 1"
  +c                # -----------------------------------------------
                    jfine = (jc-1)*lratioy + nghost + l
                    auxl(iaddaux(ma,lind+1)) = aux(ma,mitot-nghost+1,jfine)
                  else
  
positive.code:              if (auxtype(ma).eq."yleft") then
  c                # Assuming velocity at bottom-face, this fix
  c                # preserves conservation in incompressible flow:
  c                # This assumes that the coarse grid aux variable is the 
  c                # same as the fine grid variable, i.e. edges match.
                   ifine = (ic-1)*lratiox + nghost + l
                   auxl(iaddaux(ma,lind+1)) = aux(ma,ifine,mjtot-nghost+1)
                 else
  
negative.code:              if (auxtype(ma).eq."yleft") then
  c                # Assuming velocity at bottom-face, this fix
  c                # preserves conservation in incompressible flow:
  c                # -----------------------------------------------
  c                # See note for "Side 1"
  c                # -----------------------------------------------
                   auxl(iaddaux(ma,i-nghost+1)) = aux(ma,i,nghost+1)
                 else
                   auxl(iaddaux(ma,i-nghost+1)) = aux(ma,i,nghost)
  