From 36a69b3f024c9db1f83961c57a7877794ff59369 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Fri, 21 Oct 2022 09:44:43 -0400 Subject: [PATCH] Updated initial conditions to compute energy --- examples/helio_neptune_pluto/init_cond.py | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/examples/helio_neptune_pluto/init_cond.py b/examples/helio_neptune_pluto/init_cond.py index c3c6c96c7..975003c3a 100755 --- 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")