code_saturne User's Forum

Hello,

You can not use ischtp = 2 with steady flow algorithm ! You should use unsteady flow option with SIMPLEC algorithm.

You can not use ischtp = 2 with all physical models (eg. RANS models - except Rij-e). In case of use ischtp =2 with incompatible options the error messages should be proper.

You could post your setup. I will check.

It will be more convenient to set ntmabs = ntcabs when your residual is small enough (do this in cs_user_extra_operations.f90).

Regards,

The error appears only when you activate ischtp = 2 ?

Ok.

So you can use something like this to compute velocity gradient:


!===========================================
! Decl
integer inc , iprev
double precision, allocatable, dimension(:,:,:) :: gradv
!============================================
! Code
allocate(gradv(3,3,ncelet))

inc ...

Could you attach your listing here?

And 20 lines above this call stack?

There is no error message in listing ?

Hello,

You have to (re)compute it manually. You can use cs_user_extra_operation.f90 to do this.

You must use gradient function to compute velocity gradient (the gradient function name is dependant of ur Code_Saturne version).

Which version are you using ?

Regards,

Hello,

In which version of Code_Saturne?
The solution is probably in user doc.