NatsukiHosono · mikinakajima · Jul 27, 2019 · Jul 27, 2019
diff --git a/eos/ANEOS_conversion.py b/eos/ANEOS_conversion.py
@@ -134,10 +134,10 @@ def reformat(number):
     fig.savefig("Density" + str(m) + ".png")
 
 h = open(outputfilename,'w')
-h.write("%i %i %s \n" % (nr, nu , ':Grid numbers for density and internal energy'))
-h.write('Density (km/m3), Internal energy (kJ/kg), Temperature (K), Sound speed (m/s), Entropy (J/K/kg) \n')
+h.write("%s %i %i %s \n" % ('#', nr, nu , ':Grid numbers for density and internal energy'))
+h.write('# Density (km/m3), Internal energy (kJ/kg), Temperature (K), Sound speed (m/s), Entropy (J/K/kg) \n')
 
 
 for m in range(0,nr):
     for n in range(0,nu):
-        h.write("%15.8E %15.8E %15.8E %15.8E %15.8E \n" % (density[m], new_energy[n], new_temperature[m][n],  new_soundspeed[m][n], new_entropy[m][n]))
+        h.write("%15.8E %15.8E %15.8E %15.8E %15.8E %15.8E \n" % (density[m], new_energy[n], new_temperature[m][n], new_pressure[m][n],  new_soundspeed[m][n], new_entropy[m][n]))
diff --git a/eos/Density0.png b/eos/Density0.png
diff --git a/eos/Density100.png b/eos/Density100.png
diff --git a/eos/Density20.png b/eos/Density20.png
diff --git a/eos/Density40.png b/eos/Density40.png
diff --git a/eos/Density60.png b/eos/Density60.png
diff --git a/eos/Density80.png b/eos/Density80.png
diff --git a/eos/granite.rho_u.txt b/eos/granite.rho_u.txt
diff --git a/examples/input.txt b/examples/input.txt
@@ -1,4 +1,10 @@
 # input parameters
+mode = 2
+output_interval = 10
+damping = 1
+end_time = 1e5
+output_directory = results/
+impVel = 1000
 UnitMass = 5.97e+24 # target mass
 UnitRadi = 6100e+3 # target radius
 coreFracRadi = 0.54 # core fraction in radius

diff --git a/input/imp.txt b/input/imp.txt
diff --git a/input/tar.txt b/input/tar.txt
diff --git a/src/EoS.h b/src/EoS.h
@@ -93,7 +93,7 @@ namespace EoS{
 		 * A member variable that contains all the EoS data from ANEOS. The lines
 		 * represent density, the columns energy. The data values for each entry
 		 * are
-		 * Density (kg/m3) Temperature (K) Energy (J/kg) Pressure (Pa) Sound speed (m/s) Entropy (J/kg/K).
+		 * Density (kg/m3) Internal energy (J/kg) Temperature (K) Pressure (Pa) Sound speed (m/s) Entropy (J/kg/K).
 		 */
 		std::vector<std::vector<std::array<type, 6> > > eos_data;
 
@@ -232,7 +232,7 @@ namespace EoS{
 		}
 
 		inline type Pressure(const type dens, const type eng) const{
-			return get_interpolated_value(dens,eng,3);
+			return std::max(get_interpolated_value(dens,eng,3), 0.0);
 		}
 
 		inline type SoundSpeed(const type dens, const type eng) const{