Yt/master topo

Project.toml

@@ -5,6 +5,7 @@ version = "0.1.0"
 
 [deps]
 AMRVW = "233ec0c9-6b8d-4766-b9a2-a16e389fc38a"
+ArchGDAL = "c9ce4bd3-c3d5-55b8-8973-c0e20141b8c3"
 ArgParse = "c7e460c6-2fb9-53a9-8c5b-16f535851c63"
 BenchmarkTools = "6e4b80f9-dd63-53aa-95a3-0cdb28fa8baf"
 CUDA = "052768ef-5323-5732-b1bb-66c8b64840ba"
@@ -22,6 +23,9 @@ DifferentialEquations = "0c46a032-eb83-5123-abaf-570d42b7fbaa"
 DocStringExtensions = "ffbed154-4ef7-542d-bbb7-c09d3a79fcae"
 Documenter = "e30172f5-a6a5-5a46-863b-614d45cd2de4"
 ElasticArrays = "fdbdab4c-e67f-52f5-8c3f-e7b388dad3d4"
+GMT = "5752ebe1-31b9-557e-87aa-f909b540aa54"
+GeoInterface = "cf35fbd7-0cd7-5166-be24-54bfbe79505f"
+Geodesy = "0ef565a4-170c-5f04-8de2-149903a85f3d"
 Gmsh = "705231aa-382f-11e9-3f0c-b7cb4346fdeb"
 Gridap = "56d4f2e9-7ea1-5844-9cf6-b9c51ca7ce8e"
 GridapDistributed = "f9701e48-63b3-45aa-9a63-9bc6c271f355"
@@ -41,6 +45,7 @@ MPI = "da04e1cc-30fd-572f-bb4f-1f8673147195"
 MPIPreferences = "3da0fdf6-3ccc-4f1b-acd9-58baa6c99267"
 Makie = "ee78f7c6-11fb-53f2-987a-cfe4a2b5a57a"
 Metal = "dde4c033-4e86-420c-a63e-0dd931031962"
+NetCDF = "30363a11-5582-574a-97bb-aa9a979735b9"
 OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
 PProf = "e4faabce-9ead-11e9-39d9-4379958e3056"
 Parameters = "d96e819e-fc66-5662-9728-84c9c7592b0a"

meshes/gmsh_grids/hexa_TFI_DOMEX_coarse.geo

@@ -0,0 +1,65 @@
+nelemx = 25;
+nelemy = 25;
+nelemz = 10;
+
+xmin =  -25000;
+xmax =	 25000;
+ymin =  -25000;
+ymax =   25000;
+zmin =      0;
+zmax =  10000;
+gridsize = (xmax-xmin) / nelemx;
+
+Point(1) = {xmin, ymin, zmin, gridsize};
+Point(2) = {xmax, ymin, zmin, gridsize};
+Point(3) = {xmax, ymin, zmax, gridsize};
+Point(4) = {xmin, ymin, zmax, gridsize};
+
+Line(1) = {1, 2};
+Line(2) = {2, 3};
+Line(3) = {3, 4};
+Line(4) = {4, 1};
+
+
+npx = nelemx + 1;
+npy = nelemy + 1;
+npz = nelemz + 1;
+
+//Horizontal sides
+Transfinite Line {1, 3} = npx; //Ceil((xmax-xmin)/gridsize) Using Progression 1;
+//Vertical sides
+Transfinite Line {4, -2} = npz Using Progression 1.0;
+
+
+Line Loop(11) = {4, 1, 2, 3};
+Plane Surface(12) = {11};
+
+Transfinite Surface {12};
+Recombine Surface {12};
+
+surfaceVector = Extrude {0,(ymax-ymin),0} {
+  Surface{12};
+  Layers{nelemy};
+  Recombine;
+};
+//Coherence;
+
+  /* surfaceVector contains in the following order:
+     [0] - front surface (opposed to source surface)
+     [1] - extruded volume
+     [2] - bottom surface (belonging to 1st line in "Line Loop (6)")
+     [3] - right surface (belonging to 2nd line in "Line Loop (6)")
+     [4] - top surface (belonging to 3rd line in "Line Loop (6)")
+     [5] - left surface (belonging to 4th line in "Line Loop (6)")
+    */
+    Physical Surface("bottom") = {12};
+    Physical Volume("internal") = {1};
+    Physical Surface("back") = {25};
+    Physical Surface("front") = {33};
+    Physical Surface("left") = {21};
+    Physical Surface("right") = {29};
+    Physical Surface("top") = {34}; // from Plane Surface (6) ...
+  //+
+Show "*";
+//+
+Show "*";