anchor.diff:        return
         end
   
  -
  +c     #  See 'cubed_sphere_corners.ipynb'
         subroutine clawpack46_fix_corners(mx,my,mbc,meqn,q,sweep_dir,
        &             block_corner_count)
         implicit none
  
positive.code: 
                        if (use_b) then
  c                         # Transform involves B                
                            idata =  j1 + ibar(k) - jbar(k)
                            jdata = -i1 + ibar(k) + jbar(k)
                        else
  c                         # Transform involves B.transpose()             
                            idata = -j1 + ibar(k) + jbar(k)
                            jdata =  i1 - ibar(k) + jbar(k)
                        endif 
                        q(i1,j1,m) = q(idata,jdata,m)
                    end do   !! jbc
  
negative.code:             enddo
           enddo
  
           if (mcapa .gt. 0)  then
             do i = 1-mbc, mx+mbc
                 dtdx1d(i) = dtdx / aux(i,j,mcapa)
  