viscosity

Module

Description

Routines related to physical viscosity

References:

None

Owner:

Daniel Price

Runtime parameters:

• bulkvisc : magnitude of bulk viscosity

• irealvisc : physical viscosity type (0=none,1=const,2=Shakura/Sunyaev)

• shearparam : magnitude of shear viscosity (irealvisc=1) or alpha_SS (irealvisc=2)

Dependencies:

dim, eos, infile_utils, io, part, timestep

Quick access

Variables:

bulkvisc, hoverr, irealvisc, shearparam

Routines:

dt_viscosity(), read_options_viscosity(), set_defaults_viscosity(), shearfunc(), viscinfo(), write_options_viscosity()

Variables

• viscosity/bulkvisc [real,public]

• viscosity/hoverr [real,public]

• viscosity/irealvisc [integer,public]

• viscosity/shearparam [real,public]

Subroutines and functions

subroutine  viscosity/set_defaults_viscosity()

Called from:

set_default_options()

function  viscosity/shearfunc(xi, yi, zi, spsoundi)

Parameters:

• xi [real,in]

• yi [real,in]

• zi [real,in]

• spsoundi [real,in]

Return:

shearfunc [real]

Use :

part (xyzmh_ptmass()), eos (polyk(), qfacdisc())

Called from:

dt_viscosity()

function  viscosity/dt_viscosity(xi, yi, zi, hi, spsoundi)

Parameters:

• xi [real,in]

• yi [real,in]

• zi [real,in]

• hi [real,in]

• spsoundi [real,in]

Return:

dt_viscosity [real]

Use :

timestep (c_force(), bignumber())

Call to:

shearfunc()

subroutine  viscosity/viscinfo(ivisc, iprint)

Parameters:

• ivisc [integer,in]

• iprint [integer,in]

Use :

dim (maxp(), maxdvdx())

subroutine  viscosity/write_options_viscosity(iwritein)

Parameters:

iwritein [integer,in]

Use :

infile_utils (write_inopt())

subroutine  viscosity/read_options_viscosity(db, nerr)

Parameters:

• db (*) [inopts,inout]

• nerr [integer,inout]

Use :

io (error()), infile_utils (inopts(), read_inopt())