exp(850)How do I solve this problem,

You can add a constraint to your AMPL model to set a lower bound or an upper bound on a variable, or on an expression involving variables. For example,

s.t. xr_bound {i in D3 union D333, j in IP}:
   xr[i,j] + 12 <= 0;