Commit 9b6ea606 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Changed the solution script for the Sedov blast to match the ICs.


Former-commit-id: fccd3921614e1cf0f6aa7ea43e55fa863318e553
parent 08fb977c
......@@ -30,13 +30,13 @@ from numpy import *
# Parameters
rho_0 = 1. # Background Density
P_0 = 1.e-5 # Background Pressure
E_0 = 1.e5 # Energy of the explosion
E_0 = 1.e2 # Energy of the explosion
gamma = 5./3. # Gas polytropic index
t = 0.1 # Time of the solution
t = 0.15 # Time of the solution
N = 1000 # Number of radial points
R_max = 5. # Maximal radius
R_max = 3. # Maximal radius
# ---------------------------------------------------------------
# Don't touch anything after this.
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment