anchor.diff:  c
   c        # update fluxes for use in AMR:
   c        # NOTE : We update ghost cell values for subcycling
  -         do 25 m=1,meqn
  -            do 25 i=2-mbc,mx+mbc
  +         do  m=1,meqn
  +            do  i=2-mbc,mx+mbc
                  fm(i,j,m) = fm(i,j,m) + faddm(i,m)
                  fp(i,j,m) = fp(i,j,m) + faddp(i,m)
                  gm(i,j,m) = gm(i,j,m) + gaddm(i,m,1)
                  gp(i,j,m) = gp(i,j,m) + gaddp(i,m,1)
                  gm(i,j+1,m) = gm(i,j+1,m) + gaddm(i,m,2)
                  gp(i,j+1,m) = gp(i,j+1,m) + gaddp(i,m,2)
  -   25          continue
  -   50    continue
  +            end do
  +        end do
  +      end do
   c
   c
   c
   c     # perform y sweeps
   c     ==================
   c
  -c
  -      do 100 i = 2-mbc, mx+mbc-1
  +
  +c     # Set corners for an y-sweep      
  +      sweep_dir = 1
  +      call clawpack46_fix_corners(mx,my,mbc,meqn,qold,sweep_dir,
  +     &             block_corner_count)
  +
  +c
  +      do  i = 2-mbc, mx+mbc-1
   c
   c        # copy data along a slice into 1d arrays:
            do m=1,meqn
  
positive.code: c        #
  c        # update fluxes for use in AMR:
  c
           do 75 m=1,meqn
              do 75 j=2-mbc,my+mbc
                 gm(i,j,m) = gm(i,j,m) + faddm(j,m)
                 gp(i,j,m) = gp(i,j,m) + faddp(j,m)
                 fm(i,j,m) = fm(i,j,m) + gaddm(j,m,1)
                 fp(i,j,m) = fp(i,j,m) + gaddp(j,m,1)
                 fm(i+1,j,m) = fm(i+1,j,m) + gaddm(j,m,2)
                 fp(i+1,j,m) = fp(i+1,j,m) + gaddp(j,m,2)
     75          continue
    100    continue
  c
  c
        return
        end
  
negative.code:       integer i0bpadq, iused
        double precision dtdx, dtdy, cfl1d
        integer m,i,j, ma, ixy
        integer sweep_dir
  
        common /comxyt/ dtcom,dxcom,dycom,tcom,icom,jcom
  
  