# data file for c-reload.mod # AMPL coding: S. Leyffer, University of Dundee, April 2000. # ... definition of sets set I := 1 2 3 4 5 6 7 8 9 10 11 12 13 14 ; set L := 1 2 3 4 ; set M := 1 2 3 ; set T := 1 2 3 4 5 6 ; # .. scalar parameters param conspw := 380; param kfresh := 1.25; param alfa := 6e-06; param flim := 1.8; param cytime := 350; # G_ij matrix in sparse format param: G := 1 1 0.63 1 2 0.26 1 3 0.023 1 4 0.007 1 7 0.05 1 8 0.02 1 11 0.003 2 1 0.13 2 2 0.572 2 3 0.118 2 4 0.015 2 5 0.007 2 7 0.09 2 8 0.048 2 9 0.01 2 11 0.007 2 12 0.003 3 1 0.015 3 2 0.132 3 3 0.518 3 4 0.115 3 5 0.015 3 6 0.007 3 7 0.047 3 8 0.083 3 9 0.04 3 10 0.01 3 11 0.008 3 12 0.007 3 13 0.003 4 2 0.022 4 3 0.115 4 4 0.515 4 5 0.115 4 6 0.015 4 7 0.01 4 8 0.05 4 9 0.083 4 10 0.04 4 12 0.008 4 13 0.007 5 3 0.015 5 4 0.115 5 5 0.515 5 6 0.115 5 7 0.007 5 8 0.01 5 9 0.047 5 10 0.083 5 11 0.003 5 13 0.008 6 4 0.015 6 5 0.115 6 6 0.515 6 8 0.007 6 9 0.01 6 10 0.047 6 12 0.003 7 1 0.057 7 2 0.166 7 3 0.087 7 4 0.013 7 5 0.007 7 7 0.44 7 8 0.157 7 9 0.016 7 11 0.04 7 12 0.014 7 14 0.003 8 1 0.01 8 2 0.055 8 3 0.097 8 4 0.043 8 5 0.01 8 6 0.007 8 7 0.075 8 8 0.48 8 9 0.075 8 10 0.008 8 11 0.075 8 12 0.048 8 13 0.007 8 14 0.007 9 2 0.01 9 3 0.054 9 4 0.083 9 5 0.04 9 6 0.01 9 7 0.008 9 8 0.082 9 9 0.443 9 10 0.075 9 11 0.04 9 12 0.075 9 13 0.04 9 14 0.008 10 3 0.01 10 4 0.047 10 5 0.083 10 6 0.04 10 8 0.015 10 9 0.075 10 10 0.44 10 11 0.007 10 12 0.043 10 13 0.075 11 1 0.003 11 2 0.014 11 3 0.016 11 4 0.014 11 5 0.006 11 7 0.04 11 8 0.15 11 9 0.08 11 10 0.007 11 11 0.44 11 12 0.15 11 13 0.016 11 14 0.04 12 2 0.003 12 3 0.007 12 4 0.015 12 6 0.003 12 7 0.007 12 8 0.048 12 9 0.082 12 10 0.043 12 11 0.075 12 12 0.48 12 13 0.075 12 14 0.075 13 3 0.003 13 4 0.007 13 5 0.008 13 8 0.007 13 9 0.047 13 10 0.075 13 11 0.008 13 12 0.082 13 13 0.443 13 14 0.04 14 5 0.007 14 7 0.003 14 8 0.014 14 9 0.016 14 10 0.007 14 11 0.04 14 12 0.15 14 13 0.08 14 14 0.44; # ... diagonal pairs let d1[1] := 1; let d2[1] := 7; let d1[2] := 11; let d2[2] := 14; let {n in ND} v[ d1[n] ] := 0.5; let {n in ND} v[ d2[n] ] := 0.5; # ... initial Keff value let{t in T} keff[t] := 0.956145;