code_saturne User's Forum

Hello,

Can you post the whole code instead of just small pieces of it ?

Regards,

Yvan

Hello Yvan,

Please find attached the cs_user_extra_operations.f90 for a simple cases without computation of criteria.

We modified the code and program a ionic conduction module (Electrical Double Layer Effect), the scalar ipproc(iqelec) is calculated in this module and represent the space charge density at the interface solid/liquid.

icell3 is a color of only one cell close to interface (in medium liquid).

So my question is -> why when I am working with only one procs, I obtain 1.d0 in residu.dat file and 0.d0 with multi procs run ?

Thank you for your help
Regards
Paul

Hello,

Why did you remove the call to parsom ? In parallel, you need to call

parsom(residu(ntcabs))

After the loop on selected cells, and before writing the output.

Regards,

Yvan

Hello,

Thank you very much, does it works now.
I placed the line
where you suggested

Thank you again
Ragards

Paul