Updating posteriors for scanning a chemical system

Thanks @mjclarke94,

Do you have a python function that takes care of the parsing and calling GULP via subprocess??

Basically we need something like

def simulator(a, b, c):
    create input for gulp
    call gulp with the created input
    parse outuput
    return results

where a, b, c are all the necessary parameters (not necessary those names or number) for GULP and results is a numpy array. Also you will need a observed array with the same shape as results

If attaching files is giving you trouble you can create a gist or send me the files by e-mail.