anchor.diff:  !**********************************************************************
   subroutine wdump(wr,clock)
     implicit none
  +  ! Saves the w velocity at all nodes; deprecated and unused.
     integer, parameter :: lxng=$lngx$,lyng=$lngy$,lzng=$lngz$
     integer, parameter :: ngx=2**lxng,ngy=2**lyng,ngz=2**lzng
     integer i, j, k, clock
  
positive.code: !**********************************************************************
  subroutine wlindump(wlin,clock)
    implicit none
    integer, parameter :: lxng=$lngx$,lyng=$lngy$,lzng=$lngz$
    integer, parameter :: ngx=2**lxng,ngy=2**lyng,ngz=2**lzng
    integer i, clock
  
negative.code: subroutine cellcenter(klok, xfn, my_nnode, cap_i, xcenter, &
       ycenter, zcenter)
    implicit none
    ! Calculates the center of a given capsule and saves these coordinates
    integer cap_i
    double precision :: xcenter, ycenter, zcenter
    double precision :: xfn(:,:)
  