anchor.diff:  											vortFlagFn1, tol1, desc1, nParticlesBefore, nParticlesAfter)
   				elseif ( refineVariableCount == 2 ) then
   					if ( present(field2) ) then
  -						call IterateMeshRefinementTwoVariables( refine, newSphere%mesh, newSphere%relVort, vortFlagFn1, &
  +						call IterateMeshRefinementTwoVariables(refine, newSphere%mesh, newSphere%relVort, vortFlagFn1, &
   								tol1, desc1, field2, flagFn2, tol2, desc2, nParticlesBefore, nParticlesAfter)
   					else
  -						call IterateMeshRefinementTwoVariables( refine, newSphere%mesh, newSphere%relVort, vortFlagFn1, &
  +						call IterateMeshRefinementTwoVariables(refine, newSphere%mesh, newSphere%relVort, vortFlagFn1, &
   								tol1, desc1, newSphere%relVort, flagFn2, tol2, desc2, nParticlesBefore, nParticlesAfter)
   					endif
   				endif
  
positive.code: 									 newSphere%mesh%particles%z(j))
  					lat = Latitude( newSphere%mesh%particles%x(j), newSphere%mesh%particles%y(j), &
  									newSphere%mesh%particles%z(j))
  					newSphere%relVort%scalar(j) = InterpolateScalar( lon, lat, self%relVortSource, oldSphere%mesh, &
  																	 self%delTri, oldSphere%relVort)
  					newSphere%absVort%scalar(j) = InterpolateScalar( lon, lat, self%absVortSource, oldSphere%mesh, &
  
negative.code: 	real(kreal), dimension(3) :: x0, vec
  	type(RefineSetup) :: refine
  	integer(kint) :: nParticlesBefore, nParticlesAfter
  	integer(kint) :: amrVarCount
  	
  	remeshCounter = remeshCounter + 1
  	
  	!
  	!	interpolate to new uniform mesh
  	!
  	newSphere%density%N = newSphere%mesh%particles%N
  	do i = 1, size(newSphere%tracers)
  		newSphere%tracers(i)%N = newSphere%mesh%particles%N
  	enddo
  