anchor.diff:  	deallocate(reference)
   endif
   deallocate(totalMassGHills)
  +deallocate(tracerVar)
   call Delete(timekeeper)
   call Delete(remesh)
   if ( procrank == 0 ) call Delete(vtkOut)
  
positive.code: allocate(totalMassGHills(0:timesteps))
  totalMassGHills = 0.0_kreal
  mass0 = TotalMass(sphere, tracerID)
  
  
  sphereParticles => sphere%particles
  
negative.code: !
  ! test case variables
  !
  real(kreal), allocatable :: totalMassGHills(:), tracerVar(:)
  real(kreal) :: sphereL2, sphereLinf, panelsLinf, particlesLinf, phiMax, phiMin, deltaPhi, phimax0, phimin0
  real(kreal) :: mass0, var0
  
  !
  ! logging
  