Skip to content
Snippets Groups Projects
Commit 19dec53f authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Add a helper for velocity stepping plots

parent 7f737a0b
No related branches found
No related tags found
No related merge requests found
## State evolution
load ../data;
graphics_toolkit fltk;
path(path, '../../../');
plot_state
## Comparison of steady-state convergence rates, to and fro
# 600 time steps
load ../data;
graphics_toolkit fltk;
path(path, '../../../');
plot(1:144,A(121:(120+144),1)-A(121,1),1:144,-interp1(1:240, A(361:600,1), linspace(1,240,144))+A(361,1));
axis tight;
legend('decline (first)','incline (second)');
# 6000 time steps
load ../data;
graphics_toolkit fltk;
path(path, '../../../');
plot(1:1440,A(1201:(1200+1440),1)-A(1201,1),1:1440,-interp1(1:2400, A(3601:6000,1), linspace(1,2400,1440))+A(3601,1));
axis tight;
legend('decline (first)','incline (second)');
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment