Hi,
for cons3 can you try to protect out of bounds for Time with
forall(k in Freight, j in DCities, m in Mode, t in Periods )
cons3:
sum(i in CCities:
(intValue(i) in 1..nbMode) &&
(t-Time[intValue(i)][intValue(j)][intValue(m)] in Periods) )
(x[i][j][k][m][t-Time[intValue(i)][intValue(j)][intValue(m)]]) ==Demand[intValue(k)][intValue(j)][t]- e [k][j][t];
Regards
#DecisionOptimization#OPLusingCPLEXOptimizer