1 %script compares old and
new detailed simulations
4 params.coarse_factor = 16;
5 model = advection_fv_output_opt_model(params);
7 model_old = advection_fv_output_opt_model_old(params);
9 model_data = gen_model_data(model);
12 M=[0,0,0;0,0,1;0,1,0;1,0,0;0,1,1;1,0,1;1,1,1;1,1,0];
17 model = set_mu(model,M(i,:));
18 model_old = set_mu(model_old,M(i,:));
20 sim_data_new = detailed_simulation(model, model_data);
22 sim_data_old = detailed_simulation(model_old, model_data);
26 difference = sim_data_new.U - sim_data_old.U;
29 nt = size(sim_data_new.U,2);
30 real_error = zeros(1,nt);
32 real_error(n) = sqrt(difference(:,n)'*A*difference(:,n));