fix xyz

Author: simbilod

gplugins/gmsh/get_mesh.py

@@ -23,6 +23,7 @@ def get_mesh(
     type: str,
     layer_stack: LayerStack,
     layer_physical_map: dict | None = None,
+    layer_meshbool_map: dict | None = None,
     z: float | None = None,
     xsection_bounds=None,
     wafer_padding: float = 0.0,
@@ -37,6 +38,7 @@ def get_mesh(
         type: one of "xy", "uz", or "3D". Determines the type of mesh to return.
         layer_stack: LayerStack object containing layer information.
         layer_physical_map: by default, physical are tagged with layername; this dict allows you to specify custom mappings.
+        layer_meshbool_map: by default, all polygons on layer_stack layers are meshed; this dict allows you set True of False to the meshing of given layers.
         z: used to define z-slice for xy meshing.
         xsection_bounds: used to define in-plane line for uz meshing.
         wafer_padding: padding beyond bbox to add to WAFER layers.
@@ -85,6 +87,16 @@ def get_mesh(
             if layer_name not in layer_physical_map.keys():
                 layer_physical_map[layer_name] = layer_name
 
+    # Default meshing flags (all True)
+    if layer_meshbool_map is None:
+        layer_meshbool_map = {}
+        for layer_name in layer_stack.layers.keys():
+            layer_meshbool_map[layer_name] = True
+    else:
+        for layer_name in layer_stack.layers.keys():
+            if layer_name not in layer_physical_map.keys():
+                layer_meshbool_map[layer_name] = True
+
     if type == "xy":
         if z is None:
             raise ValueError(

gplugins/gmsh/xyz_mesh.py

@@ -75,6 +75,7 @@ def define_edgeport(
 def define_prisms(
     layer_polygons_dict: dict,
     layer_stack: LayerStack,
+    layer_physical_map: dict,
     model: Any,
     resolutions: dict,
     scale_factor: float = 1,
@@ -84,6 +85,7 @@ def define_prisms(
     Args:
         layer_polygons_dict: dictionary of polygons for each layer
         layer_stack: gdsfactory LayerStack to parse
+        layer_physical_map: map layer names to physical names
         model: meshwell Model object
         resolutions: Pairs {"layername": {"resolution": float, "distance": "float}} to roughly control mesh refinement.
         scale_factor: scaling factor to apply to the polygons (default: 1)
@@ -120,7 +122,9 @@ def define_prisms(
                 model=model,
                 resolution=resolutions.get(layername, None),
                 mesh_order=buffered_layer_stack.layers.get(layername).mesh_order,
-                physical_name=layername,
+                physical_name=layer_physical_map[layername]
+                if layername in layer_physical_map
+                else layername,
             )
         )
 
@@ -255,6 +259,7 @@ def xyz_mesh(
         model=model,
         scale_factor=global_scaling_premesh,
         resolutions=resolutions,
+        layer_physical_map=layer_physical_map,
     )
 
     # Add edgeports