{UTF-8}
********************************************************
	.daisyworld 1
********************************************************~
		Simulation: Peter Bützer
	|

********************************************************
	.Daisyworld 4
********************************************************~
		Simulation: Peter Bützer
	|

Einstrahlung= INTEG (
	SIN(Time*korr*0.52)/4*korr,
		0)
	~	Dmnl [0,?]
	~		|

korr=
	1*1
	~	1/Month [0,?]
	~	Faktor zur Korrektur der Einheiten bei der Sinus-Funktion
	|

vtw=
	kuw*Fläche unbewachsen*(0+Einstrahlung)
	~	m2/Month [0,?]
	~	Die Weissen wachsen bei grosser Einstrahlung besser
	|

vts=
	kus*Fläche unbewachsen*(1-Einstrahlung)
	~	m2/Month [0,?]
	~	Die Schwarzen wachsen bei kleiner Einstralung besser.
	|

Albedo=
	rw*Fläche weiss+ru*Fläche unbewachsen+rs*Fläche schwarz
	~	Dmnl [0,?]
	~		|

Fläche schwarz= INTEG (
	+vts-vst,
		0)
	~	m2 [0,?]
	~		|

Fläche unbewachsen= INTEG (
	vst+vwt-vts-vtw,
		1)
	~	m2 [0,?]
	~		|

Fläche weiss= INTEG (
	vtw-vwt,
		0)
	~	m2 [0,?]
	~		|

ksu=
	2.8
	~	1/Month [0,4]
	~	Absterben ca. 0.25 Monate -> k= 0.7/0.125 =2.8
	|

kus=
	1.4
	~	1/Month [0,1]
	~	Wachstum: ca. 1 Monat --> k=0.7/0.5=1.4
	|

kuw=
	1.4
	~	1/Month [0,2]
	~	Wachstum: ca. 1 Monat --> k=0.7/0.5=1.4
	|

kwu=
	2.8
	~	1/Month [0,4]
	~	Absterben ca. 0.25 Monate -> k= 0.7/0.125 =2.8
	|

rs=
	0.1
	~	1/m2 [0,?]
	~	Albedo für schwarze Gänseblümchen =0.1 (wie Wald)
	|

ru=
	0.3
	~	1/m2 [0,1]
	~	Reflexionsfaktor für mit Blumen unbewachsene Fläche (Wiese), Albedo=0.3 \
		(helleres grau)
	|

rw=
	0.8
	~	1/m2 [0,1]
	~	Reflexionsfaktor für weisse Fläche, Albedo=1 -> Gänseblümchen: 0.8
	|

vst=
	ksu*Fläche schwarz
	~	m2/Month [0,?]
	~		|

vwt=
	kwu*Fläche weiss
	~	m2/Month [0,?]
	~		|

********************************************************
	.Control
********************************************************~
		Simulation Control Parameters
	|

FINAL TIME  = 36
	~	Month
	~	The final time for the simulation.
	|

INITIAL TIME  = 0
	~	Month
	~	The initial time for the simulation.
	|

SAVEPER  = 
        TIME STEP 
	~	Month [0,?]
	~	The frequency with which output is stored.
	|

TIME STEP  = 0.1
	~	Month [0,?]
	~	The time step for the simulation.
	|

\\\---/// Sketch information - do not modify anything except names
V300  Do not put anything below this section - it will be ignored
*View 1
$192-192-192,0,Times New Roman|12||0-0-0|0-0-0|0-0-255|-1--1--1|-1--1--1|120,120,100,0
10,1,Fläche unbewachsen,298,193,59,30,3,131,0,4,0,0,0,0,0-0-0,192-192-192,|12||0-0-0
10,2,Fläche schwarz,502,192,47,28,3,131,0,6,0,0,0,0,0-0-0,0-0-0,|12||255-255-255
10,3,Fläche weiss,78,193,46,28,3,131,0,0,0,0,0,0
1,4,6,2,4,0,0,22,0,0,0,-1--1--1,,1|(427,166)|
1,5,6,1,100,0,0,22,0,0,0,-1--1--1,,1|(372,166)|
11,6,380,393,166,6,8,34,3,0,0,1,0,0,0
10,7,vts,393,186,14,12,40,3,0,0,-1,0,0,0
12,8,48,220,-15,10,8,0,3,0,0,-1,0,0,0
1,9,11,3,4,0,0,22,0,0,0,-1--1--1,,1|(150,171)|
1,10,11,1,100,0,0,22,0,0,0,-1--1--1,,1|(213,171)|
11,11,48,182,171,6,8,34,3,0,0,1,0,0,0
10,12,vtw,182,191,17,12,40,3,0,0,-1,0,0,0
1,13,15,1,4,0,0,22,0,0,0,-1--1--1,,1|(217,214)|
1,14,15,3,100,0,0,22,0,0,0,-1--1--1,,1|(153,214)|
11,15,252,189,214,6,8,34,3,0,0,1,0,0,0
10,16,vwt,189,234,17,12,40,3,0,0,-1,0,0,0
1,17,19,1,4,0,0,22,0,0,0,-1--1--1,,1|(374,214)|
1,18,19,2,100,0,0,22,0,0,0,-1--1--1,,1|(429,214)|
11,19,460,398,214,6,8,34,3,0,0,1,0,0,0
10,20,vst,398,234,14,12,40,3,0,0,-1,0,0,0
10,21,kuw,234,90,17,12,8,3,0,0,0,0,0,0
10,22,kwu,151,284,17,12,8,3,0,0,0,0,0,0
10,23,kus,336,89,13,12,8,3,0,0,0,0,0,0
10,24,ksu,450,283,13,12,8,3,0,0,0,0,0,0
1,25,1,11,1,0,0,0,0,64,0,-1--1--1,,1|(233,140)|
1,26,3,15,1,0,0,0,0,64,0,-1--1--1,,1|(143,241)|
1,27,22,15,1,0,0,0,0,64,0,-1--1--1,,1|(197,261)|
1,28,21,11,1,0,0,0,0,64,0,-1--1--1,,1|(187,110)|
1,29,1,6,1,0,0,0,0,64,0,-1--1--1,,1|(337,135)|
1,30,23,6,1,0,0,0,0,64,0,-1--1--1,,1|(379,106)|
1,31,2,19,1,0,0,0,0,64,0,-1--1--1,,1|(451,239)|
1,32,24,19,1,0,0,0,0,64,0,-1--1--1,,1|(402,259)|
10,33,Albedo,298,317,32,12,8,3,0,0,0,0,0,0
1,34,3,33,1,0,0,0,0,64,0,-1--1--1,,1|(93,286)|
1,35,2,33,1,0,0,0,0,64,0,-1--1--1,,1|(469,326)|
1,36,1,33,0,0,0,0,0,64,0,-1--1--1,,1|(298,257)|
10,37,rw,194,389,13,12,8,3,0,0,0,0,0,0
10,38,rs,411,391,10,12,8,3,0,0,0,0,0,0
10,39,ru,295,393,10,12,8,3,0,0,0,0,0,0
1,40,37,33,1,0,0,0,0,64,0,-1--1--1,,1|(255,379)|
1,41,38,33,1,0,0,0,0,64,0,-1--1--1,,1|(345,379)|
1,42,39,33,0,0,0,0,0,64,0,-1--1--1,,1|(295,361)|
12,43,0,891,252,332,241,3,188,0,0,1,0,0,0
Regelung
10,44,Einstrahlung,295,34,57,25,3,131,0,0,0,0,0,0
10,45,Time,102,101,34,12,8,2,0,3,-1,0,0,0,128-128-128,0-0-0,|12||128-128-128
1,46,45,44,1,0,0,0,0,64,0,-1--1--1,,1|(149,56)|
1,47,44,11,1,0,0,0,0,64,0,-1--1--1,,1|(175,77)|
1,48,44,6,1,0,0,0,0,64,0,-1--1--1,,1|(395,61)|
12,49,0,134,465,114,19,8,132,0,4,-1,0,0,0,0-0-0,0-192-255,|12||0-0-0
Gänseblümchen = Daisy
12,50,0,250,539,40,20,8,0,0,0,-1,0,0,0
10,51,korr,90,45,19,12,8,3,0,0,0,0,0,0
1,52,51,44,1,0,0,0,0,64,0,-1--1--1,,1|(170,26)|
///---\\\
:GRAPH Regelung
:TITLE Regelung
:X-AXIS Time
:X-LABEL Monate
:X-DIV 6
:X-MIN 0
:X-MAX 36
:SCALE
:VAR Einstrahlung
:Y-MIN 0
:SCALE
:VAR Albedo
:Y-MIN 0.2
:Y-MAX 0.5
:L<%^E!@
9:Current
23:0
15:0,0,0,0,0,0
19:100,0
27:2,
34:0,
5:Albedo
24:0
25:0
26:0