# Point is not inside the domain when exporting surface flux

**URL:** <https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110>\
**Category:** User Support\
**Created:** [March 24, 2025, 7:37am UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110 "2025-03-24T07:37:23Z")\
**Posts on this page:** 18\
**Page:** 1

<div class="post-metadata">

**Author:** ![Ucempvs](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/ucempvs/32/150_2.png) [@Ucempvs](https://festim.discourse.group/u/Ucempvs)\
**Post date:** [March 24, 2025, 7:37am UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/1 "2025-03-24T07:37:23Z")

</div>

Hello all,

I want to perform a simulation on a 3D disc-shaped structure that I meshed with SALOME and then export the surface flux as a derived quantity. However, I noticed that when I use a finer mesh, I get the error message

```auto
 File "/home/ucempvs/FESTIM/Disc/CAD_example.py", line 116, in <module>
    preview_plot(5.10924409e-04, 0.21747223 )
    ~~~~~~~~~~~~ ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "/home/ucempvs/FESTIM/Disc/CAD_example.py", line 108, in preview_plot
    res = TDS(*prm)
  File "/home/ucempvs/FESTIM/Disc/CAD_example.py", line 101, in TDS
    simulation.run()
    ~~~~~~~~~~~~~~ ^^
  File "/home/ucempvs/anaconda3/envs/festim-env/lib/python3.13/site-packages/festim/generic_simulation.py", line 460, in
 run
    self.run_transient()
    ~~~~~~~~~~~~~~~~~~ ^^
  File "/home/ucempvs/anaconda3/envs/festim-env/lib/python3.13/site-packages/festim/generic_simulation.py", line 483, in
 run_transient
    self.iterate()
    ~~~~~~~~~~~~ ^^
  File "/home/ucempvs/anaconda3/envs/festim-env/lib/python3.13/site-packages/festim/generic_simulation.py", line 517, in
 iterate
    self.run_post_processing()
    ~~~~~~~~~~~~~~~~~~~~~~~~ ^^
  File "/home/ucempvs/anaconda3/envs/festim-env/lib/python3.13/site-packages/festim/generic_simulation.py", line 544, in
 run_post_processing
    self.exports.write(
    ~~~~~~~~~~~~~~~~~~ ^
        self.label_to_function,
        ^^^^^^^^^^^^^^^^^^^^^^^
        self.mesh.dx,
        ^^^^^^^^^^^^^
    )
    ^
  File "/home/ucempvs/anaconda3/envs/festim-env/lib/python3.13/site-packages/festim/exports/exports.py", line 77, in write
    export.compute(self.t)
    ~~~~~~~~~~~~~~ ^^^^^^^^
  File "/home/ucempvs/anaconda3/envs/festim-env/lib/python3.13/site-packages/festim/exports/derived_quantities/derived_quantities.py", line 133, in compute
    value = quantity.compute()
  File "/home/ucempvs/anaconda3/envs/festim-env/lib/python3.13/site-packages/festim/exports/derived_quantities/surface_flux.py", line 77, in compute
    flux = f.assemble(
        self.prop * f.dot(f.grad(self.function), self.n) * self.ds(self.surface)
    )
  File "/home/ucempvs/anaconda3/envs/festim-env/lib/python3.13/site-packages/dolfin/fem/assembling.py", line 213, in assemble
    assembler.assemble(tensor, dolfin_form)
    ~~~~~~~~~~~~~~~~~~ ^^^^^^^^^^^^^^^^^^^^^
  File "/home/ucempvs/anaconda3/envs/festim-env/lib/python3.13/site-packages/dolfin/function/expression.py", line 56, in
 wrapped_eval_cell
    self.user_expression.eval_cell(values, x, cell)
    ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ^^^^^^^^^^^^^^^^^
  File "/home/ucempvs/anaconda3/envs/festim-env/lib/python3.13/site-packages/festim/materials/materials.py", line 363, in eval_cell
    value[0] = D_0 * f.exp(-E_D / k_B / self._T(x))
                                        ~~~~~~~ ^^^
  File "/home/ucempvs/anaconda3/envs/festim-env/lib/python3.13/site-packages/dolfin/function/function.py", line 337, in __call__
    self._cpp_object.eval(values, x)
    ~~~~~~~~~~~~~~~~~~~~~ ^^^^^^^^^^^
RuntimeError:

*** -------------------------------------------------------------------------
*** DOLFIN encountered an error. If you are not able to resolve this issue
*** using the information listed below, you can ask for help at
***
*** fenics-support@googlegroups.com
***
*** Remember to include the error message listed below and, if possible,
*** include a *minimal* running example to reproduce the error.
***
*** -------------------------------------------------------------------------
*** Error: Unable to evaluate function at point.
*** Reason: The point is not inside the domain. Consider calling "Function::set_allow_extrapolation(true)" on this Function to allow extrapolation.
*** Where: This error was encountered inside Function.cpp.
*** Process: 0
***
*** DOLFIN version: 2019.1.0
*** Git changeset: e7ee545889ccbd48d64c71286627e7f8ae60a690
*** -------------------------------------------------------------------------

```

Exporting to XDMF works as expected and when I remove the surface flux from the list of derived quantities, the code also runs.

Here is my code:

```auto
import matplotlib.pyplot as plt
import festim as F
import numpy as np
import warnings

warnings.filterwarnings("ignore", category=DeprecationWarning)

ni_atom_density = 9.140308E+28  
eV = 96485 # eV to J conversion factor
a_IN718 = 3.591e-10 #lattice parameter of IN718 in m

id_vol = 9

id_top = 6
id_bottom = 7
id_rim = 8

def TDS(n1, E_p1):
 
    trap_conc1 = n1 * ni_atom_density #sites/m3

    simulation = F.Simulation()
    #simulation.log_level = 20

    simulation.mesh = F.MeshFromXDMF(volume_file="Disc/Mesh/mesh_domains.xdmf", boundary_file="Disc/Mesh/mesh_boundaries.xdmf")

    IN718 = F.Material(
        id=id_vol, #ID from mesh file
        D_0=4.06e-7, # Diffusion frequency factor (m^2/s)
        E_D=48630.5 / eV, # Diffusion activation energy in eV
    )
    
    simulation.materials = IN718

    # Define trap properties
    trap_1 = F.Trap(
        k_0=IN718.D_0 /((2**0.5/2*a_IN718)**2*ni_atom_density),
        E_k=IN718.E_D,
        p_0=1.77e6, 
	    E_p=E_p1,
        density=trap_conc1,
        materials=IN718,
    )

    simulation.traps = [trap_1]

    simulation.initial_conditions = [      
        F.InitialCondition(field="1", value = trap_conc1,label=id_vol), 
        ]

    simulation.boundary_conditions = [
        F.DirichletBC(surfaces=[id_top, id_bottom, id_rim], value=0, field=0),
    ]

    ramp = 0.5 # K/s

    simulation.T = 200 + ramp * (F.t)

    simulation.dt = F.Stepsize(
        initial_value=0.1,
        stepsize_change_ratio=1.2,
        max_stepsize=lambda t: None if t < 1 else 2,
        dt_min=1e-8,
    )

    simulation.settings = F.Settings(
        absolute_tolerance=1e11,
        relative_tolerance=1e-6,
        final_time=2500,
        maximum_iterations=100,
    )

    export_folder = "TDS_Fit_Exports/Disc_No_Oxide_1_Trap"

    derived_quantities = F.DerivedQuantities(
        [
            F.SurfaceFlux(field="solute",surface=id_top),
            F.AverageVolume(field="T", volume=id_vol), #temperature
            F.TotalVolume(field="0", volume=id_vol), #mobile H   
            F.AverageVolume(field="0", volume=id_vol), #exports volume of mobile H in H/m3 for solute H 
            F.TotalVolume(field="1", volume=id_vol),  
            F.AverageVolume(field="1", volume=id_vol),          
        ],
        show_units=True,
    )

    xdmf_exports = [
        F.XDMFExport(field="T", filename=export_folder + "/T.xdmf", checkpoint=False),
        F.XDMFExport(field="0", filename=export_folder + "/Mobile H.xdmf", checkpoint=False),
        F.XDMFExport(field="1", filename=export_folder + "/Trapped H.xdmf", checkpoint=False),
    ]
     
    simulation.exports = [derived_quantities] + xdmf_exports

    simulation.initialise()
    simulation.run()

    return derived_quantities

def preview_plot(*prm):
    res = TDS(*prm)
    
    T = np.array(res.filter(fields="T").data)
    flux = -np.array(res[0].data)
    plt.plot(T, flux, linewidth=2, color="black")
    plt.show()
     

preview_plot(5.10924409e-04, 0.21747223 )

```

The mesh was created in SALOME and converted to XDMF using the code in the workshop task 8:

```auto
import meshio

def convert_med_to_xdmf(
    med_file,
    cell_file="mesh_domains.xdmf",
    facet_file="mesh_boundaries.xdmf",
    cell_type="tetra",
    facet_type="triangle",
):
    """Converts a MED mesh to XDMF
    Args:
        med_file (str): the name of the MED file
        cell_file (str, optional): the name of the file containing the
            volume markers. Defaults to "mesh_domains.xdmf".
        facet_file (str, optional): the name of the file containing the
            surface markers.. Defaults to "mesh_boundaries.xdmf".
        cell_type (str, optional): The topology of the cells. Defaults to "tetra".
        facet_type (str, optional): The topology of the facets. Defaults to "triangle".
    Returns:
        dict, dict: the correspondance dict, the cell types
    """
    msh = meshio.read(med_file)

    correspondance_dict = msh.cell_tags

    cell_data_types = msh.cell_data_dict["cell_tags"].keys()

    for mesh_block in msh.cells:
        if mesh_block.type == cell_type:

            meshio.write_points_cells(
                cell_file,
                msh.points,
                [mesh_block],
                cell_data={"f": [-1 * msh.cell_data_dict["cell_tags"][cell_type]]},
            )
        elif mesh_block.type == facet_type:
            meshio.write_points_cells(
                facet_file,
                msh.points,
                [mesh_block],
                cell_data={"f": [-1 * msh.cell_data_dict["cell_tags"][facet_type]]},
            )

    return correspondance_dict, cell_data_types

correspondance_dict, cell_data_types = convert_med_to_xdmf("Disc/Mesh/Disc.med", cell_file="mesh_domains.xdmf", facet_file="mesh_boundaries.xdmf")

print(correspondance_dict)

```

I tried using different meshing algorithms in SALOME but to no avail…  
Attached you can find one mesh where the simulation works as expected and one where it doesn’t.

Thank you for your help!

[Mesh\_Files.zip](https://festim.discourse.group/uploads/short-url/nlTvoNk9JTyuYvU0zAu0PVRww8d.zip) (435.1 KB)

---

<div class="post-metadata">

**Author:** ![remidm](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/remidm/32/4_2.png) [@remidm](https://festim.discourse.group/u/remidm)\
**Post date:** [March 24, 2025, 12:39pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/2 "2025-03-24T12:39:47Z")

</div>

Hey @Ucempvs

Are you running this in parallel by any chance?

---

<div class="post-metadata">

**Author:** ![Ucempvs](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/ucempvs/32/150_2.png) [@Ucempvs](https://festim.discourse.group/u/Ucempvs)\
**Post date:** [March 24, 2025, 2:20pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/3 "2025-03-24T14:20:55Z")

</div>

Hi Remi,

no, I’m not running in parallel.

---

<div class="post-metadata">

**Author:** ![remidm](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/remidm/32/4_2.png) [@remidm](https://festim.discourse.group/u/remidm)\
**Post date:** [March 24, 2025, 2:35pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/4 "2025-03-24T14:35:56Z")

</div>

ok thanks, I’ll try to reproduce your error. In the meantime, I see that your file isn’t really a MWE (minimal reproducible example). For example, the XDMF exports aren’t needed, the TDS plot isn’t needed, etc. Also the mesh filenames don’t match the ones in your script, making it harder for me to reproduce your issue

---

<div class="post-metadata">

**Author:** ![remidm](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/remidm/32/4_2.png) [@remidm](https://festim.discourse.group/u/remidm)\
**Post date:** [March 24, 2025, 2:52pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/5 "2025-03-24T14:52:30Z")

</div>

Been able to reproduce the issue with:

```python
import festim as F

id_vol = 9

id_top = 6
id_bottom = 7
id_rim = 8

simulation = F.Simulation()
# simulation.log_level = 20

simulation.mesh = F.MeshFromXDMF(
    volume_file="mesh_domains.xdmf",
    boundary_file="mesh_boundaries.xdmf",
)

simulation.materials = F.Material(id=id_vol, D_0=1, E_D=0)

simulation.T = 300

simulation.dt = 0.1

simulation.settings = F.Settings(
    absolute_tolerance=1e11,
    relative_tolerance=1e-10,
    final_time=1,
)

derived_quantities = F.DerivedQuantities(
    [F.SurfaceFlux(field="solute", surface=id_top)],
    show_units=True,
)

simulation.exports = [derived_quantities]

simulation.initialise()
simulation.run()

```

---

<div class="post-metadata">

**Author:** ![remidm](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/remidm/32/4_2.png) [@remidm](https://festim.discourse.group/u/remidm)\
**Post date:** [March 24, 2025, 3:05pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/6 "2025-03-24T15:05:44Z")

</div>

@Ucempvs looks to me like it’s a mesh-precision error, which would explain why you don’t have this bug on a finer mesh.

Been able to reproduce the error with a pure FEniCS example:

```python
import fenics as f

mesh = f.Mesh()
f.XDMFFile("mesh_domains.xdmf").read(mesh)
V = f.FunctionSpace(mesh, "CG", 1)
u = f.Function(V)

class Diff(f.UserExpression):
    def eval_cell(self, value, x, ufc_cell):
        D_0 = 1
        E_D = 0.1
        k_B = 8.6173303e-5
        value[0] = D_0 * f.exp(-E_D / k_B / u(x))

    def value_shape(self):
        return ()

my_coeff = Diff()

n = f.FacetNormal(mesh)

flux = f.assemble(my_coeff * f.dot(f.grad(u), n) * f.ds)

```

---

<div class="post-metadata">

**Author:** ![remidm](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/remidm/32/4_2.png) [@remidm](https://festim.discourse.group/u/remidm)\
**Post date:** [March 24, 2025, 3:14pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/7 "2025-03-24T15:14:14Z")

</div>

In the above example, setting

```auto
u.set_allow_extrapolation(True)

```

Seems to fix the issue.

Which means the MWE I provided above can be fixed with:

```python
import festim as F

id_vol = 9

id_top = 6
id_bottom = 7
id_rim = 8

simulation = F.Simulation()
# simulation.log_level = 20

simulation.mesh = F.MeshFromXDMF(
    volume_file="mesh_domains.xdmf",
    boundary_file="mesh_boundaries.xdmf",
)

simulation.materials = F.Material(id=id_vol, D_0=1, E_D=0)

simulation.T = 300

simulation.dt = 0.1

simulation.settings = F.Settings(
    absolute_tolerance=1e11,
    relative_tolerance=1e-10,
    final_time=1,
)

derived_quantities = F.DerivedQuantities(
    [F.SurfaceFlux(field="solute", surface=id_bottom)],
    show_units=True,
)

simulation.exports = [derived_quantities]

simulation.initialise()

# after initialise, allow extrapolation of temeperature
simulation.T.T.set_allow_extrapolation(True)

simulation.run()

```

This makes this MWE run fine. However this looks a bit hacky so you may want to investigate these mesh inaccuracies. Have you considered using GMSH instead of SALOME?

---

<div class="post-metadata">

**Author:** ![Ucempvs](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/ucempvs/32/150_2.png) [@Ucempvs](https://festim.discourse.group/u/Ucempvs)\
**Post date:** [March 24, 2025, 3:45pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/8 "2025-03-24T15:45:47Z")

</div>

Amazing, thank you for going through the effort of reproducing the issue and sorry about leaving some unnecessary stuff in the code.  
I haven’t tried GMSH yet but can give it a go if you think it might prevent this issue!

---

<div class="post-metadata">

**Author:** ![remidm](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/remidm/32/4_2.png) [@remidm](https://festim.discourse.group/u/remidm)\
**Post date:** [March 24, 2025, 3:53pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/9 "2025-03-24T15:53:28Z")

</div>

It’s definitely something to do with mesh precision. I must admit this doesn’t occur a lot… if ever… maybe there is a way to “clean” the mesh in SALOME? is it made of several volumes?

We have a [tutorial with GMSH](https://festim.readthedocs.io/en/latest/userguide/mesh.html#gmsh-example) in the documentation written by @celyn

---

<div class="post-metadata">

**Author:** ![remidm](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/remidm/32/4_2.png) [@remidm](https://festim.discourse.group/u/remidm)\
**Post date:** [March 24, 2025, 3:57pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/10 "2025-03-24T15:57:47Z")

</div>

Also, maybe moving to FESTIM2 could help? I believe all the functionality is already there to simulate this problem.

---

<div class="post-metadata">

**Author:** ![Ucempvs](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/ucempvs/32/150_2.png) [@Ucempvs](https://festim.discourse.group/u/Ucempvs)\
**Post date:** [March 25, 2025, 7:02am UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/11 "2025-03-25T07:02:36Z")

</div>

oh is it? I was holding back from moving to FESTIM 2 but I might give it a go now!

---

<div class="post-metadata">

**Author:** ![VVKulagin](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/vvkulagin/32/51_2.png) [@VVKulagin](https://festim.discourse.group/u/VVKulagin)\
**Post date:** [March 25, 2025, 8:00am UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/12 "2025-03-25T08:00:40Z")

</div>

Hi @Ucempvs !

At first glance, your problem is axisymmetric. Did you consider solving your problem in cylindrical coordinates?

---

<div class="post-metadata">

**Author:** ![remidm](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/remidm/32/4_2.png) [@remidm](https://festim.discourse.group/u/remidm)\
**Post date:** [March 25, 2025, 12:29pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/13 "2025-03-25T12:29:10Z")

</div>

Yes it can do trapping although the way we handle it changed slightly. It seems this is all you need.

To @VVKulagin 's comment, if your system is indeed axisymmetric you may want to run it in 2D-cylindrical to save up some computational resources.

FESTIM2 doesn’t handle cylindrical yet though 😅

---

<div class="post-metadata">

**Author:** ![Ucempvs](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/ucempvs/32/150_2.png) [@Ucempvs](https://festim.discourse.group/u/Ucempvs)\
**Post date:** [March 25, 2025, 1:11pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/14 "2025-03-25T13:11:52Z")

</div>

@VVKulagin, thank you for your suggestion, it is indeed axissymmetric. The only reason I went for the 3D simulation is that I need to consider H desorption from three sides (top, sides, bottom) of the disc as the thickness is quite large.  
My understanding was that I can’t do that with cylindrical coordinates.

To be honest, for finding the initial trapping parameters, I am quite happy with computation time when I use a relatively coarse mesh (~3 minutes per iteration on my laptop).

---

<div class="post-metadata">

**Author:** ![remidm](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/remidm/32/4_2.png) [@remidm](https://festim.discourse.group/u/remidm)\
**Post date:** [March 25, 2025, 1:24pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/15 "2025-03-25T13:24:50Z")

</div>

> [@Ucempvs](#):
>
> thank you for your suggestion, it is indeed axissymmetric. The only reason I went for the 3D simulation is that I need to consider H desorption from three sides (top, sides, bottom) of the disc as the thickness is quite large.  
> My understanding was that I can’t do that with cylindrical coordinates.

This is totally doable in 2D-cylindrical. On the following diagram, you simply set surfaces 1, 2, and 3 with c=0 and then set three different derived quantities.

 ![image](https://global.discourse-cdn.com/free1/uploads/festim/original/1X/b96d7bc116d9b5608b1cd2ac7bbbf0915937a47f.png)

---

<div class="post-metadata">

**Author:** ![Ucempvs](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/ucempvs/32/150_2.png) [@Ucempvs](https://festim.discourse.group/u/Ucempvs)\
**Post date:** [March 25, 2025, 7:06pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/16 "2025-03-25T19:06:45Z")

</div>

Gave it a go and looks quite close to the results from the 3D simulation  
**2D:**

 ![image](https://global.discourse-cdn.com/free1/uploads/festim/original/1X/95004cd15f8459ae9d932d4015ac2477f2064391.png)  
 ![image](https://global.discourse-cdn.com/free1/uploads/festim/original/1X/1ffef4796a18633f4970a4e7b20132eeba2b33e9.jpeg)

**3D:**

 ![image](https://global.discourse-cdn.com/free1/uploads/festim/original/1X/7b42a83f9d73a933621ae4b5a8643f9b384115df.png)  
 ![image](https://global.discourse-cdn.com/free1/uploads/festim/original/1X/8da3938f0b3b240115475b4de02e8e4ccc49ce62.png)

There’s a discrepancy of ~1e18 in the H flux between 2D and 3D but I’m guessing that this could be a consequence of the coarser mesh for 3D simulations? Refining it seems to narrow the gap between the two.

If that’s the case then I’ll definitely stick with 2D

---

<div class="post-metadata">

**Author:** ![remidm](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/remidm/32/4_2.png) [@remidm](https://festim.discourse.group/u/remidm)\
**Post date:** [March 25, 2025, 7:31pm UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/17 "2025-03-25T19:31:49Z")

</div>

Careful with the units here. Are you sure it’s H/m2/s and not H/s?

I would try and have either some local refinement near the surfaces or a finer mesh overall to avoid these oscillations (ref. 2D field).

---

<div class="post-metadata">

**Author:** ![Ucempvs](https://yyz2.discourse-cdn.com/free1/user_avatar/festim.discourse.group/ucempvs/32/150_2.png) [@Ucempvs](https://festim.discourse.group/u/Ucempvs)\
**Post date:** [March 26, 2025, 6:35am UTC](https://festim.discourse.group/t/point-is-not-inside-the-domain-when-exporting-surface-flux/110/18 "2025-03-26T06:35:13Z")

</div>

good shout!  
The units should be right, I divided the flux by the respective surface areas to be able to compare it to experimental data