diff --git a/examples/symba_energy_momentum/hitandrun.in b/examples/symba_energy_momentum/hitandrun.in
new file mode 100644
index 000000000..45f3654f8
--- /dev/null
+++ b/examples/symba_energy_momentum/hitandrun.in
@@ -0,0 +1,13 @@
+2
+1 1e-07 0.0009
+7e-06
+1.0 -4.20E-05 0.0
+0.00 6.28 0.0
+0.4 0.4 0.4 !Ip
+0.0 0.0 6.0e4 !rot
+2 7e-10 0.0004
+3.25e-06
+1.0 4.20E-05 0.0
+-1.50 -6.28 0.0
+0.4 0.4 0.4 !Ip
+0.0 0.0 1.0e5 !rot
diff --git a/examples/symba_energy_momentum/param.hitandrun.in b/examples/symba_energy_momentum/param.hitandrun.in
new file mode 100644
index 000000000..c4c29a05b
--- /dev/null
+++ b/examples/symba_energy_momentum/param.hitandrun.in
@@ -0,0 +1,32 @@
+T0             0.0e0
+TSTOP          0.000100           ! simulation length in seconds = 100 years
+DT             0.0000001              ! stepsize in seconds
+CB_IN          cb.in
+PL_IN          hitandrun.in
+TP_IN          tp.in
+IN_TYPE        ASCII
+ISTEP_OUT      1                 ! output cadence every year
+BIN_OUT        bin.hitandrun.nc
+PARTICLE_OUT   particle.hitandrun.dat
+OUT_TYPE       NETCDF_DOUBLE     ! double precision real output
+OUT_FORM       XVEL              ! osculating element output
+OUT_STAT       REPLACE
+ISTEP_DUMP     1                 ! system dump cadence
+CHK_CLOSE      yes               ! check for planetary close encounters
+CHK_RMIN       0.005
+CHK_RMAX       1e6
+CHK_EJECT      -1.0              ! ignore this check
+CHK_QMIN       -1.0              ! ignore this check
+ENC_OUT        enc.hitandrun.dat
+DISCARD_OUT    discard.hitandrun.out
+EXTRA_FORCE    no                ! no extra user-defined forces
+BIG_DISCARD    no                ! output all planets if anything discarded
+RHILL_PRESENT  yes               ! Hill's sphere radii in input file
+GMTINY          1.0e-16
+FRAGMENTATION  yes
+MU2KG          1.98908e30
+TU2S           3.1556925e7
+DU2M           1.49598e11
+ENERGY         yes
+ROTATION       yes
+SEED           8 933097 -220886113 -118730874 233084005 32111237 -823335422 524551114 -61162322