global N N_p N_fp N_faces NODETOL;
global P E N_v K r s x y;
global V invV M invM D_r D_s D_r_weak D_s_weak LIFT;
global r_x r_y s_x s_y J;
global n_x n_y Fmask Fx Fy surfJ Fscale; % surfJ contains lenghts for all the faces of elements
global EtoE EtoF;
global vmapM mapM vmapP mapP vmapB mapB;
global alpha_flux;
global dt N_steps;
global rk4a rk4b rk4c;

global mesh_size_parametr;

% Advection2D problem related variables
global a_advec b_advec RHS IC BC;
global vmapI mapI vmapO mapO;

% LEE2D problem related variables
global poisson_const rho_0 u_0 v_0 p_0 a_0 M_0 A_1 A_2;
global eps_1 eps_2 eps_3 kappa_1 kappa_2 kappa_3;
global rho_IC u_IC v_IC p_IC;
global rho_RHS u_RHS v_RHS p_RHS;
global C_flux; % Lax-Friedrichs flux
global TDplusTinv_flux_cell TDminusTinv_flux_cell; % upwind flux
global x_l x_u y_l y_u;
global vmapW mapW vmapF mapF;
global elements_Euler;
global sigma_sponge D_sponge elements_sponge; % for sponge zone implemented as NRBC
global sigma_x_PML sigma_y_PML D_PML coef_PML; % for PML implemented as NRBC
global sigma_x_PML_plus_sigma_y_PML sigma_x_PML_times_sigma_y_PML sigma_y_PML_times_n_x sigma_x_PML_times_n_y;
global elements_PML vmapPML_M mapPML_M vmapPML_P mapPML_P vmapW_PML mapW_PML;
global rho_states u_states v_states p_states; % intermidiate results in time
global rho_wall u_wall v_wall p_wall; % olution on wall to test wall BC
global rho_boundary_R u_boundary_R v_boundary_R p_boundary_R rho_boundary_T u_boundary_T v_boundary_T p_boundary_T rho_boundary_L u_boundary_L v_boundary_L p_boundary_L;