Dear all,
In SISYPHE there is a algorithm for restraining the max slope in the mesh (maxslope.f). It adjusts the elevation of each node when the average slope exceeds a threshold.
In Telemac-3D, I had a problem near the river banks, where excessive sedimentation is always found and the bottom elevation increases over time. So I tried to migrate the algorithm in SISYPHE to Telemac-3D.
Attached is a working example. I used PRIVE2D array defined in Telemac-3D to store the adjusted elevations so I don't have to declare a new one. The attached code is meant for a simple case, one bed-layer and one sediment class, although the original algorithm in SISYPHE can deal with multiple layers and sediment classes.
At this moment the code seems working fine for me. I saw the adjustments of the mesh over time steps. But I didn't find any documentation about this algorithm so I'm not sure if it's ok to use it in 3D model?
Kind regards,
Qilong