Problem generator for streaming instability test in non-stratified disks. More...

Functions
double	ran2 (long int *idum)
void	OutputModeAmplitude (Grid pGrid, Domain pDomain, Output *pOut)
	Output the perturbation amplitude function.
static Real	ShearingBoxPot (const Real x1, const Real x2, const Real x3)
	static Real ShearingBoxPot(const Real x1, const Real x2, const Real x3)
static Real	pert_even (Real fR, Real fI, Real x, Real z, Real t)
	even perturbation mode
static Real	pert_odd (Real fR, Real fI, Real x, Real z, Real t)
	odd perturbation mode
static int	property_mybin (const Grain gr, const GrainAux grsub)
	user defined particle selection function (1: true; 0: false)
Real	expr_V1par (const Grid *pG, const int i, const int j, const int k)
	Wrapper for particle 1-velocity.
Real	expr_V2par (const Grid *pG, const int i, const int j, const int k)
	Wrapper for particle 2-velocity.
Real	expr_V3par (const Grid *pG, const int i, const int j, const int k)
	Wrapper for particle 3-velocity.
Real	expr_V2 (const Grid *pG, const int i, const int j, const int k)
void	problem (Grid pGrid, Domain pDomain)
void	problem_write_restart (Grid pG, Domain pD, FILE *fp)
void	problem_read_restart (Grid pG, Domain pD, FILE *fp)
static Real	expr_rhodif (const Grid *pG, const int i, const int j, const int k)
	difd
static Real	expr_dVx (const Grid *pG, const int i, const int j, const int k)
	dVx
static Real	expr_dVy (const Grid *pG, const int i, const int j, const int k)
	dVy
static Real	expr_rhopardif (const Grid *pG, const int i, const int j, const int k)
	difdpar
static Real	expr_dVxpar (const Grid *pG, const int i, const int j, const int k)
	dVxpar
static Real	expr_dVypar (const Grid *pG, const int i, const int j, const int k)
	dVypar
Gasfun_t	get_usr_expr (const char *expr)
	Sets pointer to expression for special output data.
VGFunout_t	get_usr_out_fun (const char *name)
	Returns a user defined output function pointer.
PropFun_t	get_usr_par_prop (const char *name)
	Returns a user defined particle selection function.
void	gasvshift (const Real x1, const Real x2, const Real x3, Real u1, Real u2, Real *u3)
	Gas velocity shift.
void	Userforce_particle (Vector *ft, const Real x1, const Real x2, const Real x3, const Real v1, const Real v2, const Real v3)
void	Userwork_in_loop (Grid pGrid, Domain pDomain)
void	Userwork_after_loop (Grid pGrid, Domain pDomain)
Variables
Real	rho0
Real	mratio
Real	etavk
Real	uxNSH
Real	uyNSH
Real	wxNSH
Real	wyNSH
int	Npar
int	Npar2
int	downsamp
int	Nx
Real	Reux
Real	Imux
Real	Reuy
Real	Imuy
Real	Reuz
Real	Imuz
Real	Rewx
Real	Imwx
Real	Rewy
Real	Imwy
Real	Rewz
Real	Imwz
Real	Rerho
Real	Imrho
Real	omg
Real	s
Real	amp
Real	kx
Real	kz
int	ipert
int	nwave
char	name [50]

Detailed Description

Problem generator for streaming instability test in non-stratified disks.

PURPOSE: Problem generator for streaming instability test in non-stratified disks. This code works in 2D, with single particle species ONLY. Isothermal eos is assumed, and the value etavk/iso_sound is fixed. In the case of growth rate test, since one requires an integer number of wavelengths across the grid, the isothermal sound speed is re-evaluated.

Perturbation modes: