anchor.diff:        subroutine setaux_manifold(mbc,mx,my,
  -     &      xlower,ylower,dx,dy,
  +     &      blockno,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
  +      integer mbc, mx,my, meqn, maux, blockno
         double precision dx,dy, xlower, ylower, z
         double precision  aux(1-mbc:mx+mbc,1-mbc:my+mbc, maux)
   
  
positive.code:             xd2(1) = xlower + (i-1)*dx
              xd2(2) = ylower + (j-1)*dy
  
              call get_vel_psi_comp(xd1,xd2,dy,vn,t)
              aux(i,j,2) = vn
           enddo
        enddo
  
negative.code: 
  
        subroutine compute_velocity_psi_comp(mx,my,mbc,
       &      dx,dy,t,blockno,xlower,ylower,aux,maux)
        implicit none
  
        integer mx,my,mbc,maux,blockno
        double precision dx,dy, t, xlower,ylower
  
        double precision xd1(2),xd2(2)
  