anchor.diff:  !
   ! program-specific variables
   !
  -
  +real(kreal) :: test1method1Error(4), test1method2Error(4), test1method3Error(4)
   
   !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   !	PROGRAM START
  
positive.code: 	real(kreal) :: A(31,10), AT(10,31), ATA(10,10), ATAInvAT(10,31), coeffs(10), scalarData(31), datalocs(3,31)
  	integer(kint) :: adjPanels(8), nAdj, edgeList(8), vertList(8), nVerts, nearbyParticles(20), nNear
  	logical(klog) :: isNewPoint
  
  	aParticles => aMesh%particles
  	aPanels => aMesh%panels
  
  	do j = 1, aPanels%N
  		if ( .NOT. aPanels%hasChildren(j) ) then
  
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
  
  !
  