open TELEMAC-MASCARET The mathematically superior suite of solvers

Dear telemac users,

I am looking for an effective method to model spatial diffusion on an unstructured grid in 2D, for example to simulate lateral expansion of vegetation.

So far I used the telemac programming guide to incorporate a cellular automata approach modeling spatial expansion by adjusting the bottom friction variable in the subroutine CORSTR(which is called at every time step). I basically use the connectivity table(IKLE) to find neighboring nodes which I give an increasing value every iteration(in my case increase in the friction coefficient). I thought about of correcting the unequal node distances by calculating simple median distances(between nodes) using trigonometry.

This method works, but is not very effective which brings me to my question. Is there a better way to implement a spatial diffusion mechanism into TELEMAC by for example using the MATDIF function which is applied for viscosity calculations??

My second question is, is there another method to change that interval the subroutine is called(e.g. not at every timestep), except manipulating the "AT" variable within the fortran script?



Regards,
Christian