import numpy as np
import model3
calc = {	'num_vars':6,
		'names': ['volume','uvalue','insideEnv','outside','cond','n'],
		'bounds': [[20,500],[0.8,1],[19,21],[-5,-3],[4.7,4.9],[0.9,1.1]]
	}
from SALib.sample import saltelli
inputs = saltelli.sample(calc, 1000)
np.savetxt("inputs.csv", inputs)
outputs = model3.evaluate(inputs)
np.savetxt("outputs.csv", outputs)
from SALib.analyze import sobol
Si = sobol.analyze(calc,outputs)
print(Si['S1'])
print(Si['ST'])