prob/par_shwave1d.c File Reference

Problem generator radial epicyclic wave test with Lagrangian particles. More...

Go to the source code of this file.

Functions
int	GetPosition (Grain *gr)
	get particle status (grid/crossing)
static Real	ShearingBoxPot (const Real x1, const Real x2, const Real x3)
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	ScalarDen (const Grid *pG, const int i, const int j, const int k)
	Scalar density.
static Real	dratio (const Grid *pG, const int i, const int j, const int k)
	Density ratio.
static Real	expr_dV3 (const Grid *pG, const int i, const int j, const int k)
	3-component of velocity
static Real	expr_dV3par (const Grid *pG, const int i, const int j, const int k)
	3-component of particle velocity
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)

---

Detailed Description

Problem generator radial epicyclic wave test with Lagrangian particles.

PURPOSE: Problem generator radial epicyclic wave test with Lagrangian particles.

Only works for 2D, and the wavevector is in x1 direction. Particles can be initialized either with uniform density or with density profile same as the gas. SHEARING_BOX must be turned on, FARGO is optional. ipert: 1 for linear wave; 2 for non-linear wave from Fromang & Papaloizou.

Reference: Fromang & Papaloizou, A&A, 468, 1-18 (2007).

Definition in file par_shwave1d.c.

---

Function Documentation

static Real dratio	(	const Grid *	pG,
		const int	i,
		const int	j,
		const int	k
	)			[static]

Density ratio.

Definition at line 284 of file par_shwave1d.c.

static Real expr_dV3	(	const Grid *	pG,
		const int	i,
		const int	j,
		const int	k
	)			[static]

3-component of velocity

Definition at line 296 of file par_shwave1d.c.

References cc_pos(), Omega_0, and qshear.

Here is the call graph for this function:

static Real expr_dV3par	(	const Grid *	pG,
		const int	i,
		const int	j,
		const int	k
	)			[static]

3-component of particle velocity

Definition at line 311 of file par_shwave1d.c.

References cc_pos(), expr_V3par(), Omega_0, and qshear.

Here is the call graph for this function:

void gasvshift	(	const Real	x1,
		const Real	x2,
		const Real	x3,
		Real *	u1,
		Real *	u2,
		Real *	u3
	)

Gas velocity shift.

Definition at line 346 of file par_shwave1d.c.

Gasfun_t get_usr_expr

(

const char *

expr

)

Sets pointer to expression for special output data.

Definition at line 323 of file par_shwave1d.c.

VGFunout_t get_usr_out_fun

(

const char *

name

)

Returns a user defined output function pointer.

Definition at line 336 of file par_shwave1d.c.

PropFun_t get_usr_par_prop

(

const char *

name

)

Returns a user defined particle selection function.

Definition at line 341 of file par_shwave1d.c.

int GetPosition

(

Grain *

gr

)

get particle status (grid/crossing)

Get particle status (grid/crossing).

Definition at line 457 of file par_shwave1d.c.

References x1lpar, x1upar, x2lpar, x2upar, x3lpar, and x3upar.

void problem	(	Grid *	pGrid,
		Domain *	pDomain
	)

Definition at line 48 of file par_shwave1d.c.

References amp, ath_error(), cc_pos(), GetPosition(), ipert, Iso_csound, kx, n, Npar, nwave, Omega_0, omg, par_getd(), par_getd_def(), par_geti(), par_geti_def(), particle_realloc(), qshear, ShearingBoxPot, StaticGravPot, tstop0, x1max, x1min, x2max, and x2min.

Here is the call graph for this function:

void problem_read_restart	(	Grid *	pG,
		Domain *	pD,
		FILE *	fp
	)

Definition at line 262 of file par_shwave1d.c.

References Omega_0, par_getd_def(), qshear, ShearingBoxPot, and StaticGravPot.

Here is the call graph for this function:

void problem_write_restart	(	Grid *	pG,
		Domain *	pD,
		FILE *	fp
	)

Definition at line 257 of file par_shwave1d.c.

static Real ScalarDen	(	const Grid *	pG,
		const int	i,
		const int	j,
		const int	k
	)			[static]

Scalar density.

Definition at line 275 of file par_shwave1d.c.

static Real ShearingBoxPot	(	const Real	x1,
		const Real	x2,
		const Real	x3
	)			[static]

Definition at line 472 of file par_shwave1d.c.

References Omega_0, and qshear.

void Userforce_particle	(	Vector *	ft,
		const Real	x1,
		const Real	x2,
		const Real	x3,
		const Real	v1,
		const Real	v2,
		const Real	v3
	)

Definition at line 352 of file par_shwave1d.c.

void Userwork_after_loop	(	Grid *	pGrid,
		Domain *	pDomain
	)

Definition at line 370 of file par_shwave1d.c.

References amp, ath_fname(), cc_pos(), Iso_csound, kx, n, nwave, Omega_0, omg, par_getd(), par_getd_def(), par_geti(), particle_to_grid(), property_all(), qshear, x1max, x1min, x2max, and x2min.

Here is the call graph for this function:

void Userwork_in_loop	(	Grid *	pGrid,
		Domain *	pDomain
	)

Definition at line 359 of file par_shwave1d.c.