prob/cshock1d.c File Reference

Problem generator for 1-D standing C-type shock test. More...

Go to the source code of this file.

Functions
Real	RK4 (Real D, Real A, Real M, Real theta, Real h)
	4th order Runge-Kutta integrator
Real	Dprime (Real D, Real A, Real M, Real theta)
	dD/dt
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	get_eta_user (GridS *pG, int i, int j, int k, Real *eta_O, Real *eta_H, Real *eta_A)
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 ConsS *	RootSoln = NULL

---

Detailed Description

Problem generator for 1-D standing C-type shock test.

PURPOSE: Problem generator for 1-D standing C-type shock test. Only works for grid-aligned shock, the shock is along the x1 direction.

REFERENCE: Mac Low, M-M et al., "Incorporation of Ambipolar Diffusion into the Zeus Magnetohydrodynamics Code", 1995, ApJ, 442, 726

Definition in file cshock1d.c.

---

Function Documentation

Real Dprime	(	Real	D,
		Real	A,
		Real	M,
		Real	theta
	)

dD/dt

Definition at line 275 of file cshock1d.c.

Referenced by RK4().

Here is the caller graph for this function:

void get_eta_user	(	GridS *	pG,
		int	i,
		int	j,
		int	k,
		Real *	eta_O,
		Real *	eta_H,
		Real *	eta_A
	)

Definition at line 230 of file cshock1d.c.

ConsFun_t get_usr_expr

(

const char *

expr

)

Sets pointer to expression for special output data.

Definition at line 220 of file cshock1d.c.

VOutFun_t get_usr_out_fun

(

const char *

name

)

Returns a user defined output function pointer.

Definition at line 225 of file cshock1d.c.

void problem

(

DomainS *

pDomain

)

Problem initial conditions, edited for each problem.

Definition at line 46 of file cshock1d.c.

References ath_error(), B0, ConsS::B1c, GridS::B1i, ConsS::B2c, GridS::B2i, ConsS::B3c, GridS::B3i, calloc_1d_array(), cc_pos(), ConsS::d, d0, d_ind, DomainS::Grid, GridS::ie, GridS::is, Iso_csound, GridS::je, GridS::js, GridS::ke, GridS::ks, DomainS::Level, Lx, ConsS::M1, ConsS::M2, ConsS::M3, Mach, GridS::Nx, nx1, par_getd(), par_getd_def(), Q_AD, RK4(), DomainS::RootMaxX, DomainS::RootMinX, Soln, theta, GridS::U, x01, x02, x1max, and x1min.

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 215 of file cshock1d.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 210 of file cshock1d.c.

Real RK4	(	Real	D,
		Real	A,
		Real	M,
		Real	theta,
		Real	h
	)

4th order Runge-Kutta integrator

Definition at line 260 of file cshock1d.c.

References Dprime().

Referenced by problem().

Here is the call graph for this function:

Here is the caller graph for this function:

void Userwork_after_loop

(

MeshS *

pM

)

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

Definition at line 251 of file cshock1d.c.

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 241 of file cshock1d.c.

---

Variable Documentation

ConsS* RootSoln = NULL [static]

Definition at line 41 of file cshock1d.c.