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: 		call LogMessage(exelog, TRACE_LOGGING_LEVEL, 'day = ', t/ONE_DAY)
  
  		write(vtkFile, '(A,I0.4,A)') trim(vtkRoot), frameCounter, '.vtk'
  		call UpdateFilename(vtkOut, vtkFile)
  		call VTKOutput(vtkOut, sphere)
  
  		frameCounter = frameCounter + 1
  	endif
  
negative.code: 
  	write(vtkRoot,'(A,A,A,A,A)') trim(outputDir), '/vtkOut/',trim(jobPrefix),trim(amrString),'_'
  	write(vtkFile,'(A,I0.4,A)') trim(vtkRoot),0,'.vtk'
  	
  	write(vtkMeshFile'(A,A,I0.4,A)') trim(vtkRoot), '_mesh_', 0, '.vtk'
  	
  	write(summaryFile,'(A,A,A,A)') trim(outputDir), trim(jobPrefix), trim(amrString), '_summary.txt'
  	write(datafile,'(A,A,A,A)') trim(outputDir), trim(jobPrefix), trim(amrstring), '_calculatedData.m'
  	call New(vtkOut, sphere, vtkFile, 'moving vortices')
  	call New(vtkMeshOut, sphere, vtkMeshFile, 'moving vortices')
  	call VTKOutput(vtkOut, sphere)
  	call VTKOutputMidpointRule(vtkMeshOut, sphere)
  endif
  
  !
  