prob/cyladvect.c File Reference

A simple density-pulse advection test in cylindrical coordinates with no pressure or tension forces. More...

Go to the source code of this file.

Functions
static Real	grav_pot (const Real x1, const Real x2, const Real x3)
	Gravitational potential.
static Real	grav_acc (const Real x1, const Real x2, const Real x3)
	Gravitational acceleration.
Real	d (const Real x1, const Real x2, const Real x3)
	Density.
Real	M2 (const Real x1, const Real x2, const Real x3)
	2-component of momentum
void	cyladvect_ix1 (GridS *pG)
	Inner-R boundary conditions.
void	cyladvect_ox1 (GridS *pG)
	B_R = B_0/R boundary conditions, Outer x1 boundary.
void	problem (DomainS *pDomain)
	Problem initial conditions, edited for each problem.
void	problem_write_restart (MeshS *pM, FILE *fp)
	Problem specific subroutine for use with writing restart files, Edited for each problem.
void	problem_read_restart (MeshS *pM, FILE *fp)
	Problem specific subroutine for use with reading restart files, Edited for each problem.
ConsFun_t	get_usr_expr (const char *expr)
	Sets pointer to expression for special output data.
VOutFun_t	get_usr_out_fun (const char *name)
	Returns a user defined output function pointer.
void	Userwork_in_loop (MeshS *pM)
	Problem work executed during each step in the main Athena loop, see main.c.
void	Userwork_after_loop (MeshS *pM)
	Problem work executed after the main Athena loop, see main.c Edited for each problem.
Variables
static Real	omega0
static Real	rho0
static Real	bz0
static Real	vz0
static Real	Pgas0
static Real	amp
static Real	R0
static Real	phi0
static Real	z0
static Real	rad
static Real	alpha
static Real	x2min
static Real	x2max
static int	iprob
static ConsS ***	RootSoln = NULL

---

Detailed Description

A simple density-pulse advection test in cylindrical coordinates with no pressure or tension forces.

Definition in file cyladvect.c.

---

Function Documentation

void cyladvect_ix1

(

GridS *

pG

)

Inner-R boundary conditions.

d, M2, B1, B1i, and P are all functions of R, phi, and t.

Definition at line 237 of file cyladvect.c.

References alpha, avg2d(), ConsS::B3c, bz0, d(), ConsS::d, ConsS::E, Gamma_1, GridS::is, GridS::je, GridS::js, GridS::ke, GridS::ks, ConsS::M1, M2(), ConsS::M2, ConsS::M3, omega0, Pgas0, phi0, GridS::time, GridS::U, vz0, and x2min.

Referenced by problem().

Here is the call graph for this function:

Here is the caller graph for this function:

void cyladvect_ox1

(

GridS *

pG

)

B_R = B_0/R boundary conditions, Outer x1 boundary.

Definition at line 279 of file cyladvect.c.

References alpha, avg2d(), ConsS::B3c, bz0, d(), ConsS::d, ConsS::E, Gamma_1, GridS::ie, GridS::je, GridS::js, GridS::ke, GridS::ks, ConsS::M1, M2(), ConsS::M2, ConsS::M3, omega0, Pgas0, phi0, GridS::time, GridS::U, vz0, and x2min.

Referenced by problem().

Here is the call graph for this function:

Here is the caller graph for this function:

Real d	(	const Real	x1,
		const Real	x2,
		const Real	x3
	)

Density.

Definition at line 197 of file cyladvect.c.

References alpha, amp, iprob, phi0, r, R0, rad, rho0, x0, and z0.

Referenced by cyladvect_ix1(), cyladvect_ox1(), cylbr_ix1(), cylbr_ox1(), dmrbv_ijb(), dump_vtk(), integrate_1d_ctu(), integrate_1d_vl(), integrate_2d_ctu(), integrate_2d_vl(), integrate_3d_ctu(), integrate_3d_vl(), integrate_emf1_corner(), integrate_emf2_corner(), integrate_emf3_corner(), lr_states(), M2(), problem(), ProCon(), reflect_ix2(), reflect_ix3(), reflect_ox2(), and reflect_ox3().

Here is the caller graph for this function:

ConsFun_t get_usr_expr

(

const char *

expr

)

Sets pointer to expression for special output data.

Definition at line 172 of file cyladvect.c.

VOutFun_t get_usr_out_fun

(

const char *

name

)

Returns a user defined output function pointer.

Definition at line 177 of file cyladvect.c.

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

Gravitational acceleration.

Definition at line 227 of file cyladvect.c.

References omega0.

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

Gravitational potential.

Definition at line 221 of file cyladvect.c.

References omega0.

Real M2	(	const Real	x1,
		const Real	x2,
		const Real	x3
	)

2-component of momentum

Definition at line 215 of file cyladvect.c.

References d(), and omega0.

Referenced by cyladvect_ix1(), cyladvect_ox1(), cylbr_ix1(), cylbr_ox1(), dmrbv_ijb(), problem(), and ProCon().

Here is the call graph for this function:

Here is the caller graph for this function:

void problem

(

DomainS *

pDomain

)

Problem initial conditions, edited for each problem.

Definition at line 37 of file cyladvect.c.

References alpha, amp, ath_error(), avg2d(), ConsS::B3c, GridS::B3i, bvals_mhd_fun(), bz0, calloc_3d_array(), cc_pos(), cyladvect_ix1(), cyladvect_ox1(), d(), ConsS::d, do_nothing_bc(), ConsS::E, Gamma_1, DomainS::Grid, GridS::ie, il, iprob, GridS::is, iu, GridS::je, jl, GridS::js, ju, GridS::ke, kl, GridS::ks, ku, left_x1, ConsS::M1, M2(), ConsS::M2, ConsS::M3, nx1, nx2, nx3, omega0, par_getd(), par_geti(), Pgas0, phi0, R0, rad, rho0, right_x1, StaticGravPot, GridS::U, vz0, x1GravAcc, x1vc(), x2max, x2min, and z0.

Here is the call graph for this function:

void problem_read_restart	(	MeshS *	pM,
		FILE *	fp
	)

Problem specific subroutine for use with reading restart files, Edited for each problem.

Definition at line 167 of file cyladvect.c.

void problem_write_restart	(	MeshS *	pM,
		FILE *	fp
	)

Problem specific subroutine for use with writing restart files, Edited for each problem.

Definition at line 162 of file cyladvect.c.

void Userwork_after_loop

(

MeshS *

pM

)

Problem work executed after the main Athena loop, see main.c Edited for each problem.

Definition at line 188 of file cyladvect.c.

References compute_l1_error().

Here is the call graph for this function:

void Userwork_in_loop

(

MeshS *

pM

)

Problem work executed during each step in the main Athena loop, see main.c.

Edited for each problem.

Definition at line 181 of file cyladvect.c.

References compute_div_b(), MeshS::Domain, and DomainS::Grid.

Here is the call graph for this function:

---

Variable Documentation

Real alpha [static]

Definition at line 24 of file cyladvect.c.

Referenced by cyladvect_ix1(), cyladvect_ox1(), d(), and problem().

Real amp [static]

Definition at line 24 of file cyladvect.c.

Referenced by A3(), d(), ParticlePosition(), ParticleVelocity(), problem(), problem_read_restart(), and Userwork_after_loop().

Real bz0 [static]

Definition at line 24 of file cyladvect.c.

Referenced by cyladvect_ix1(), cyladvect_ox1(), and problem().

int iprob [static]

Definition at line 25 of file cyladvect.c.

Referenced by d(), and problem().

Real omega0 [static]

Definition at line 24 of file cyladvect.c.

Referenced by cyladvect_ix1(), cyladvect_ox1(), grav_acc(), grav_pot(), M2(), and problem().

Real Pgas0 [static]

Definition at line 24 of file cyladvect.c.

Referenced by cyladvect_ix1(), cyladvect_ox1(), and problem().

Real phi0 [static]

Definition at line 24 of file cyladvect.c.

Referenced by cyladvect_ix1(), cyladvect_ox1(), d(), and problem().

Real R0 [static]

Definition at line 24 of file cyladvect.c.

Referenced by d(), and problem().

Real rad [static]

Definition at line 24 of file cyladvect.c.

Referenced by d(), density(), grav_acc(), grav_pot(), jet_iib(), PlanetPot(), problem(), problem_read_restart(), and selfg_multig_3d().

Real rho0 [static]

Definition at line 24 of file cyladvect.c.

Referenced by BpNet(), d(), density(), expr_rhodif(), expr_rhopardif(), grav_acc(), grav_pot(), M2(), Pgas(), problem(), problem_read_restart(), problem_write_restart(), and Userwork_in_loop().

ConsS*** RootSoln = NULL [static]

Definition at line 32 of file cyladvect.c.

Real vz0 [static]

Definition at line 24 of file cyladvect.c.

Referenced by cyladvect_ix1(), cyladvect_ox1(), and problem().

Real x2max [static]

Definition at line 24 of file cyladvect.c.

Real x2min [static]

Definition at line 24 of file cyladvect.c.

Real z0 [static]

Definition at line 24 of file cyladvect.c.

Referenced by d(), problem(), and ShearingBoxPot().