Add ability to constrain model using observations in any box

[mrghg]

At present, I think it's only possible to use surface measurements. Should be able to use any number of boxes.