vtf-logo

weno/applications/euler/3d/Homogeneous/src/cles_transport.f

      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_CST
      kappa_mode = CLES_TRANSP_KAPPA_CST
      diff_mode = CLES_TRANSP_DIFF_SCHMIDT

      ierr = 0

      return
      end

      subroutine cles_transport(ux, ncomps, vx, nvars, visc, kappa,
     $     diff, nscal, n, ierr)

      implicit none
      
      include 'cuser.i'

      integer nvars, ncomps, n, nscal, ierr
      double precision ux(ncomps,*), vx(nvars,*)
      double precision visc, kappa, diff(*)
      integer i
      
      visc = 1.0d0/ Re
      kappa = 1.0d0 / alpha

      do i = 1, nscal
         diff(i) = 1.0d0
      enddo

      ierr = 0

      return
      end