K-epsilon Realizable: modifying routine turbke.f90

[Bruno Furieri]

I modified some parameters in the turbke.f90 routine to obtain the same system of equations
proposed by the Realizable model. In the comparison of these two models
(standard and realizable) some terms must be deleted and also some terms must
be inserted.

My routine compiles and begins the simulation without any problem. However, after a certain number
of iterations it diverges (always in the same iteration). For information, I launched
all the simulations without restart file.

Somebody has already tried to modify this routine to obtain the k-epsilon Realizable or
even somebody knows other way to implement this modified turbulence model?

Thank you all!

Bruno FURIERI

[Mickael Hassanaly]

Hello Bruno,
You forgot terms array in the arguments of the subroutine
Best regards
Mickaël