param b{i in D, p in P, (i,j) in A} = min(d[i,p: p = E[k]],u[i,j]);
The problem here is that the first line defining the "Upper bound" is not mathematically well-formed. Thus any attempt to transcribe it into AMPL will not work. Probably you just want to write
param b {(i,j) in A, p in P} = min (d[i,p], u[i,j]);
param E {k in K} symbolic in P;
and then your constraint can be x[i,j,k] <= b[i,j,E[k]] * z[i,j,k].
Bob Fourer
am...@googlegroups.com
s.t. MaxFlow1{k in K, (i,j) in A}: x[i,j,k] <= b[i,j,E[k]] * z[i,j,k];
s.t. MaxFlow1{k in K: E[k] = p, (i,j) in A}: x[i,j,k] <= b[i,j,E[k]] * z[i,j,k];