diff --git a/examples/G2-gassphere/analyze_snapshots.py b/examples/G2-gassphere/analyze_snapshots.py
new file mode 100644
index 0000000000000000000000000000000000000000..5d8792e8998d63625740de2b166ad1127a297fd5
--- /dev/null
+++ b/examples/G2-gassphere/analyze_snapshots.py
@@ -0,0 +1,166 @@
+#!/usr/bin/env python3
+"""
+Code that plots different radial profiles for the Evrard collapse
+
+"""
+# load libraries
+import sys # load sys; needed for exit codes
+import numpy as np # load numpy
+import h5py # load h5py; needed to read snapshots
+import matplotlib
+import matplotlib.pyplot as plt ## needs to be active for plotting!
+import csv
+
+matplotlib.rc_file_defaults()
+FloatType = np.float64
+
+#loads the gadget snapshot name
+def load_snapshot(name):
+ gamma = 5./3.
+
+ try:
+ data = h5py.File(name, "r")
+ except:
+ print("could not open file : "+ name+" !")
+ exit(1)
+
+ time = FloatType(data["Header"].attrs["Time"])
+ Pos = np.array(data["PartType0"]["Coordinates"], dtype = FloatType)
+ Density = np.array(data["PartType0"]["Density"], dtype = FloatType)
+ Velocity = np.array(data["PartType0"]["Velocities"], dtype = FloatType)
+ Uthermal = np.array(data["PartType0"]["InternalEnergy"], dtype = FloatType)
+
+ radius = np.sqrt(Pos[:,0]**2 + Pos[:,1]**2 + Pos[:,2]**2)
+ vr = (Velocity[:,0] * Pos[:,0] + Velocity[:,1] * Pos[:,1] + Velocity[:,2] * Pos[:,2]) / radius[:]
+ origin_particles = np.argwhere(radius == 0.0)
+ vr[origin_particles] = 0
+
+ print(np.max(Density))
+ press = (gamma-1)*Uthermal*Density
+
+ A = press / Density**gamma
+
+ return radius,Density, vr, A, time
+
+#bins the particle properties
+def get_data_bins(radius,Density, vr, A):
+ number_bins = 100
+ min_r = 0.01
+ max_r = 1.0
+
+ r_bin = np.zeros(number_bins)
+ vr_bin = np.zeros(number_bins)
+
+ rho_bin = np.zeros(number_bins)
+ count_bin = np.zeros(number_bins)
+ entropy_bin = np.zeros(number_bins)
+
+ for i in range(radius.size):
+ if(radius[i] < min_r or radius[i] > max_r):
+ continue
+ bin_value = int((np.log10(radius[i] / min_r))/(np.log10(max_r/min_r)) * number_bins)
+ count_bin[bin_value] = count_bin[bin_value] + 1
+ vr_bin[bin_value] = vr_bin[bin_value] + vr[i]
+ rho_bin[bin_value] = rho_bin[bin_value] + Density[i]
+ entropy_bin[bin_value] = entropy_bin[bin_value] + A[i]
+
+
+ vr_bin /= count_bin
+ rho_bin /= count_bin
+ entropy_bin /= count_bin
+
+ for i in range(number_bins):
+ r_bin[i] = (i+0.5)* (np.log10(max_r/min_r)/number_bins) + np.log10(min_r)
+ r_bin[i] = 10**r_bin[i]
+ print(count_bin[i])
+
+ return r_bin,rho_bin, vr_bin, entropy_bin
+
+#loads 1d ppm results
+def load_ppm_result():
+ gamma = 5./3.
+ rost = 3./4./np.pi
+ est = 1.054811e-1 / 1.05
+ pst = rost*est
+ vst = np.sqrt(est)
+ rst = 1.257607
+ time = 0
+
+ radius = np.zeros(350)
+ rho = np.zeros(350)
+ vr = np.zeros(350)
+ press = np.zeros(350)
+
+ with open('ppm_profile/ppm1oaf') as csvfile:
+ readCSV = csv.reader(csvfile)
+ line = 0
+ for row in readCSV:
+ line = line+1
+ values = row[0].split()
+ if(line == 1):
+ time = values[1]
+ continue
+ if(line == 352):
+ break
+
+ radius[line -2] = float(values[1]) /rst*1e-11
+ rho[line -2] = float(values[2]) /rost
+ vr[line -2] = float(values[4]) /vst*1e-8
+ press[line -2] = float(values[3])/pst*1e-16
+
+ rho=rho*(3.0/(4*np.pi))
+ press = press*(3.0/(4*np.pi))
+
+ entropy = press / rho**gamma
+
+ return radius, rho, vr, entropy
+
+
+i_file = int(sys.argv[1])
+
+filename = 'output/snapshot_%03d.hdf5' % i_file
+
+radius,Density, vr, A, time = load_snapshot(filename)
+
+radius,Density, vr, A = get_data_bins(radius,Density, vr, A)
+
+fig, ax = plt.subplots(1,3,figsize=(18,6))
+
+
+ax[0].scatter(radius, Density, facecolors='none', edgecolors='b',alpha = 0.5)
+ax[1].scatter(radius, vr, facecolors='none', edgecolors='b',alpha = 0.5)
+ax[2].scatter(radius, A, facecolors='none', edgecolors='b',alpha = 0.5)
+
+radius_ppm, rho_ppm, vr_ppm, entropy_ppm = load_ppm_result()
+
+ax[0].plot(radius_ppm, rho_ppm, color="k")
+ax[1].plot(radius_ppm, vr_ppm, color="k")
+ax[2].plot(radius_ppm, entropy_ppm, color="k")
+
+ax[0].set_xscale('log')
+ax[1].set_xscale('log')
+ax[2].set_xscale('log')
+ax[0].set_yscale('log')
+ax[0].set_xlim(0.01,1)
+ax[1].set_xlim(0.01,1)
+ax[2].set_xlim(0.01,1)
+
+ax[0].set_ylim(0.004,700)
+ax[1].set_ylim(-1.8,0)
+ax[2].set_ylim(0.0,0.2)
+
+ax[0].set_xlabel("R")
+ax[1].set_xlabel("R")
+ax[2].set_xlabel("R")
+
+ax[0].set_ylabel(r"$\rho$")
+ax[1].set_ylabel(r"$V_R$")
+ax[2].set_ylabel(r"$P/\rho^\gamma$")
+plt.tight_layout()
+
+
+
+
+plt.savefig("evrard_"+str(i_file)+".eps")
+plt.show()
+
diff --git a/examples/G2-gassphere/ppm_profile/ppm1oaf b/examples/G2-gassphere/ppm_profile/ppm1oaf
new file mode 100644
index 0000000000000000000000000000000000000000..91de8b65047d69f7a40b5a71f19bf084bf7ae31d
--- /dev/null
+++ b/examples/G2-gassphere/ppm_profile/ppm1oaf
@@ -0,0 +1,351 @@
+ 985 3.01231E+03
+ 1 2.09742E+08 6.08883E+02 2.32980E+17 -4.51975E+05
+ 2 6.29225E+08 4.39246E+02 2.17950E+17 -6.36015E+04
+ 3 1.04871E+09 3.40706E+02 2.00443E+17 -1.59249E+05
+ 4 1.46819E+09 2.73573E+02 1.84939E+17 -3.45453E+05
+ 5 1.88767E+09 2.33603E+02 1.70469E+17 -5.03304E+05
+ 6 2.30716E+09 2.00144E+02 1.55148E+17 -4.72734E+05
+ 7 2.72664E+09 1.74907E+02 1.40944E+17 -5.76558E+05
+ 8 3.14612E+09 1.55563E+02 1.28594E+17 -7.81126E+05
+ 9 3.56561E+09 1.38856E+02 1.17701E+17 -9.39298E+05
+ 10 3.98509E+09 1.25568E+02 1.07223E+17 -9.69142E+05
+ 11 4.40457E+09 1.12794E+02 9.72219E+16 -9.82451E+05
+ 12 4.82406E+09 1.02851E+02 8.83261E+16 -1.18091E+06
+ 13 5.24354E+09 9.31950E+01 8.01524E+16 -1.46879E+06
+ 14 5.66302E+09 8.48749E+01 7.24827E+16 -1.62407E+06
+ 15 6.08251E+09 7.62484E+01 6.56727E+16 -1.63254E+06
+ 16 6.50199E+09 6.82103E+01 5.97061E+16 -1.70207E+06
+ 17 6.92147E+09 6.08448E+01 5.46841E+16 -1.90813E+06
+ 18 7.34095E+09 5.48109E+01 4.98511E+16 -2.15424E+06
+ 19 7.76044E+09 4.92200E+01 4.52647E+16 -2.33907E+06
+ 20 8.17992E+09 4.46141E+01 4.13352E+16 -2.49095E+06
+ 21 8.59940E+09 4.08728E+01 3.81248E+16 -2.69366E+06
+ 22 9.01889E+09 3.74649E+01 3.53740E+16 -2.86898E+06
+ 23 9.43837E+09 3.43386E+01 3.27234E+16 -2.95118E+06
+ 24 9.85785E+09 3.16545E+01 3.00964E+16 -3.07925E+06
+ 25 1.02773E+10 2.89887E+01 2.77112E+16 -3.38759E+06
+ 26 1.06968E+10 2.68896E+01 2.56969E+16 -3.82324E+06
+ 27 1.11163E+10 2.50261E+01 2.39036E+16 -4.14251E+06
+ 28 1.15358E+10 2.32115E+01 2.22001E+16 -4.22695E+06
+ 29 1.19553E+10 2.17133E+01 2.06823E+16 -4.23208E+06
+ 30 1.23748E+10 2.03094E+01 1.93821E+16 -4.39251E+06
+ 31 1.27942E+10 1.90937E+01 1.81672E+16 -4.80535E+06
+ 32 1.32137E+10 1.78770E+01 1.68878E+16 -5.25835E+06
+ 33 1.36332E+10 1.66777E+01 1.55941E+16 -5.58162E+06
+ 34 1.40527E+10 1.56397E+01 1.44580E+16 -5.85496E+06
+ 35 1.44722E+10 1.47773E+01 1.35899E+16 -6.18214E+06
+ 36 1.48917E+10 1.40399E+01 1.28437E+16 -6.57158E+06
+ 37 1.53111E+10 1.32903E+01 1.20618E+16 -6.84091E+06
+ 38 1.57306E+10 1.25319E+01 1.11725E+16 -7.06526E+06
+ 39 1.61501E+10 1.13578E+01 9.63474E+15 -9.43537E+06
+ 40 1.65696E+10 6.60165E+00 4.04442E+15 -2.66858E+07
+ 41 1.69891E+10 3.43397E+00 6.48042E+14 -4.89991E+07
+ 42 1.74086E+10 3.26692E+00 5.49214E+14 -4.94270E+07
+ 43 1.78280E+10 3.13475E+00 5.18155E+14 -4.89533E+07
+ 44 1.82475E+10 3.03538E+00 4.91510E+14 -4.84845E+07
+ 45 1.86670E+10 2.92910E+00 4.64872E+14 -4.80780E+07
+ 46 1.90865E+10 2.83326E+00 4.41255E+14 -4.76644E+07
+ 47 1.95060E+10 2.74176E+00 4.19138E+14 -4.72712E+07
+ 48 1.99254E+10 2.65508E+00 3.98619E+14 -4.68918E+07
+ 49 2.03449E+10 2.57305E+00 3.79564E+14 -4.65253E+07
+ 50 2.07644E+10 2.49516E+00 3.61813E+14 -4.61716E+07
+ 51 2.11839E+10 2.42120E+00 3.45264E+14 -4.58299E+07
+ 52 2.16034E+10 2.35088E+00 3.29811E+14 -4.54995E+07
+ 53 2.20229E+10 2.28396E+00 3.15361E+14 -4.51798E+07
+ 54 2.24423E+10 2.22021E+00 3.01832E+14 -4.48704E+07
+ 55 2.28618E+10 2.15943E+00 2.89147E+14 -4.45707E+07
+ 56 2.32813E+10 2.10144E+00 2.77241E+14 -4.42803E+07
+ 57 2.37008E+10 2.04604E+00 2.66052E+14 -4.39988E+07
+ 58 2.41203E+10 1.99310E+00 2.55524E+14 -4.37257E+07
+ 59 2.45398E+10 1.94245E+00 2.45608E+14 -4.34606E+07
+ 60 2.49592E+10 1.89397E+00 2.36258E+14 -4.32032E+07
+ 61 2.53787E+10 1.84753E+00 2.27433E+14 -4.29532E+07
+ 62 2.57982E+10 1.80300E+00 2.19094E+14 -4.27102E+07
+ 63 2.62177E+10 1.76028E+00 2.11207E+14 -4.24739E+07
+ 64 2.66372E+10 1.71928E+00 2.03741E+14 -4.22441E+07
+ 65 2.70567E+10 1.67989E+00 1.96667E+14 -4.20205E+07
+ 66 2.74761E+10 1.64203E+00 1.89958E+14 -4.18029E+07
+ 67 2.78956E+10 1.60562E+00 1.83590E+14 -4.15909E+07
+ 68 2.83151E+10 1.57057E+00 1.77541E+14 -4.13844E+07
+ 69 2.87346E+10 1.53683E+00 1.71790E+14 -4.11831E+07
+ 70 2.91541E+10 1.50433E+00 1.66318E+14 -4.09870E+07
+ 71 2.95736E+10 1.47299E+00 1.61107E+14 -4.07956E+07
+ 72 2.99930E+10 1.44277E+00 1.56142E+14 -4.06090E+07
+ 73 3.04125E+10 1.41362E+00 1.51407E+14 -4.04269E+07
+ 74 3.08320E+10 1.38546E+00 1.46889E+14 -4.02491E+07
+ 75 3.12515E+10 1.35827E+00 1.42575E+14 -4.00755E+07
+ 76 3.16710E+10 1.33200E+00 1.38452E+14 -3.99060E+07
+ 77 3.20905E+10 1.30660E+00 1.34510E+14 -3.97403E+07
+ 78 3.25099E+10 1.28203E+00 1.30738E+14 -3.95784E+07
+ 79 3.29294E+10 1.25826E+00 1.27127E+14 -3.94202E+07
+ 80 3.33489E+10 1.23524E+00 1.23669E+14 -3.92655E+07
+ 81 3.37684E+10 1.21296E+00 1.20353E+14 -3.91141E+07
+ 82 3.41879E+10 1.19136E+00 1.17174E+14 -3.89660E+07
+ 83 3.46074E+10 1.17044E+00 1.14123E+14 -3.88212E+07
+ 84 3.50268E+10 1.15014E+00 1.11194E+14 -3.86794E+07
+ 85 3.54463E+10 1.13046E+00 1.08380E+14 -3.85405E+07
+ 86 3.58658E+10 1.11137E+00 1.05676E+14 -3.84046E+07
+ 87 3.62853E+10 1.09283E+00 1.03076E+14 -3.82714E+07
+ 88 3.67048E+10 1.07483E+00 1.00574E+14 -3.81409E+07
+ 89 3.71243E+10 1.05735E+00 9.81665E+13 -3.80131E+07
+ 90 3.75437E+10 1.04036E+00 9.58480E+13 -3.78878E+07
+ 91 3.79632E+10 1.02385E+00 9.36142E+13 -3.77650E+07
+ 92 3.83827E+10 1.00780E+00 9.14612E+13 -3.76445E+07
+ 93 3.88022E+10 9.92181E-01 8.93852E+13 -3.75264E+07
+ 94 3.92217E+10 9.76991E-01 8.73825E+13 -3.74106E+07
+ 95 3.96412E+10 9.62207E-01 8.54498E+13 -3.72969E+07
+ 96 4.00606E+10 9.47816E-01 8.35837E+13 -3.71854E+07
+ 97 4.04801E+10 9.33803E-01 8.17814E+13 -3.70759E+07
+ 98 4.08996E+10 9.20153E-01 8.00398E+13 -3.69685E+07
+ 99 4.13191E+10 9.06854E-01 7.83564E+13 -3.68630E+07
+ 100 4.17386E+10 8.93894E-01 7.67286E+13 -3.67594E+07
+ 101 4.21581E+10 8.81260E-01 7.51539E+13 -3.66576E+07
+ 102 4.25775E+10 8.68940E-01 7.36300E+13 -3.65577E+07
+ 103 4.29970E+10 8.56924E-01 7.21548E+13 -3.64595E+07
+ 104 4.34165E+10 8.45202E-01 7.07263E+13 -3.63630E+07
+ 105 4.38360E+10 8.33764E-01 6.93425E+13 -3.62682E+07
+ 106 4.42555E+10 8.22599E-01 6.80015E+13 -3.61750E+07
+ 107 4.46750E+10 8.11699E-01 6.67015E+13 -3.60834E+07
+ 108 4.50944E+10 8.01055E-01 6.54410E+13 -3.59933E+07
+ 109 4.55139E+10 7.90658E-01 6.42183E+13 -3.59048E+07
+ 110 4.59334E+10 7.80501E-01 6.30319E+13 -3.58176E+07
+ 111 4.63529E+10 7.70576E-01 6.18805E+13 -3.57320E+07
+ 112 4.67724E+10 7.60875E-01 6.07625E+13 -3.56477E+07
+ 113 4.71919E+10 7.51392E-01 5.96768E+13 -3.55647E+07
+ 114 4.76113E+10 7.42119E-01 5.86221E+13 -3.54831E+07
+ 115 4.80308E+10 7.33051E-01 5.75972E+13 -3.54028E+07
+ 116 4.84503E+10 7.24180E-01 5.66011E+13 -3.53237E+07
+ 117 4.88698E+10 7.15500E-01 5.56325E+13 -3.52459E+07
+ 118 4.92893E+10 7.07007E-01 5.46906E+13 -3.51692E+07
+ 119 4.97088E+10 6.98693E-01 5.37743E+13 -3.50938E+07
+ 120 5.01282E+10 6.90555E-01 5.28827E+13 -3.50195E+07
+ 121 5.05477E+10 6.82586E-01 5.20149E+13 -3.49463E+07
+ 122 5.09672E+10 6.74783E-01 5.11701E+13 -3.48742E+07
+ 123 5.13867E+10 6.67140E-01 5.03476E+13 -3.48032E+07
+ 124 5.18062E+10 6.59653E-01 4.95465E+13 -3.47333E+07
+ 125 5.22257E+10 6.52317E-01 4.87660E+13 -3.46644E+07
+ 126 5.26451E+10 6.45125E-01 4.80051E+13 -3.45964E+07
+ 127 5.30646E+10 6.38069E-01 4.72626E+13 -3.45293E+07
+ 128 5.34841E+10 6.31144E-01 4.65379E+13 -3.44630E+07
+ 129 5.39036E+10 6.24353E-01 4.58311E+13 -3.43977E+07
+ 130 5.43231E+10 6.17714E-01 4.51442E+13 -3.43337E+07
+ 131 5.47425E+10 6.11254E-01 4.44803E+13 -3.42714E+07
+ 132 5.51620E+10 6.04994E-01 4.38419E+13 -3.42115E+07
+ 133 5.55815E+10 5.98918E-01 4.32267E+13 -3.41535E+07
+ 134 5.60010E+10 5.92934E-01 4.26233E+13 -3.40960E+07
+ 135 5.64205E+10 5.86860E-01 4.20100E+13 -3.40355E+07
+ 136 5.68400E+10 5.80455E-01 4.13583E+13 -3.39675E+07
+ 137 5.72594E+10 5.73483E-01 4.06408E+13 -3.38874E+07
+ 138 5.76789E+10 5.65785E-01 3.98401E+13 -3.37921E+07
+ 139 5.80984E+10 5.57331E-01 3.89543E+13 -3.36807E+07
+ 140 5.85179E+10 5.48212E-01 3.79959E+13 -3.35545E+07
+ 141 5.89374E+10 5.38593E-01 3.69859E+13 -3.34165E+07
+ 142 5.93569E+10 5.28659E-01 3.59468E+13 -3.32700E+07
+ 143 5.97763E+10 5.18568E-01 3.48979E+13 -3.31179E+07
+ 144 6.01958E+10 5.08438E-01 3.38527E+13 -3.29626E+07
+ 145 6.06153E+10 4.98342E-01 3.28197E+13 -3.28053E+07
+ 146 6.10348E+10 4.88323E-01 3.18039E+13 -3.26469E+07
+ 147 6.14543E+10 4.78408E-01 3.08080E+13 -3.24879E+07
+ 148 6.18738E+10 4.68611E-01 2.98335E+13 -3.23286E+07
+ 149 6.22932E+10 4.58945E-01 2.88815E+13 -3.21691E+07
+ 150 6.27127E+10 4.49418E-01 2.79527E+13 -3.20098E+07
+ 151 6.31322E+10 4.40036E-01 2.70475E+13 -3.18506E+07
+ 152 6.35517E+10 4.30802E-01 2.61660E+13 -3.16916E+07
+ 153 6.39712E+10 4.21717E-01 2.53079E+13 -3.15329E+07
+ 154 6.43907E+10 4.12784E-01 2.44732E+13 -3.13745E+07
+ 155 6.48101E+10 4.04001E-01 2.36615E+13 -3.12165E+07
+ 156 6.52296E+10 3.95367E-01 2.28724E+13 -3.10586E+07
+ 157 6.56491E+10 3.86884E-01 2.21056E+13 -3.09011E+07
+ 158 6.60686E+10 3.78549E-01 2.13607E+13 -3.07439E+07
+ 159 6.64881E+10 3.70361E-01 2.06373E+13 -3.05869E+07
+ 160 6.69076E+10 3.62320E-01 1.99350E+13 -3.04303E+07
+ 161 6.73270E+10 3.54424E-01 1.92532E+13 -3.02738E+07
+ 162 6.77465E+10 3.46670E-01 1.85917E+13 -3.01176E+07
+ 163 6.81660E+10 3.39059E-01 1.79499E+13 -2.99616E+07
+ 164 6.85855E+10 3.31587E-01 1.73275E+13 -2.98058E+07
+ 165 6.90050E+10 3.24255E-01 1.67239E+13 -2.96502E+07
+ 166 6.94245E+10 3.17059E-01 1.61387E+13 -2.94948E+07
+ 167 6.98439E+10 3.09999E-01 1.55716E+13 -2.93396E+07
+ 168 7.02634E+10 3.03072E-01 1.50221E+13 -2.91846E+07
+ 169 7.06829E+10 2.96278E-01 1.44898E+13 -2.90297E+07
+ 170 7.11024E+10 2.89615E-01 1.39742E+13 -2.88750E+07
+ 171 7.15219E+10 2.83081E-01 1.34750E+13 -2.87204E+07
+ 172 7.19414E+10 2.76674E-01 1.29916E+13 -2.85660E+07
+ 173 7.23608E+10 2.70393E-01 1.25238E+13 -2.84118E+07
+ 174 7.27803E+10 2.64235E-01 1.20711E+13 -2.82576E+07
+ 175 7.31998E+10 2.58200E-01 1.16331E+13 -2.81036E+07
+ 176 7.36193E+10 2.52285E-01 1.12094E+13 -2.79497E+07
+ 177 7.40388E+10 2.46488E-01 1.07996E+13 -2.77959E+07
+ 178 7.44583E+10 2.40807E-01 1.04033E+13 -2.76422E+07
+ 179 7.48777E+10 2.35241E-01 1.00202E+13 -2.74886E+07
+ 180 7.52972E+10 2.29787E-01 9.64978E+12 -2.73350E+07
+ 181 7.57167E+10 2.24445E-01 9.29178E+12 -2.71815E+07
+ 182 7.61362E+10 2.19211E-01 8.94583E+12 -2.70281E+07
+ 183 7.65557E+10 2.14085E-01 8.61156E+12 -2.68748E+07
+ 184 7.69752E+10 2.09064E-01 8.28864E+12 -2.67215E+07
+ 185 7.73946E+10 2.04147E-01 7.97673E+12 -2.65682E+07
+ 186 7.78141E+10 1.99332E-01 7.67551E+12 -2.64150E+07
+ 187 7.82336E+10 1.94617E-01 7.38466E+12 -2.62618E+07
+ 188 7.86531E+10 1.90000E-01 7.10385E+12 -2.61087E+07
+ 189 7.90726E+10 1.85480E-01 6.83279E+12 -2.59555E+07
+ 190 7.94921E+10 1.81055E-01 6.57116E+12 -2.58024E+07
+ 191 7.99115E+10 1.76723E-01 6.31868E+12 -2.56493E+07
+ 192 8.03310E+10 1.72482E-01 6.07505E+12 -2.54963E+07
+ 193 8.07505E+10 1.68331E-01 5.84001E+12 -2.53432E+07
+ 194 8.11700E+10 1.64269E-01 5.61327E+12 -2.51902E+07
+ 195 8.15895E+10 1.60292E-01 5.39458E+12 -2.50371E+07
+ 196 8.20090E+10 1.56401E-01 5.18368E+12 -2.48841E+07
+ 197 8.24284E+10 1.52593E-01 4.98034E+12 -2.47311E+07
+ 198 8.28479E+10 1.48867E-01 4.78429E+12 -2.45780E+07
+ 199 8.32674E+10 1.45221E-01 4.59531E+12 -2.44250E+07
+ 200 8.36869E+10 1.41654E-01 4.41317E+12 -2.42719E+07
+ 201 8.41064E+10 1.38165E-01 4.23765E+12 -2.41189E+07
+ 202 8.45259E+10 1.34751E-01 4.06852E+12 -2.39658E+07
+ 203 8.49453E+10 1.31412E-01 3.90558E+12 -2.38127E+07
+ 204 8.53648E+10 1.28145E-01 3.74862E+12 -2.36596E+07
+ 205 8.57843E+10 1.24950E-01 3.59744E+12 -2.35065E+07
+ 206 8.62038E+10 1.21826E-01 3.45186E+12 -2.33533E+07
+ 207 8.66233E+10 1.18770E-01 3.31169E+12 -2.32002E+07
+ 208 8.70428E+10 1.15782E-01 3.17675E+12 -2.30470E+07
+ 209 8.74622E+10 1.12861E-01 3.04686E+12 -2.28939E+07
+ 210 8.78817E+10 1.10005E-01 2.92186E+12 -2.27407E+07
+ 211 8.83012E+10 1.07213E-01 2.80159E+12 -2.25875E+07
+ 212 8.87207E+10 1.04483E-01 2.68588E+12 -2.24343E+07
+ 213 8.91402E+10 1.01816E-01 2.57459E+12 -2.22812E+07
+ 214 8.95597E+10 9.92094E-02 2.46757E+12 -2.21281E+07
+ 215 8.99791E+10 9.66625E-02 2.36467E+12 -2.19750E+07
+ 216 9.03986E+10 9.41743E-02 2.26577E+12 -2.18220E+07
+ 217 9.08181E+10 9.17438E-02 2.17072E+12 -2.16691E+07
+ 218 9.12376E+10 8.93702E-02 2.07941E+12 -2.15162E+07
+ 219 9.16571E+10 8.70525E-02 1.99170E+12 -2.13635E+07
+ 220 9.20765E+10 8.47899E-02 1.90747E+12 -2.12110E+07
+ 221 9.24960E+10 8.25815E-02 1.82660E+12 -2.10585E+07
+ 222 9.29155E+10 8.04263E-02 1.74899E+12 -2.09063E+07
+ 223 9.33350E+10 7.83232E-02 1.67450E+12 -2.07542E+07
+ 224 9.37545E+10 7.62712E-02 1.60302E+12 -2.06023E+07
+ 225 9.41740E+10 7.42691E-02 1.53444E+12 -2.04506E+07
+ 226 9.45934E+10 7.23156E-02 1.46864E+12 -2.02991E+07
+ 227 9.50129E+10 7.04093E-02 1.40551E+12 -2.01477E+07
+ 228 9.54324E+10 6.85487E-02 1.34492E+12 -1.99964E+07
+ 229 9.58519E+10 6.67322E-02 1.28677E+12 -1.98451E+07
+ 230 9.62714E+10 6.49581E-02 1.23093E+12 -1.96939E+07
+ 231 9.66909E+10 6.32249E-02 1.17731E+12 -1.95426E+07
+ 232 9.71103E+10 6.15308E-02 1.12579E+12 -1.93911E+07
+ 233 9.75298E+10 5.98745E-02 1.07628E+12 -1.92394E+07
+ 234 9.79493E+10 5.82546E-02 1.02869E+12 -1.90875E+07
+ 235 9.83688E+10 5.66702E-02 9.82952E+11 -1.89353E+07
+ 236 9.87883E+10 5.51204E-02 9.39001E+11 -1.87828E+07
+ 237 9.92078E+10 5.36051E-02 8.96783E+11 -1.86301E+07
+ 238 9.96272E+10 5.21242E-02 8.56255E+11 -1.84772E+07
+ 239 1.00047E+11 5.06780E-02 8.17382E+11 -1.83242E+07
+ 240 1.00466E+11 4.92670E-02 7.80134E+11 -1.81713E+07
+ 241 1.00886E+11 4.78916E-02 7.44481E+11 -1.80185E+07
+ 242 1.01305E+11 4.65525E-02 7.10392E+11 -1.78661E+07
+ 243 1.01725E+11 4.52497E-02 6.77828E+11 -1.77142E+07
+ 244 1.02144E+11 4.39830E-02 6.46742E+11 -1.75627E+07
+ 245 1.02564E+11 4.27518E-02 6.17073E+11 -1.74118E+07
+ 246 1.02983E+11 4.15547E-02 5.88749E+11 -1.72614E+07
+ 247 1.03403E+11 4.03897E-02 5.61688E+11 -1.71114E+07
+ 248 1.03822E+11 3.92545E-02 5.35798E+11 -1.69615E+07
+ 249 1.04242E+11 3.81464E-02 5.10989E+11 -1.68116E+07
+ 250 1.04661E+11 3.70627E-02 4.87174E+11 -1.66614E+07
+ 251 1.05081E+11 3.60014E-02 4.64283E+11 -1.65106E+07
+ 252 1.05500E+11 3.49609E-02 4.42260E+11 -1.63591E+07
+ 253 1.05919E+11 3.39408E-02 4.21076E+11 -1.62069E+07
+ 254 1.06339E+11 3.29413E-02 4.00719E+11 -1.60541E+07
+ 255 1.06758E+11 3.19639E-02 3.81194E+11 -1.59009E+07
+ 256 1.07178E+11 3.10100E-02 3.62512E+11 -1.57476E+07
+ 257 1.07597E+11 3.00817E-02 3.44685E+11 -1.55945E+07
+ 258 1.08017E+11 2.91805E-02 3.27721E+11 -1.54420E+07
+ 259 1.08436E+11 2.83077E-02 3.11610E+11 -1.52904E+07
+ 260 1.08856E+11 2.74636E-02 2.96328E+11 -1.51398E+07
+ 261 1.09275E+11 2.66473E-02 2.81831E+11 -1.49901E+07
+ 262 1.09695E+11 2.58571E-02 2.68059E+11 -1.48412E+07
+ 263 1.10114E+11 2.50898E-02 2.54940E+11 -1.46927E+07
+ 264 1.10534E+11 2.43418E-02 2.42399E+11 -1.45442E+07
+ 265 1.10953E+11 2.36094E-02 2.30369E+11 -1.43953E+07
+ 266 1.11373E+11 2.28893E-02 2.18797E+11 -1.42454E+07
+ 267 1.11792E+11 2.21797E-02 2.07649E+11 -1.40944E+07
+ 268 1.12212E+11 2.14804E-02 1.96915E+11 -1.39423E+07
+ 269 1.12631E+11 2.07932E-02 1.86606E+11 -1.37892E+07
+ 270 1.13051E+11 2.01210E-02 1.76742E+11 -1.36356E+07
+ 271 1.13470E+11 1.94669E-02 1.67341E+11 -1.34820E+07
+ 272 1.13890E+11 1.88331E-02 1.58416E+11 -1.33292E+07
+ 273 1.14309E+11 1.82219E-02 1.49979E+11 -1.31774E+07
+ 274 1.14729E+11 1.76357E-02 1.42041E+11 -1.30270E+07
+ 275 1.15148E+11 1.70762E-02 1.34594E+11 -1.28776E+07
+ 276 1.15568E+11 1.65418E-02 1.27581E+11 -1.27293E+07
+ 277 1.15987E+11 1.60268E-02 1.20901E+11 -1.25826E+07
+ 278 1.16407E+11 1.55258E-02 1.14479E+11 -1.24370E+07
+ 279 1.16826E+11 1.50379E-02 1.08330E+11 -1.22908E+07
+ 280 1.17246E+11 1.45644E-02 1.02503E+11 -1.21413E+07
+ 281 1.17665E+11 1.40991E-02 9.69490E+10 -1.19884E+07
+ 282 1.18085E+11 1.36288E-02 9.15155E+10 -1.18350E+07
+ 283 1.18504E+11 1.31485E-02 8.61333E+10 -1.16837E+07
+ 284 1.18923E+11 1.26718E-02 8.09398E+10 -1.15330E+07
+ 285 1.19343E+11 1.22163E-02 7.61240E+10 -1.13797E+07
+ 286 1.19762E+11 1.17839E-02 7.17099E+10 -1.12240E+07
+ 287 1.20182E+11 1.13656E-02 6.75985E+10 -1.10694E+07
+ 288 1.20601E+11 1.09590E-02 6.37323E+10 -1.09181E+07
+ 289 1.21021E+11 1.05675E-02 6.00915E+10 -1.07704E+07
+ 290 1.21440E+11 1.01901E-02 5.66298E+10 -1.06267E+07
+ 291 1.21860E+11 9.82391E-03 5.33275E+10 -1.04857E+07
+ 292 1.22279E+11 9.46989E-03 5.02148E+10 -1.03441E+07
+ 293 1.22699E+11 9.12839E-03 4.72987E+10 -1.01989E+07
+ 294 1.23118E+11 8.80032E-03 4.45619E+10 -1.00483E+07
+ 295 1.23538E+11 8.48516E-03 4.19670E+10 -9.89227E+06
+ 296 1.23957E+11 8.17224E-03 3.94112E+10 -9.73502E+06
+ 297 1.24377E+11 7.85138E-03 3.68251E+10 -9.58123E+06
+ 298 1.24796E+11 7.52908E-03 3.42816E+10 -9.43056E+06
+ 299 1.25216E+11 7.21994E-03 3.19050E+10 -9.27964E+06
+ 300 1.25635E+11 6.93413E-03 2.97691E+10 -9.12636E+06
+ 301 1.26076E+11 6.66011E-03 2.77768E+10 -8.96417E+06
+ 302 1.26560E+11 6.38054E-03 2.57953E+10 -8.78921E+06
+ 303 1.27093E+11 6.09145E-03 2.37959E+10 -8.60395E+06
+ 304 1.27679E+11 5.78592E-03 2.17409E+10 -8.40815E+06
+ 305 1.28324E+11 5.45709E-03 1.96251E+10 -8.19414E+06
+ 306 1.29034E+11 5.09522E-03 1.74480E+10 -7.94755E+06
+ 307 1.29814E+11 4.68767E-03 1.51915E+10 -7.65664E+06
+ 308 1.30672E+11 4.23909E-03 1.28979E+10 -7.32190E+06
+ 309 1.31616E+11 3.77363E-03 1.06729E+10 -6.96757E+06
+ 310 1.32655E+11 3.35419E-03 8.84545E+09 -6.61119E+06
+ 311 1.33797E+11 2.95771E-03 7.29746E+09 -6.24865E+06
+ 312 1.35054E+11 2.50423E-03 5.83764E+09 -5.79988E+06
+ 313 1.36436E+11 2.00209E-03 4.26338E+09 -5.20782E+06
+ 314 1.37957E+11 1.65718E-03 3.14724E+09 -4.66074E+06
+ 315 1.39630E+11 1.28499E-03 2.48345E+09 -4.11704E+06
+ 316 1.41470E+11 9.17612E-04 1.51142E+09 -3.34812E+06
+ 317 1.43493E+11 6.60525E-04 1.09040E+09 -2.66972E+06
+ 318 1.45720E+11 4.47073E-04 7.08716E+08 -1.78571E+06
+ 319 1.48169E+11 2.99549E-04 4.90797E+08 -7.97456E+05
+ 320 1.50862E+11 1.93898E-04 3.44200E+08 3.65393E+05
+ 321 1.53826E+11 1.18569E-04 2.40493E+08 1.70686E+06
+ 322 1.57085E+11 6.54233E-05 1.59985E+08 3.25139E+06
+ 323 1.60671E+11 3.06864E-05 9.76580E+07 5.06501E+06
+ 324 1.64614E+11 1.11861E-05 5.14685E+07 7.23895E+06
+ 325 1.68953E+11 2.81074E-06 1.86659E+07 9.55535E+06
+ 326 1.73725E+11 5.59510E-07 3.32772E+06 1.11952E+07
+ 327 1.78975E+11 8.79840E-08 4.23759E+05 1.23519E+07
+ 328 1.84749E+11 1.03500E-08 5.05319E+04 1.29185E+07
+ 329 1.91101E+11 8.93073E-10 2.66109E+04 1.28303E+07
+ 330 1.98088E+11 5.89332E-11 2.90515E+04 1.26619E+07
+ 331 2.05774E+11 3.12172E-12 3.00473E+04 1.19373E+07
+ 332 2.14228E+11 2.36918E-13 3.00429E+04 1.12451E+07
+ 333 2.23528E+11 1.27890E-13 3.00518E+04 1.03364E+07
+ 334 2.33757E+11 1.24662E-13 3.00584E+04 9.51410E+06
+ 335 2.45010E+11 1.24458E-13 3.00647E+04 8.57112E+06
+ 336 2.57388E+11 1.24377E-13 3.00752E+04 7.84800E+06
+ 337 2.71004E+11 1.24234E-13 3.00791E+04 6.94648E+06
+ 338 2.85981E+11 1.23843E-13 3.00909E+04 6.27334E+06
+ 339 3.02456E+11 1.20834E-13 2.91413E+04 4.23227E+06
+ 340 3.20578E+11 1.10508E-13 2.50406E+04 1.10077E+06
+ 341 3.40513E+11 1.03101E-13 2.15873E+04 -1.48394E+06
+ 342 3.62441E+11 1.01046E-13 2.05438E+04 -2.48993E+06
+ 343 3.86562E+11 1.00339E-13 2.01755E+04 -2.56440E+06
+ 344 4.13096E+11 1.00000E-13 2.00000E+04 -2.32120E+06
+ 345 4.42282E+11 1.00009E-13 2.00045E+04 -2.03176E+06
+ 346 4.74387E+11 1.00006E-13 2.00031E+04 -1.76922E+06
+ 347 5.09703E+11 1.00004E-13 2.00022E+04 -1.53483E+06
+ 348 5.48550E+11 1.00000E-13 2.00000E+04 -1.32666E+06
+ 349 5.91282E+11 1.00000E-13 2.00000E+04 -1.14209E+06
+ 350 6.38287E+11 1.00000E-13 2.00000E+04 -1.02375E+06