Refactor `create_model` function

[ojdo]

Right now, this function does too much. While there is probably no harm in defining all sets, parameters, variables, objective and constraints in one long function, the pre-processing steps, most notably deriving m.peak at the top, should be moved to a separate function, similar to the prepare_edge function in file runhg15.py.