[Gmsh] Error for (.pos) file with scalar points (SP)

[G Devaraj]

Hello Gmsh team,

I am using Gmsh version 2.16.0 to create unstructured grids for tsunami
simulations.
For that, I use a rectangular grid (via .pos file) to specify the
characteristic length.

*Q)* Is it possible to specify a view (.pos) file with scalar points (SP)?

If yes, then I request your help with the following error (pasted below) :
it seems there is a problem with the boundaries.
I get no errors when the grid is specified via SQ elements in the .pos file.
(I would like to use SP instead of SQ since I have a a very large domain
and the SQ mesh becomes quite bulky.)

My file for Gmsh is :
//////////////////////////////////////////////////////////
Mesh.CharacteristicLengthMin = 100;
Mesh.CharacteristicLengthMax = 100000;
Merge "Mesh_Geometry.geo1";
Merge "Mesh_Geometry.pos";
Background Mesh View[0];
//////////////////////////////////////////////////////////

My *Mesh_Geometry.geo1* file is :
//////////////////////////////////////////////////////////
// A. POINTS
Point(1) = {-328428, -222175, 0, 444350}; // Bottom-Left   corner of
rectangular domain
Point(2) = {328428, -222175, 0, 444350};  // Bottom-Right corner of
rectangular domain
Point(3) = {328428, 222175, 0, 444350};   // Top-Right      corner of
rectangular domain
Point(4) = {-328428, 222175, 0, 444350};  // Top-Left        corner of
rectangular domain

// B. LINES
Line(1) = {1,2}; // Bottom of rectangular domain
Line(2) = {2,3}; // Right    of rectangular domain
Line(3) = {3,4}; // Top      of rectangular domain
Line(4) = {4,1}; // Left      of rectangular domain

// C. SURFACES
Line Loop(5) = {1,2,3,4};
Plane Surface(6) = {5}; // The rectangle
//////////////////////////////////////////////////////////

My *Mesh_Geometry.pos* file for specifying characteristic length is of the
form:
//////////////////////////////////////////////////////////
View "BackgroundMesh" {
SP(-328428,-222175,0){7152};
SP(-328428,-221866,0){7156};
SP(-328428,-221557,0){7161};
 ...<etc..>...
SP(328428,221247,0){100000};
SP(328428,221557,0){100000};
SP(328428,221866,0){100000};
SP(328428,222175,0){100000};
};
//////////////////////////////////////////////////////////

I get errors of the following nature :
//////////////////////////////////////////////////////////
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]
Info    : Started on Thu Mar 30 17:54:03 2017
Info    : Reading '.././Input_Folder/Mesh_Geometry.geo'...
Info    : Reading '.././Input_Folder/Mesh_Geometry.geo1'...
Info    : Done reading '.././Input_Folder/Mesh_Geometry.geo1'
Info    : Reading '.././Input_Folder/Mesh_Geometry.pos'...
Info    : Done reading '.././Input_Folder/Mesh_Geometry.pos'
Info    : Done reading '.././Input_Folder/Mesh_Geometry.geo'
Info    : Finalized high order topology of periodic connections
Info    : Meshing 1D...
Info    : *Meshing curve 1* (Line)
Info    : No scalar element found containing point (0,*-222175*,0)
Info    : No scalar element found containing point (-164214,*-222175*,0)
 ...<etc..>...
Info    : No scalar element found containing point (310428,*-222175*,0)
Info    :* Meshing curve 2 *(Line)
Info    : No scalar element found containing point (*328428*,0,0)
Info    : No scalar element found containing point (*328428*,-111088,0)
 ...<etc..>...
Info    : No scalar element found containing point (*328428*,74510.8,0)
Info    : *Meshing curve 3* (Line)
Info    : No scalar element found containing point (0,*222175*,0)
Info    : No scalar element found containing point (164214,*222175*,0)
 ...<etc..>...
Info    : No scalar element found containing point (-314434,*222175*,0)
Info    : *Meshing curve 4* (Line)
Info    : No scalar element found containing point (*-328428*,0,0)
Info    : No scalar element found containing point (*-328428*,111088,0)
 ...<etc..>...
Info    : No scalar element found containing point (*-328428*,-200755,0)
Info    : Done meshing 1D (2.44 s)
Info    : Meshing 2D...
Info    : Meshing surface 6 (Plane, Delaunay)
Info    : No scalar element found containing point (-294207,-222175,0)
 ...<etc..>...
Info    : No scalar element found containing point (-311318,-204813,0)
Info    : Done meshing 2D (0.00637889 s)
Info    : 125 vertices 252 elements
//////////////////////////////////////////////////////////

Thanks and regards,
Devaraj.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/gmsh/attachments/20170330/dc4f3c44/attachment.html>