mus_turb_calcVisc Subroutine

public subroutine mus_turb_calcVisc(turbData, turbConfig, calcTurbVisc, state, neigh, auxField, gradData, densPos, velPos, nSize, nSolve, nScalars, nAuxScalars, layout, dxL, dtL, viscKine)

This routine compute turbulence viscosity and stores in turbulence data type


type(mus_turbulence_data_type), intent(inout) :: turbData

turbulence data type

type(mus_turbulence_config_type), intent(in) :: turbConfig

turbulence configuration

type(mus_turbulence_visc_proc_type), intent(in) :: calcTurbVisc

turbulence function

real(kind=rk), intent(in) :: state(:)

state array

integer, intent(in) :: neigh(:)

neigh array to obtain precollision pdf

real(kind=rk), intent(in) :: auxField(:)

Auxiliary field variable array

type(mus_gradData_type), intent(in) :: gradData

gradient data

integer, intent(in) :: densPos

position of density in auxField

integer, intent(in) :: velPos(3)

position of velocity components in auxField

integer, intent(in) :: nSize

number of elements in state array

integer, intent(in) :: nSolve

Number of element to solve in this level

integer, intent(in) :: nScalars

number of scalars in state array

integer, intent(in) :: nAuxScalars

number of scalars in auxField array

type(mus_scheme_layout_type), intent(in) :: layout

scheme layout

real(kind=rk), intent(in) :: dxL

lattice element size in current level

real(kind=rk), intent(in) :: dtL

current level lattice time step size

real(kind=rk), intent(in) :: viscKine(:)

Background kinematic viscosity divided by dtL

Called by

proc~~mus_turb_calcvisc~~CalledByGraph proc~mus_turb_calcvisc mus_turb_calcVisc proc~mus_update_visckine mus_update_viscKine proc~mus_update_visckine->proc~mus_turb_calcvisc proc~mus_update_relaxparams mus_update_relaxParams proc~mus_update_relaxparams->proc~mus_update_visckine proc~do_recursive_multilevel do_recursive_multiLevel proc~do_recursive_multilevel->proc~mus_update_relaxparams proc~do_fast_singlelevel do_fast_singleLevel proc~do_fast_singlelevel->proc~mus_update_relaxparams proc~do_benchmark do_benchmark proc~do_benchmark->proc~mus_update_relaxparams