arbitrary phase shift now added to model
model output now plotted on top of data
pymc script now working with new advection routine module
arbitrary phase shift added to velocity and bouyancy functions