1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
|
#papersize(20.5,61.0)#
#pagelength(59.0)#
#type("picac")#
NOTE 25.04.1985 14.30
note forstbetriebsmodell
note
note stehender Holzvorrat
l holz.k=holz.j+dt*(zuwachs.jk-ernte.jk)
n holz=startholz
c startholz=15000 Angabe in Festmetern(fm)
r zuwachs.kl=holz.k*vorratsfaktor*pflegefaktor.k
c vorratsfaktor=0.035 Zuwachs in % von holz
a pflegefaktor.k=tabhl(tpffaktor,pfzustand.k,0,1,0.1)
t tpffaktor=.7/.7/.7/.7/.7/.8/.9/.95/1/1.05/1.2
r ernte.kl=ernte1.k
a ernte1.k=holz.k*ernteproz.k*ernteprozfak.k
a ernteproz.k=tabhl(ternproz,ratio1.k,0.5,1.2,0.1)
t ternproz=0.02/0.025/0.03/0.03/0.03/0.035/0.06/0.08
a ratio1.k=holz.k/maxvorrat
c maxvorrat=16000
note
note resourcen
l resourcen.k=resourcen.j+dt*(einnahme.jk-eig.jk-pfausgaben.jk)
n resourcen=startresourcen
c startresourcen=100000 Geldeinheiten
c preis=190 Geldeinheiten pro fm
r einnahme.kl=ernte1.k*preis
r pfausgaben.kl=resourcen.k-eigenent.k
r eig.kl=eigenent.k
a anpassungsfaktor.k=tabhl(tanpass,ratio2.k,0.5,1.5,0.1)
t tanpass=.5/.55/.6/.7/.9/1/1/1/1.1/1.2/1.3
l eigenent.k=min(eigenent.j*anpass.jk,resourcen.j)
r anpass.kl=anpassungsfaktor.k
n eigenent=eigenentstart
c eigenentstart=20000
note
note arbeitseinheiten
note
l arbeit.k=arbeit.j+dt*(pfausgaben.jk/preisae-arbeitsverbrauch.jk)
n arbeit=startarbeit
c startarbeit=800
c preisae=100 ( preis pro arbeitseinheit in geldeinheiten )
r arbeitsverbrauch.kl=min(arbeit.k,notwarbeit.k)
a notwarbeit.k=tabhl(tnotwarb,pfzustand.k,0.0,1.0,0.1)
t tnotwarb=1600/1550/1500/1450/1400/1300/1150/950/800/700/600
a ratio2.k=arbeit.k/notwarbeit.k
a ernteprozfak.k=tabhl(ternteprozfak,ratio2.k,0.2,1.6,0.2)
t ternteprozfak=.4/.5/1/2/1.05/1/.9/.7
note
note Pflegezustand (pfzustand; dimensionslose Größe)
l pfzustand.k=pfzustand.j+dt*(pfverbess.jk-pfversch.jk)
n pfzustand=startpfzustand
c startpfzustand=0.8
r pfverbess.kl=smooth(pfx1.k,pfverzoeg)
a pfx1.k=pfverbfaktor.k*(1-pfzustand.k)
a pfverbfaktor.k=tabhl(tpfverb,ratio2.k,.8,1.4,.2)
t tpfverb=0/0/.1/0.2
r pfversch.kl=smooth(pfverschfaktor.k*pfzustand.k,pfverzoeg)
c pfverzoeg=2
a pfverschfaktor.k=tabhl(tpfversch,ratio2.k,0,.8,.2)
t tpfversch=.4/.2/.1/.05/0
note
note
note print ratio1,ratio2,eigenent,arbeit,pfzustand
plot holz=h(1e4,2e4)/eigenent=e(0,2e5)/pfzustand=P(0,1)/ratio2=2(0,5)
c dt=1
c length=50
note prtper=1
c pltper=1
run dyn.forst7
|
