Restart with a change in rho algorithm option (idilat)

Questions and remarks about code_saturne usage
Forum rules
Please read the forum usage recommendations before posting.
Post Reply
Cbarnaud
Posts: 36
Joined: Tue Oct 22, 2019 10:28 am

Restart with a change in rho algorithm option (idilat)

Post by Cbarnaud »

Hello

I'm using CS6.0.0 and I achieved simulations with the parameter idilat=2 in cs_user_parameters.f90 .
I would like to use this simulation as a restart, this time with the parameter idilat=3 (low-mach algorithm) with LES and exactly the same input file (xml).
But the job fails at the first iteration of the restart, I have this error :

Call stack:
1: 0x2ad0c762e340 <+0x36340> (libc.so.6)
2: 0x2ad0c15767f7 <pthrbm_+0x130b> (libsaturne-6.0.so)
3: 0x2ad0c15bace2 <tridim_+0x54e6> (libsaturne-6.0.so)
4: 0x2ad0c12ce928 <caltri_+0x40ab> (libsaturne-6.0.so)
5: 0x2ad0b78cbb23 <cs_run+0x559> (libcs_solver-6.0.so)
6: 0x2ad0b78cbe59 <main+0x19f> (libcs_solver-6.0.so)
7: 0x2ad0c761a495 <__libc_start_main+0xf5> (libc.so.6)
8: 0x40a51f <> (cs_solver)
End of stack

And I checked, it seems I have an error with the line 262 of base/pthrbm.f90 :
if (irangp.ge.0 .or. iperio.eq.1) then
call synsca(crom)
endif

But I don't know where it cames from.
Does it exist a way to do a restart of this kind or do I have to start from zero a new simulation?

Thank you a lot,
Camille
Yvan Fournier
Posts: 4070
Joined: Mon Feb 20, 2012 3:25 pm

Re: Restart with a change in rho algorithm option (idilat)

Post by Yvan Fournier »

Hello,

This looks like a bug.

Could you check with version 6.0.6 to see if the bug is still there ? Otherwise do you have a small test case so I can check/fix this ?

Regards,

Yvan
Post Reply