$ontext
    Adhya4 pooling problem data.
    Author: Mohammed Alfaki, Wed Nov  9 15:36:07 2011
$offtext

$eolcom #

# Declare sets
    set i    / 1*15 /;
    set s(i) / 1*8  /;
    set t(i) /11*15 /;
    set k    / 1*4  /;

alias (i,j);

# The arc unit cost c_{ij}
table c(i,j)
          9      10      11      12      13      14      15
  1   15.00    0.00    0.00    0.00    0.00    0.00    0.00
  2    7.00    0.00    0.00    0.00    0.00    0.00    0.00
  3    4.00    0.00    0.00    0.00    0.00    0.00    0.00
  4    5.00    0.00    0.00    0.00    0.00    0.00    0.00
  5    0.00    6.00    0.00    0.00    0.00    0.00    0.00
  6    0.00    3.00    0.00    0.00    0.00    0.00    0.00
  7    0.00    5.00    0.00    0.00    0.00    0.00    0.00
  8    0.00    5.00    0.00    0.00    0.00    0.00    0.00
  9    0.00    0.00  -10.00  -25.00  -30.00   -6.00  -10.00
 10    0.00    0.00  -10.00  -25.00  -30.00   -6.00  -10.00 ;

# The adjacency matrix (the arcs set A)
table a(i,j)
      9  10  11  12  13  14  15
  1   1   0   0   0   0   0   0
  2   1   0   0   0   0   0   0
  3   1   0   0   0   0   0   0
  4   1   0   0   0   0   0   0
  5   0   1   0   0   0   0   0
  6   0   1   0   0   0   0   0
  7   0   1   0   0   0   0   0
  8   0   1   0   0   0   0   0
  9   0   0   1   1   1   1   1
 10   0   0   1   1   1   1   1 ;

# Source qualities/terminal quality upper bounds
table q(i,k)
          1       2       3       4
  1    0.50    1.90    1.30    1.00
  2    1.40    1.80    1.70    1.60
  3    1.20    1.90    1.40    1.40
  4    1.50    1.20    1.70    1.30
  5    1.60    1.80    1.60    2.00
  6    1.20    1.10    1.40    2.00
  7    1.50    1.50    1.50    1.50
  8    1.40    1.60    1.20    1.60
 11    1.20    1.70    1.40    1.70
 12    1.40    1.30    1.80    1.40
 13    1.30    1.30    1.90    1.90
 14    1.20    1.10    1.70    1.60
 15    1.60    1.90    2.00    2.50 ;

# Node capacity lower bound
parameter bl(i) /  1     0.00
                   2     0.00
                   3     0.00
                   4     0.00
                   5     0.00
                   6     0.00
                   7     0.00
                   8     0.00
                  11     0.00
                  12     0.00
                  13     0.00
                  14     0.00
                  15     0.00 / ;

# Node capacity upper bound
parameter bu(i) /  1    85.00
                   2    85.00
                   3    85.00
                   4    85.00
                   5    85.00
                   6    85.00
                   7    85.00
                   8    85.00
                   9    85.00
                  10    85.00
                  11    15.00
                  12    25.00
                  13    10.00
                  14    20.00
                  15    15.00 / ;

$include xmodel.gms