init.c File Reference

#include "pluto.h"

Include dependency graph for init.c:

Go to the source code of this file.

Functions
void	Init (double *us, double x, double y, double z)
void	Analysis (const Data *d, Grid *grid)
void	UserDefBoundary (const Data *d, RBox *box, int side, Grid *grid)

Function Documentation

void Analysis	(	const Data *	d,
		Grid *	grid
	)

Perform runtime data analysis.

Parameters

[in]	d	the PLUTO Data structure
[in]	grid	pointer to array of Grid structures

Compute volume-integrated magnetic pressure, Maxwell and Reynolds stresses. Save them to "averages.dat"

Definition at line 33 of file init.c.

{

}

void Init	(	double *	us,
		double	x,
		double	y,
		double	z
	)

The Init() function can be used to assign initial conditions as as a function of spatial position.

Parameters

[out]	v	a pointer to a vector of primitive variables
[in]	x1	coordinate point in the 1st dimension
[in]	x2	coordinate point in the 2nd dimension
[in]	x3	coordinate point in the 3rdt dimension

The meaning of x1, x2 and x3 depends on the geometry:

Variable names are accessed by means of an index v[nv], where nv = RHO is density, nv = PRS is pressure, nv = (VX1, VX2, VX3) are the three components of velocity, and so forth.

Definition at line 4 of file init.c.

{
  double r, x0=1,y0=1,scrh;

  g_gamma = 1.4;

  r = sqrt( (x-x0)*(x-x0) + (y-y0)*(y-y0));
  us[RHO] = 1.0;
  us[VX1] = 0.0;
  us[VX2] = 0.0;
  us[VX3] = 0.0;
  us[TRC] = 0.0;

  if (x<0.1){
    us[PRS] = 1.e3;
  }else if(x<0.9) {
    us[PRS] = 1.e-2;
  }else {
    us[PRS] = 1.e2;
  }

}

void UserDefBoundary	(	const Data *	d,
		RBox *	box,
		int	side,
		Grid *	grid
	)

Assign user-defined boundary conditions.

Parameters

[in,out]	d	pointer to the PLUTO data structure containing cell-centered primitive quantities (d->Vc) and staggered magnetic fields (d->Vs, when used) to be filled.
[in]	box	pointer to a RBox structure containing the lower and upper indices of the ghost zone-centers/nodes or edges at which data values should be assigned.
[in]	side	specifies the boundary side where ghost zones need to be filled. It can assume the following pre-definite values: X1_BEG, X1_END, X2_BEG, X2_END, X3_BEG, X3_END. The special value side == 0 is used to control a region inside the computational domain.
[in]	grid	pointer to an array of Grid structures.

Definition at line 43 of file init.c.

{ }