Pons Michaël

Hi everyone, is it fine if I leave a comment here for a new publication entry ? thank you in advance for your help ! @article{https://doi.org/10.1029/2023TC007765, author = {Pons, Michaël...

Hi Timo, thanks for your reply, the shrinking still happens when the global level of refinement is set to 0 (video : https://youtu.be/rXHt5qS6Vww). The model crashed after 20Myr. ``` subsection...

Oh yes, sorry, here is the initial model with global refinement 1 and time steps between mesh refinement 0 (https://youtu.be/xMQs2K0hdZ8). Shrinkage still occurs, but is much less -60 km mean...

@MFraters advised me to switch to the normal projection, I tried the normal projection instead and it seems to work https://youtu.be/PfKEuJMkuV8. The temperature also matches the initial boundary condition. I...

> @Minerallo I've been reading through some parts of the code and found that there are a number of member variables that are either not used, or only ever set....

yes I would think that n_used_vertices() make more sense to use to set the vectors using the patch. | Quantity | Value | |----------------------------------|-------| | `surface_mesh.n_active_cells()` | 100 | |...

If we set n_grid_nodes = surface_mesh.n_active_cells(); during initialization, it should work for both Box and SphericalShell geometries. That means I might be able to remove the line: `n_grid_nodes = spherical_vertex_index_map.size();`...

Hi @gassmoeller, Bernhard has recalculated a new corrected lookup table (cc) and the purple profile is resulting from the model run and was extracted via paraview (Bernhard could you confirm...

This sounds great! Just to clarify, does this mean it will allow interpolation of values at remote points, for example on a higher-resolution mesh used by an external surface model?...