anchor.diff:          xcenterold(i) = xcenter(i)
           ycenterold(i) = ycenter(i)
           zcenterold(i) = zcenter(i)
  +        ! Save the n vectors, the six vectors from the center of the capsule
  +        ! to the six tracer points.
           call dumpnvec(klok, xfn(:,cap_n_start(i):cap_n_end(i)), xcenter(i), &
                ycenter(i), zcenter(i), nvec_i(:, i), i)
           call shape(h64,klok, i, xfn(1:3, cap_n_start(i): &
  
positive.code:   nfsize2 = 0
  
    ncap=$ncap$
    if (ncap > 0) then
       allocate(rad(ncap), xcenter(ncap), ycenter(ncap), zcenter(ncap), &
            xcenterold(ncap), ycenterold(ncap), zcenterold(ncap))
       allocate(cap_n_start(ncap), cap_n_end(ncap), cap_e_start(ncap), &
            cap_e_end(ncap))
  
negative.code:         xcenterold(i) = xcenter(i)
          ycenterold(i) = ycenter(i)
          zcenterold(i) = zcenter(i)
          ! Save the n vectors, the six vectors from the center of the capsule
          ! to the six tracer points.
          call dumpnvec(klok, xfn(:,cap_n_start(i):cap_n_end(i)), xcenter(i), &
               ycenter(i), zcenter(i), nvec_i(:, i), i)
          call shape(h64,klok, i, xfn(1:3, cap_n_start(i): &
  