distances = pd.read_csv("distances.csv", header=None).T.to_dict() D = {} for j in range(0, 32): for i in range(0, 13): D[(i, j)] = distances[j][i] model.D = Param(model.i, model.j, initialize=D, doc='Distances') demands = pd.read_csv("demands.csv", header=None).T.to_dict() CD = {} for j in range(0, 32): for k in ['A', 'B']: if k == 'A': CD[(j, k)] = demands[j][2]*1000 if k == 'B': CD[(j, k)] = demands[j][5]*1000 model.CD = Param(model.j, model.k, initialize=CD, doc='Demands') MC = { 0: 1404000, 1: 1452000, 2: 1476000, 3: 1482000, 4: 1500000, 5: 1488000, 6: 1434000, 7: 1356000, 8: 1464000, 9: 1470000, 10: 1458000, 11: 1392000, 12: 1392000 }