Updated initial conditions to compute energy

MintonGroup · Oct 21, 2022 · 36a69b3 · 36a69b3
1 parent eb770cd
commit 36a69b3
Showing 1 changed file with 5 additions and 4 deletions.
diff --git a/examples/helio_neptune_pluto/init_cond.py b/examples/helio_neptune_pluto/init_cond.py
@@ -1,7 +1,7 @@
 #!/usr/bin/env python3
 import swiftest
 
-sim = swiftest.Simulation()
+sim = swiftest.Simulation(readbin=False)
 sim.param['IN_TYPE'] = "NETCDF_DOUBLE"
 sim.param['NC_IN'] = "init_cond.nc"
 sim.param['IN_FORM'] = "XV"
@@ -15,8 +15,8 @@
 sim.param['T0'] = 0.0
 sim.param['DT'] = 5.0 
 sim.param['TSTOP'] = 1.e5
-sim.param['ISTEP_OUT']  = 100
-sim.param['ISTEP_DUMP'] = 100
+sim.param['ISTEP_OUT']  = 20 
+sim.param['ISTEP_DUMP'] = 20
 sim.param['CHK_QMIN_COORD'] = "HELIO"
 sim.param['CHK_QMIN'] = swiftest.RSun / swiftest.AU2M
 sim.param['CHK_QMIN_RANGE'] = f"{swiftest.RSun / swiftest.AU2M} 1000.0"
@@ -25,6 +25,7 @@
 sim.param['CHK_EJECT'] = 1000.0
 sim.param['OUT_STAT'] = "REPLACE"
 sim.param['RHILL_PRESENT'] = "NO"
+sim.param['ENERGY'] = "YES"
 sim.param['GR'] = 'NO'
 sim.param['CHK_CLOSE'] = "NO"
 
@@ -36,7 +37,7 @@
 
 for name, id in bodyid.items():
    sim.add(name, idval=id, date="2027-04-30")
-   
+
 sim.save("param.in")