anchor.diff:               if (auxtype(ma).eq."xleft") then
   c                # Assuming velocity at left-face, this fix
   c                # preserves conservation in incompressible flow:
  +c                # -----------------------------------------------
  +c                # This branch only gets referenced if 'use_fix' = false.
  +c                # If 'use_fix = true', then all data is assumed
  +c                # to be cell-centered, since each cell stores information
  +c                # about all four faces.
  +c                # -----------------------------------------------
                    auxl(iaddaux(ma,j-nghost+1)) = aux(ma,nghost+1,j)
                  else
   c                # Normal case -- we set the aux arrays 
  
positive.code:              if (auxtype(ma).eq."yleft") then
  c                # Assuming velocity at bottom-face, this fix
  c                # preserves conservation in incompressible flow:
                   auxl(iaddaux(ma,i-nghost+1)) = aux(ma,i,nghost+1)
                 else
                   auxl(iaddaux(ma,i-nghost+1)) = aux(ma,i,nghost)
  
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                # -----------------------------------------------
                   ifine = (ic-1)*lratiox + nghost + l
                   auxl(iaddaux(ma,lind+1)) = aux(ma,ifine,mjtot-nghost+1)
                 else
  