subroutine usvist                                                 &
 ( nvar   , nscal  ,                                              &
   mbrom  )

!===============================================================================

!===============================================================================
! Module files
!===============================================================================

use paramx
use pointe
use numvar
use optcal
use cstphy
use cstnum
use entsor
use parall
use period
use ppppar
use ppthch
use ppincl
use field
use mesh
use lagran
use cs_c_bindings

!===============================================================================

implicit none

! Arguments

integer          nvar   , nscal

integer          mbrom

integer          iel

double precision, dimension(:), pointer :: dens
double precision, dimension(:), pointer :: visc
double precision, dimension(:), pointer :: visct
double precision, dimension(:,:), pointer :: cvara_vel
double precision, dimension(:), pointer :: mu_t
double precision, dimension(:), pointer :: rey_z

double precision dh_mut

dh_mut = 10.

 call field_get_val_prev_s(iprpfl(irom), dens)
 call field_get_val_prev_s(iprpfl(iviscl), visc)
 call field_get_val_prev_s(iprpfl(ivisct), visct)
 call field_get_val_prev_v(ivarfl(iu), cvara_vel)
 
 do iel = 1, ncel
    rey_z(iel) = cvara_vel(3,iel)*dens(iel)*dh_mut/visc(iel)
    visct(iel) = 1.d0 !visc(iel)*1.d-3*rey_z(iel)
 enddo
 
!  call field_get_val_s_by_name('mut', mu_t)
 
!    do iel = 1, ncel
!      mu_t(iel) = visct(iel)
!    enddo

return
end subroutine usvist