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

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

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

korr=
	1*1
	~	1/Month [0,?]
	~	Korrektur, damit das Argument in der Sinusfunktion keine Einheit hat.
	|

Albedo total=
	(Albedo weiss+Albedo unbewachsen+Albedo schwarz)
	~	m*m [0,?]
	~		|

Einheitkorr=
	50
	~	W [40,60]
	~	Korrektur, damit die Temperatur für die verwendete Fläche stimmt.
	|

Tu=
	SQRT(SQRT((2+Einstrahlung)*Einheitkorr/(SB Konstante*Albedo total)))-273
	~	K [0,?]
	~		|

Albedo schwarz=
	rs*Fläche schwarz
	~	m*m [0,?]
	~		|

Albedo unbewachsen=
	ru*Fläche unbewachsen
	~	m*m [0,?]
	~		|

Albedo weiss=
	rw*Fläche weiss
	~	m*m [0,?]
	~		|

SB Konstante=
	5.67e-008
	~	W/( K*K*K*K*m*m)
	~		|

vtw=
	kuw*Fläche unbewachsen*(0.2+Einstrahlung)
	~	m*m/Month [0,?]
	~	Je grösser die Einstrahlung, desto rascher das Wachstum
	|

vts=
	kus*Fläche unbewachsen*(1.2-Einstrahlung)
	~	m*m/Month [0,?]
	~	Je kleiner die Einstrahlung, desto grösser das Wachstum
	|

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

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

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

ksu=
	3.5
	~	1/Month [0,4]
	~	Absterben ca. 0.2 Monate -> k= 0.7/0.2 =3.5, die Schwarzen sterben weniger \
		rasch als die Weissen
	|

kus=
	1.4
	~	1/Month [0,1]
	~	Wachstum: ca. 2 Monate --> k=1.4
	|

kuw=
	0.7
	~	1/Month [0,2]
	~	Wachstum: ca. 1 Monate --> k=0.7, die weissen wachsen rascher als die \
		schwarzen (mehr Licht)
	|

kwu=
	7
	~	1/Month [0,4]
	~	Absterben ca. 0.1 Monate -> k= 0.7/0.1 =7, die Weissen sterben rascher, \
		als die Schwarzen
	|

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

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

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

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

vwt=
	kwu*Fläche weiss
	~	m*m/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.01
	~	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,310,246,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,514,245,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,90,246,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|(442,219)|
1,5,6,1,100,0,0,22,0,0,0,-1--1--1,,1|(387,219)|
11,6,508,411,219,6,8,34,3,0,0,1,0,0,0
10,7,vts,411,239,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|(162,224)|
1,10,11,1,100,0,0,22,0,0,0,-1--1--1,,1|(225,224)|
11,11,48,194,224,6,8,34,3,0,0,1,0,0,0
10,12,vtw,194,244,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|(229,267)|
1,14,15,3,100,0,0,22,0,0,0,-1--1--1,,1|(165,267)|
11,15,524,201,267,6,8,34,3,0,0,1,0,0,0
10,16,vwt,201,287,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|(386,267)|
1,18,19,2,100,0,0,22,0,0,0,-1--1--1,,1|(441,267)|
11,19,604,410,267,6,8,34,3,0,0,1,0,0,0
10,20,vst,410,287,14,12,40,3,0,0,-1,0,0,0
10,21,kuw,246,143,17,12,8,3,0,0,0,0,0,0
10,22,kwu,163,337,17,12,8,3,0,0,0,0,0,0
10,23,kus,348,142,13,12,8,3,0,0,0,0,0,0
10,24,ksu,462,336,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|(245,193)|
1,26,3,15,1,0,0,0,0,64,0,-1--1--1,,1|(155,294)|
1,27,22,15,1,0,0,0,0,64,0,-1--1--1,,1|(209,314)|
1,28,21,11,1,0,0,0,0,64,0,-1--1--1,,1|(199,163)|
1,29,1,6,1,0,0,0,0,64,0,-1--1--1,,1|(354,177)|
1,30,23,6,1,0,0,0,0,64,0,-1--1--1,,1|(397,155)|
1,31,2,19,1,0,0,0,0,64,0,-1--1--1,,1|(463,292)|
1,32,24,19,1,0,0,0,0,64,0,-1--1--1,,1|(414,312)|
10,33,Albedo total,314,448,39,20,8,131,0,0,0,0,0,0
10,34,rw,60,295,13,12,8,3,0,0,0,0,0,0
10,35,rs,550,291,10,12,8,3,0,0,0,0,0,0
10,36,ru,258,299,10,12,8,3,0,0,0,0,0,0
10,37,Einstrahlung,307,87,57,25,3,131,0,0,0,0,0,0
10,38,Time,119,164,34,12,8,2,0,3,-1,0,0,0,128-128-128,0-0-0,|12||128-128-128
1,39,37,11,1,0,0,0,0,64,0,-1--1--1,,1|(196,135)|
1,40,37,6,1,0,0,0,0,64,0,-1--1--1,,1|(408,122)|
12,41,0,262,592,40,20,8,0,0,0,-1,0,0,0
10,42,Tu,313,534,13,12,8,3,0,0,0,0,0,0
1,43,37,42,1,0,0,0,0,64,0,-1--1--1,,1|(220,348)|
10,44,SB Konstante,463,517,56,12,8,3,0,0,0,0,0,0
1,45,44,42,1,0,0,0,0,64,0,-1--1--1,,1|(365,553)|
10,46,Albedo weiss,91,405,41,24,8,131,0,0,0,0,0,0
10,47,Albedo schwarz,516,404,35,22,8,3,0,0,0,0,0,0
10,48,Albedo unbewachsen,312,359,56,22,8,3,0,0,0,0,0,0
1,49,3,46,1,0,0,0,0,64,0,-1--1--1,,1|(121,344)|
1,50,34,46,1,0,0,0,0,64,0,-1--1--1,,1|(66,350)|
1,51,1,48,1,0,0,0,0,64,0,-1--1--1,,1|(316,299)|
1,52,36,48,1,0,0,0,0,64,0,-1--1--1,,1|(288,317)|
1,53,35,47,1,0,0,0,0,64,0,-1--1--1,,1|(555,362)|
1,54,2,47,1,0,0,0,0,64,0,-1--1--1,,1|(497,327)|
1,55,46,33,1,0,0,0,0,64,0,-1--1--1,,1|(215,444)|
1,56,48,33,0,0,0,0,0,64,0,-1--1--1,,1|(312,397)|
1,57,47,33,1,0,0,0,0,64,0,-1--1--1,,1|(423,447)|
1,58,33,42,1,0,0,0,0,64,0,-1--1--1,,1|(315,487)|
12,59,0,888,596,293,190,3,188,0,0,1,0,0,0
Albedo
12,60,0,887,208,293,186,3,188,0,0,1,0,0,0
Temperatur
10,61,Einheitkorr,114,495,47,12,8,3,0,0,0,0,0,0
1,62,61,42,1,0,0,0,0,64,0,-1--1--1,,1|(221,544)|
10,63,korr,104,121,19,12,8,3,0,0,0,0,0,0
1,64,63,37,1,0,0,0,0,64,0,-1--1--1,,1|(148,87)|
1,65,38,37,1,0,0,0,0,64,0,-1--1--1,,1|(181,97)|
///---\\\
:GRAPH Temperatur
:TITLE Einstrahlung-Temperatur
:X-AXIS Time
:X-LABEL Monate
:X-DIV 6
:X-MIN 0
:X-MAX 36
:SCALE
:VAR Einstrahlung
:Y-MIN 0
:Y-MAX 1
:SCALE
:VAR Tu
:UNITS °C
:Y-MIN 10
:Y-MAX 30

:GRAPH Albedo
:TITLE Albedo
:X-AXIS Time
:SCALE
:VAR Albedo weiss
:Y-MIN 0
:VAR Albedo unbewachsen
:Y-MIN 0
:VAR Albedo schwarz
:Y-MIN 0
:SCALE
:VAR Einstrahlung
:L<%^E!@
9:Current
23:0
15:0,0,0,0,0,0
19:100,0
27:2,
34:0,
5:Albedo schwarz
24:0
25:0
26:0