anchor.diff:        enddo
   
         call compute_velocity_psi_comp(mx,my,mbc,dx,dy,
  -     &      t,xlower,ylower,aux,maux)
  +     &      t,blockno,xlower,ylower,aux,maux)
   
   c     # Compute the velocity in Cartesian coordinates;
   c     # project onto edge normals and scale. This
  
positive.code:       subroutine setaux_manifold(mbc,mx,my,
       &      xlower,ylower,dx,dy,
       &      maux,aux,xp,yp,zp,xd,yd,zd,
       &      xnormals,ynormals,edge_lengths,area)
        implicit none
  
        integer mbc, mx,my, meqn, maux
        double precision dx,dy, xlower, ylower, z
        double precision  aux(1-mbc:mx+mbc,1-mbc:my+mbc, maux)
  
  
negative.code:             xd2(1) = xlower + (i-1)*dx
              xd2(2) = ylower + (j-1)*dy
  
              call get_vel_psi_comp(blockno,xd1,xd2,dy,vn,t)
              aux(i,j,2) = vn
           enddo
        enddo
  