anchor.diff:        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
   
  
positive.code:             enddo
           enddo
  
           if (i .eq. 0) then
              do m = 1,meqn
                 if (block_corner_count(0) .eq. 3) then
                    do jbc = 1,mbc
                       q1d(1-jbc,m) = qold(0,jbc,m)
                    enddo
                 endif
                 if (block_corner_count(2) .eq. 3) then
                    do jbc = 1,mbc
                       q1d(my+jbc,m) = qold(0,my-jbc+1,m)
                    enddo
                 endif
              enddo
           else if (i .eq. mx+1) then
              do m = 1,meqn
                 if (block_corner_count(1) .eq. 3) then
                    do jbc = 1,mbc
                       q1d(1-jbc,m) = qold(mx+1,jbc,m)
                    enddo
                 endif
                 if (block_corner_count(3) .eq. 3) then
                    do jbc= 1,mbc
                       q1d(my+jbc,m) = qold(mx+1,my-jbc+1,m)
                    enddo
                 endif
              enddo
           endif
  
  c
           if (mcapa .gt. 0)  then
             do 71 j = 1-mbc, my+mbc
                 dtdy1d(j) = dtdy / aux(i,j,mcapa)
     71          continue
             endif
  c
           if (maux .gt. 0)  then
               do 72 ma=1,maux
                 do 72 j = 1-mbc, my+mbc
                   aux1(j,ma) = aux(i-1,j,ma)
                   aux2(j,ma) = aux(i,  j,ma)
                   aux3(j,ma) = aux(i+1,j,ma)
     72            continue
             endif
  c
  c
  c        # Store the value of i along this slice in the common block
  
negative.code:             enddo
           enddo
  
  c         if (j .le. 0) then
  c            do m = 1,meqn
  c               if (block_corner_count(0) .eq. 3) then
  c                  do ibc = 1,mbc
  c                     q1d(1-ibc,m) = qold(ibc,0,m)
  c                  enddo
  c               endif
  c               if (block_corner_count(1) .eq. 3) then
  c                  do ibc = 1,mbc
  c                     q1d(mx+ibc,m) = qold(mx-ibc+1,0,m)
  c                  enddo
  c               endif
  c            enddo
  c         else if (j .eq. my+1) then
  c            do m = 1,meqn
  c               if (block_corner_count(2) .eq. 3) then
  c                  do ibc = 1,mbc
  c                     q1d(1-ibc,m) = qold(ibc,my+1,m)
  c                  enddo
  c               endif
  c               if (block_corner_count(3) .eq. 3) then
  c                  do ibc = 1,mbc
  c                     q1d(mx+ibc,m) = qold(mx-ibc+1,my+1,m)
  c                  enddo
  c               endif
  c            enddo
  c         endif
  c
           if (mcapa .gt. 0)  then
             do i = 1-mbc, mx+mbc
                 dtdx1d(i) = dtdx / aux(i,j,mcapa)
             end do
           endif
  c
           if (maux .gt. 0)  then
               do  ma=1,maux
                 do  i = 1-mbc, mx+mbc
                   aux1(i,ma) = aux(i,j-1,ma)
                   aux2(i,ma) = aux(i,j  ,ma)
                   aux3(i,ma) = aux(i,j+1,ma)
                 end do
              end do
          endif
  c
  c
  c        # Store the value of j along this slice in the common block
  