anchor.diff:  allocate(totalMassGHills(0:timesteps))
   totalMassGHills = 0.0_kreal
   mass0 = TotalMass(sphere, tracerID)
  +allocate(tracerVar(0:timesteps))
  +tracerVar = 0.0_kreal
  +var0 = TracerVariation(sphere, tracerID)
   
   
   sphereParticles => sphere%particles
  
positive.code: 				write(WRITE_UNIT_1,'(F24.15,A)') totalMassGHills(j), ' ; ...'
  			enddo
  			write(WRITE_UNIT_1,'(F24.15,A)') totalMassGHills(timesteps), ' ] ;'
  
  		endif
  		close(WRITE_UNIT_1)
  
  
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
  