Add source terms before the FARGO advection step. More...

#include "pluto.h"

Include dependency graph for fargo_source.c:

Functions
void	FARGO_Source (Data_Arr UU, double dt, Grid *grid)

Detailed Description

Add source terms before the FARGO advection step.

This function is called prior to the FARGO advection algorithm to add source terms. At present, we use it only for the energy equation in the ShearingBox module:

$\pd{E'}{t} + w\pd{E'}{y} = (B_xB_y - \rho v_xv'_y)\pd{w}{x}$

where $w = -q\Omega x$ . The discretization follows the algorithm of [GS10], see Eq. (51) and (63) of that paper.

Reference

[GS10] "Implementation of the shearing box approximation in Athena", Stone & Gardiner, ApJS (2010) 189, 142.

Author: A. Mignone (migno.nosp@m.ne@p.nosp@m.h.uni.nosp@m.to.i.nosp@m.t)

Date: Aug 26, 2015

Definition in file fargo_source.c.

Function Documentation

void FARGO_Source	(	Data_Arr	UU,
		double	dt,
		Grid *	grid
	)

Parameters

[in,out]	UU	an array of conserved variables
[in]	dt	the current time increment
[in]	grid	pointer to an array of Grid structures

Definition at line 28 of file fargo_source.c.

 {
 #if (HAVE_ENERGY) && (defined SHEARINGBOX)
   int    i,j,k;
   double scrh, rho, mx, my, Bx, By;
 
   scrh = SB_Q*SB_OMEGA;
   DOM_LOOP(k,j,i){
     rho = UU[k][j][i][RHO];
     mx  = UU[k][j][i][MX1];
     my  = UU[k][j][i][MX2];
     Bx  = UU[k][j][i][BX1];
     By  = UU[k][j][i][BX2];
 
     UU[k][j][i][ENG] += - dt*scrh*Bx*(By - 0.5*dt*Bx*scrh)
                         + dt*scrh*mx*my/rho;
   }
 #endif
 }

Here is the call graph for this function:

Here is the caller graph for this function:

Functions

Detailed Description

Function Documentation