anchor.diff:  !
   ! test case variables
   !
  -real(kreal), allocatable :: totalMassGHills(:)
  +real(kreal), allocatable :: totalMassGHills(:), tracerVar(:)
   real(kreal) :: sphereL2, sphereLinf, panelsLinf, particlesLinf, phiMax, phiMin, deltaPhi, phimax0, phimin0
  -real(kreal) :: mass0
  +real(kreal) :: mass0, var0
   
   !
   ! logging
  
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: 	deallocate(reference)
  endif
  deallocate(totalMassGHills)
  deallocate(tracerVar)
  call Delete(timekeeper)
  call Delete(remesh)
  if ( procrank == 0 ) call Delete(vtkOut)
  