Initialization for cases with variable density and idilat>1

[mathrack]

Dear users, developers,

I had some issues to initialize a case with idilat=2, here is my workaround:

• set idilat=2 inside the subroutine usipsu (cs_user_parameters.f90)
• in cs_user_initialization.f90, add use ihmpre and the following

  
  ! Apply equation of state from graphical interface
  if (iihmpr.eq.1) then
  call uiphyv(iviscv, itempk, visls0, viscv0)
  endif

! Apply equation of state from user subroutine call usphyv(nvar, nscal, 0, dt)

! Put it into previous density call field_current_to_previous(icrom)

Best,
Cédric