code_saturne User's Forum

Hello,

Yes, this is what I meant in my previous post. I'll have to check this, but in any case, if your density is constant, you can use ro0 (or a local variable with the density) in your user subroutine.

Regards,

Yvan

Hello,

that's what I'm doing now, but I still have a problem with the time, as it still returns an error if I'm using dt() (in a steady or an unsteady case), then I have to use dtref which isn't the used time step for all cells.
I have especially to use it as I'm computing in a steady case with SIMPLEC algorithm, that means I had to modify your code to get the time step for each modified cell.

Anyway, thanks a lot for your help!