import numpy as np
import sys
import os
import time

# my modules
sys.path.insert(0,'/home/user/home/advection-streams/')
import hydromixing.mixer as mixer
import hydromixing.run as run

sys.path.insert(0,'/home/user/home/PyPPM/')
from ppmpy import ppm

# runid
runid = 'N16'

# dictionary keys
shellKeys = mixer.Stream.advectionData
shellHKeys = mixer.Stream.advectionHeader
centralKeys = ("X", "(1-X)")

# path to where I save everything
substr = 'standard'
repeats = 100
subdir = '{0:s}-{1:s}'.format(runid, substr)
initbase = os.path.join(subdir, '{0:s}-{1:s}'.format(runid, substr))
dataDir = '{0:s}-{1:s}/'.format(runid, substr)

# create the run object
simulation = run.Simulation(initbase, 'advection', burnFunc='PPMStarBurn', dataDir=dataDir, initdump=100)

# time simulation
start_time = time.time()
simulation.run( repeat0=repeats)

print('The simulation took {0:.0f} seconds to run {1:d} dumps'.format(time.time() - start_time, len(simulation.allDumps)))