1 #ifndef HYDRO_HYDRO_HPP_
2 #define HYDRO_HYDRO_HPP_
16 #include "../athena.hpp"
17 #include "../athena_arrays.hpp"
18 #include "../bvals/cc/hydro/bvals_hydro.hpp"
19 #include "hydro_diffusion/hydro_diffusion.hpp"
83 #if !MAGNETIC_FIELDS_ENABLED
85 const int k,
const int j,
const int il,
const int iu,
91 const int k,
const int j,
const int il,
const int iu,
Real(*)(MeshBlock *pmb) TimeStepFunc
AthenaArray< Real > lambdas_m_l_
AthenaArray< Real > gflx[3]
AthenaArray< Real > cell_volume_
AthenaArray< Real > coarse_prim_
AthenaArray< Real > scr1_nkji_
AthenaArray< Real > x1face_area_
AthenaArray< Real > cons_
void AddGravityFluxWithGflx()
AthenaArray< Real > lambdas_p_l_
AthenaArray< Real > x3face_area_
AthenaArray< Real > bb_normal_
AthenaArray< Real > scr2_nkji_
void CalculateGravityFlux(AthenaArray< Real > &phi_in)
HydroBoundaryVariable hbvar
AthenaArray< Real > lambdas_p_r_
AthenaArray< Real > x3face_area_p1_
void CalculateFluxes(AthenaArray< Real > &w, FaceField &b, AthenaArray< Real > &bcc, const int order)
void CalculateFluxes_STS()
Hydro(MeshBlock *pmb, ParameterInput *pin)
void AddFluxDivergence(const Real wght, AthenaArray< Real > &u_out)
void RiemannSolver(const int k, const int j, const int il, const int iu, const int ivx, AthenaArray< Real > &wl, AthenaArray< Real > &wr, AthenaArray< Real > &flx, const AthenaArray< Real > &dxw)
The HLLC Riemann solver for adiabatic hydrodynamics (use HLLE for isothermal)
AthenaArray< Real > wl3d_
void AddDiffusionFluxes()
AthenaArray< Real > laplacian_l_fc_
AthenaArray< Real > wr3d_
Real GetWeightForCT(Real dflx, Real rhol, Real rhor, Real dx, Real dt)
AthenaArray< Real > laplacian_r_fc_
AthenaArray< Real > x2face_area_
AthenaArray< Real > flux[3]
AthenaArray< Real > coarse_cons_
AthenaArray< Real > x2face_area_p1_
TimeStepFunc UserTimeStep_
AthenaArray< Real > dflx_
AthenaArray< Real > gflx_old[3]
void VaporConcentrationLimiter(AthenaArray< Real > &u)
AthenaArray< Real > lambdas_m_r_