add new geo files

Meshes/dmr_hybrid.geo

@@ -66,5 +66,4 @@ Recombine Surface(surf4);
 //Recombine Surface (surf3);
 
 
-Physical Surface(0) = {surf1,surf2,surf3,surf4};
-Mesh.MshFileVersion =2;
+Physical Surface(0) = {surf1,surf2,surf3,surf4};
\ No newline at end of file

Meshes/grav_quad.geo

+// square
+
+lc = 0.1; // Mesh parameter
+R=0.5;
+
+// Points
+p1 = newp;
+Point(p1) = {-R, -R, 0.0, lc};
+p2 = newp;
+Point(p2) = {R, -R, 0.0, lc};
+p3 = newp;
+Point(p3) = {R, R, 0.0, lc};
+p4=newp;
+Point(p4)={-R,R,0.0, lc};
+
+
+// lines
+c1 = newreg;
+Line(c1) = {p1, p2};
+c2 = newreg;
+Line(c2) = {p2, p3};
+c3 = newreg;
+Line(c3) = {p3, p4};
+c4=newreg;
+Line(c4)={p4,p1};
+
+
+// Line loop
+ll1 = newreg;
+Line Loop(ll1) = {c1,c2,c3,c4};
+
+// Plane surfaces
+surf1 = news;
+Plane Surface(surf1) = {ll1}; 
+
+
+// Set boundary tags
+//Physical Line("inflow") = {c4};          // inflow
+//Physical Line("outflow") = {c2};          // outflow
+Physical Line("outflow") = {c1,c2,c4,c3}; // wall
+Physical Surface(0) = {surf1};
+
+n1=128;
+n2=128;
+
+// Make one square structured.
+Transfinite Line {c1,c3}=n1 Using Progression 1;
+Transfinite Line {c4,c2}=n2 Using Progression 1;
+Transfinite Surface {surf1}={1,2,3,4};
+
+
+
+//Recombine Surface {surf1};
+
+

Meshes/shockVortex.geo

+// square
+
+lc = 0.01; // Mesh parameter
+R=1.;R2=1.; R3=0.;
+
+// Points
+p1 = newp;
+Point(p1) = {-0., 0, 0.0, lc};
+p2 = newp;
+Point(p2) = {R, 0, 0.0, lc};
+p3 = newp;
+Point(p3) = {R, R2, 0.0, lc};
+p4=newp;
+Point(p4)={-0.,R2,0.0, lc};
+p5=newp;
+Point(p5)={0.5,0.,0.,lc};
+p6=newp;
+Point(p6)= {0.5,R2,0.,lc};
+
+
+// lines
+c1 = newreg;
+Line(c1) = {p1, p5};
+c2 = newreg;
+Line(c2) = {p5, p6};
+c3 = newreg;
+Line(c3) = {p6, p4};
+c4=newreg;
+Line(c4)={p4,p1};
+c5=newreg;
+Line(c5)={p5,p2};
+c6=newreg;
+Line(c6)={p2,p3};
+c7=newreg;
+Line(c7)={p3,p6};
+
+
+// Line loop
+ll1 = newreg;
+Line Loop(ll1) = {c1,c2,c3,c4};
+ll2=newreg;
+Line Loop(ll2)={c5,c6,c7,-c2};
+
+// Plane surfaces
+surf1 = news;
+Plane Surface(surf1) = {ll1};
+surf2=news;
+Plane Surface(surf2)={ll2};
+
+ny=100;
+nx=50;
+Transfinite Line{c4,c2}=ny Using Progression 1;
+Transfinite Line{c1,c3}=nx Using Progression 1;
+Transfinite Surface{surf1};
+Recombine Surface(surf1);
+
+// Set boundary tags
+Physical Line("inflow") = {c4};          // inflow
+Physical Line("outflow") = {c6};          // outflow
+Physical Line("Wall") = {c1,c5,c3,c7}; // wall
+Physical Surface(0) = {surf1,surf2};
+

Meshes/shockVortexLong.geo

+// square
+
+lc = 0.01; // Mesh parameter
+R=1.;R2=1.; R3=0.;
+
+// Points
+p1 = newp;
+Point(p1) = {-0.5, 0, 0.0, lc};
+p2 = newp;
+Point(p2) = {R, 0, 0.0, lc};
+p3 = newp;
+Point(p3) = {R, R2, 0.0, lc};
+p4=newp;
+Point(p4)={-0.5,R2,0.0, lc};
+p5=newp;
+Point(p5)={0.5,0.,0.,lc};
+p6=newp;
+Point(p6)= {0.5,R2,0.,lc};
+
+
+// lines
+c1 = newreg;
+Line(c1) = {p1, p5};
+c2 = newreg;
+Line(c2) = {p5, p6};
+c3 = newreg;
+Line(c3) = {p6, p4};
+c4=newreg;
+Line(c4)={p4,p1};
+c5=newreg;
+Line(c5)={p5,p2};
+c6=newreg;
+Line(c6)={p2,p3};
+c7=newreg;
+Line(c7)={p3,p6};
+
+
+// Line loop
+ll1 = newreg;
+Line Loop(ll1) = {c1,c2,c3,c4};
+ll2=newreg;
+Line Loop(ll2)={c5,c6,c7,-c2};
+
+// Plane surfaces
+surf1 = news;
+Plane Surface(surf1) = {ll1};
+surf2=news;
+Plane Surface(surf2)={ll2};
+
+n1=100;
+n2=100;
+Transfinite Line{c4,c2}=n1 Using Progression 1;
+Transfinite Line{c1,c3}=n2 Using Progression 1;
+Transfinite Surface{surf1};
+Recombine Surface(surf1);
+
+// Set boundary tags
+Physical Line("inflow") = {c4};          // inflow
+Physical Line("outflow") = {c6};          // outflow
+Physical Line("Wall") = {c1,c5,c3,c7}; // wall
+Physical Surface(0) = {surf1,surf2};
+

Meshes/shockVortexQ.geo

+// square
+
+lc = 0.01; // Mesh parameter
+R=1.;R2=1.;
+
+// Points
+p1 = newp;
+Point(p1) = {0, 0, 0.0, lc};
+p2 = newp;
+Point(p2) = {R, 0, 0.0, lc};
+p3 = newp;
+Point(p3) = {R, R2, 0.0, lc};
+p4=newp;
+Point(p4)={0,R2,0.0, lc};
+p5=newp;
+Point(p5)={0.5,0.,0.,lc};
+p6=newp;
+Point(p6)= {0.5,R2,0.,lc};
+
+
+// lines
+c1 = newreg;
+Line(c1) = {p1, p5};
+c2 = newreg;
+Line(c2) = {p5, p6};
+c3 = newreg;
+Line(c3) = {p6, p4};
+c4=newreg;
+Line(c4)={p4,p1};
+c5=newreg;
+Line(c5)={p5,p2};
+c6=newreg;
+Line(c6)={p2,p3};
+c7=newreg;
+Line(c7)={p3,p6};
+
+
+// Line loop
+ll1 = newreg;
+Line Loop(ll1) = {c1,c2,c3,c4};
+ll2=newreg;
+Line Loop(ll2)={c5,c6,c7,-c2};
+
+// Plane surfaces
+surf1 = news;
+Plane Surface(surf1) = {ll1};
+surf2=news;
+Plane Surface(surf2)={ll2};
+
+n1=100;
+n2=50;
+Transfinite Line{c4,c2}=n1 Using Progression 1;
+Transfinite Line{c1,c3}=n2 Using Progression 1;
+Transfinite Surface{surf1};
+Recombine Surface(surf1);
+
+Transfinite Line{c6,c2}=n1 Using Progression 1;
+Transfinite Line{c5,c7}=n2 Using Progression 1;
+Transfinite Surface{surf2};
+Recombine Surface(surf2);
+
+// Set boundary tags
+Physical Line("inflow") = {c4};          // inflow
+Physical Line("outflow") = {c6};          // outflow
+Physical Line("Wall") = {c1,c5,c3,c7}; // wall
+Physical Surface(0) = {surf1,surf2};
+

Meshes/shockVortex_sym.geo

+// square
+
+lc = 0.01; // Mesh parameter
+R=1.;R2=1.; R3=0.;
+
+// Points
+p1 = newp;
+Point(p1) = {-0., 0, 0.0, lc};
+p2 = newp;
+Point(p2) = {R, 0, 0.0, lc};
+p3 = newp;
+Point(p3) = {R, R2, 0.0, lc};
+p4=newp;
+Point(p4)={-0.,R2,0.0, lc};
+p5=newp;
+Point(p5)={0.5,0.,0.,lc};
+p6=newp;
+Point(p6)= {0.5,R2,0.,lc};
+
+
+// lines
+c1 = newreg;
+Line(c1) = {p1, p5};
+c2 = newreg;
+Line(c2) = {p5, p6};
+c3 = newreg;
+Line(c3) = {p6, p4};
+c4=newreg;
+Line(c4)={p4,p1};
+c5=newreg;
+Line(c5)={p5,p2};
+c6=newreg;
+Line(c6)={p2,p3};
+c7=newreg;
+Line(c7)={p3,p6};
+
+
+// Line loop
+ll1 = newreg;
+Line Loop(ll1) = {c1,c2,c3,c4};
+ll2=newreg;
+Line Loop(ll2)={c5,c6,c7,-c2};
+
+// Plane surfaces
+surf1 = news;
+Plane Surface(surf1) = {ll1};
+surf2=news;
+Plane Surface(surf2)={ll2};
+
+ny=100;
+nx=50;
+Transfinite Line{c2,c6}=ny Using Progression 1;
+Transfinite Line{c5,c7}=nx Using Progression 1;
+Transfinite Surface{surf2};
+Recombine Surface(surf2);
+
+// Set boundary tags
+Physical Line("inflow") = {c4};          // inflow
+Physical Line("outflow") = {c6};          // outflow
+Physical Line("Wall") = {c1,c5,c3,c7}; // wall
+Physical Surface(0) = {surf1,surf2};
+