Set n / n1*n150 /;
Set s / a1*a3 /;

yupp('n1','n1','a1')=0;

Parameter fcost(m,n,s) /
$ondelim
   n2      n6      a1              1
   n6      n5      a1              1
   n6      n5      a2              1
   n2      n7      a1              1
   n7      n5      a1              1
   n2      n8      a1              1
   n8      n5      a1              1
   n2      n9      a1              1
   n2      n9      a2              1
   n9      n5      a1              1
   n5      n11     a1              1
   n11     n10     a1              1
   n10     n4      a1              1
   n10     n4      a2              1
   n5      n13     a1              1
   n5      n13     a2              1
   n13     n12     a1              1
   n12     n4      a1              1
   n12     n4      a2              1
   n2      n16     a1              1
   n16     n15     a1              1
   n2      n15     a1              1
   n15     n14     a1              1
   n2      n17     a1              1
   n17     n14     a1              1
   n17     n14     a2              1
   n2      n18     a1              1
   n18     n14     a1              1
   n14     n20     a1              1
   n20     n19     a1              1
   n20     n19     a2              1
   n14     n19     a1              1
   n19     n21     a1              1
   n21     n4      a1              1
   n19     n4      a1              1
   n19     n4      a2              1
   n4      n23     a1              1
   n4      n23     a2              1
   n23     n24     a1              1
   n23     n24     a2              1
   n24     n22     a1              1
   n4      n25     a1              1
   n25     n22     a1              1
   n4      n26     a1              1
   n26     n22     a1              1
   n4      n22     a1              1
   n22     n28     a1              1
   n22     n28     a2              1
   n28     n27     a1              1
   n22     n27     a1              1
   n27     n3      a1              1
   n27     n29     a1              1
   n29     n3      a1              1
   n27     n30     a1              1
   n30     n3      a1              1
   n4      n32     a1              1
   n32     n31     a1              1
   n32     n31     a2              1
   n4      n33     a1              1
   n33     n31     a1              1
   n31     n3      a1              1
   n31     n34     a1              1
   n34     n3      a1              1
   n31     n35     a1              1
   n35     n3      a1              1
   n4      n36     a1              1
   n4      n37     a1              1
   n37     n36     a1              1
   n37     n36     a2              1
   n36     n39     a1              1
   n36     n39     a2              1
   n39     n38     a1              1
   n38     n3      a1              1
   n38     n3      a2              1
   n2      n41     a1              1
   n41     n40     a1              1
   n40     n3      a1              1
   n40     n3      a2              1
   n2      n43     a1              1
   n2      n43     a2              1
   n43     n42     a1              1
   n42     n3      a1              1
   n42     n3      a2              1
   n2      n44     a1              1
   n44     n3      a1              1
   n2      n45     a1              1
   n45     n3      a1              1
   n2      n3      a1              1
   n2      n46     a1              1
   n46     n3      a1              1
   n2      n48     a1              1
   n48     n47     a1              1
   n47     n3      a1              1
   n2      n51     a1              1
   n2      n51     a2              1
   n51     n50     a1              1
   n51     n52     a1              1
   n52     n50     a1              1
   n52     n50     a2              1
   n50     n53     a1              1
   n50     n53     a2              1
   n53     n54     a1              1
   n53     n54     a2              1
   n54     n49     a1              1
   n49     n3      a1              1
   n49     n55     a1              1
   n55     n3      a1              1
   n49     n56     a1              1
   n56     n3      a1              1
   n49     n57     a1              1
   n57     n3      a1              1
   n49     n58     a1              1
   n58     n3      a1              1
   n2      n62     a1              1
   n62     n61     a1              1
   n61     n60     a1              1
   n61     n60     a2              1
   n2      n64     a1              1
   n64     n63     a1              1
   n2      n63     a1              1
   n63     n60     a1              1
   n63     n60     a2              1
   n60     n65     a1              1
   n65     n59     a1              1
   n60     n66     a1              1
   n60     n66     a2              1
   n66     n59     a1              1
   n60     n67     a1              1
   n60     n67     a2              1
   n67     n68     a1              1
   n68     n59     a1              1
   n67     n59     a1              1
   n59     n70     a1              1
   n70     n69     a1              1
   n59     n71     a1              1
   n59     n71     a2              1
   n71     n69     a1              1
   n59     n72     a1              1
   n59     n72     a2              1
   n72     n69     a1              1
   n69     n73     a1              1
   n73     n3      a1              1
   n69     n74     a1              1
   n74     n3      a1              1
   n69     n3      a1              1
   n69     n75     a1              1
   n75     n3      a1              1
   n69     n76     a1              1
   n69     n76     a2              1
   n76     n3      a1              1
   n3      n79     a1              1
   n79     n80     a1              1
   n80     n78     a1              1
   n80     n78     a2              1
   n3      n81     a1              1
   n3      n81     a2              1
   n81     n78     a1              1
   n81     n82     a1              1
   n82     n78     a1              1
   n3      n84     a1              1
   n84     n83     a1              1
   n3      n85     a1              1
   n85     n83     a1              1
   n3      n83     a1              1
   n83     n86     a1              1
   n86     n78     a1              1
   n83     n87     a1              1
   n83     n87     a2              1
   n87     n78     a1              1
   n78     n88     a1              1
   n88     n77     a1              1
   n78     n89     a1              1
   n89     n77     a1              1
   n78     n77     a1              1
   n78     n90     a1              1
   n78     n90     a2              1
   n90     n91     a1              1
   n91     n77     a1              1
   n78     n92     a1              1
   n78     n93     a1              1
   n93     n92     a1              1
   n78     n94     a1              1
   n94     n92     a1              1
   n92     n77     a1              1
   n92     n77     a2              1
   n92     n95     a1              1
   n95     n77     a1              1
   n77     n97     a1              1
   n97     n96     a1              1
   n96     n1      a1              1
   n96     n98     a1              1
   n98     n1      a1              1
   n77     n99     a1              1
   n99     n1      a1              1
   n77     n100    a1              1
   n100    n1      a1              1
   n100    n1      a2              1
   n77     n102    a1              1
   n77     n102    a2              1
   n102    n103    a1              1
   n102    n103    a2              1
   n103    n101    a1              1
   n101    n104    a1              1
   n101    n104    a2              1
   n104    n1      a1              1
   n104    n1      a2              1
   n77     n105    a1              1
   n77     n106    a1              1
   n106    n105    a1              1
   n105    n107    a1              1
   n105    n107    a2              1
   n107    n1      a1              1
   n77     n108    a1              1
   n77     n108    a2              1
   n108    n109    a1              1
   n109    n1      a1              1
   n77     n111    a1              1
   n111    n110    a1              1
   n111    n110    a2              1
   n77     n112    a1              1
   n112    n110    a1              1
   n110    n1      a1              1
   n110    n1      a2              1
   n110    n113    a1              1
   n113    n1      a1              1
   n3      n117    a1              1
   n117    n116    a1              1
   n117    n116    a2              1
   n3      n118    a1              1
   n118    n116    a1              1
   n3      n119    a1              1
   n119    n116    a1              1
   n3      n120    a1              1
   n120    n116    a1              1
   n120    n116    a2              1
   n116    n121    a1              1
   n121    n122    a1              1
   n122    n115    a1              1
   n122    n115    a2              1
   n116    n124    a1              1
   n116    n124    a2              1
   n124    n123    a1              1
   n123    n115    a1              1
   n123    n115    a2              1
   n115    n127    a1              1
   n127    n126    a1              1
   n115    n126    a1              1
   n115    n126    a2              1
   n115    n126    a3              1
   n126    n129    a1              1
   n126    n129    a2              1
   n129    n128    a1              1
   n128    n125    a1              1
   n128    n125    a2              1
   n125    n130    a1              1
   n125    n130    a2              1
   n130    n114    a1              1
   n130    n114    a2              1
   n125    n131    a1              1
   n131    n132    a1              1
   n132    n114    a1              1
   n132    n114    a2              1
   n114    n133    a1              1
   n114    n134    a1              1
   n134    n133    a1              1
   n133    n135    a1              1
   n133    n135    a2              1
   n135    n1      a1              1
   n114    n136    a1              1
   n114    n136    a2              1
   n136    n137    a1              1
   n137    n1      a1              1
   n114    n139    a1              1
   n114    n139    a2              1
   n139    n138    a1              1
   n139    n138    a2              1
   n138    n140    a1              1
   n138    n141    a1              1
   n141    n140    a1              1
   n140    n1      a1              1
   n140    n1      a2              1
   n114    n143    a1              1
   n114    n143    a2              1
   n143    n142    a1              1
   n142    n1      a1              1
   n142    n1      a2              1
   n114    n144    a1              1
   n114    n145    a1              1
   n145    n144    a1              1
   n144    n1      a1              1
   n114    n147    a1              1
   n114    n147    a2              1
   n147    n148    a1              1
   n148    n146    a1              1
   n148    n146    a2              1
   n146    n149    a1              1
   n149    n1      a1              1
   n146    n1      a1              1
   n146    n150    a1              1
   n150    n1      a1              1
$offdelim
/;

Parameter vcost(m,n,s) /
$ondelim
   n2      n6      a1              0
   n6      n5      a1              0
   n6      n5      a2              0
   n2      n7      a1              0
   n7      n5      a1              0
   n2      n8      a1              0
   n8      n5      a1              0
   n2      n9      a1              0
   n2      n9      a2              0
   n9      n5      a1              0
   n5      n11     a1              0
   n11     n10     a1              0
   n10     n4      a1              0
   n10     n4      a2              0
   n5      n13     a1              0
   n5      n13     a2              0
   n13     n12     a1              0
   n12     n4      a1              0
   n12     n4      a2              0
   n2      n16     a1              0
   n16     n15     a1              0
   n2      n15     a1              0
   n15     n14     a1              0
   n2      n17     a1              0
   n17     n14     a1              0
   n17     n14     a2              0
   n2      n18     a1              0
   n18     n14     a1              0
   n14     n20     a1              0
   n20     n19     a1              0
   n20     n19     a2              0
   n14     n19     a1              0
   n19     n21     a1              0
   n21     n4      a1              0
   n19     n4      a1              0
   n19     n4      a2              0
   n4      n23     a1              0
   n4      n23     a2              0
   n23     n24     a1              0
   n23     n24     a2              0
   n24     n22     a1              0
   n4      n25     a1              0
   n25     n22     a1              0
   n4      n26     a1              0
   n26     n22     a1              0
   n4      n22     a1              0
   n22     n28     a1              0
   n22     n28     a2              0
   n28     n27     a1              0
   n22     n27     a1              0
   n27     n3      a1              0
   n27     n29     a1              0
   n29     n3      a1              0
   n27     n30     a1              0
   n30     n3      a1              0
   n4      n32     a1              0
   n32     n31     a1              0
   n32     n31     a2              0
   n4      n33     a1              0
   n33     n31     a1              0
   n31     n3      a1              0
   n31     n34     a1              0
   n34     n3      a1              0
   n31     n35     a1              0
   n35     n3      a1              0
   n4      n36     a1              0
   n4      n37     a1              0
   n37     n36     a1              0
   n37     n36     a2              0
   n36     n39     a1              0
   n36     n39     a2              0
   n39     n38     a1              0
   n38     n3      a1              0
   n38     n3      a2              0
   n2      n41     a1              0
   n41     n40     a1              0
   n40     n3      a1              0
   n40     n3      a2              0
   n2      n43     a1              0
   n2      n43     a2              0
   n43     n42     a1              0
   n42     n3      a1              0
   n42     n3      a2              0
   n2      n44     a1              0
   n44     n3      a1              0
   n2      n45     a1              0
   n45     n3      a1              0
   n2      n3      a1              0
   n2      n46     a1              0
   n46     n3      a1              0
   n2      n48     a1              0
   n48     n47     a1              0
   n47     n3      a1              0
   n2      n51     a1              0
   n2      n51     a2              0
   n51     n50     a1              0
   n51     n52     a1              0
   n52     n50     a1              0
   n52     n50     a2              0
   n50     n53     a1              0
   n50     n53     a2              0
   n53     n54     a1              0
   n53     n54     a2              0
   n54     n49     a1              0
   n49     n3      a1              0
   n49     n55     a1              0
   n55     n3      a1              0
   n49     n56     a1              0
   n56     n3      a1              0
   n49     n57     a1              0
   n57     n3      a1              0
   n49     n58     a1              0
   n58     n3      a1              0
   n2      n62     a1              0
   n62     n61     a1              0
   n61     n60     a1              0
   n61     n60     a2              0
   n2      n64     a1              0
   n64     n63     a1              0
   n2      n63     a1              0
   n63     n60     a1              0
   n63     n60     a2              0
   n60     n65     a1              0
   n65     n59     a1              0
   n60     n66     a1              0
   n60     n66     a2              0
   n66     n59     a1              0
   n60     n67     a1              0
   n60     n67     a2              0
   n67     n68     a1              0
   n68     n59     a1              0
   n67     n59     a1              0
   n59     n70     a1              0
   n70     n69     a1              0
   n59     n71     a1              0
   n59     n71     a2              0
   n71     n69     a1              0
   n59     n72     a1              0
   n59     n72     a2              0
   n72     n69     a1              0
   n69     n73     a1              0
   n73     n3      a1              0
   n69     n74     a1              0
   n74     n3      a1              0
   n69     n3      a1              0
   n69     n75     a1              0
   n75     n3      a1              0
   n69     n76     a1              0
   n69     n76     a2              0
   n76     n3      a1              0
   n3      n79     a1              0
   n79     n80     a1              0
   n80     n78     a1              0
   n80     n78     a2              0
   n3      n81     a1              0
   n3      n81     a2              0
   n81     n78     a1              0
   n81     n82     a1              0
   n82     n78     a1              0
   n3      n84     a1              0
   n84     n83     a1              0
   n3      n85     a1              0
   n85     n83     a1              0
   n3      n83     a1              0
   n83     n86     a1              0
   n86     n78     a1              0
   n83     n87     a1              0
   n83     n87     a2              0
   n87     n78     a1              0
   n78     n88     a1              0
   n88     n77     a1              0
   n78     n89     a1              0
   n89     n77     a1              0
   n78     n77     a1              0
   n78     n90     a1              0
   n78     n90     a2              0
   n90     n91     a1              0
   n91     n77     a1              0
   n78     n92     a1              0
   n78     n93     a1              0
   n93     n92     a1              0
   n78     n94     a1              0
   n94     n92     a1              0
   n92     n77     a1              0
   n92     n77     a2              0
   n92     n95     a1              0
   n95     n77     a1              0
   n77     n97     a1              0
   n97     n96     a1              0
   n96     n1      a1              0
   n96     n98     a1              0
   n98     n1      a1              0
   n77     n99     a1              0
   n99     n1      a1              0
   n77     n100    a1              0
   n100    n1      a1              0
   n100    n1      a2              0
   n77     n102    a1              0
   n77     n102    a2              0
   n102    n103    a1              0
   n102    n103    a2              0
   n103    n101    a1              0
   n101    n104    a1              0
   n101    n104    a2              0
   n104    n1      a1              0
   n104    n1      a2              0
   n77     n105    a1              0
   n77     n106    a1              0
   n106    n105    a1              0
   n105    n107    a1              0
   n105    n107    a2              0
   n107    n1      a1              0
   n77     n108    a1              0
   n77     n108    a2              0
   n108    n109    a1              0
   n109    n1      a1              0
   n77     n111    a1              0
   n111    n110    a1              0
   n111    n110    a2              0
   n77     n112    a1              0
   n112    n110    a1              0
   n110    n1      a1              0
   n110    n1      a2              0
   n110    n113    a1              0
   n113    n1      a1              0
   n3      n117    a1              0
   n117    n116    a1              0
   n117    n116    a2              0
   n3      n118    a1              0
   n118    n116    a1              0
   n3      n119    a1              0
   n119    n116    a1              0
   n3      n120    a1              0
   n120    n116    a1              0
   n120    n116    a2              0
   n116    n121    a1              0
   n121    n122    a1              0
   n122    n115    a1              0
   n122    n115    a2              0
   n116    n124    a1              0
   n116    n124    a2              0
   n124    n123    a1              0
   n123    n115    a1              0
   n123    n115    a2              0
   n115    n127    a1              0
   n127    n126    a1              0
   n115    n126    a1              0
   n115    n126    a2              0
   n115    n126    a3              0
   n126    n129    a1              0
   n126    n129    a2              0
   n129    n128    a1              0
   n128    n125    a1              0
   n128    n125    a2              0
   n125    n130    a1              0
   n125    n130    a2              0
   n130    n114    a1              0
   n130    n114    a2              0
   n125    n131    a1              0
   n131    n132    a1              0
   n132    n114    a1              0
   n132    n114    a2              0
   n114    n133    a1              0
   n114    n134    a1              0
   n134    n133    a1              0
   n133    n135    a1              0
   n133    n135    a2              0
   n135    n1      a1              0
   n114    n136    a1              0
   n114    n136    a2              0
   n136    n137    a1              0
   n137    n1      a1              0
   n114    n139    a1              0
   n114    n139    a2              0
   n139    n138    a1              0
   n139    n138    a2              0
   n138    n140    a1              0
   n138    n141    a1              0
   n141    n140    a1              0
   n140    n1      a1              0
   n140    n1      a2              0
   n114    n143    a1              0
   n114    n143    a2              0
   n143    n142    a1              0
   n142    n1      a1              0
   n142    n1      a2              0
   n114    n144    a1              0
   n114    n145    a1              0
   n145    n144    a1              0
   n144    n1      a1              0
   n114    n147    a1              0
   n114    n147    a2              0
   n147    n148    a1              0
   n148    n146    a1              0
   n148    n146    a2              0
   n146    n149    a1              0
   n149    n1      a1              0
   n146    n1      a1              0
   n146    n150    a1              0
   n150    n1      a1              0
$offdelim
/;

Parameter demand(n) /
   n1           1126
   n3            184
   n5             -2
   n15           -21
   n18            -7
   n47            -3
   n49             1
   n53           -49
   n55             1
   n59            21
   n60             7
   n61             2
   n62          -139
   n63             3
   n64          -529
   n65             1
   n75            -9
   n77            46
   n80             1
   n82            -1
   n88          -202
   n94             1
   n96            20
   n97            12
   n98             4
   n99             1
   n100            1
   n115            1
   n116           -1
   n121            1
   n124           -2
   n137          -21
   n138            1
   n139         -398
   n141            1
   n146          -55
   n149            2
   n150            1
/;

* Upper bound for flow on any arc
u = -sum(n$(demand(n)<0), demand(n));
xupp(m,n,s)$(fcost(m,n,s) or vcost(m,n,s)) = u;