# Face has no cross edge

**URL:** <https://community.geodynamics.org/t/face-has-no-cross-edge/1300>\
**Category:** PyLith\
**Created:** [June 22, 2020, 6:01pm UTC](https://community.geodynamics.org/t/face-has-no-cross-edge/1300 "2020-06-22T18:01:13Z")\
**Posts on this page:** 4\
**Page:** 1

<div class="post-metadata">

**Author:** ![albarodriguez](https://avatars.discourse-cdn.com/v4/letter/a/f05b48/32.png) [@albarodriguez](https://community.geodynamics.org/u/albarodriguez)\
**Post date:** [June 22, 2020, 6:01pm UTC](https://community.geodynamics.org/t/face-has-no-cross-edge/1300/1 "2020-06-22T18:01:13Z")

</div>

Hi Pylith team,

My mesh is yielding the following error: " RuntimeError: Error occurred while adjusting topology to create cohesive cells for fault ‘LCRF’.Internal error while creating fault mesh. Face 48877827 has no cross edge."

I am using a version of a mesh with three faults that I have used extensively before, but the faults are now built using nonplanar geometries. Let me know if more information is required to address this issue.

Thanks,  
Alba

---

<div class="post-metadata">

**Author:** ![baagaard](https://yyz2.discourse-cdn.com/flex036/user_avatar/community.geodynamics.org/baagaard/32/615_2.png) [@baagaard](https://community.geodynamics.org/u/baagaard)\
**Post date:** [June 22, 2020, 8:19pm UTC](https://community.geodynamics.org/t/face-has-no-cross-edge/1300/2 "2020-06-22T20:19:04Z")

</div>

Please provide an image of your fault geometry (for example, a screenshot of the surfaces in CUBIT/Trelis) so that we can see how nonplanar it is and how the location of the faults relative to each other.

---

<div class="post-metadata">

**Author:** ![albarodriguez](https://avatars.discourse-cdn.com/v4/letter/a/f05b48/32.png) [@albarodriguez](https://community.geodynamics.org/u/albarodriguez)\
**Post date:** [June 22, 2020, 9:40pm UTC](https://community.geodynamics.org/t/face-has-no-cross-edge/1300/3 "2020-06-22T21:40:08Z")

</div>

The fault in the middle is the problematic one. The two large ones are meshed based on USGS geometry. The middle one is meshed based on field data.

 ![Screen Shot 2020-06-22 at 2.39.51 PM](https://canada1.discourse-cdn.com/flex036/uploads/geodynamics/original/1X/5382ea9433ddbb6d3d1a69b45dc91f0a60706c66.png)  
 ![Screen Shot 2020-06-22 at 2.38.12 PM](https://canada1.discourse-cdn.com/flex036/uploads/geodynamics/original/1X/c92c4c006d6730e43f223a48259bac1924b1a321.png) ![Screen Shot 2020-06-22 at 2.28.27 PM](https://canada1.discourse-cdn.com/flex036/uploads/geodynamics/original/1X/7b91c698863597110402d75d123ea2228145690b.png) ![Screen Shot 2020-06-22 at 2.37.47 PM](https://canada1.discourse-cdn.com/flex036/uploads/geodynamics/original/1X/674a2b6e3281b13a9cadd3fc34aef411ef1164bc.png) ![Screen Shot 2020-06-22 at 2.36.50 PM](https://canada1.discourse-cdn.com/flex036/uploads/geodynamics/original/1X/f5fd8b3a66444d091b40c34ebf791e05a231f8b1.png)

---

<div class="post-metadata">

**Author:** ![baagaard](https://yyz2.discourse-cdn.com/flex036/user_avatar/community.geodynamics.org/baagaard/32/615_2.png) [@baagaard](https://community.geodynamics.org/u/baagaard)\
**Post date:** [June 23, 2020, 2:39pm UTC](https://community.geodynamics.org/t/face-has-no-cross-edge/1300/4 "2020-06-23T14:39:40Z")

</div>

The order you specify the faults in the `interfaces` array is important as it influences how the cohesive cells are inserted into the mesh and what happens where the faults intersect.

You should specify the two major (non-intersecting) faults first and then the connecting fault.

```auto
interfaces = [fault_major_1, fault_major_2, connecting_fault]

```

You also need to make sure the buried edges are marked correctly and included in each fault specification.

If you haven’t already done so, I would run the following simulations to check to make sure the individual and two major faults are fine:

1. A static simulation each fault separately (i.e., three different static simulations) run properly.
2. A static simulation using the two major faults runs properly.