source: coopr.pysp/stable/2.1/examples/pysp/forestry/chile18-nb-yr/ForestChile2.dat @ 2068

Last change on this file since 2068 was 2068, checked in by wehart, 10 years ago

Merged revisions 1952-2067 via svnmerge from
https://software.sandia.gov/svn/public/coopr/coopr.pysp/trunk

........

r1956 | jwatson | 2009-12-02 17:56:53 -0700 (Wed, 02 Dec 2009) | 3 lines


Added --scenario-solver-options and --ef-solver-options options to the "runph" script.

........

r1957 | dlwoodr | 2009-12-03 14:17:35 -0700 (Thu, 03 Dec 2009) | 2 lines


Documentation updates for pysp

........

r1974 | wehart | 2009-12-06 17:20:56 -0700 (Sun, 06 Dec 2009) | 2 lines


Updating PyPI categories

........

r1978 | jwatson | 2009-12-10 21:29:33 -0700 (Thu, 10 Dec 2009) | 3 lines


Eliminated exception-handling logic when loading user-defined extension modules in PH. The range of exceptions is too large, and for debugging purposes, it is more useful to see the raw trace output.

........

r1979 | jwatson | 2009-12-10 22:23:17 -0700 (Thu, 10 Dec 2009) | 5 lines


Biggest enhancement: The efwriter command-line script now has the option to output a CVaR-weighted objective term. Not extensively tested, but behaves sane on a number of small test cases.


Improved exception handling and error diagnostics in both the runph and efwriter scripts.

........

r1985 | jwatson | 2009-12-12 10:45:17 -0700 (Sat, 12 Dec 2009) | 3 lines


Modified PH to only use warm-starts if a solver has the capability!

........

r1998 | jwatson | 2009-12-13 15:17:58 -0700 (Sun, 13 Dec 2009) | 3 lines


Changed references to _component to active_component.

........

r2026 | wehart | 2009-12-21 23:27:06 -0700 (Mon, 21 Dec 2009) | 2 lines


Attempting to update PH. I'm not sure if this works, since I don't know how to test PH.

........

r2029 | jwatson | 2009-12-22 09:52:21 -0700 (Tue, 22 Dec 2009) | 3 lines


Some fixes to the ef writer based on Bill's recent changes to _initialize_constraint.

........

r2035 | jwatson | 2009-12-22 21:10:32 -0700 (Tue, 22 Dec 2009) | 3 lines


Added --scenario-mipgap option to PH script. Added _mipgap attribute to PH object. This attribute is mirrored to the solver plugin at the initiation of each iteration, after any PH extensions have the opportunity to provide a new value to the attribute. It is currently made use of by the WW PH extension.

........

r2037 | dlwoodr | 2009-12-23 14:38:43 -0700 (Wed, 23 Dec 2009) | 2 lines


add this example from Fernando

........

r2038 | dlwoodr | 2009-12-23 14:46:56 -0700 (Wed, 23 Dec 2009) | 3 lines


finish the job: we can now duplicate Fernando's example

........

r2039 | jwatson | 2009-12-23 15:13:54 -0700 (Wed, 23 Dec 2009) | 3 lines


Missed fix with new constraint initialization syntax in PH linearization.

........

r2058 | jwatson | 2009-12-29 10:57:58 -0700 (Tue, 29 Dec 2009) | 3 lines


Missed some _initialize_constraint function calls within the PySP EF writer during the recent switch to the corresponding "add" method.

........

r2059 | jwatson | 2009-12-29 10:58:34 -0700 (Tue, 29 Dec 2009) | 3 lines


Enabling garbage collection by default in PH.

........

r2060 | jwatson | 2009-12-29 10:59:05 -0700 (Tue, 29 Dec 2009) | 3 lines


Elimnating some debug output.

........

r2061 | jwatson | 2009-12-29 11:07:47 -0700 (Tue, 29 Dec 2009) | 3 lines


Fixing some option documentation in PH.

........

r2062 | jwatson | 2009-12-29 12:00:37 -0700 (Tue, 29 Dec 2009) | 3 lines


Added ef-mipgap option to PH scripts.

........

  • Property svn:executable set to *
File size: 8.7 KB
Line 
1# fbv
2# mnr May 3, 2009
3
4# general nodes
5set Nodes := E1
6            I1 I2 I3
7            C01 C02 C03 C04 C05 C06 C07 C08 C09 ;
8
9# following the notation in the paper (note: nodes \equiv origins)
10
11#
12# Sets
13#
14
15set Times := Ano1 Ano2 Ano3 Ano4;
16
17set OriginNodes := C01 C02 C03 C04 C05 C06 C07 C08 C09 ;
18
19set IntersectionNodes := I1 I2 I3;
20
21set ExitNodes := E1;
22
23set HarvestCells := U1 U2 U3 U4 U5 U6 U7 U8 U9 U10 U11 U12 U13
24                        U14 U15 U16 U17 U18 U19 U20 U21 U22 U23 U24 U25;
25
26set HCellsForOrigin[C01] := U1 U2 U4 ;
27set HCellsForOrigin[C02] := U6 U7 U9 U10 ;
28set HCellsForOrigin[C03] := U12 U13 ;
29set HCellsForOrigin[C04] := U11 U15 U16 ;
30set HCellsForOrigin[C05] := U14 U23 U24 ;
31set HCellsForOrigin[C06] := U17 U18 U19 ;
32set HCellsForOrigin[C07] := U22 U25 ;
33set HCellsForOrigin[C08] := U20 U21 ;
34set HCellsForOrigin[C09] := U3 U5 U8 ;
35
36set ExistingRoads :=
37C02 C03
38C03 E1
39C05 I1
40I1 I2
41I2 I3
42I3 E1
43;
44
45set PotentialRoads :=
46C01 C09
47C01 C02
48C02 C09
49C09 C02
50C09 E1
51C09 C03
52C04 C03
53C04 I3
54C06 I3
55C08 C06
56C08 I2
57C07 C08
58C07 I1
59C07 C05
60;
61
62set AllRoads :=
63C02 C03
64C03 E1
65C05 I1
66I1 I2
67I2 I3
68I3 E1
69C01 C09
70C01 C02
71C02 C09
72C09 C02
73C09 E1
74C09 C03
75C04 C03
76C04 I3
77C06 I3
78C08 C06
79C08 I2
80C07 C08
81C07 I1
82C07 C05
83;
84
85#
86# Deterministic Parameters
87#
88
89# productivity of cell h if it is harvested in period t
90param a :=
91U1      Ano1    362
92U1      Ano2    362
93U1      Ano3    362
94U1      Ano4    362
95U2      Ano1    362
96U2      Ano2    362
97U2      Ano3    362
98U2      Ano4    362
99U3      Ano1    362
100U3      Ano2    362
101U3      Ano3    362
102U3      Ano4    362
103U4      Ano1    642
104U4      Ano2    642
105U4      Ano3    642
106U4      Ano4    642
107U5      Ano1    362
108U5      Ano2    362
109U5      Ano3    362
110U5      Ano4    362
111U6      Ano1    642
112U6      Ano2    642
113U6      Ano3    642
114U6      Ano4    642
115U7      Ano1    642
116U7      Ano2    642
117U7      Ano3    642
118U7      Ano4    642
119U8      Ano1    362
120U8      Ano2    362
121U8      Ano3    362
122U8      Ano4    362
123U9      Ano1    642
124U9      Ano2    642
125U9      Ano3    642
126U9      Ano4    642
127U10     Ano1    642
128U10     Ano2    642
129U10     Ano3    642
130U10     Ano4    642
131U11     Ano1    600
132U11     Ano2    600
133U11     Ano3    600
134U11     Ano4    600
135U12     Ano1    518
136U12     Ano2    518
137U12     Ano3    518
138U12     Ano4    518
139U13     Ano1    518
140U13     Ano2    518
141U13     Ano3    518
142U13     Ano4    518
143U14     Ano1    567
144U14     Ano2    567
145U14     Ano3    567
146U14     Ano4    567
147U15     Ano1    600
148U15     Ano2    600
149U15     Ano3    600
150U15     Ano4    600
151U16     Ano1    600
152U16     Ano2    600
153U16     Ano3    600
154U16     Ano4    600
155U17     Ano1    518
156U17     Ano2    518
157U17     Ano3    518
158U17     Ano4    518
159U18     Ano1    401
160U18     Ano2    401
161U18     Ano3    401
162U18     Ano4    401
163U19     Ano1    401
164U19     Ano2    401
165U19     Ano3    401
166U19     Ano4    401
167U20     Ano1    401
168U20     Ano2    401
169U20     Ano3    401
170U20     Ano4    401
171U21     Ano1    373
172U21     Ano2    373
173U21     Ano3    373
174U21     Ano4    373
175U22     Ano1    373
176U22     Ano2    373
177U22     Ano3    373
178U22     Ano4    373
179U23     Ano1    567
180U23     Ano2    567
181U23     Ano3    567
182U23     Ano4    567
183U24     Ano1    567
184U24     Ano2    567
185U24     Ano3    567
186U24     Ano4    567
187U25     Ano1    373
188U25     Ano2    373
189U25     Ano3    373
190U25     Ano4    373
191;
192
193# Area of cell h to be harvested
194param A :=
195U1      10.1
196U2      10.1
197U3      10.1
198U4      13.4
199U5      10.3
200U6      10.2
201U7      10.1
202U8      10.4
203U9      16.7
204U10     12.8
205U11     12.2
206U12     15.2
207U13     11.3
208U14     12.7
209U15     12.6
210U16     12.6
211U17     10
212U18     14.4
213U19     10.8
214U20     10.4
215U21     16.4
216U22     10
217U23     11.4
218U24     14.8
219U25     10.1
220;
221
222# harvesting cost of one hectare of cell h in time t
223param P :=
224U1      Ano1    8
225U1      Ano2    8
226U1      Ano3    8
227U1      Ano4    8
228U2      Ano1    8
229U2      Ano2    8
230U2      Ano3    8
231U2      Ano4    8
232U3      Ano1    8
233U3      Ano2    8
234U3      Ano3    8
235U3      Ano4    8
236U4      Ano1    8
237U4      Ano2    8
238U4      Ano3    8
239U4      Ano4    8
240U5      Ano1    8
241U5      Ano2    8
242U5      Ano3    8
243U5      Ano4    8
244U6      Ano1    8
245U6      Ano2    8
246U6      Ano3    8
247U6      Ano4    8
248U7      Ano1    8
249U7      Ano2    8
250U7      Ano3    8
251U7      Ano4    8
252U8      Ano1    8
253U8      Ano2    8
254U8      Ano3    8
255U8      Ano4    8
256U9      Ano1    8
257U9      Ano2    8
258U9      Ano3    8
259U9      Ano4    8
260U10     Ano1    8
261U10     Ano2    8
262U10     Ano3    8
263U10     Ano4    8
264U11     Ano1    8
265U11     Ano2    8
266U11     Ano3    8
267U11     Ano4    8
268U12     Ano1    8
269U12     Ano2    8
270U12     Ano3    8
271U12     Ano4    8
272U13     Ano1    8
273U13     Ano2    8
274U13     Ano3    8
275U13     Ano4    8
276U14     Ano1    8
277U14     Ano2    8
278U14     Ano3    8
279U14     Ano4    8
280U15     Ano1    8
281U15     Ano2    8
282U15     Ano3    8
283U15     Ano4    8
284U16     Ano1    8
285U16     Ano2    8
286U16     Ano3    8
287U16     Ano4    8
288U17     Ano1    8
289U17     Ano2    8
290U17     Ano3    8
291U17     Ano4    8
292U18     Ano1    8
293U18     Ano2    8
294U18     Ano3    8
295U18     Ano4    8
296U19     Ano1    8
297U19     Ano2    8
298U19     Ano3    8
299U19     Ano4    8
300U20     Ano1    8
301U20     Ano2    8
302U20     Ano3    8
303U20     Ano4    8
304U21     Ano1    8
305U21     Ano2    8
306U21     Ano3    8
307U21     Ano4    8
308U22     Ano1    8
309U22     Ano2    8
310U22     Ano3    8
311U22     Ano4    8
312U23     Ano1    8
313U23     Ano2    8
314U23     Ano3    8
315U23     Ano4    8
316U24     Ano1    8
317U24     Ano2    8
318U24     Ano3    8
319U24     Ano4    8
320U25     Ano1    8
321U25     Ano2    8
322U25     Ano3    8
323U25     Ano4    8
324;
325
326
327# prod cost
328param Q :=
329C01     Ano1    0.1
330C01     Ano2    0.1
331C01     Ano3    0.1
332C01     Ano4    0.1
333C02     Ano1    0.1
334C02     Ano2    0.1
335C02     Ano3    0.1
336C02     Ano4    0.1
337C03     Ano1    0.1
338C03     Ano2    0.1
339C03     Ano3    0.1
340C03     Ano4    0.1
341C04     Ano1    0.1
342C04     Ano2    0.1
343C04     Ano3    0.1
344C04     Ano4    0.1
345C05     Ano1    0.1
346C05     Ano2    0.1
347C05     Ano3    0.1
348C05     Ano4    0.1
349C06     Ano1    0.1
350C06     Ano2    0.1
351C06     Ano3    0.1
352C06     Ano4    0.1
353C07     Ano1    0.1
354C07     Ano2    0.1
355C07     Ano3    0.1
356C07     Ano4    0.1
357C08     Ano1    0.1
358C08     Ano2    0.1
359C08     Ano3    0.1
360C08     Ano4    0.1
361C09     Ano1    0.1
362C09     Ano2    0.1
363C09     Ano3    0.1
364C09     Ano4    0.1
365;
366
367# construction cost of one road in arc (k,l) in time t
368param C :=
369C01     C09     Ano1    1440
370C01     C09     Ano2    1309.09
371C01     C09     Ano3    1190.08
372C01     C09     Ano4    1081.89
373C01     C02     Ano1    1950
374C01     C02     Ano2    1772.73
375C01     C02     Ano3    1611.57
376C01     C02     Ano4    1465.06
377C02     C09     Ano1    1690
378C02     C09     Ano2    1536.36
379C02     C09     Ano3    1396.69
380C02     C09     Ano4    1269.72
381C09     C02     Ano1    1690
382C09     C02     Ano2    1536.36
383C09     C02     Ano3    1396.69
384C09     C02     Ano4    1269.72
385C09     E1      Ano1    2840
386C09     E1      Ano2    2581.82
387C09     E1      Ano3    2347.11
388C09     E1      Ano4    2133.73
389C09     C03     Ano1    2190
390C09     C03     Ano2    1990.91
391C09     C03     Ano3    1809.92
392C09     C03     Ano4    1645.38
393C04     C03     Ano1    1170
394C04     C03     Ano2    1063.64
395C04     C03     Ano3    966.94
396C04     C03     Ano4    879.04
397C04     I3      Ano1    1320
398C04     I3      Ano2    1200
399C04     I3      Ano3    1090.91
400C04     I3      Ano4    991.74
401C06     I3      Ano1    150
402C06     I3      Ano2    136.36
403C06     I3      Ano3    123.97
404C06     I3      Ano4    112.7
405C08     C06     Ano1    1950
406C08     C06     Ano2    1772.73
407C08     C06     Ano3    1611.57
408C08     C06     Ano4    1465.06
409C08     I2      Ano1    1110
410C08     I2      Ano2    1009.09
411C08     I2      Ano3    917.36
412C08     I2      Ano4    833.96
413C07     C08     Ano1    930
414C07     C08     Ano2    845.45
415C07     C08     Ano3    768.6
416C07     C08     Ano4    698.72
417C07     I1      Ano1    1950
418C07     I1      Ano2    1772.73
419C07     I1      Ano3    1611.57
420C07     I1      Ano4    1465.06
421C07     C05     Ano1    2130
422C07     C05     Ano2    1936.36
423C07     C05     Ano3    1760.33
424C07     C05     Ano4    1600.3
425;
426
427# unit transport cost through arc (k,l) in time t
428param D :=
429C01     C09     Ano1    1.9
430C01     C02     Ano1    3.25
431C02     C09     Ano1    2.35
432C09     C02     Ano1    2.35
433C09     E1      Ano1    4.1
434C09     C03     Ano1    3.65
435C04     C03     Ano1    1.95
436C04     I3      Ano1    2.2
437C06     I3      Ano1    0.25
438C08     C06     Ano1    3.25
439C08     I2      Ano1    1.85
440C07     C08     Ano1    1.55
441C07     I1      Ano1    3.15
442C07     C05     Ano1    3.25
443C02     C03     Ano1    3.6
444C03     E1      Ano1    1.4
445C05     I1      Ano1    0.95
446I1      I2      Ano1    1.1
447I2      I3      Ano1    1.45
448I3      E1      Ano1    1.95
449C01     C09     Ano2    1.9
450C01     C02     Ano2    3.25
451C02     C09     Ano2    2.35
452C09     C02     Ano2    2.35
453C09     E1      Ano2    4.1
454C09     C03     Ano2    3.65
455C04     C03     Ano2    1.95
456C04     I3      Ano2    2.2
457C06     I3      Ano2    0.25
458C08     C06     Ano2    3.25
459C08     I2      Ano2    1.85
460C07     C08     Ano2    1.55
461C07     I1      Ano2    3.15
462C07     C05     Ano2    3.25
463C02     C03     Ano2    3.6
464C03     E1      Ano2    1.4
465C05     I1      Ano2    0.95
466I1      I2      Ano2    1.1
467I2      I3      Ano2    1.45
468I3      E1      Ano2    1.95
469C01     C09     Ano3    1.9
470C01     C02     Ano3    3.25
471C02     C09     Ano3    2.35
472C09     C02     Ano3    2.35
473C09     E1      Ano3    4.1
474C09     C03     Ano3    3.65
475C04     C03     Ano3    1.95
476C04     I3      Ano3    2.2
477C06     I3      Ano3    0.25
478C08     C06     Ano3    3.25
479C08     I2      Ano3    1.85
480C07     C08     Ano3    1.55
481C07     I1      Ano3    3.15
482C07     C05     Ano3    3.25
483C02     C03     Ano3    3.6
484C03     E1      Ano3    1.4
485C05     I1      Ano3    0.95
486I1      I2      Ano3    1.1
487I2      I3      Ano3    1.45
488I3      E1      Ano3    1.95
489C01     C09     Ano4    1.9
490C01     C02     Ano4    3.25
491C02     C09     Ano4    2.35
492C09     C02     Ano4    2.35
493C09     E1      Ano4    4.1
494C09     C03     Ano4    3.65
495C04     C03     Ano4    1.95
496C04     I3      Ano4    2.2
497C06     I3      Ano4    0.25
498C08     C06     Ano4    3.25
499C08     I2      Ano4    1.85
500C07     C08     Ano4    1.55
501C07     I1      Ano4    3.15
502C07     C05     Ano4    3.25
503C02     C03     Ano4    3.6
504C03     E1      Ano4    1.4
505C05     I1      Ano4    0.95
506I1      I2      Ano4    1.1
507I2      I3      Ano4    1.45
508I3      E1      Ano4    1.95
509;
510                       
511#
512# Stochastic Parameters
513#
514
515# sale price at exit e in time period t
516param R := E1 Ano1 45
517           E1 Ano2 60
518           E1 Ano3 65
519           E1 Ano4 57 ;
520
521# upper and lower bounds on wood supplied
522param Zlb := Ano1 30000
523             Ano2 27000
524             Ano3 28000
525             Ano4 20000 ;
526param Zub := Ano1 40000
527             Ano2 50000
528             Ano3 52000
529             Ano4 51000 ;
530param yr:=
531Ano1    1
532Ano2    1
533Ano3    1
534Ano4    1
535;
536
537set COriginNodeForCell[U1]:= C01;
538set COriginNodeForCell[U2]:= C01;
539set COriginNodeForCell[U3]:= C09;
540set COriginNodeForCell[U4]:= C01;
541set COriginNodeForCell[U5]:= C09;
542set COriginNodeForCell[U6]:= C02;
543set COriginNodeForCell[U7]:= C02;
544set COriginNodeForCell[U8]:= C09;
545set COriginNodeForCell[U9]:= C02;
546set COriginNodeForCell[U10]:= C02;
547set COriginNodeForCell[U11]:= C04;
548set COriginNodeForCell[U12]:= C03;
549set COriginNodeForCell[U13]:= C03;
550set COriginNodeForCell[U14]:= C05;
551set COriginNodeForCell[U15]:= C04;
552set COriginNodeForCell[U16]:= C04;
553set COriginNodeForCell[U17]:= C06;
554set COriginNodeForCell[U18]:= C06;
555set COriginNodeForCell[U19]:= C06;
556set COriginNodeForCell[U20]:= C08;
557set COriginNodeForCell[U21]:= C08;
558set COriginNodeForCell[U22]:= C07;
559set COriginNodeForCell[U23]:= C05;
560set COriginNodeForCell[U24]:= C05;
561set COriginNodeForCell[U25]:= C07;
562
Note: See TracBrowser for help on using the repository browser.