anchor.diff:  		! delete objects associated with old mesh
   		!
   		call Delete(timekeeper)
  -		if ( procrank == 0 ) call Delete(vtkOUt)
  +		if ( procrank == 0 ) then
  +			call Delete(vtkOUt)
  +			call Delete(vtkMeshOut)
  +		endif
   		!
   		! create new associated objects for new mesh
   		!
   		call New(timekeeper, sphere, numProcs)
  -		if ( procRank == 0 ) call New(vtkOut, sphere, vtkFile, 'Gaussian hills advection')
  +		if ( procRank == 0 ) then
  +			call New(vtkOut, sphere, vtkFile, 'moving vortices')
  +			call New(vtkMeshOut, sphere, vtkMeshFile, 'moving vortices')
  +		endif
   		sphereParticles => sphere%particles
   		spherePanels => sphere%panels
   	endif ! remesh
  
positive.code: deallocate(phiMin)
  call Delete(timekeeper)
  call Delete(remesh)
  if ( procrank == 0 ) call Delete(vtkOut)
  call Delete(sphere)
  call Delete(testCaseTracer)
  call Delete(exeLog)
  
negative.code: 		call LogMessage(exelog, TRACE_LOGGING_LEVEL, 'day = ', t/ONE_DAY)
  
  		write(vtkFile, '(A,I0.4,A)') trim(vtkRoot), frameCounter, '.vtk'
  		write(vtkMeshFile, '(A,A, I0.4,A)') trim(vtkRoot), '_mesh_', frameCounter, '.vtk'
  
  		call UpdateFilename(vtkOut, vtkFile)
  		call UpdateFilename(vtkMeshOut, vtkMeshFile)
  		
  		call VTKOutput(vtkOut, sphere)
  		call VTKOutputMidpointRule(vtkMeshOut, sphere)
  
  		frameCounter = frameCounter + 1
  	endif
  