import numpy as np
import bungee

calc = {'num_vars': 3,
	'names': ['H', 'M', 'S'],
	'bounds': [[40, 60],[67, 74],[20, 40]]
	}

from SALib.sample import saltelli

inputs = saltelli.sample(calc, 1000)
np.savetxt("inputs.csv", inputs)
outputs = bungee.evaluate(inputs)
np.savetxt("outputs.csv", outputs)

from SALib.analyze import sobol

Si = sobol.analyze(calc,outputs)
print(Si['S1'])
print(Si['ST'])