[makaveli_lcf]

Hi all!

I'm trying to add some source terms to the equations of turbulent model variables (say k and epsilon). Is there any approach to add some sink terms to equations for k and epsilon rather then rewriting turbulent model source code?
A need just simple source terms such as k1*k and k2*epsilon in corresponding equations. The problem is, that equation system matrix is deffined in correct() method and solved at place. How to acces that matrices? Also k1 and k2 are based on external fields, not described in turbulence model.

So, looking forward to hear some advices from you! Thank you, best regards!