subroutine cles_setup_transport(visc_mode, kappa_mode,
$ diff_mode, ierr)
implicit none
include 'cles.i'
integer visc_mode, kappa_mode, diff_mode, ierr
visc_mode = CLES_TRANSP_VISC_VAR
kappa_mode = CLES_TRANSP_KAPPA_VAR
diff_mode = CLES_TRANSP_DIFF_LEWIS
ierr = 0
return
end
subroutine cles_transport(ux, ncomps, vx, nvars, visc, kappa,
$ diff, nscal, n, ierr)
implicit none
integer nvars, ncomps, n, nscal, ierr
double precision ux(ncomps,*), vx(nvars,*)
double precision visc(*), kappa(*), diff(*)
integer i
double precision tmp
do i = 1, n
tmp = (ux(nvars+1,i)/300.0d0)**(0.7d0)
visc(i) = tmp * 1.8d-5
kappa(i) = tmp * (2.65d-2+5.0d-2*vx(6,i))
enddo
do i = 1, nscal
diff(i) = 1.0d0
enddo
ierr = 0
return
end