anchor.diff:  !**********************************************************************
   subroutine wprofile(wr,clock)
     implicit none
  +  ! Gives the w velocity profile along i=ngx/2, y = ngz/2-1
     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 volume_area(KLOK,XFN,elmnew, cap_i)
    IMPLICIT NONE
    ! Calculates and outputs the volume and surface area of a given capsule.
    INTEGER NFSIZE,NFSIZE2
    INTEGER i,KLOK, cap_i
    integer j1,j2,j3
  