I want to limit my choice set to certain index combinations and am therefore trying to declare equations over a dynamic multidimensional set.

My model ran until I tried limiting the set so my feeling is my syntax is wrong. Could you please hep me on how to get the syntax right? As my model is pretty lengthy and all equations are only to be calculated for the members of the dynamic set I am only giving you the code for one equation. I hope it's sufficient.

Code: Select all

```
sets
i Zone (district) /suburb1, city, suburb2/,
m mode /car, UAM/
;
alias (i,j,k,l,ii,jj,kk,ll);
alias (m,mm,mmm);
set
ijm(i,j,m) /(suburb1,city,suburb2).(suburb1,city,suburb2).car, suburb1.(suburb1,city).UAM, city. (suburb1,city,suburb2). UAM, suburb2.(city,suburb2).UAM/
ijkm(i,j,k,m) /(suburb1,city,suburb2).(suburb1,city,suburb2).(suburb1,city,suburb2).car, suburb1.(suburb1,city).(suburb1,city,suburb2).UAM, city. (suburb1,city,suburb2).(suburb1,city,suburb2).UAM, suburb2.(city,suburb2).(suburb1,city,suburb2).UAM /;
variable
q(i,j,m) Land demand low skilled Household Type ij (m^2);
equations
Eq_Land_demand_q(i,j,m) land demand -> q,
Eq_Land_demand_q(ijm(i,j,m)).. q(ijm) =e= phi*(FullInc(i,j,m))/r(i) ;
option mcp=path;
option iterlim=900000;
model UrbanCGE
/
Eq_Land_demand_q.q;
solve UrbanCGE using mcp;
```

Thanks for taking your time to look at this!

Best,

Anna