anchor.diff:  subroutine calculateDF(clock, cap_i, xfn, my_nnode, lambda1, lambda2, &
        xcenter, ycenter, zcenter)
     implicit none
  +  ! This subroutine calculates the DF based on the points nearest to and
  +  ! farthest from the centroid.
  +  ! It also does other interesting things, like giving the stretch ratios
  +  ! at those points, and recording the directions to those points as mvec
     integer my_nnode
     double precision :: xfn(:, :)
     integer cap_i, i
  
positive.code: subroutine cellcenter(klok, xfn, my_nnode, cap_i, xcenter, &
       ycenter, zcenter)
    implicit none
    integer cap_i
    double precision :: xcenter, ycenter, zcenter
    double precision :: xfn(:,:)
  
negative.code:   !     Alex Szatmary
    !     8-24-06
    !     Records the coordinates of each solid node
    ! It's since been stretched to output the force on each solid node.
    !     Format: x, y, z
  
    implicit none
  