physics_euler_mod Module Reference

physics module for 1D,2D and 3D non-isothermal Euler equations More...

Data Types
interface	cons2prim
type	physics_euler
interface	prim2cons
interface	setflux
interface	statevector_euler

Functions/Subroutines
subroutine	initphysics (this, Mesh, config, IO)
	constructor of physics_euler class More...
subroutine	new_statevector (this, new_sv, flavour, num)
	allocate and initialize new non-isothermal state vector More...
pure subroutine	convert2primitive_all (this, cvar, pvar)
	Converts conservative to primitive variables on the whole mesh. More...
pure subroutine	convert2primitive_subset (this, i1, i2, j1, j2, k1, k2, cvar, pvar)
	Converts conservative to primitive variables on a subset of the data. More...
pure subroutine	convert2conservative_all (this, pvar, cvar)
	Converts primitive to conservative variables on the whole mesh. More...
pure subroutine	convert2conservative_subset (this, i1, i2, j1, j2, k1, k2, pvar, cvar)
	Converts primitive to conservative variables on a subset of the data. More...
pure subroutine	calcfluxesx (this, Mesh, nmin, nmax, prim, cons, xfluxes)
	Calculate Fluxes in x-direction. More...
pure subroutine	calcfluxesy (this, Mesh, nmin, nmax, prim, cons, yfluxes)
	Calculate Fluxes in y-direction. More...
pure subroutine	calcfluxesz (this, Mesh, nmin, nmax, prim, cons, zfluxes)
	Calculate Fluxes in z-direction. More...
pure subroutine	calculatecharsystemx (this, Mesh, i1, i2, pvar, lambda, xvar)
pure subroutine	calculatecharsystemy (this, Mesh, j1, j2, pvar, lambda, xvar)
pure subroutine	calculatecharsystemz (this, Mesh, k1, k2, pvar, lambda, xvar)
pure subroutine	calculateboundarydatax (this, Mesh, i1, i2, xvar, pvar)
	extrapolate pvar using characteristic pseudo variables (absorbing boundaries) More...
pure subroutine	calculateboundarydatay (this, Mesh, j1, j2, xvar, pvar)
pure subroutine	calculateboundarydataz (this, Mesh, k1, k2, xvar, pvar)
pure subroutine	calculateprim2riemannx (this, Mesh, i, pvar, lambda, Rinv)
	Conversion from primitive to riemann invariants for farfield boundaries. More...
pure subroutine	calculateprim2riemanny (this, Mesh, j, pvar, lambda, Rinv)
	Conversion from primitive to riemann invariants for farfield boundaries. More...
pure subroutine	calculateprim2riemannz (this, Mesh, k, pvar, lambda, Rinv)
	Conversion from primitive to riemann invariants for farfield boundaries. More...
pure subroutine	calculateriemann2primx (this, Mesh, i, Rinv, pvar)
	Convert Riemann invariants to primitives for farfield boundaries. More...
pure subroutine	calculateriemann2primy (this, Mesh, j, Rinv, pvar)
	Convert Riemann invariants to primitives for farfield boundaries. More...
pure subroutine	calculateriemann2primz (this, Mesh, k, Rinv, pvar)
	Convert Riemann invariants to primitives for farfield boundaries. More...
pure subroutine	geometricalsources (this, Mesh, pvar, cvar, sterm)
	Calculates geometrical sources. More...
subroutine	externalsources (this, accel, pvar, cvar, sterm)
	compute momentum and energy sources given an external force More...
subroutine	viscositysources (this, Mesh, pvar, btxx, btxy, btxz, btyy, btyz, btzz, sterm)
pure subroutine	addbackgroundvelocityx (this, Mesh, w, pvar, cvar)
	Adds a background velocity field for fargo routines. More...
pure subroutine	addbackgroundvelocityy (this, Mesh, w, pvar, cvar)
	Adds a background velocity field for fargo routines. More...
pure subroutine	addbackgroundvelocityz (this, Mesh, w, pvar, cvar)
	Adds a background velocity field for fargo routines. More...
pure subroutine	subtractbackgroundvelocityx (this, Mesh, w, pvar, cvar)
	Substracts a background velocity field for fargo routines. More...
pure subroutine	subtractbackgroundvelocityy (this, Mesh, w, pvar, cvar)
	Substracts a background velocity field for fargo routines. More...
pure subroutine	subtractbackgroundvelocityz (this, Mesh, w, pvar, cvar)
	Substracts a background velocity field for fargo routines. More...
pure subroutine	addfargosourcesx (this, Mesh, w, pvar, cvar, sterm)
	sources terms for fargo advection along x-direction More...
pure subroutine	addfargosourcesy (this, Mesh, w, pvar, cvar, sterm)
	sources terms for fargo advection along y-direction More...
pure subroutine	addfargosourcesz (this, Mesh, w, pvar, cvar, sterm)
	sources terms for fargo advection along z-direction More...
pure subroutine	updatesoundspeed (this, pvar)
subroutine	finalize (this)
	Destructor of the physics_euler class. More...
type(statevector_euler) function	createstatevector (Physics, flavour, num)
	Constructor of statevector_euler. More...
subroutine	assignmarray_0 (this, ma)
	assigns one state vector to another state vector More...
subroutine	finalize_statevector (this)
	actual destructor of the statevector_eulerisotherm type More...
elemental real function	getsoundspeed (gamma, density, pressure)
elemental subroutine	setroeaverages (gamma, rhoL, rhoR, ul, uR, vL, vR, pL, pR, eL, eR, u, cs)
elemental subroutine	setwavespeeds (cs, v, minwav, maxwav)
	set minimal and maximal wave speeds More...
elemental subroutine	seteigenvalues1d (gamma, rho, v, P, l1, l2, l3)
	set all eigenvalues for 1D transport (used in absorbing boundary conditions) More...
elemental subroutine	seteigenvalues2d (gamma, rho, v, P, l1, l2, l3, l4)
	set all eigenvalues for 2D transport (used in absorbing boundary conditions) More...
elemental subroutine	seteigenvalues3d (gamma, rho, v, P, l1, l2, l3, l4, l5)
	set all eigenvalues for 3D transport (used in absorbing boundary conditions) More...
elemental subroutine	setintermediatestate (rhoL, rhoR, uL, uR, vl, vR, pL, pR, eL, eR, amin, amax, rho, mu, mv, e, a)
elemental subroutine	setcharvars1d (gamma, rho1, rho2, u1, u2, P1, P2, l1, l3, xvar1, xvar2, xvar3)
	compute characteristic variables for 1D transport More...
elemental subroutine	setcharvars2d (gamma, rho1, rho2, u1, u2, v1, v2, P1, P2, l1, l4, xvar1, xvar2, xvar3, xvar4)
	compute characteristic variables for 2D transport More...
elemental subroutine	setcharvars3d (gamma, rho1, rho2, u1, u2, v1, v2, w1, w2, P1, P2, l1, l5, xvar1, xvar2, xvar3, xvar4, xvar5)
	compute characteristic variables for 3D transport More...
elemental subroutine	setboundarydata1d (delta, gamma, rho1, u1, P1, xvar1, xvar2, xvar3, rho2, u2, P2)
	extrapolate primitive variables using characteristic pseudo pevariables 1D transport More...
elemental subroutine	setboundarydata2d (delta, gamma, rho1, u1, v1, P1, xvar1, xvar2, xvar3, xvar4, rho2, u2, v2, P2)
	extrapolate primitive variables using characteristic pseudo pevariables 2D transport More...
elemental subroutine	setboundarydata3d (delta, gamma, rho1, u1, v1, w1, P1, xvar1, xvar2, xvar3, xvar4, xvar5, rho2, u2, v2, w2, P2)
	extrapolate primitive variables using characteristic pseudo pevariables 2D transport More...
elemental subroutine	prim2riemann1d (gamma, rho, vx, p, l1, l2, l3, Rminus, Rs, Rplus)
elemental subroutine	prim2riemann2d (gamma, rho, vx, vy, p, l1, l2, l3, l4, Rminus, Rs, Rvt, Rplus)
elemental subroutine	prim2riemann3d (gamma, rho, vx, vy, vz, p, l1, l2, l3, l4, l5, Rminus, Rs, Rvt, Rwt, Rplus)
elemental subroutine	riemann2prim1d (gamma, Rminus, Rs, Rplus, rho, vx, p)
elemental subroutine	riemann2prim2d (gamma, Rminus, Rs, Rvt, Rplus, rho, vx, vy, p)
elemental subroutine	riemann2prim3d (gamma, Rminus, Rs, Rvt, Rwt, Rplus, rho, vx, vy, vz, p)
elemental subroutine	setflux1d (rho, u, P, mu, E, f1, f2, f3)
	set mass, 1D momentum and energy flux for transport along the 1st dimension More...
elemental subroutine	setflux2d (rho, u, P, mu, mv, E, f1, f2, f3, f4)
	set mass, 2D momentum and energy flux for transport along the 1st dimension More...
elemental subroutine	setflux3d (rho, u, P, mu, mv, mw, E, f1, f2, f3, f4, f5)
	set mass, 3D momentum and energy flux for transport along the 1st dimension More...
elemental subroutine	cons2prim1d (gamma, rho_in, mu, E, rho_out, u, P)
	Convert from 1D conservative to primitive variables. More...
elemental subroutine	cons2prim2d (gamma, rho_in, mu, mv, E, rho_out, u, v, P)
	Convert from 2D conservative to primitive variables. More...
elemental subroutine	cons2prim3d (gamma, rho_in, mu, mv, mw, E, rho_out, u, v, w, P)
	Convert from 3D conservative to primitive variables. More...
elemental subroutine	prim2cons1d (gamma, rho_in, u, P, rho_out, mu, E)
	Convert from 1D primitive to conservative variables. More...
elemental subroutine	prim2cons2d (gamma, rho_in, u, v, P, rho_out, mu, mv, E)
	Convert from 2D primitive to conservative variables. More...
elemental subroutine	prim2cons3d (gamma, rho_in, u, v, w, P, rho_out, mu, mv, mw, E)
	Convert from 3D primitive to conservative variables. More...
elemental real function	getgeometricalsourcex (cxyx, cxzx, cyxy, czxz, vx, vy, vz, P, my, mz)
	geometrical momentum source terms P is the either isothermal pressure rho*cs**2 or the real pressure. More...
elemental real function	getgeometricalsourcey (cxyx, cyxy, cyzy, czyz, vx, vy, vz, P, mx, mz)
	y-momentum geometrical source term More...
elemental real function	getgeometricalsourcez (cxzx, cyzy, czxz, czyz, vx, vy, vz, P, mx, my)
	z-momentum geometrical source term More...

Variables
character(len=32), parameter	problem_name = "Euler"

Detailed Description

physics module for 1D,2D and 3D non-isothermal Euler equations

Author

Tobias Illenseer

Björn Sperling

Jannes Klee

Function/Subroutine Documentation

addbackgroundvelocityx()

pure subroutine physics_euler_mod::addbackgroundvelocityx ( class(physics_euler), intent(inout)  this, class(mesh_base), intent(in)  Mesh, real, dimension(mesh%jgmin:mesh%jgmax,mesh%kgmin:mesh%kgmax), intent(in)  w, class(marray_compound), intent(inout)  pvar, class(marray_compound), intent(inout)  cvar  )

private

Adds a background velocity field for fargo routines.

Calculates

\begin{eqnarray*} E &=& E' + m_x' w + \frac{1}{2}\varrho w^2 \\ v_x &=& v_x' + w \\ m_x &=& m_x' + \varrho w, \end{eqnarray*}

with \( E, v_y, m_y \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift.

Definition at line 1955 of file physics_euler.f90.

addbackgroundvelocityy()

pure subroutine physics_euler_mod::addbackgroundvelocityy ( class(physics_euler), intent(inout)  this, class(mesh_base), intent(in)  Mesh, real, dimension(mesh%igmin:mesh%igmax,mesh%kgmin:mesh%kgmax), intent(in)  w, class(marray_compound), intent(inout)  pvar, class(marray_compound), intent(inout)  cvar  )

private

Adds a background velocity field for fargo routines.

Calculates

\begin{eqnarray*} E &=& E' + m_y' w + \frac{1}{2}\varrho w^2 \\ v_y &=& v_y' + w \\ m_y &=& m_y' + \varrho w, \end{eqnarray*}

with \( E, v_y, m_y \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift.

Definition at line 2003 of file physics_euler.f90.

addbackgroundvelocityz()

pure subroutine physics_euler_mod::addbackgroundvelocityz ( class(physics_euler), intent(inout)  this, class(mesh_base), intent(in)  Mesh, real, dimension(mesh%igmin:mesh%igmax,mesh%jgmin:mesh%jgmax), intent(in)  w, class(marray_compound), intent(inout)  pvar, class(marray_compound), intent(inout)  cvar  )

private

Adds a background velocity field for fargo routines.

Calculates

\begin{eqnarray*} E &=& E' + m_z' w + \frac{1}{2}\varrho w^2 \\ v_z &=& v_z' + w \\ m_z &=& m_z' + \varrho w, \end{eqnarray*}

with \( E, v_z, m_z \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift.

Definition at line 2059 of file physics_euler.f90.

addfargosourcesx()

pure subroutine physics_euler_mod::addfargosourcesx	(	class(physics_euler), intent(inout)	this,
		class(mesh_base), intent(in)	Mesh,
		real, dimension(mesh%jgmin:mesh%jgmax,mesh%kgmin:mesh%kgmax), intent(in)	w,
		class(marray_compound), intent(inout)	pvar,
		class(marray_compound), intent(inout)	cvar,
		class(marray_compound), intent(inout)	sterm
	)

sources terms for fargo advection along x-direction

If the background velocity \(\vec{w}=w\,\hat{e}_\xi\) with \(w\) independent of \(\xi\) and \(t\) is subtracted from the overall velocity of the flow, an additional source term occurs in the \(\xi\)-momentum equation:

\[ S_\mathrm{Fargo} = -\varrho \vec{v} \cdot \nabla \vec{w} \]

where \(\vec{v}\) is the real velocity (including the background velocity \(\vec{w}\) ).

Definition at line 2267 of file physics_euler.f90.

Here is the call graph for this function:

addfargosourcesy()

pure subroutine physics_euler_mod::addfargosourcesy	(	class(physics_euler), intent(inout)	this,
		class(mesh_base), intent(in)	Mesh,
		real, dimension(mesh%igmin:mesh%igmax,mesh%kgmin:mesh%kgmax), intent(in)	w,
		class(marray_compound), intent(inout)	pvar,
		class(marray_compound), intent(inout)	cvar,
		class(marray_compound), intent(inout)	sterm
	)

sources terms for fargo advection along y-direction

If the background velocity \(\vec{w}=w\,\hat{e}_\eta\) with \(w\) independent of \(\eta\) and \(t\) is subtracted from the overall velocity of the flow, an additional source term occurs in the \(\eta\)-momentum equation:

\[ S_\mathrm{Fargo} = -\varrho \vec{v} \cdot \nabla \vec{w} \]

where \(\vec{v}\) is the real velocity (including the background velocity \(\vec{w}\) ).

Definition at line 2368 of file physics_euler.f90.

Here is the call graph for this function:

addfargosourcesz()

pure subroutine physics_euler_mod::addfargosourcesz	(	class(physics_euler), intent(inout)	this,
		class(mesh_base), intent(in)	Mesh,
		real, dimension(mesh%igmin:mesh%igmax,mesh%jgmin:mesh%jgmax), intent(in)	w,
		class(marray_compound), intent(inout)	pvar,
		class(marray_compound), intent(inout)	cvar,
		class(marray_compound), intent(inout)	sterm
	)

sources terms for fargo advection along z-direction

If the background velocity \(\vec{w}=w\,\hat{e}_phi\) with \(w\) independent of \(\phi\) and \(t\) is subtracted from the overall velocity of the flow, an additional source term occurs in the \(\phi\)-momentum equation:

\[ S_\mathrm{Fargo} = -\varrho \vec{v} \cdot \nabla \vec{w} \]

where \(\vec{v}\) is the real velocity (including the background velocity \(\vec{w}\) ).

Definition at line 2468 of file physics_euler.f90.

Here is the call graph for this function:

assignmarray_0()

subroutine physics_euler_mod::assignmarray_0 ( class(statevector_euler), intent(inout)  this, class(marray_base), intent(in)  ma  )

private

assigns one state vector to another state vector

Definition at line 2632 of file physics_euler.f90.

calcfluxesx()

pure subroutine physics_euler_mod::calcfluxesx ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  nmin, integer, intent(in)  nmax, class(marray_compound), intent(inout)  prim, class(marray_compound), intent(inout)  cons, class(marray_compound), intent(inout)  xfluxes  )

private

Calculate Fluxes in x-direction.

Definition at line 492 of file physics_euler.f90.

calcfluxesy()

pure subroutine physics_euler_mod::calcfluxesy ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  nmin, integer, intent(in)  nmax, class(marray_compound), intent(inout)  prim, class(marray_compound), intent(inout)  cons, class(marray_compound), intent(inout)  yfluxes  )

private

Calculate Fluxes in y-direction.

Definition at line 548 of file physics_euler.f90.

calcfluxesz()

pure subroutine physics_euler_mod::calcfluxesz ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  nmin, integer, intent(in)  nmax, class(marray_compound), intent(inout)  prim, class(marray_compound), intent(inout)  cons, class(marray_compound), intent(inout)  zfluxes  )

private

Calculate Fluxes in z-direction.

Definition at line 604 of file physics_euler.f90.

calculateboundarydatax()

pure subroutine physics_euler_mod::calculateboundarydatax ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  i1, integer, intent(in)  i2, real, dimension(mesh%jmin:mesh%jmax,mesh%kmin:mesh%kmax,this%vnum), intent(in)  xvar, class(marray_compound), intent(inout)  pvar  )

private

extrapolate pvar using characteristic pseudo variables (absorbing boundaries)

Definition at line 1072 of file physics_euler.f90.

Here is the call graph for this function:

calculateboundarydatay()

pure subroutine physics_euler_mod::calculateboundarydatay ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  j1, integer, intent(in)  j2, real, dimension(mesh%imin:mesh%imax,mesh%kmin:mesh%kmax,this%vnum), intent(in)  xvar, class(marray_compound), intent(inout)  pvar  )

private

Definition at line 1131 of file physics_euler.f90.

Here is the call graph for this function:

calculateboundarydataz()

pure subroutine physics_euler_mod::calculateboundarydataz ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  k1, integer, intent(in)  k2, real, dimension(mesh%imin:mesh%imax,mesh%jmin:mesh%jmax,this%vnum), intent(in)  xvar, class(marray_compound), intent(inout)  pvar  )

private

Definition at line 1204 of file physics_euler.f90.

Here is the call graph for this function:

calculatecharsystemx()

pure subroutine physics_euler_mod::calculatecharsystemx ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  i1, integer, intent(in)  i2, class(marray_compound), intent(inout)  pvar, real, dimension(mesh%jmin:mesh%jmax,mesh%kmin:mesh%kmax,this%vnum), intent(out)  lambda, real, dimension(mesh%jmin:mesh%jmax,mesh%kmin:mesh%kmax,this%vnum), intent(out)  xvar  )

private

Definition at line 743 of file physics_euler.f90.

Here is the call graph for this function:

calculatecharsystemy()

pure subroutine physics_euler_mod::calculatecharsystemy ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  j1, integer, intent(in)  j2, class(marray_compound), intent(inout)  pvar, real, dimension(mesh%imin:mesh%imax,mesh%kmin:mesh%kmax,this%vnum), intent(out)  lambda, real, dimension(mesh%imin:mesh%imax,mesh%kmin:mesh%kmax,this%vnum), intent(out)  xvar  )

private

Definition at line 848 of file physics_euler.f90.

Here is the call graph for this function:

calculatecharsystemz()

pure subroutine physics_euler_mod::calculatecharsystemz ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  k1, integer, intent(in)  k2, class(marray_compound), intent(inout)  pvar, real, dimension(mesh%imin:mesh%imax,mesh%jmin:mesh%jmax,this%vnum), intent(out)  lambda, real, dimension(mesh%imin:mesh%imax,mesh%jmin:mesh%jmax,this%vnum), intent(out)  xvar  )

private

Definition at line 966 of file physics_euler.f90.

Here is the call graph for this function:

calculateprim2riemannx()

pure subroutine physics_euler_mod::calculateprim2riemannx ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  i, class(marray_compound), intent(in)  pvar, real, dimension(mesh%jmin:mesh%jmax,mesh%kmin:mesh%kmax,this%vnum), intent(out)  lambda, real, dimension(mesh%jmin:mesh%jmax,mesh%kmin:mesh%kmax,this%vnum), intent(out)  Rinv  )

private

Conversion from primitive to riemann invariants for farfield boundaries.

Definition at line 1266 of file physics_euler.f90.

Here is the call graph for this function:

calculateprim2riemanny()

pure subroutine physics_euler_mod::calculateprim2riemanny ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  j, class(marray_compound), intent(in)  pvar, real, dimension(mesh%imin:mesh%imax,mesh%kmin:mesh%kmax,this%vnum), intent(out)  lambda, real, dimension(mesh%imin:mesh%imax,mesh%kmin:mesh%kmax,this%vnum), intent(out)  Rinv  )

private

Conversion from primitive to riemann invariants for farfield boundaries.

Definition at line 1360 of file physics_euler.f90.

Here is the call graph for this function:

calculateprim2riemannz()

pure subroutine physics_euler_mod::calculateprim2riemannz ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  k, class(marray_compound), intent(in)  pvar, real, dimension(mesh%imin:mesh%imax,mesh%jmin:mesh%jmax,this%vnum), intent(out)  lambda, real, dimension(mesh%imin:mesh%imax,mesh%jmin:mesh%jmax,this%vnum), intent(out)  Rinv  )

private

Conversion from primitive to riemann invariants for farfield boundaries.

Definition at line 1456 of file physics_euler.f90.

Here is the call graph for this function:

calculateriemann2primx()

pure subroutine physics_euler_mod::calculateriemann2primx ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  i, real, dimension(mesh%jmin:mesh%jmax,mesh%kmin:mesh%kmax,this%vnum), intent(in)  Rinv, class(marray_compound), intent(inout)  pvar  )

private

Convert Riemann invariants to primitives for farfield boundaries.

Definition at line 1555 of file physics_euler.f90.

Here is the call graph for this function:

calculateriemann2primy()

pure subroutine physics_euler_mod::calculateriemann2primy ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  j, real, dimension(mesh%imin:mesh%imax,mesh%kmin:mesh%kmax,this%vnum), intent(in)  Rinv, class(marray_compound), intent(inout)  pvar  )

private

Convert Riemann invariants to primitives for farfield boundaries.

Definition at line 1606 of file physics_euler.f90.

Here is the call graph for this function:

calculateriemann2primz()

pure subroutine physics_euler_mod::calculateriemann2primz ( class(physics_euler), intent(in)  this, class(mesh_base), intent(in)  Mesh, integer, intent(in)  k, real, dimension(mesh%imin:mesh%imax,mesh%jmin:mesh%jmax,this%vnum), intent(in)  Rinv, class(marray_compound), intent(inout)  pvar  )

private

Convert Riemann invariants to primitives for farfield boundaries.

Definition at line 1657 of file physics_euler.f90.

Here is the call graph for this function:

cons2prim1d()

elemental subroutine physics_euler_mod::cons2prim1d ( real, intent(in)  gamma, real, intent(in)  rho_in, real, intent(in)  mu, real, intent(in)  E, real, intent(out)  rho_out, real, intent(out)  u, real, intent(out)  P  )

private

Convert from 1D conservative to primitive variables.

Definition at line 3033 of file physics_euler.f90.

cons2prim2d()

elemental subroutine physics_euler_mod::cons2prim2d ( real, intent(in)  gamma, real, intent(in)  rho_in, real, intent(in)  mu, real, intent(in)  mv, real, intent(in)  E, real, intent(out)  rho_out, real, intent(out)  u, real, intent(out)  v, real, intent(out)  P  )

private

Convert from 2D conservative to primitive variables.

Definition at line 3048 of file physics_euler.f90.

cons2prim3d()

elemental subroutine physics_euler_mod::cons2prim3d ( real, intent(in)  gamma, real, intent(in)  rho_in, real, intent(in)  mu, real, intent(in)  mv, real, intent(in)  mw, real, intent(in)  E, real, intent(out)  rho_out, real, intent(out)  u, real, intent(out)  v, real, intent(out)  w, real, intent(out)  P  )

private

Convert from 3D conservative to primitive variables.

Definition at line 3064 of file physics_euler.f90.

convert2conservative_all()

pure subroutine physics_euler_mod::convert2conservative_all ( class(physics_euler), intent(in)  this, class(marray_compound), intent(inout)  pvar, class(marray_compound), intent(inout)  cvar  )

private

Converts primitive to conservative variables on the whole mesh.

Definition at line 367 of file physics_euler.f90.

convert2conservative_subset()

pure subroutine physics_euler_mod::convert2conservative_subset ( class(physics_euler), intent(in)  this, integer, intent(in)  i1, integer, intent(in)  i2, integer, intent(in)  j1, integer, intent(in)  j2, integer, intent(in)  k1, integer, intent(in)  k2, class(marray_compound), intent(inout)  pvar, class(marray_compound), intent(inout)  cvar  )

private

Converts primitive to conservative variables on a subset of the data.

Definition at line 405 of file physics_euler.f90.

convert2primitive_all()

pure subroutine physics_euler_mod::convert2primitive_all ( class(physics_euler), intent(in)  this, class(marray_compound), intent(inout)  cvar, class(marray_compound), intent(inout)  pvar  )

private

Converts conservative to primitive variables on the whole mesh.

Definition at line 243 of file physics_euler.f90.

convert2primitive_subset()

pure subroutine physics_euler_mod::convert2primitive_subset ( class(physics_euler), intent(in)  this, integer, intent(in)  i1, integer, intent(in)  i2, integer, intent(in)  j1, integer, intent(in)  j2, integer, intent(in)  k1, integer, intent(in)  k2, class(marray_compound), intent(inout)  cvar, class(marray_compound), intent(inout)  pvar  )

private

Converts conservative to primitive variables on a subset of the data.

Definition at line 281 of file physics_euler.f90.

createstatevector()

type(statevector_euler) function physics_euler_mod::createstatevector	(	class(physics_euler), intent(in)	Physics,
		integer, intent(in), optional	flavour,
		integer, intent(in), optional	num
	)

Constructor of statevector_euler.

methods of class statevector_euler

This is not a class member itself, instead its an ordinary module procedure. The function name is overloaded with the class name.

Definition at line 2593 of file physics_euler.f90.

Here is the caller graph for this function:

externalsources()

subroutine physics_euler_mod::externalsources	(	class(physics_euler), intent(in)	this,
		class(marray_base), intent(in)	accel,
		class(marray_compound), intent(in)	pvar,
		class(marray_compound), intent(in)	cvar,
		class(marray_compound), intent(inout)	sterm
	)

compute momentum and energy sources given an external force

Definition at line 1849 of file physics_euler.f90.

finalize()

subroutine physics_euler_mod::finalize

(

class(physics_euler), intent(inout)

this

)

Destructor of the physics_euler class.

Definition at line 2577 of file physics_euler.f90.

finalize_statevector()

subroutine physics_euler_mod::finalize_statevector ( type(statevector_euler), intent(inout)  this )

private

actual destructor of the statevector_eulerisotherm type

Definition at line 2666 of file physics_euler.f90.

Here is the caller graph for this function:

geometricalsources()

pure subroutine physics_euler_mod::geometricalsources ( class(physics_euler), intent(inout)  this, class(mesh_base), intent(in)  Mesh, class(marray_compound), intent(inout)  pvar, class(marray_compound), intent(inout)  cvar, class(marray_compound), intent(inout)  sterm  )

private

Calculates geometrical sources.

Definition at line 1706 of file physics_euler.f90.

Here is the call graph for this function:

getgeometricalsourcex()

elemental real function physics_euler_mod::getgeometricalsourcex ( real, intent(in)  cxyx, real, intent(in)  cxzx, real, intent(in)  cyxy, real, intent(in)  czxz, real, intent(in)  vx, real, intent(in)  vy, real, intent(in)  vz, real, intent(in)  P, real, intent(in)  my, real, intent(in)  mz  )

private

geometrical momentum source terms P is the either isothermal pressure rho*cs**2 or the real pressure.

Attention

These elemental functions exist multiple times for performance reasons (inlining). Please keep this in mind for changes. Other modules with this function:

• physics_eulerisotherm_mod

x-momentum geometrical source term

Definition at line 3128 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

getgeometricalsourcey()

elemental real function physics_euler_mod::getgeometricalsourcey ( real, intent(in)  cxyx, real, intent(in)  cyxy, real, intent(in)  cyzy, real, intent(in)  czyz, real, intent(in)  vx, real, intent(in)  vy, real, intent(in)  vz, real, intent(in)  P, real, intent(in)  mx, real, intent(in)  mz  )

private

y-momentum geometrical source term

Definition at line 3139 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

getgeometricalsourcez()

elemental real function physics_euler_mod::getgeometricalsourcez ( real, intent(in)  cxzx, real, intent(in)  cyzy, real, intent(in)  czxz, real, intent(in)  czyz, real, intent(in)  vx, real, intent(in)  vy, real, intent(in)  vz, real, intent(in)  P, real, intent(in)  mx, real, intent(in)  my  )

private

z-momentum geometrical source term

Definition at line 3150 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

getsoundspeed()

elemental real function physics_euler_mod::getsoundspeed ( real, intent(in)  gamma, real, intent(in)  density, real, intent(in)  pressure  )

private

elemental non-class subroutines / functions

Definition at line 2682 of file physics_euler.f90.

Here is the caller graph for this function:

initphysics()

subroutine physics_euler_mod::initphysics ( class(physics_euler), intent(inout)  this, class(mesh_base), intent(in)  Mesh, type(dict_typ), pointer  config, type(dict_typ), pointer  IO  )

private

constructor of physics_euler class

Todo:

remove / improve in future version set array indices for 1st,2nd,3rd non-vanishing velocities this may actually not coincide with the x,y and z-velocities

Definition at line 135 of file physics_euler.f90.

new_statevector()

subroutine physics_euler_mod::new_statevector	(	class(physics_euler), intent(in)	this,
		class(marray_compound), pointer	new_sv,
		integer, intent(in), optional	flavour,
		integer, intent(in), optional	num
	)

allocate and initialize new non-isothermal state vector

Definition at line 225 of file physics_euler.f90.

prim2cons1d()

elemental subroutine physics_euler_mod::prim2cons1d ( real, intent(in)  gamma, real, intent(in)  rho_in, real, intent(in)  u, real, intent(in)  P, real, intent(out)  rho_out, real, intent(out)  mu, real, intent(out)  E  )

private

Convert from 1D primitive to conservative variables.

Definition at line 3081 of file physics_euler.f90.

prim2cons2d()

elemental subroutine physics_euler_mod::prim2cons2d ( real, intent(in)  gamma, real, intent(in)  rho_in, real, intent(in)  u, real, intent(in)  v, real, intent(in)  P, real, intent(out)  rho_out, real, intent(out)  mu, real, intent(out)  mv, real, intent(out)  E  )

private

Convert from 2D primitive to conservative variables.

Definition at line 3093 of file physics_euler.f90.

prim2cons3d()

elemental subroutine physics_euler_mod::prim2cons3d ( real, intent(in)  gamma, real, intent(in)  rho_in, real, intent(in)  u, real, intent(in)  v, real, intent(in)  w, real, intent(in)  P, real, intent(out)  rho_out, real, intent(out)  mu, real, intent(out)  mv, real, intent(out)  mw, real, intent(out)  E  )

private

Convert from 3D primitive to conservative variables.

Definition at line 3106 of file physics_euler.f90.

prim2riemann1d()

elemental subroutine physics_euler_mod::prim2riemann1d ( real, intent(in)  gamma, real, intent(in)  rho, real, intent(in)  vx, real, intent(in)  p, real, intent(in)  l1, real, intent(in)  l2, real, intent(in)  l3, real, intent(out)  Rminus, real, intent(out)  Rs, real, intent(out)  Rplus  )

private

Definition at line 2893 of file physics_euler.f90.

Here is the caller graph for this function:

prim2riemann2d()

elemental subroutine physics_euler_mod::prim2riemann2d ( real, intent(in)  gamma, real, intent(in)  rho, real, intent(in)  vx, real, intent(in)  vy, real, intent(in)  p, real, intent(in)  l1, real, intent(in)  l2, real, intent(in)  l3, real, intent(in)  l4, real, intent(out)  Rminus, real, intent(out)  Rs, real, intent(out)  Rvt, real, intent(out)  Rplus  )

private

Definition at line 2913 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

prim2riemann3d()

elemental subroutine physics_euler_mod::prim2riemann3d ( real, intent(in)  gamma, real, intent(in)  rho, real, intent(in)  vx, real, intent(in)  vy, real, intent(in)  vz, real, intent(in)  p, real, intent(in)  l1, real, intent(in)  l2, real, intent(in)  l3, real, intent(in)  l4, real, intent(in)  l5, real, intent(out)  Rminus, real, intent(out)  Rs, real, intent(out)  Rvt, real, intent(out)  Rwt, real, intent(out)  Rplus  )

private

Definition at line 2930 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

riemann2prim1d()

elemental subroutine physics_euler_mod::riemann2prim1d ( real, intent(in)  gamma, real, intent(in)  Rminus, real, intent(in)  Rs, real, intent(in)  Rplus, real, intent(out)  rho, real, intent(out)  vx, real, intent(out)  p  )

private

Definition at line 2945 of file physics_euler.f90.

Here is the caller graph for this function:

riemann2prim2d()

elemental subroutine physics_euler_mod::riemann2prim2d ( real, intent(in)  gamma, real, intent(in)  Rminus, real, intent(in)  Rs, real, intent(in)  Rvt, real, intent(in)  Rplus, real, intent(out)  rho, real, intent(out)  vx, real, intent(out)  vy, real, intent(out)  p  )

private

Definition at line 2965 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

riemann2prim3d()

elemental subroutine physics_euler_mod::riemann2prim3d ( real, intent(in)  gamma, real, intent(in)  Rminus, real, intent(in)  Rs, real, intent(in)  Rvt, real, intent(in)  Rwt, real, intent(in)  Rplus, real, intent(out)  rho, real, intent(out)  vx, real, intent(out)  vy, real, intent(out)  vz, real, intent(out)  p  )

private

Definition at line 2983 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

setboundarydata1d()

elemental subroutine physics_euler_mod::setboundarydata1d ( integer, intent(in)  delta, real, intent(in)  gamma, real, intent(in)  rho1, real, intent(in)  u1, real, intent(in)  P1, real, intent(in)  xvar1, real, intent(in)  xvar2, real, intent(in)  xvar3, real, intent(out)  rho2, real, intent(out)  u2, real, intent(out)  P2  )

private

extrapolate primitive variables using characteristic pseudo pevariables 1D transport

Definition at line 2843 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

setboundarydata2d()

elemental subroutine physics_euler_mod::setboundarydata2d ( integer, intent(in)  delta, real, intent(in)  gamma, real, intent(in)  rho1, real, intent(in)  u1, real, intent(in)  v1, real, intent(in)  P1, real, intent(in)  xvar1, real, intent(in)  xvar2, real, intent(in)  xvar3, real, intent(in)  xvar4, real, intent(out)  rho2, real, intent(out)  u2, real, intent(out)  v2, real, intent(out)  P2  )

private

extrapolate primitive variables using characteristic pseudo pevariables 2D transport

Definition at line 2865 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

setboundarydata3d()

elemental subroutine physics_euler_mod::setboundarydata3d ( integer, intent(in)  delta, real, intent(in)  gamma, real, intent(in)  rho1, real, intent(in)  u1, real, intent(in)  v1, real, intent(in)  w1, real, intent(in)  P1, real, intent(in)  xvar1, real, intent(in)  xvar2, real, intent(in)  xvar3, real, intent(in)  xvar4, real, intent(in)  xvar5, real, intent(out)  rho2, real, intent(out)  u2, real, intent(out)  v2, real, intent(out)  w2, real, intent(out)  P2  )

private

extrapolate primitive variables using characteristic pseudo pevariables 2D transport

Definition at line 2879 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

setcharvars1d()

elemental subroutine physics_euler_mod::setcharvars1d ( real, intent(in)  gamma, real, intent(in)  rho1, real, intent(in)  rho2, real, intent(in)  u1, real, intent(in)  u2, real, intent(in)  P1, real, intent(in)  P2, real, intent(in)  l1, real, intent(in)  l3, real, intent(out)  xvar1, real, intent(out)  xvar2, real, intent(out)  xvar3  )

private

compute characteristic variables for 1D transport

Definition at line 2799 of file physics_euler.f90.

Here is the caller graph for this function:

setcharvars2d()

elemental subroutine physics_euler_mod::setcharvars2d ( real, intent(in)  gamma, real, intent(in)  rho1, real, intent(in)  rho2, real, intent(in)  u1, real, intent(in)  u2, real, intent(in)  v1, real, intent(in)  v2, real, intent(in)  P1, real, intent(in)  P2, real, intent(in)  l1, real, intent(in)  l4, real, intent(out)  xvar1, real, intent(out)  xvar2, real, intent(out)  xvar3, real, intent(out)  xvar4  )

private

compute characteristic variables for 2D transport

Definition at line 2818 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

setcharvars3d()

elemental subroutine physics_euler_mod::setcharvars3d ( real, intent(in)  gamma, real, intent(in)  rho1, real, intent(in)  rho2, real, intent(in)  u1, real, intent(in)  u2, real, intent(in)  v1, real, intent(in)  v2, real, intent(in)  w1, real, intent(in)  w2, real, intent(in)  P1, real, intent(in)  P2, real, intent(in)  l1, real, intent(in)  l5, real, intent(out)  xvar1, real, intent(out)  xvar2, real, intent(out)  xvar3, real, intent(out)  xvar4, real, intent(out)  xvar5  )

private

compute characteristic variables for 3D transport

Definition at line 2830 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

seteigenvalues1d()

elemental subroutine physics_euler_mod::seteigenvalues1d ( real, intent(in)  gamma, real, intent(in)  rho, real, intent(in)  v, real, intent(in)  P, real, intent(out)  l1, real, intent(out)  l2, real, intent(out)  l3  )

private

set all eigenvalues for 1D transport (used in absorbing boundary conditions)

Definition at line 2729 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

seteigenvalues2d()

elemental subroutine physics_euler_mod::seteigenvalues2d ( real, intent(in)  gamma, real, intent(in)  rho, real, intent(in)  v, real, intent(in)  P, real, intent(out)  l1, real, intent(out)  l2, real, intent(out)  l3, real, intent(out)  l4  )

private

set all eigenvalues for 2D transport (used in absorbing boundary conditions)

Definition at line 2746 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

seteigenvalues3d()

elemental subroutine physics_euler_mod::seteigenvalues3d ( real, intent(in)  gamma, real, intent(in)  rho, real, intent(in)  v, real, intent(in)  P, real, intent(out)  l1, real, intent(out)  l2, real, intent(out)  l3, real, intent(out)  l4, real, intent(out)  l5  )

private

set all eigenvalues for 3D transport (used in absorbing boundary conditions)

Definition at line 2757 of file physics_euler.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

setflux1d()

elemental subroutine physics_euler_mod::setflux1d ( real, intent(in)  rho, real, intent(in)  u, real, intent(in)  P, real, intent(in)  mu, real, intent(in)  E, real, intent(out)  f1, real, intent(out)  f2, real, intent(out)  f3  )

private

set mass, 1D momentum and energy flux for transport along the 1st dimension

Definition at line 2999 of file physics_euler.f90.

Here is the caller graph for this function:

setflux2d()

elemental subroutine physics_euler_mod::setflux2d ( real, intent(in)  rho, real, intent(in)  u, real, intent(in)  P, real, intent(in)  mu, real, intent(in)  mv, real, intent(in)  E, real, intent(out)  f1, real, intent(out)  f2, real, intent(out)  f3, real, intent(out)  f4  )

private

set mass, 2D momentum and energy flux for transport along the 1st dimension

Definition at line 3011 of file physics_euler.f90.

Here is the caller graph for this function:

setflux3d()

elemental subroutine physics_euler_mod::setflux3d ( real, intent(in)  rho, real, intent(in)  u, real, intent(in)  P, real, intent(in)  mu, real, intent(in)  mv, real, intent(in)  mw, real, intent(in)  E, real, intent(out)  f1, real, intent(out)  f2, real, intent(out)  f3, real, intent(out)  f4, real, intent(out)  f5  )

private

set mass, 3D momentum and energy flux for transport along the 1st dimension

Definition at line 3022 of file physics_euler.f90.

setintermediatestate()

elemental subroutine physics_euler_mod::setintermediatestate ( real, intent(in)  rhoL, real, intent(in)  rhoR, real, intent(in)  uL, real, intent(in)  uR, real, intent(in)  vl, real, intent(in)  vR, real, intent(in)  pL, real, intent(in)  pR, real, intent(in)  eL, real, intent(in)  eR, real, intent(in)  amin, real, intent(in)  amax, real, intent(out)  rho, real, intent(out)  mu, real, intent(out)  mv, real, intent(out)  e, real, intent(out)  a  )

private

Definition at line 2769 of file physics_euler.f90.

setroeaverages()

elemental subroutine physics_euler_mod::setroeaverages ( real, intent(in)  gamma, real, intent(in)  rhoL, real, intent(in)  rhoR, real, intent(in)  ul, real, intent(in)  uR, real, intent(in)  vL, real, intent(in)  vR, real, intent(in)  pL, real, intent(in)  pR, real, intent(in)  eL, real, intent(in)  eR, real, intent(out)  u, real, intent(out)  cs  )

private

Todo:

NOT VERIFIED only for advanced wavespeeds

Definition at line 2693 of file physics_euler.f90.

setwavespeeds()

elemental subroutine physics_euler_mod::setwavespeeds ( real, intent(in)  cs, real, intent(in)  v, real, intent(out)  minwav, real, intent(out)  maxwav  )

private

set minimal and maximal wave speeds

Definition at line 2717 of file physics_euler.f90.

subtractbackgroundvelocityx()

pure subroutine physics_euler_mod::subtractbackgroundvelocityx ( class(physics_euler), intent(inout)  this, class(mesh_base), intent(in)  Mesh, real, dimension(mesh%jgmin:mesh%jgmax,mesh%kgmin:mesh%kgmax), intent(in)  w, class(marray_compound), intent(inout)  pvar, class(marray_compound), intent(inout)  cvar  )

private

Substracts a background velocity field for fargo routines.

Calculates

\begin{eqnarray*} E' &=& E - m_x w + \frac{1}{2}\varrho w^2 \\ v_x' &=& v_x - w \\ m_x' &=& m_x - \varrho w, \end{eqnarray*}

with \( E, v_x, m_x \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift. ATTENTION: the "+" before the 3rd term in the energy transformation is correct!

Definition at line 2109 of file physics_euler.f90.

subtractbackgroundvelocityy()

pure subroutine physics_euler_mod::subtractbackgroundvelocityy ( class(physics_euler), intent(inout)  this, class(mesh_base), intent(in)  Mesh, real, dimension(mesh%igmin:mesh%igmax,mesh%kgmin:mesh%kgmax), intent(in)  w, class(marray_compound), intent(inout)  pvar, class(marray_compound), intent(inout)  cvar  )

private

Substracts a background velocity field for fargo routines.

Calculates

\begin{eqnarray*} E' &=& E - m_y w + \frac{1}{2}\varrho w^2 \\ v_y' &=& v_y - w \\ m_y' &=& m_y - \varrho w, \end{eqnarray*}

with \( E, v_y, m_y \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift. ATTENTION: the "+" before the 3rd term in the energy transformation is correct!

Definition at line 2159 of file physics_euler.f90.

subtractbackgroundvelocityz()

pure subroutine physics_euler_mod::subtractbackgroundvelocityz ( class(physics_euler), intent(inout)  this, class(mesh_base), intent(in)  Mesh, real, dimension(mesh%igmin:mesh%igmax,mesh%jgmin:mesh%jgmax), intent(in)  w, class(marray_compound), intent(inout)  pvar, class(marray_compound), intent(inout)  cvar  )

private

Substracts a background velocity field for fargo routines.

Calculates

\begin{eqnarray*} E' &=& E - m_z w + \frac{1}{2}\varrho w^2 \\ v_z' &=& v_z - w \\ m_z' &=& m_z - \varrho w, \end{eqnarray*}

with \( E, v_z, m_z \) the total energy, velocity and momentum. The \( ' \) denotes the residual part. \( w \) is the velocity shift. ATTENTION: the "+" before the 3rd term in the energy transformation is correct!

Definition at line 2217 of file physics_euler.f90.

updatesoundspeed()

pure subroutine physics_euler_mod::updatesoundspeed ( class(physics_euler), intent(inout)  this, class(statevector_euler), intent(inout)  pvar  )

private

Definition at line 2555 of file physics_euler.f90.

Here is the call graph for this function:

viscositysources()

subroutine physics_euler_mod::viscositysources ( class(physics_euler), intent(inout)  this, class(mesh_base), intent(in)  Mesh, class(marray_compound), intent(inout)  pvar, real, dimension(mesh%igmin:mesh%igmax,mesh%jgmin:mesh%jgmax,mesh%kgmin:mesh%kgmax), intent(in)  btxx, real, dimension(mesh%igmin:mesh%igmax,mesh%jgmin:mesh%jgmax,mesh%kgmin:mesh%kgmax), intent(in)  btxy, real, dimension(mesh%igmin:mesh%igmax,mesh%jgmin:mesh%jgmax,mesh%kgmin:mesh%kgmax), intent(in)  btxz, real, dimension(mesh%igmin:mesh%igmax,mesh%jgmin:mesh%jgmax,mesh%kgmin:mesh%kgmax), intent(in)  btyy, real, dimension(mesh%igmin:mesh%igmax,mesh%jgmin:mesh%jgmax,mesh%kgmin:mesh%kgmax), intent(in)  btyz, real, dimension(mesh%igmin:mesh%igmax,mesh%jgmin:mesh%jgmax,mesh%kgmin:mesh%kgmax), intent(in)  btzz, class(marray_compound), intent(inout)  sterm  )

private

Definition at line 1887 of file physics_euler.f90.

Variable Documentation

problem_name

character(len=32), parameter physics_euler_mod::problem_name = "Euler"

private

Definition at line 52 of file physics_euler.f90.