Okay, figured it out: For 2D, you don't want "Physical Volume" but "Physical Surface".<div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Nov 2, 2012 at 3:51 PM, Nico Schlömer <span dir="ltr"><<a href="mailto:nico.schloemer@gmail.com" target="_blank">nico.schloemer@gmail.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi all,<div><br></div><div>I'm trying to define a 2D mesh with regions. To this end</div><div><br></div><div>=================== *snip* ===================<br>

</div><div><div>Function Rectangle</div><div>  // Points.</div>
<div>  cp1 = newp;</div><div>  Point(cp1) = {xmin,ymin,z,lcar};</div><div>  cp2 = newp;</div><div>  Point(cp2) = {xmax,ymin,z,lcar};</div><div>  cp3 = newp;</div><div>  Point(cp3) = {xmax,ymax,z,lcar};</div><div>  cp4 = newp;</div>


<div>  Point(cp4) = {xmin,ymax,z,lcar};</div><div>  </div><div>  // Lines.</div><div>  cl1 = newreg;</div><div>  Line(cl1) = {cp1,cp2};</div><div>  cl2 = newreg;</div><div>  Line(cl2) = {cp2,cp3};</div><div>  cl3 = newreg;</div>


<div>  Line(cl3) = {cp3,cp4};</div><div>  cl4 = newreg;</div><div>  Line(cl4) = {cp4,cp1};</div><div><br></div><div>  loop = newreg;</div><div>  Line Loop(loop) = {cl1,cl2,cl3,cl4};</div><div>  surf = newreg;</div><div>  Plane Surface(surf) = {loop};</div>


<div><br></div><div>  Physical Volume(Sprintf("rectangle %g", t)) = surf;</div><div><br></div><div>Return</div><div><br></div><div>// Build two rectangles.</div><div>lcar = 0.1;</div><div>z = 0;</div><div>ymin = -1.0;</div>


<div>ymax= 1.0;</div><div><br></div><div>t = 0;</div><div>xmin = -1.0;</div><div>xmax = 0.0;</div><div>Call Rectangle;</div><div><br></div><div>t = 1;</div><div>xmin = 0.0;</div><div>xmax = 1.0;</div><div>Call Rectangle;</div>


<div><br></div><div>Coherence;</div></div><div>=================== *snap* ===================</div><div><br></div><div>Meshing this works well, and when saving the data as MSH-file, I get</div><div><br></div><div>=================== *snip* ===================<br>


</div><div><div>$MeshFormat</div><div>2.2 0 8</div><div>$EndMeshFormat</div><div>$PhysicalNames</div><div>2</div><div>3 1 "rectangle 0"</div><div>3 2 "rectangle 1"</div><div>$EndPhysicalNames</div><div>


$Nodes</div><div>504</div><div>1 -1 -1 0</div><div>2 0 -1 0</div><div>3 0 1 0</div></div><div>[...]</div><div><div>$EndNodes</div><div>$Elements</div><div>1031</div><div>1 15 2 0 1 1</div></div><div>[...]</div><div><div>

1023 2 2 0 12 344 455 494</div>
<div>1024 2 2 0 12 355 476 464</div><div>1025 2 2 0 12 344 441 501</div><div>1026 2 2 0 12 344 494 441</div></div><div>[...]</div><div>=================== *snap* ===================<br></div><div><br></div><div>I'm wondering about the tags for the elements here. It seems that the tags are "0" and "12", one of which I would expect to correspond with a PhysicalName.</div>


<div>Can anyone enlighten me about the meaning of the tags here and how they correspond with physical volumes?</div><div><br></div><div>Cheers,</div><div>Nico</div>
</blockquote></div><br></div>