Numerical Advection Schemes in Two Dimensions

REFERENCES
REFERENCES
(1/2)*cv*(f(k,pp)-f(k,p));
fiminusm(k,p) = (1/2)*(f(km,p)+f(k,p)) -
(1/2)*cu*(f(k,p)-f(km,p));
fjminusm(k,p) = (1/2)*(f(k,pm)+f(k,p)) -
(1/2)*cv*(f(k,p)-f(k,pm));
end
f(k,p) = f(k,p) - cu*(fiplusm(k,p)-fiminusm(k,p)) -
cv*(fjplusm(k,p)-fjminusm(k,p)) +
cu*cv*(f(k,p)-f(k,pm)-f(km,p)+f(km,pm));
% Plot f at this timestep
contourf(x,y,f)
pause(0.1)
end
err = abs(f - h);
rss = sum(sum(err.^2));
rss3 = sum(sum(err.^3));
error2 = (rss/(n^2))^(1/2);
error3 = (rss3/(n^2))^(1/3);
errorinf = max(max(err));
errorstuff = [error2;error3;errorinf]
contourf(x,y,f);
xlabel(’x’)
ylabel(’y’)
19