clc; clear; close all;

Globals2D;

tic;
Problem;

N = 6; % degree of polynomial representing solution within an element (local solution representation)

StartUp2D; % computing local operators and mesh information

NonreflectingBCStartup;

%LEE2DNumericalFluxStartup_Problem6;

[rho,u,v,p] = Timestepping_PML(T); % time-stepping for the ODE system
%[rho,u,v,p] = Timestepping_sponge(T); % time-stepping for the ODE system
computation_time = toc;