open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello,

I want to calculate the critical shields parameter AC as a grain-Reynolds-number dependent value.
In the current release the crit. shields parameter AC is calculated in init_sediment.f dependening on the nondimensional grain diameter.

Now my idea: calculate AC in bedload_formula.f which is called every time step and where I have all the necessary variables.
My problem is, that AC is in doubleprecision format AC(NSICLM).
Can I convert AC to a bief object, so that it has the same vector structure like e.g. TETAP (=nondimensional shields number)?

Or is there generally a better way to do that?

Thanks for any hints!
Clemens