function err = err_fit_wtd_zeroed(x,D,G,N)

%used in var_spher.m

%Weighted sum-of-squares:

%           
%err = sum( (Ni/(Vi)^2)*(Gi - Vi)^2 )
%

if size(x,1) > 1
    
    x = x';
    
end

p = [x,0]; %add zero nugget

v = spher(p,D);

err = sum( (N./(v.^2)).*((G - v).^2) ) ; %Weighted Sum of Squares