The whole model comprises 14335 continuous variables, 3372 binary variables, 4421 equality and 7524 inequality constraints. The GAMS package (Brooke et al.,
1998) was used in order to implement the MILP optimization model and generate their solutions. The OSL (IBM, 1991) solver within GAMS was used to solve the
problem.