anchor.diff:  		newSphere%mesh%particles%y0(i) = x0(2)
   		newSphere%mesh%particles%z0(i) = x0(3)
   		
  +		newSphere%density%scalar(i) = InterpolateScalar( lon, lat, self%densitySource, oldSphere%mesh, &
  +											self%delTri, oldSphere%density)
  +		
   		do j = 1, size(newSphere%tracers)
   			if ( newSphere%tracers(j)%nDim == 1 ) then
   				newSphere%tracers(j)%scalar(i) = InterpolateScalar(lon, lat, self%tracerSource(j), oldSphere%mesh, &
   					self%delTri, oldSphere%tracers(j) )
   			else
  -				vec = InterpolateVector(lon, lat, self%tracerSource(j), oldSphere%mesh, self%delTri, oldSphere%tracers(j))
  +				vec = InterpolateVector(lon, lat, self%tracerSource(j), oldSphere%mesh, self%delTri,&
  +					 oldSphere%tracers(j))
   				newSphere%tracers(j)%xComp(i) = vec(1)
   				newSphere%tracers(j)%yComp(i) = vec(2)
   				newSphere%tracers(j)%zComp(i) = vec(3)
  
positive.code: 	endif!AMR
  	
  	call SetVelocityOnMesh( newSphere, velFn, t )
  	call SetDivergenceOnMesh( newSphere, divFn, t )
  end subroutine
  
  !> @brief Performs a remesh/remap of an LPM @ref SphereBVE simulation using indirect interpolation for vorticity variables and up to two scalar tracer variables in a BVE mesh.  Additional tracers are directly interpolated.
  
negative.code: 									 newSphere%mesh%particles%z(j))
  					lat = Latitude( newSphere%mesh%particles%x(j), newSphere%mesh%particles%y(j), &
  									newSphere%mesh%particles%z(j))
  					
  					x0 = InterpolateLagParam( lon, lat, self%lagParamSource, oldSphere%mesh, self%delTri)
  					newSphere%mesh%particles%x0(j) = x0(1)
  					newSphere%mesh%particles%y0(j) = x0(2)
  					newSphere%mesh%particles%z0(j) = x0(3)				
  									
  					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, &
  