Peano
Loading...
Searching...
No Matches
plotSolution.py File Reference

Go to the source code of this file.

Namespaces

namespace  plotSolution
 

Variables

float plotSolution.gas_gamma = 5.0 / 3.0
 This file is part of SWIFT.
 
float plotSolution.rho0 = 1.0
 
float plotSolution.P0 = 1.0e-6
 
int plotSolution.v0 = 1
 
float plotSolution.time = 0.2
 
str plotSolution.pvdfile = "particles.pvd"
 
 plotSolution.reader = ParticleVTUReader(pvdfile=pvdfile, snapshot_time=time)
 
 plotSolution.partData = reader.load()
 
 plotSolution.x = partData.x[:, 0]
 
 plotSolution.y = partData.x[:, 1]
 
 plotSolution.vx = partData.v[:, 0]
 
 plotSolution.vy = partData.v[:, 1]
 
 plotSolution.u = partData.u
 
 plotSolution.P = partData.pressure
 
 plotSolution.rho = partData.density
 
 plotSolution.r = np.sqrt((x - 1.0) ** 2 + (y - 1.0) ** 2)
 
 plotSolution.v = -np.sqrt(vx**2 + vy**2)
 
 plotSolution.r_bin_edge = np.arange(0.0, 1.0, 0.02)
 
float plotSolution.r_bin = 0.5 * (r_bin_edge[1:] + r_bin_edge[:-1])
 
 plotSolution.rho_bin
 
 plotSolution._
 
 plotSolution.statistic
 
 plotSolution.bins
 
 plotSolution.v_bin
 
 plotSolution.P_bin
 
 plotSolution.u_bin
 
 plotSolution.rho2_bin
 
 plotSolution.v2_bin
 
 plotSolution.P2_bin
 
 plotSolution.u2_bin
 
 plotSolution.rho_sigma_bin = np.sqrt(rho2_bin - rho_bin**2)
 
 plotSolution.v_sigma_bin = np.sqrt(v2_bin - v_bin**2)
 
 plotSolution.P_sigma_bin = np.sqrt(P2_bin - P_bin**2)
 
 plotSolution.u_sigma_bin = np.sqrt(u2_bin - u_bin**2)
 
int plotSolution.N = 1000
 
float plotSolution.x_s = np.arange(0, 2.0, 2.0 / N) - 1.0
 
float plotSolution.rho_s = np.ones(N) * rho0
 
float plotSolution.P_s = np.ones(N) * rho0
 
int plotSolution.v_s = np.ones(N) * v0
 
float plotSolution.u0 = rho0 * P0 * (gas_gamma - 1)
 
float plotSolution.us = 0.5 * (gas_gamma - 1) * v0
 
float plotSolution.rs = us * time
 
float plotSolution.u_s = P_s / (rho_s * (gas_gamma - 1.0))
 
float plotSolution.s_s = P_s / rho_s**gas_gamma
 
 plotSolution.figsize
 
str plotSolution.line_color = "C4"
 
str plotSolution.binned_color = "C2"
 
int plotSolution.binned_marker_size = 4
 
 plotSolution.scatter_props
 
 plotSolution.errorbar_props = dict(color=binned_color, ms=binned_marker_size, fmt=".", lw=1.2)
 
int plotSolution.rows = 2
 
int plotSolution.cols = 2
 
 plotSolution.color
 
 plotSolution.alpha
 
 plotSolution.lw
 
 plotSolution.yerr
 
 plotSolution.labelpad
 
str plotSolution.fname = f"NohProfile{time:.4f}.png"
 
 plotSolution.dpi