lemon-project-template-glpk

annotate deps/glpk/examples/dist.mod @ 9:33de93886c88

find changesets by author, revision, files, or words in the commit message
Import GLPK 4.47
author Alpar Juttner <alpar@cs.elte.hu>
date Sun, 06 Nov 2011 20:59:10 +0100
parents
children
rev   line source
alpar@9 1 # DIST, a product distribution model
alpar@9 2 #
alpar@9 3 # References:
alpar@9 4 # Robert Fourer, David M. Gay and Brian W. Kernighan, "A Modeling Language
alpar@9 5 # for Mathematical Programming." Management Science 36 (1990) 519-554.
alpar@9 6
alpar@9 7 ### SHIPPING SETS AND PARAMETERS ###
alpar@9 8
alpar@9 9 set whse 'warehouses'; # Locations from which demand is satisfied
alpar@9 10
alpar@9 11 set dctr 'distribution centers' within whse;
alpar@9 12
alpar@9 13 # Locations from which product may be shipped
alpar@9 14
alpar@9 15 param sc 'shipping cost' {dctr,whse} >= 0;
alpar@9 16
alpar@9 17 # Shipping costs, to whse from dctr, in $ / 100 lb
alpar@9 18
alpar@9 19 param huge 'largest shipping cost' > 0;
alpar@9 20
alpar@9 21 # Largest cost allowed for a usable shipping route
alpar@9 22
alpar@9 23 param msr 'minimum size restriction' {dctr,whse} logical;
alpar@9 24
alpar@9 25 # True indicates a minimum-size restriction on
alpar@9 26 # direct shipments using this dctr --> whse route
alpar@9 27
alpar@9 28 param dsr 'direct shipment requirement' {dctr} >= 0;
alpar@9 29
alpar@9 30 # Minimum total demand, in pallets, needed to
alpar@9 31 # allow shipment on routes subject to the
alpar@9 32 # minimum size restriction
alpar@9 33
alpar@9 34 ### PLANT SETS AND PARAMETERS ###
alpar@9 35
alpar@9 36 set fact 'factories' within dctr;
alpar@9 37
alpar@9 38 # Locations where product is manufactured
alpar@9 39
alpar@9 40 param rtmin 'regular-time total minimum' >= 0;
alpar@9 41
alpar@9 42 # Lower limit on (average) total regular-time
alpar@9 43 # crews employed at all factories
alpar@9 44
alpar@9 45 param rtmax 'regular-time total maximum' >= rtmin;
alpar@9 46
alpar@9 47 # Upper limit on (average) total regular-time
alpar@9 48 # crews employed at all factories
alpar@9 49
alpar@9 50 param otmin 'overtime total minimum' >= 0;
alpar@9 51
alpar@9 52 # Lower limit on total overtime hours at all factories
alpar@9 53
alpar@9 54 param otmax 'overtime total maximum' >= otmin;
alpar@9 55
alpar@9 56 # Upper limit on total overtime hours at all factories
alpar@9 57
alpar@9 58 param rmin 'regular-time minimums' {fact} >= 0;
alpar@9 59
alpar@9 60 # Lower limits on (average) regular-time crews
alpar@9 61
alpar@9 62 param rmax 'regular-time maximums' {f in fact} >= rmin[f];
alpar@9 63
alpar@9 64 # Upper limits on (average) regular-time crews
alpar@9 65
alpar@9 66 param omin 'overtime minimums' {fact} >= 0;
alpar@9 67
alpar@9 68 # Lower limits on overtime hours
alpar@9 69
alpar@9 70 param omax 'overtime maximums' {f in fact} >= omin[f];
alpar@9 71
alpar@9 72 # Upper limits on overtime hours
alpar@9 73
alpar@9 74 param hd 'hours per day' {fact} >= 0;
alpar@9 75
alpar@9 76 # Regular-time hours per working day
alpar@9 77
alpar@9 78 param dp 'days in period' {fact} > 0;
alpar@9 79
alpar@9 80 # Working days in the current planning period
alpar@9 81
alpar@9 82 ### PRODUCT SETS AND PARAMETERS ###
alpar@9 83
alpar@9 84 set prd 'products'; # Elements of the product group
alpar@9 85
alpar@9 86 param wt 'weight' {prd} > 0;
alpar@9 87
alpar@9 88 # Weight in 100 lb / 1000 cases
alpar@9 89
alpar@9 90 param cpp 'cases per pallet' {prd} > 0;
alpar@9 91
alpar@9 92 # Cases of product per shipping pallet
alpar@9 93
alpar@9 94 param tc 'transshipment cost' {prd} >= 0;
alpar@9 95
alpar@9 96 # Transshipment cost in $ / 1000 cases
alpar@9 97
alpar@9 98 param pt 'production time' {prd,fact} >= 0;
alpar@9 99
alpar@9 100 # Crew-hours to produce 1000 cases
alpar@9 101
alpar@9 102 param rpc 'regular-time production cost' {prd,fact} >= 0;
alpar@9 103
alpar@9 104 # Cost of production on regular time,
alpar@9 105 # in $ / 1000 cases
alpar@9 106
alpar@9 107 param opc 'overtime production cost' {prd,fact} >= 0;
alpar@9 108
alpar@9 109 # Cost of production on overtime, in $ / 1000 cases
alpar@9 110
alpar@9 111 ### DEMAND SETS AND PARAMETERS ###
alpar@9 112
alpar@9 113 param dt 'total demand' {prd} >= 0;
alpar@9 114
alpar@9 115 # Total demands for products, in 1000s
alpar@9 116
alpar@9 117 param ds 'demand shares' {prd,whse} >= 0.0, <= 1.0;
alpar@9 118
alpar@9 119 # Historical demand data, from which each
alpar@9 120 # warehouse's share of total demand is deduced
alpar@9 121
alpar@9 122 param dstot {p in prd} := sum {w in whse} ds[p,w];
alpar@9 123
alpar@9 124 # Total of demand shares; should be 1, but often isn't
alpar@9 125
alpar@9 126 param dem 'demand' {p in prd, w in whse} := dt[p] * ds[p,w] / dstot[p];
alpar@9 127
alpar@9 128 # Projected demands to be satisfied, in 1000s
alpar@9 129
alpar@9 130 set rt 'shipping routes available' :=
alpar@9 131
alpar@9 132 {d in dctr, w in whse:
alpar@9 133 d <> w and sc[d,w] < huge and
alpar@9 134 (w in dctr or sum {p in prd} dem[p,w] > 0) and
alpar@9 135 not (msr[d,w] and sum {p in prd} 1000*dem[p,w]/cpp[p] < dsr[d]) };
alpar@9 136
alpar@9 137 # List of ordered pairs that represent routes
alpar@9 138 # on which shipments are allowed
alpar@9 139
alpar@9 140 ### VARIABLES ###
alpar@9 141
alpar@9 142 var Rprd 'regular-time production' {prd,fact} >= 0;
alpar@9 143
alpar@9 144 # Regular-time production of each product
alpar@9 145 # at each factory, in 1000s of cases
alpar@9 146
alpar@9 147 var Oprd 'overtime production' {prd,fact} >= 0;
alpar@9 148
alpar@9 149 # Overtime production of each product
alpar@9 150 # at each factory, in 1000s of cases
alpar@9 151
alpar@9 152 var Ship 'shipments' {prd,rt} >= 0;
alpar@9 153
alpar@9 154 # Shipments of each product on each allowed route,
alpar@9 155 # in 1000s of cases
alpar@9 156
alpar@9 157 var Trans 'transshipments' {prd,dctr} >= 0;
alpar@9 158
alpar@9 159 # Transshipments of each product at each
alpar@9 160 # distribution center, in 1000s of cases
alpar@9 161
alpar@9 162 ### OBJECTIVE ###
alpar@9 163
alpar@9 164 minimize cost: sum {p in prd, f in fact} rpc[p,f] * Rprd[p,f] +
alpar@9 165 sum {p in prd, f in fact} opc[p,f] * Oprd[p,f] +
alpar@9 166 sum {p in prd, (d,w) in rt} sc[d,w] * wt[p] * Ship[p,d,w] +
alpar@9 167 sum {p in prd, d in dctr} tc[p] * Trans[p,d];
alpar@9 168
alpar@9 169 # Total cost: regular production, overtime
alpar@9 170 # production, shipping, and transshipment
alpar@9 171
alpar@9 172 ### CONSTRAINTS ###
alpar@9 173
alpar@9 174 rtlim 'regular-time total limits':
alpar@9 175
alpar@9 176 rtmin <= sum {p in prd, f in fact}
alpar@9 177 (pt[p,f] * Rprd[p,f]) / (dp[f] * hd[f]) <= rtmax;
alpar@9 178
alpar@9 179 # Total crews must lie between limits
alpar@9 180
alpar@9 181 otlim 'overtime total limits':
alpar@9 182
alpar@9 183 otmin <= sum {p in prd, f in fact} pt[p,f] * Oprd[p,f] <= otmax;
alpar@9 184
alpar@9 185 # Total overtime must lie between limits
alpar@9 186
alpar@9 187 rlim 'regular-time limits' {f in fact}:
alpar@9 188
alpar@9 189 rmin[f] <= sum {p in prd}
alpar@9 190 (pt[p,f] * Rprd[p,f]) / (dp[f] * hd[f]) <= rmax[f];
alpar@9 191
alpar@9 192 # Crews at each factory must lie between limits
alpar@9 193
alpar@9 194 olim 'overtime limits' {f in fact}:
alpar@9 195
alpar@9 196 omin[f] <= sum {p in prd} pt[p,f] * Oprd[p,f] <= omax[f];
alpar@9 197
alpar@9 198 # Overtime at each factory must lie between limits
alpar@9 199
alpar@9 200 noRprd 'no regular production' {p in prd, f in fact: rpc[p,f] = 0}:
alpar@9 201
alpar@9 202 Rprd[p,f] = 0;
alpar@9 203
alpar@9 204 noOprd 'no overtime production' {p in prd, f in fact: opc[p,f] = 0}:
alpar@9 205
alpar@9 206 Oprd[p,f] = 0; # Do not produce where specified cost is zero
alpar@9 207
alpar@9 208 bal 'material balance' {p in prd, w in whse}:
alpar@9 209
alpar@9 210 sum {(v,w) in rt}
alpar@9 211 Ship [p,v,w] + (if w in fact then Rprd[p,w] + Oprd[p,w]) =
alpar@9 212
alpar@9 213 dem[p,w] + (if w in dctr then sum {(w,v) in rt} Ship[p,w,v]);
alpar@9 214
alpar@9 215 # Demand is satisfied by shipment into warehouse
alpar@9 216 # plus production (if it is a factory)
alpar@9 217 # minus shipment out (if it is a distn. center)
alpar@9 218
alpar@9 219 trdef 'transshipment definition' {p in prd, d in dctr}:
alpar@9 220
alpar@9 221 Trans[p,d] >= sum {(d,w) in rt} Ship [p,d,w] -
alpar@9 222 (if d in fact then Rprd[p,d] + Oprd[p,d]);
alpar@9 223
alpar@9 224 # Transshipment at a distribution center is
alpar@9 225 # shipments out less production (if any)
alpar@9 226
alpar@9 227 ### DATA -- 3 PRODUCTS ###
alpar@9 228
alpar@9 229 data;
alpar@9 230
alpar@9 231 set prd := 18REG 24REG 24PRO ;
alpar@9 232
alpar@9 233 set whse := w01 w02 w03 w04 w05 w06 w08 w09 w12 w14 w15 w17
alpar@9 234 w18 w19 w20 w21 w24 w25 w26 w27 w28 w29 w30 w31
alpar@9 235 w32 w33 w34 w35 w36 w37 w38 w39 w40 w41 w42 w43
alpar@9 236 w44 w45 w46 w47 w48 w49 w50 w51 w53 w54 w55 w56
alpar@9 237 w57 w59 w60 w61 w62 w63 w64 w65 w66 w68 w69 w71
alpar@9 238 w72 w73 w74 w75 w76 w77 w78 w79 w80 w81 w82 w83
alpar@9 239 w84 w85 w86 w87 w89 w90 w91 w92 w93 w94 w95 w96
alpar@9 240 w98 x22 x23 ;
alpar@9 241
alpar@9 242 set dctr := w01 w02 w03 w04 w05 w62 w76 w96 ;
alpar@9 243
alpar@9 244 set fact := w01 w05 w96 ;
alpar@9 245
alpar@9 246 param huge := 99. ;
alpar@9 247
alpar@9 248 param rtmin := 0.0 ;
alpar@9 249 param rtmax := 8.0 ;
alpar@9 250
alpar@9 251 param otmin := 0.0 ;
alpar@9 252 param otmax := 96.0 ;
alpar@9 253
alpar@9 254 param rmin := w01 0.00 w05 0.00 w96 0.00 ;
alpar@9 255 param rmax := w01 3.00 w05 2.00 w96 3.00 ;
alpar@9 256
alpar@9 257 param omin := w01 0.0 w05 0.0 w96 0.0 ;
alpar@9 258 param omax := w01 48.0 w05 0.0 w96 48.0 ;
alpar@9 259
alpar@9 260 param hd := w01 8.0 w05 8.0 w96 8.0 ;
alpar@9 261
alpar@9 262 param dp := w01 19.0 w05 19.0 w96 19.0 ;
alpar@9 263
alpar@9 264 param wt := 18REG 47.3 24REG 63.0 24PRO 63.0 ;
alpar@9 265
alpar@9 266 param tc := 18REG 40.00 24REG 45.00 24PRO 45.00 ;
alpar@9 267
alpar@9 268 param dt := 18REG 376.0 24REG 172.4 24PRO 316.3 ;
alpar@9 269
alpar@9 270 param cpp := 18REG 102. 24REG 91. 24PRO 91. ;
alpar@9 271
alpar@9 272 param dsr := w01 96. w02 96. w03 96. w04 96. w05 96.
alpar@9 273 w62 96. w76 96. w96 96. ;
alpar@9 274
alpar@9 275 param pt (tr) :
alpar@9 276
alpar@9 277 18REG 24REG 24PRO :=
alpar@9 278
alpar@9 279 w01 1.194 1.429 1.429
alpar@9 280 w05 1.194 1.509 1.509
alpar@9 281 w96 0.000 1.600 1.600 ;
alpar@9 282
alpar@9 283 param rpc (tr) :
alpar@9 284
alpar@9 285 18REG 24REG 24PRO :=
alpar@9 286
alpar@9 287 w01 2119. 2653. 2617.
alpar@9 288 w05 2489. 3182. 3176.
alpar@9 289 w96 0. 2925. 2918. ;
alpar@9 290
alpar@9 291 param opc (tr) :
alpar@9 292
alpar@9 293 18REG 24REG 24PRO :=
alpar@9 294
alpar@9 295 w01 2903. 3585. 3579.
alpar@9 296 w05 0. 0. 0.
alpar@9 297 w96 0. 3629. 3622. ;
alpar@9 298
alpar@9 299 param sc default 99.99 (tr) :
alpar@9 300
alpar@9 301 w01 w02 w03 w04 w05 w62 w76 w96 :=
alpar@9 302
alpar@9 303 w01 . 2.97 1.14 2.08 2.37 1.26 2.42 1.43
alpar@9 304 w02 4.74 . 4.17 6.12 7.41 3.78 7.04 5.21
alpar@9 305 w03 2.45 4.74 . 3.67 2.84 0.90 2.41 2.55
alpar@9 306 w04 1.74 5.03 2.43 . 3.19 2.45 2.69 0.58
alpar@9 307 w05 2.70 5.16 2.84 2.85 . 3.26 3.34 2.71
alpar@9 308 w06 1.99 4.17 2.13 2.19 2.52 2.06 2.00 1.51
alpar@9 309 w08 0.21 2.92 1.24 2.07 2.29 1.25 2.32 1.55
alpar@9 310 w09 0.66 3.76 1.41 2.47 1.82 1.66 . 1.87
alpar@9 311 w12 1.38 3.83 1.68 2.53 2.39 . 1.96 1.94
alpar@9 312 w14 2.47 1.58 2.40 3.59 3.85 2.25 . 3.05
alpar@9 313 w15 1.06 4.95 2.48 1.39 3.41 1.96 . 1.02
alpar@9 314 w17 0.88 3.39 1.46 2.00 2.67 1.45 . 1.46
alpar@9 315 w18 7.90 6.57 7.79 9.59 10.81 . . 6.70
alpar@9 316 w19 1.42 4.12 1.96 1.99 3.52 1.88 . 1.26
alpar@9 317 w20 3.03 1.59 2.34 4.76 3.98 1.88 . 3.73
alpar@9 318 w24 1.58 2.80 2.27 2.87 3.19 1.31 . 2.05
alpar@9 319 w25 1.51 5.05 2.74 0.57 2.98 . 2.95 0.27
alpar@9 320 w26 1.75 3.61 2.70 1.54 4.07 3.52 . 1.03
alpar@9 321 w27 2.48 6.87 3.17 1.59 2.08 3.45 . 0.99
alpar@9 322 w28 2.05 6.83 2.97 1.13 2.91 . . 1.26
alpar@9 323 w29 4.03 3.68 4.46 3.20 5.50 . . 3.20
alpar@9 324 w30 2.48 5.78 2.99 2.24 1.79 3.10 . 1.39
alpar@9 325 w31 2.34 5.41 2.87 1.67 1.66 . . 1.39
alpar@9 326 w32 14.36 . . . . . . .
alpar@9 327 w33 3.87 4.27 5.11 3.48 5.66 4.03 . 3.05
alpar@9 328 w34 3.26 4.80 3.21 2.70 4.14 . . 1.77
alpar@9 329 w35 2.34 2.84 2.89 3.35 3.78 2.68 . 2.52
alpar@9 330 w36 2.43 5.69 2.96 2.95 1.02 2.61 1.07 2.54
alpar@9 331 w37 2.23 4.64 2.41 1.99 4.30 2.61 . 1.44
alpar@9 332 w38 4.66 4.36 5.23 3.04 4.46 . . 3.82
alpar@9 333 w39 1.11 3.51 1.10 2.53 3.07 1.12 . 2.23
alpar@9 334 w40 2.99 4.78 4.23 1.57 3.92 . . 1.80
alpar@9 335 w41 4.93 4.00 5.43 4.45 6.31 . . 3.81
alpar@9 336 w42 3.86 6.55 5.03 2.11 4.41 . . 2.63
alpar@9 337 w43 4.61 4.45 3.77 1.22 4.31 . . 2.35
alpar@9 338 w44 2.05 4.48 1.06 3.70 3.46 1.10 . 3.21
alpar@9 339 w45 0.92 3.42 1.58 3.04 1.82 1.94 . 2.52
alpar@9 340 w46 1.36 2.44 0.95 3.08 2.78 0.39 2.16 2.37
alpar@9 341 w47 1.30 3.39 1.60 2.49 4.29 2.04 . 1.68
alpar@9 342 w48 1.65 3.78 1.03 2.97 2.21 1.31 . 2.74
alpar@9 343 w49 1.96 3.00 1.50 3.24 3.68 1.00 . 2.99
alpar@9 344 w50 0.90 4.14 1.60 1.95 3.61 1.61 . 1.52
alpar@9 345 w51 1.59 3.95 0.25 2.96 2.58 1.00 2.41 2.71
alpar@9 346 w53 1.59 3.79 1.28 3.12 3.10 0.89 . 2.98
alpar@9 347 w54 1.72 4.36 1.61 2.92 2.34 1.91 1.97 3.05
alpar@9 348 w55 2.45 2.73 2.21 4.47 4.30 2.57 . 4.48
alpar@9 349 w56 1.10 3.73 1.59 2.74 2.33 1.45 . 2.44
alpar@9 350 w57 0.95 3.39 1.37 2.30 2.47 1.15 . 1.95
alpar@9 351 w59 3.29 5.35 3.32 3.81 1.52 3.38 1.34 4.08
alpar@9 352 w60 2.41 6.12 2.46 3.65 2.35 . 1.37 4.06
alpar@9 353 w61 3.32 5.50 3.41 3.38 1.23 . 0.99 4.28
alpar@9 354 w62 1.12 3.00 0.82 3.22 2.95 . 3.33 2.53
alpar@9 355 w63 3.59 6.36 3.25 4.12 1.84 3.59 1.46 4.03
alpar@9 356 w64 1.85 4.45 2.17 3.43 2.13 2.03 . 4.02
alpar@9 357 w65 2.78 4.79 2.81 2.94 1.54 2.90 1.07 2.94
alpar@9 358 w66 3.90 5.79 3.05 3.65 1.36 3.39 1.22 3.57
alpar@9 359 w68 2.61 5.20 2.90 2.34 1.68 3.19 1.48 2.31
alpar@9 360 w69 2.94 5.21 2.78 3.43 0.21 3.26 0.68 2.54
alpar@9 361 w71 2.06 4.98 2.38 2.44 1.59 2.97 1.05 2.55
alpar@9 362 w72 2.61 5.50 2.83 3.12 1.35 3.23 0.88 2.99
alpar@9 363 w73 8.52 6.16 8.03 8.83 10.44 7.38 10.26 .
alpar@9 364 w74 6.11 5.46 9.07 9.38 10.80 . . 8.25
alpar@9 365 w75 2.66 4.94 2.87 3.69 1.52 3.15 1.24 4.00
alpar@9 366 w76 1.99 5.26 2.23 3.36 0.58 3.17 . 2.50
alpar@9 367 w77 4.32 3.07 5.05 3.88 6.04 . . 4.15
alpar@9 368 w78 5.60 2.59 5.78 5.56 7.10 . . 5.60
alpar@9 369 w79 4.25 2.32 4.93 4.57 6.04 . . 4.58
alpar@9 370 w80 5.94 4.00 5.60 7.02 9.46 . . 7.51
alpar@9 371 w81 5.39 2.21 5.10 6.22 6.46 . . 6.58
alpar@9 372 w82 8.80 5.69 9.29 9.88 11.69 8.63 11.52 .
alpar@9 373 w83 4.40 . 5.24 5.21 5.81 3.91 7.04 5.33
alpar@9 374 w84 5.87 5.43 6.17 5.70 7.63 . . 5.70
alpar@9 375 w85 3.90 3.65 3.38 4.57 5.64 3.05 . 5.04
alpar@9 376 w86 5.48 2.10 5.70 6.37 7.33 . . 6.19
alpar@9 377 w87 8.88 5.54 9.50 9.71 11.64 8.85 11.68 .
alpar@9 378 w89 4.62 4.01 4.03 6.30 6.30 3.81 . 7.77
alpar@9 379 w90 4.35 2.72 4.61 4.01 5.60 . . 3.20
alpar@9 380 w91 7.61 4.42 7.83 6.85 8.79 . . 7.66
alpar@9 381 w92 7.15 2.69 6.91 7.20 . . . 7.06
alpar@9 382 w93 3.17 3.95 4.37 3.74 5.05 . . 2.40
alpar@9 383 w94 1.21 3.07 0.90 2.74 3.17 . 2.63 2.39
alpar@9 384 w95 5.82 3.29 6.55 7.06 11.47 . . 7.83
alpar@9 385 w96 1.77 5.20 2.72 0.59 3.47 2.48 . .
alpar@9 386 w98 3.04 1.92 3.64 3.70 4.90 3.05 . 3.88
alpar@9 387 x22 4.08 6.25 4.15 4.30 1.77 . 1.77 .
alpar@9 388 x23 3.39 5.74 3.55 4.08 1.69 . 1.47 . ;
alpar@9 389
alpar@9 390 param msr (tr) :
alpar@9 391
alpar@9 392 w01 w02 w03 w04 w05 w62 w76 w96 :=
alpar@9 393
alpar@9 394 w01 0 0 0 0 0 0 1 0
alpar@9 395 w02 0 0 0 0 0 0 1 0
alpar@9 396 w03 0 0 0 0 0 0 1 0
alpar@9 397 w04 0 0 0 0 0 0 1 0
alpar@9 398 w05 0 0 0 0 0 0 0 0
alpar@9 399 w06 0 1 1 1 1 1 1 1
alpar@9 400 w08 0 1 1 1 1 1 1 1
alpar@9 401 w09 0 1 1 1 1 1 0 1
alpar@9 402 w12 0 1 1 1 1 0 1 1
alpar@9 403 w14 1 1 1 1 1 0 0 1
alpar@9 404 w15 0 1 1 1 1 1 0 1
alpar@9 405 w17 0 1 1 1 1 1 0 1
alpar@9 406 w18 0 1 1 1 1 0 0 1
alpar@9 407 w19 0 1 1 1 1 0 0 1
alpar@9 408 w20 1 1 1 1 1 0 0 1
alpar@9 409 w24 0 1 1 1 1 0 0 1
alpar@9 410 w25 0 1 1 1 1 0 1 0
alpar@9 411 w26 1 1 1 0 1 1 0 1
alpar@9 412 w27 1 1 1 0 1 1 0 1
alpar@9 413 w28 1 1 1 0 1 0 0 1
alpar@9 414 w29 0 1 1 1 1 0 0 1
alpar@9 415 w30 1 1 1 0 1 1 0 1
alpar@9 416 w31 1 1 1 0 1 0 0 1
alpar@9 417 w32 0 0 0 0 0 0 0 0
alpar@9 418 w33 1 0 1 1 1 1 0 1
alpar@9 419 w34 1 1 1 0 1 0 0 1
alpar@9 420 w35 1 1 1 1 1 0 0 1
alpar@9 421 w36 0 1 1 1 0 1 1 1
alpar@9 422 w37 1 1 1 0 1 1 0 1
alpar@9 423 w38 1 1 1 0 1 0 0 1
alpar@9 424 w39 0 1 1 1 1 1 0 1
alpar@9 425 w40 1 1 1 0 1 0 0 1
alpar@9 426 w41 1 0 1 1 1 0 0 1
alpar@9 427 w42 1 1 1 0 1 0 0 1
alpar@9 428 w43 1 1 1 0 1 0 0 1
alpar@9 429 w44 1 1 1 1 1 0 0 1
alpar@9 430 w45 0 1 1 1 1 1 0 1
alpar@9 431 w46 0 1 1 1 1 0 1 1
alpar@9 432 w47 0 1 1 1 1 1 0 1
alpar@9 433 w48 0 1 1 1 1 0 0 1
alpar@9 434 w49 1 1 1 1 1 0 0 1
alpar@9 435 w50 0 1 1 1 1 1 0 1
alpar@9 436 w51 0 1 1 1 1 0 1 1
alpar@9 437 w53 1 1 1 1 1 0 0 1
alpar@9 438 w54 0 1 1 1 1 1 1 1
alpar@9 439 w55 0 1 1 1 1 0 0 1
alpar@9 440 w56 0 1 1 1 1 1 0 1
alpar@9 441 w57 0 1 1 1 1 1 0 1
alpar@9 442 w59 0 1 1 1 0 1 1 1
alpar@9 443 w60 0 1 1 1 1 0 1 1
alpar@9 444 w61 0 1 1 1 0 0 1 1
alpar@9 445 w62 0 0 0 0 0 0 1 0
alpar@9 446 w63 0 1 1 1 0 1 1 1
alpar@9 447 w64 0 1 1 1 1 1 0 1
alpar@9 448 w65 0 1 1 1 0 1 1 1
alpar@9 449 w66 0 1 1 1 0 1 1 1
alpar@9 450 w68 0 1 1 1 0 1 1 1
alpar@9 451 w69 0 1 1 1 0 1 1 1
alpar@9 452 w71 0 1 1 1 0 1 1 1
alpar@9 453 w72 0 1 1 1 0 1 1 1
alpar@9 454 w73 0 1 1 1 0 1 1 0
alpar@9 455 w74 0 1 1 1 0 0 0 1
alpar@9 456 w75 0 1 1 1 0 1 1 1
alpar@9 457 w76 0 0 0 0 0 0 0 0
alpar@9 458 w77 1 0 1 1 1 0 0 1
alpar@9 459 w78 1 0 1 1 1 0 0 1
alpar@9 460 w79 1 0 1 1 1 0 0 1
alpar@9 461 w80 1 0 1 1 1 0 0 1
alpar@9 462 w81 1 0 1 1 1 0 0 1
alpar@9 463 w82 1 0 1 1 1 1 1 0
alpar@9 464 w83 1 0 1 1 1 0 1 1
alpar@9 465 w84 1 0 1 1 1 0 0 1
alpar@9 466 w85 1 1 1 1 1 0 0 1
alpar@9 467 w86 1 0 1 1 1 0 0 1
alpar@9 468 w87 1 0 1 1 1 1 1 0
alpar@9 469 w89 1 0 1 1 1 1 0 1
alpar@9 470 w90 0 1 1 1 1 0 0 1
alpar@9 471 w91 1 0 1 1 1 0 0 1
alpar@9 472 w92 1 0 1 1 1 0 0 1
alpar@9 473 w93 1 1 1 0 1 0 0 1
alpar@9 474 w94 0 0 1 1 1 0 1 1
alpar@9 475 w95 1 0 1 1 1 0 0 1
alpar@9 476 w96 0 0 0 0 0 0 0 0
alpar@9 477 w98 1 0 1 1 1 1 0 1
alpar@9 478 x22 1 1 1 1 0 0 1 0
alpar@9 479 x23 1 1 1 1 0 0 1 0 ;
alpar@9 480
alpar@9 481 param ds default 0.000 (tr) :
alpar@9 482
alpar@9 483 18REG 24REG 24PRO :=
alpar@9 484
alpar@9 485 w01 0.000 0.000 0.008
alpar@9 486 w02 0.004 0.000 0.000
alpar@9 487 w03 0.000 0.000 0.000
alpar@9 488 w04 0.010 0.002 0.000
alpar@9 489 w05 0.000 0.000 0.000
alpar@9 490 w06 0.010 0.008 0.008
alpar@9 491 w08 0.030 0.024 0.024
alpar@9 492 w09 0.014 0.018 0.020
alpar@9 493 w12 0.014 0.012 0.010
alpar@9 494 w14 0.007 0.007 0.012
alpar@9 495 w15 0.010 0.019 0.018
alpar@9 496 w17 0.013 0.010 0.011
alpar@9 497 w19 0.015 0.012 0.009
alpar@9 498 w20 0.012 0.021 0.022
alpar@9 499 w21 0.000 0.000 0.000
alpar@9 500 w24 0.012 0.022 0.018
alpar@9 501 w25 0.019 0.025 0.020
alpar@9 502 w26 0.006 0.015 0.021
alpar@9 503 w27 0.008 0.010 0.015
alpar@9 504 w28 0.011 0.016 0.019
alpar@9 505 w29 0.008 0.020 0.013
alpar@9 506 w30 0.011 0.013 0.015
alpar@9 507 w31 0.011 0.013 0.017
alpar@9 508 w32 0.006 0.000 0.000
alpar@9 509 w33 0.000 0.015 0.014
alpar@9 510 w34 0.008 0.007 0.005
alpar@9 511 w35 0.002 0.006 0.014
alpar@9 512 w36 0.015 0.013 0.005
alpar@9 513 w37 0.017 0.016 0.015
alpar@9 514 w38 0.015 0.009 0.012
alpar@9 515 w39 0.007 0.017 0.022
alpar@9 516 w40 0.009 0.014 0.020
alpar@9 517 w41 0.003 0.014 0.011
alpar@9 518 w42 0.017 0.011 0.012
alpar@9 519 w43 0.009 0.013 0.011
alpar@9 520 w44 0.002 0.012 0.012
alpar@9 521 w45 0.016 0.025 0.028
alpar@9 522 w46 0.038 0.062 0.040
alpar@9 523 w47 0.007 0.010 0.010
alpar@9 524 w48 0.003 0.015 0.016
alpar@9 525 w49 0.005 0.016 0.017
alpar@9 526 w50 0.011 0.008 0.007
alpar@9 527 w51 0.010 0.022 0.021
alpar@9 528 w53 0.004 0.026 0.020
alpar@9 529 w54 0.020 0.017 0.025
alpar@9 530 w55 0.004 0.019 0.028
alpar@9 531 w56 0.004 0.010 0.008
alpar@9 532 w57 0.014 0.020 0.018
alpar@9 533 w59 0.012 0.006 0.007
alpar@9 534 w60 0.019 0.010 0.009
alpar@9 535 w61 0.028 0.010 0.012
alpar@9 536 w62 0.000 0.000 0.000
alpar@9 537 w63 0.070 0.027 0.037
alpar@9 538 w64 0.009 0.004 0.005
alpar@9 539 w65 0.022 0.015 0.016
alpar@9 540 w66 0.046 0.017 0.020
alpar@9 541 w68 0.005 0.012 0.016
alpar@9 542 w69 0.085 0.036 0.039
alpar@9 543 w71 0.011 0.013 0.010
alpar@9 544 w72 0.089 0.031 0.034
alpar@9 545 w75 0.026 0.012 0.010
alpar@9 546 w77 0.001 0.004 0.002
alpar@9 547 w78 0.002 0.004 0.002
alpar@9 548 w79 0.001 0.004 0.002
alpar@9 549 w80 0.001 0.001 0.002
alpar@9 550 w81 0.001 0.003 0.002
alpar@9 551 w83 0.009 0.010 0.008
alpar@9 552 w84 0.001 0.002 0.002
alpar@9 553 w85 0.001 0.004 0.005
alpar@9 554 w86 0.001 0.002 0.002
alpar@9 555 w87 0.002 0.003 0.000
alpar@9 556 w89 0.001 0.001 0.002
alpar@9 557 w90 0.006 0.017 0.013
alpar@9 558 w91 0.002 0.010 0.013
alpar@9 559 w92 0.000 0.003 0.002
alpar@9 560 w93 0.002 0.006 0.007
alpar@9 561 w95 0.001 0.007 0.007
alpar@9 562 w96 0.000 0.000 0.000
alpar@9 563 w98 0.006 0.005 0.002 ;
alpar@9 564
alpar@9 565 end;