import numpy as np
import model2
calc = {	'num_vars':6,
		'names': ['volume','uvalue','insideEnv','outside','cond','n'],
		'bounds': [[999, 1000],[0.5,5],[18,22],[-5,5],[4.7,4.9],[0.5,1.5]]
	}
from SALib.sample import saltelli
inputs = saltelli.sample(calc, 1000)
np.savetxt("inputs.csv", inputs)
outputs = model2.evaluate(inputs)
np.savetxt("outputs.csv", outputs)
from SALib.analyze import sobol
Si = sobol.analyze(calc,outputs)
print(Si['S1'])
print(Si['ST'])