<div dir="ltr"><div><div><div>Hello Gmsh team,<br><br></div><div></div>I am using Gmsh version 2.16.0 to create unstructured grids for tsunami simulations.<br></div><div>For that, I use a rectangular grid (via .pos file) to specify the characteristic length.<br></div><div><br></div><div><b>Q)</b> Is it possible to specify a view (.pos) file with scalar points (SP)?<br><br></div><div>If yes, then I request your help with the following error (pasted below) : it seems there is a problem with the boundaries.<br></div><div>I get no errors when the grid is specified via SQ elements in the .pos file.<br></div><div></div><div>(I would like to use SP instead of SQ since I have a a very large domain and the SQ mesh becomes quite bulky.)<br></div><div><br></div><div>My file for Gmsh is :<br>//////////////////////////////////////////////////////////<br>Mesh.CharacteristicLengthMin = 100; <br>Mesh.CharacteristicLengthMax = 100000; <br>Merge "Mesh_Geometry.geo1"; <br>Merge "Mesh_Geometry.pos"; <br>Background Mesh View[0]; <br></div><div>//////////////////////////////////////////////////////////<br><br></div>My <b>Mesh_Geometry.geo1</b> file is :<br>//////////////////////////////////////////////////////////<br>// A. POINTS <br>Point(1) = {-328428, -222175, 0, 444350}; // Bottom-Left   corner of rectangular domain <br>Point(2) = {328428, -222175, 0, 444350};  // Bottom-Right corner of rectangular domain <br>Point(3) = {328428, 222175, 0, 444350};   // Top-Right      corner of rectangular domain <br>Point(4) = {-328428, 222175, 0, 444350};  // Top-Left        corner of rectangular domain <br><br>// B. LINES <br>Line(1) = {1,2}; // Bottom of rectangular domain <br>Line(2) = {2,3}; // Right    of rectangular domain <br>Line(3) = {3,4}; // Top      of rectangular domain <br>Line(4) = {4,1}; // Left      of rectangular domain <br><br>// C. SURFACES <br>Line Loop(5) = {1,2,3,4}; <br>Plane Surface(6) = {5}; // The rectangle <br>//////////////////////////////////////////////////////////<br><br>My <b>Mesh_Geometry.pos</b> file for specifying characteristic length is of the form:<br>//////////////////////////////////////////////////////////<br>View "BackgroundMesh" {<br>SP(-328428,-222175,0){7152}; <br>SP(-328428,-221866,0){7156}; <br>SP(-328428,-221557,0){7161}; <br></div><div> ...<etc..>...<br></div><div>SP(328428,221247,0){100000}; <br>SP(328428,221557,0){100000}; <br>SP(328428,221866,0){100000}; <br>SP(328428,222175,0){100000}; <br>}; <br>//////////////////////////////////////////////////////////<br><br></div>I get errors of the following nature :<br>//////////////////////////////////////////////////////////<div>Info    : Running './gmsh .././Input_Folder/Mesh_Geometry.geo -2 -o .././Input_Folder_MakranNERC/Mesh_ASCII.msh' [Gmsh 2.16.0, 1 node, max. 1 thread]<br>Info    : Started on Thu Mar 30 17:54:03 2017<br>Info    : Reading '.././Input_Folder/Mesh_Geometry.geo'...<br>Info    : Reading '.././Input_Folder/Mesh_Geometry.geo1'...<br>Info    : Done reading '.././Input_Folder/Mesh_Geometry.geo1'<br>Info    : Reading '.././Input_Folder/Mesh_Geometry.pos'...<br>Info    : Done reading '.././Input_Folder/Mesh_Geometry.pos'<br>Info    : Done reading '.././Input_Folder/Mesh_Geometry.geo'<br>Info    : Finalized high order topology of periodic connections<br>Info    : Meshing 1D...<br>Info    : <b>Meshing curve 1</b> (Line)<br>Info    : No scalar element found containing point (0,<b>-222175</b>,0)<br>Info    : No scalar element found containing point (-164214,<b>-222175</b>,0)<br> ...<etc..>...<br>Info    : No scalar element found containing point (310428,<b>-222175</b>,0)<br>Info    :<b> Meshing curve 2 </b>(Line)<br>Info    : No scalar element found containing point (<b>328428</b>,0,0)<br>Info    : No scalar element found containing point (<b>328428</b>,-111088,0)<br> ...<etc..>...<br>Info    : No scalar element found containing point (<b>328428</b>,74510.8,0)<br>Info    : <b>Meshing curve 3</b> (Line)<br>Info    : No scalar element found containing point (0,<b>222175</b>,0)<br>Info    : No scalar element found containing point (164214,<b>222175</b>,0)<br> ...<etc..>...<br>Info    : No scalar element found containing point (-314434,<b>222175</b>,0)<br>Info    : <b>Meshing curve 4</b> (Line)<br>Info    : No scalar element found containing point (<b>-328428</b>,0,0)<br>Info    : No scalar element found containing point (<b>-328428</b>,111088,0)<br> ...<etc..>...<br>Info    : No scalar element found containing point (<b>-328428</b>,-200755,0)<br>Info    : Done meshing 1D (2.44 s)<br>Info    : Meshing 2D...<br>Info    : Meshing surface 6 (Plane, Delaunay)<br>Info    : No scalar element found containing point (-294207,-222175,0)<br> ...<etc..>...<br>Info    : No scalar element found containing point (-311318,-204813,0)<br>Info    : Done meshing 2D (0.00637889 s)<br>Info    : 125 vertices 252 elements<br>//////////////////////////////////////////////////////////<br><br></div><div>Thanks and regards,<br></div><div>Devaraj.<br></div></div>