anchor.diff:  	write(logstring,'(A,F15.4,A)') 'singular kernel elapsed time = ', MPI_WTIME() - testClock, ' seconds.'
   	call LogMessage(exeLog, TRACE_LOGGING_LEVEL,'test 1, method 3 : ', trim(logString) )
   
  -	call VTKOutput(vtkOut, sphere)
  +	if ( procRank == 0 ) then
  +		write(vtkFile,'(A,A,I0.4,A)') trim(vtkRoot), 'test1', 2, '.vtk'
  +		call UpdateFilename(vtkOut, vtkFile)
  +		call UpdateTitle(vtkOut, 'singular kernels')
  +		call VTKOutput(vtkOut, sphere)
  +	endif
   call EndSection(exeLog)
   
   !
   !	test function 2
   !
  -call TestSphericalHarmonic54(sphere)
  -call UpdateTitle(vtkOut, 'spherical harmonic 5 4')
  +!call TestSphericalHarmonic54(sphere)
  +call TestSphericalHarmonic52(sphere)
  +call UpdateTitle(vtkOut, 'spherical harmonic 5 2')
   
   
   !
  
positive.code: 			!	data collected. ready for least squares
  			!
  
  
  		endif
  	enddo
  end subroutine
  
  
  end program
  
negative.code: spherePanels => sphere%panels
  
  !	TO DO : AMR
  write(amrString,'(A,I1)') 'quadNest', initNest
  
  if ( procRank == 0 ) then
  	write(vtkRoot, '(5A)') trim(outputDir), '/vtkOut/', trim(jobPrefix), trim(amrString), '_'
  	write(vtkFile,'(A,A,I0.4,A)') trim(vtkRoot),'test1', 0, '.vtk'
  	call New(vtkOut, sphere, vtkFile, 'quadratic polynomials in r3')
  endif
  
  !
  