
      include("/MDVRP/src/run.jl")
      args = ["/MDVRP/data/toy.vrp", "-r"]
      src_folder = "/MDVRP/src"
      app = parse_commandline(args, src_folder);
      
      data = readMDVRPData(app);
      data.nb_depots
      data.depot_ids
      customers(data)
      edges(data)
       
      Add the parameter depot_id to the function build_graph and adapt this function to use the depot_id instead of 0 (you should modify the construction of the VrpGraph, and adjust the mapping of the x variables).
      Create multiple graphs (one per depot) using a for loop over data.depot_ids
      Replace lines 65-67 by lines 69-73. This step adjusts the definition of packing sets, capacity cuts separator and distance matrices to take into account multiple graphs.

      include("/MDVRP/src/model.jl") 
      (model, x) = build_model(data, app);
      print(model.formulation) 'check that formulation is ok'
      print(model.graphs) 'two graphs'

      main(args)   
