anchor.diff:        !   matrix (transposed) as the columns of A, and zero out
         !   the lower triangle (upper triangle of A) with Givens
         !   rotations.
  +
         do  i = 1, 6
            ed = eds (i)
            !Apply the rotation for this point
  -         call aplyr (mesh%edhx(ed)%c%p(1), mesh%edhx(ed)%c%p(2), &
  -              mesh%edhx(ed)%c%p(3), cx, sx, cy, sy, xp, yp, zp)
  +         if (var%pos == 3)then
  +            call aplyr (mesh%edhx(ed)%c%p(1), mesh%edhx(ed)%c%p(2), &
  +                 mesh%edhx(ed)%c%p(3), cx, sx, cy, sy, xp, yp, zp)
  +         else
  +            call aplyr (mesh%ed(ed)%c%p(1), mesh%ed(ed)%c%p(2), &
  +                 mesh%ed(ed)%c%p(3), cx, sx, cy, sy, xp, yp, zp)
  +         endif
   
            !Define the weight of the edge
            wt = 1. / (1. - zp) - rinv
  
positive.code:     do i=7, n
         ed = eds (i)
         !Apply the rotation for this point
         call aplyr (mesh%edhx(ed)%c%p(1), mesh%edhx(ed)%c%p(2), &
              mesh%edhx(ed)%c%p(3), cx, sx, cy, sy, xp, yp, zp)
  
         !Define the weight of the edge
         wt = 1. / (1. - zp) - rinv
  
negative.code:           vector_reconstruct=vecrecon_pered(pindtmp, var, mesh)
         case("dtred") !Dual triangle recons for edges (tangent vector reconstruction)
            vector_reconstruct=vecrecon_dtred(pindtmp, var, mesh)
         case("lsq") !Least Square fit
            vector_reconstruct=vecrecon_lsq (p, var, stencil, mesh, pindtmp)
         case default
            print*, "VECTOR RECON ERROR: Don't know how to do this interpolation:", &
                 kindinterpol, var%pos
  