> On 22 Feb 2018, at 12:56, Sathyanarayan Rao <sathyanarayan....@uclouvain.be> > wrote: > > Dear gmsh team, > > Is it not possible to use both line loop type surface definition and open > cascade rectangle. I always have problem with using both.
Sure you can: but you need to put SetFactory("OpenCASCADE") at the begining of your file: otherwise you create the first entities with the built-in CAD kernel ; and the rest (after SetFactory) with open cascade. > > For example in the below code, surface 105 is unrecognized ! > > // mesh size description > > cl_1 = 0.3; > cl_2 = 0.02; > cl3 = 0.02; > Mesh.CharacteristicLengthMax = 40; > extend = 15; > > // define inversion region > > Point(1) = {-0.05, 0.005, 0, cl_2}; > Point(2) = {7.25, 0.005, 0, cl_2}; > Point(3) = {7.25, -1.45, 0, cl_2}; > Point(4) = {-0.05, -1.45, 0, cl_2}; > > // define bounding box > > Point(5) = {-0.05-extend, 0.005, 0, cl_1}; > Point(6) = {7.25+extend, 0.005, 0, cl_1}; > Point(7) = {7.25+extend, -1.45-extend, 0, cl_1}; > Point(8) = {-0.05-extend, -1.45-extend, 0, cl_1}; > > Line(1) = {1, 2}; > Line(2) = {2, 6}; > Line(3) = {6, 7}; > Line(4) = {7, 8}; > Line(5) = {8, 5}; > Line(6) = {5, 1}; > Line(7) = {1, 4}; > Line(8) = {4, 3}; > Line(9) = {3, 2}; > > // Definition of regions > Line Loop(101) = {2, 3, 4, 5, 6, 7, 8, 9}; > Plane Surface(102) = {101}; // Outer box > Line Loop(105) = {1, -9, -8, -7}; > > > // geological body > > SetFactory("OpenCASCADE"); > Rectangle(1) = {0.1, -0.2, 0, 0.5, 0.1, 0}; > Rectangle(2) = {0.1, -1.0, 0, 0.3, 0.7, 0}; > Rectangle(3) = {3.1, -0.6, 0, 1, 0.5, 0}; > Rectangle(4) = {4.2, -1.2, 0, 1, 0.5, 0}; > Rectangle(5) = {5.6, -0.9, 0, 1, 0.5, 0}; > Plane Surface(110) = {105,1,2,3,4,5}; // Inner box without geological body > > > // Definition of boundaries > Physical Line(1) = {6, 2, 1}; // Free surface > Physical Line(2) = {3,4, 5}; // Mixed boundary conditions > > // Definition of physical regions > Physical Surface(1) = {102}; // No-inversion region > Physical Surface(2) = {110}; // Inversion region > Physical Surface(3) = {1}; // Geological body 1 > Physical Surface(4) = {2}; // Geological body 2 > Physical Surface(5) = {3}; // Geological body 3 > Physical Surface(6) = {4}; // Geological body 4 > Physical Surface(7) = {5}; // Geological body 5 > > > > Best Regards, > > Sathyanarayan Rao, PhD student > Earth and Life Institute/Environmental Sciences (ELI-e) > Université catholique de Louvain > c.037, Croix du Sud 2, Louvain-la-Neuve, Belgium > Phone: 010473827 ( intercom 73827) > _______________________________________________ > gmsh mailing list > gmsh@onelab.info <mailto:gmsh@onelab.info> > http://onelab.info/mailman/listinfo/gmsh > <http://onelab.info/mailman/listinfo/gmsh> — Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info
_______________________________________________ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh