{UTF-8}
AA mothies= INTEG (
	birth 1-death 1,
		250)
	~	
	~		|

AA mothies freq=
	AA mothies/total moths
	~	
	~		|

Aa mothlets= INTEG (
	birth 2-death 2,
		500)
	~	
	~		|

Aa mothlets freq=
	Aa mothlets/total moths
	~	
	~		|

aa moths= INTEG (
	birth 3-death 3,
		250)
	~	
	~		|

aa moths freq=
	aa moths/total moths
	~	
	~		|

big A allele=
	((2*AA mothies)+Aa mothlets)/(2*total moths)
	~	
	~		|

bird pred light=
	pollution*bird pred rate
	~	
	~		|

bird pred rate=
	0.04
	~	 [0,0.1,0.01]
	~		|

bird predation dark=
	bird pred rate*(1-pollution)
	~	
	~		|

birth 1=
	repro rate*total moths*(big A allele^2)
	~	
	~		|

birth 2=
	repro rate*total moths*(2*big A allele*little a allele)
	~	
	~		|

birth 3=
	repro rate*total moths*(little a allele^2)
	~	
	~		|

dark freq=
	(AA mothies+Aa mothlets)/total moths
	~	
	~		|

death 1=
	AA mothies*(natural mortality+bird predation dark)
	~	
	~		|

death 2=
	Aa mothlets*(bird predation dark+natural mortality)
	~	
	~		|

death 3=
	aa moths*(bird pred light+natural mortality)
	~	
	~		|

light freq=
	aa moths/total moths
	~	
	~		|

little a allele=
	((2*aa moths)+Aa mothlets)/(2*total moths)
	~	
	~		|

natural mortality=
	0.045
	~	 [0,0.1,0.005]
	~		|

pollution=
	0
	~	 [0,1,1]
	~		|

repro rate=
	0.055
	~	 [0,0.1,0.005]
	~		|

total moths=
	AA mothies+Aa mothlets+aa moths
	~	
	~		|

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

FINAL TIME  = 200
	~	Year
	~	The final time for the simulation.
	|

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

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

TIME STEP  = 0.125
	~	Year [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,Comic Sans MS|12||0-0-0|0-0-0|0-0-255|-1--1--1|-1--1--1|120,120,100,0
10,1,AA mothies,489,101,78,34,3,131,0,0,0,0,0,0
10,2,Aa mothlets,486,315,79,33,3,131,0,0,0,0,0,0
10,3,aa moths,487,522,74,32,3,131,0,0,0,0,0,0
12,4,48,172,95,10,8,0,3,0,0,-1,0,0,0
1,6,8,1,4,0,0,22,0,0,0,-1--1--1,,1|(353,95)|
1,7,8,4,100,0,0,22,0,0,0,-1--1--1,,1|(233,95)|
11,8,0,290,95,6,8,34,3,0,0,1,0,0,0
10,9,birth 1,290,118,34,15,40,3,0,0,-1,0,0,0
12,10,48,782,97,10,8,0,3,0,0,-1,0,0,0
1,12,14,10,4,0,0,22,0,0,0,-1--1--1,,1|(722,97)|
1,13,14,1,100,0,0,22,0,0,0,-1--1--1,,1|(613,97)|
11,14,0,666,97,6,8,34,3,0,0,1,0,0,0
10,15,death 1,666,120,36,15,40,3,0,0,-1,0,0,0
12,16,48,184,319,10,8,0,3,0,0,-1,0,0,0
1,18,20,2,4,0,0,22,0,0,0,-1--1--1,,1|(356,319)|
1,19,20,16,100,0,0,22,0,0,0,-1--1--1,,1|(244,319)|
11,20,0,300,319,6,8,34,3,0,0,1,0,0,0
10,21,birth 2,300,342,35,15,40,3,0,0,-1,0,0,0
12,22,48,801,315,10,8,0,3,0,0,-1,0,0,0
1,24,26,22,4,0,0,22,0,0,0,-1--1--1,,1|(737,315)|
1,25,26,2,100,0,0,22,0,0,0,-1--1--1,,1|(618,315)|
11,26,0,678,315,6,8,34,3,0,0,1,0,0,0
10,27,death 2,678,338,38,15,40,3,0,0,-1,0,0,0
12,28,48,191,517,10,8,0,3,0,0,-1,0,0,0
1,30,32,3,4,0,0,22,0,0,0,-1--1--1,,1|(363,517)|
1,31,32,28,100,0,0,22,0,0,0,-1--1--1,,1|(251,517)|
11,32,0,307,517,6,8,34,3,0,0,1,0,0,0
10,33,birth 3,307,540,35,15,40,3,0,0,-1,0,0,0
12,34,48,802,522,10,8,0,3,0,0,-1,0,0,0
1,36,38,34,4,0,0,22,0,0,0,-1--1--1,,1|(737,522)|
1,37,38,3,100,0,0,22,0,0,0,-1--1--1,,1|(615,522)|
11,38,0,676,522,6,8,34,3,0,0,1,0,0,0
10,39,death 3,676,545,38,15,40,3,0,0,-1,0,0,0
10,40,big A allele,347,200,53,15,8,3,0,0,0,0,0,0
10,41,natural mortality,578,193,44,28,8,3,0,0,0,0,0,0
10,42,bird predation dark,800,224,71,28,8,3,0,0,0,0,0,0
10,43,pollution,904,127,41,15,8,3,0,0,0,0,0,0
10,44,repro rate,109,250,51,15,8,3,0,0,0,0,0,0
10,45,little a allele,333,419,58,15,8,3,0,0,0,0,0,0
10,46,bird pred rate,828,425,71,14,8,3,0,0,0,0,0,0
10,47,bird pred light,734,642,73,14,8,3,0,0,0,0,0,0
10,48,total moths,475,734,54,15,8,3,0,0,0,0,0,0
10,49,aa moths,744,865,51,15,8,2,0,3,-1,0,0,0,128-128-128,0-0-0,|0||128-128-128
10,50,Aa mothlets,419,846,67,14,8,2,0,3,-1,0,0,0,128-128-128,0-0-0,|0||128-128-128
10,51,AA mothies,161,848,65,14,8,2,0,3,-1,0,0,0,128-128-128,0-0-0,|0||128-128-128
10,52,AA mothies freq,199,745,57,28,8,3,0,0,0,0,0,0
10,53,aa moths freq,780,743,68,14,8,3,0,0,0,0,0,0
10,54,Aa mothlets freq,463,945,58,28,8,3,0,0,0,0,0,0
10,55,dark freq,187,959,43,43,2,131,0,1,0,3,0,0,0-0-255,0-0-0,|0||0-0-0
10,56,light freq,645,960,45,45,2,131,0,1,0,3,0,0,255-0-0,0-0-0,|0||0-0-0
1,57,51,48,0,0,0,0,0,64,0,-1--1--1,,1|(309,793)|
1,58,50,48,0,0,0,0,0,64,0,-1--1--1,,1|(443,796)|
1,59,49,48,0,0,0,0,0,64,0,-1--1--1,,1|(615,802)|
1,60,51,52,1,0,0,0,0,64,0,-1--1--1,,1|(147,802)|
1,61,48,52,1,0,0,0,0,64,0,-1--1--1,,1|(344,712)|
1,62,49,53,1,0,0,0,0,64,0,-1--1--1,,1|(797,811)|
1,63,48,53,1,0,0,0,0,64,0,-1--1--1,,1|(609,714)|
1,64,50,54,1,0,0,0,0,64,0,-1--1--1,,1|(416,887)|
1,65,48,54,1,0,0,0,0,64,0,-1--1--1,,1|(527,840)|
1,66,49,56,1,0,0,0,0,64,0,-1--1--1,,1|(714,910)|
1,67,48,56,1,0,0,0,0,64,0,-1--1--1,,1|(574,824)|
1,68,50,55,1,0,0,0,0,64,0,-1--1--1,,1|(319,917)|
1,69,51,55,1,0,0,0,0,64,0,-1--1--1,,1|(148,883)|
1,70,48,55,1,0,0,0,0,64,0,-1--1--1,,1|(333,823)|
10,71,pollution,885,591,49,15,8,2,0,3,-1,0,0,0,128-128-128,0-0-0,|0||128-128-128
10,72,natural mortality,558,614,48,28,8,2,0,3,-1,0,0,0,128-128-128,0-0-0,|0||128-128-128
10,73,total moths,148,607,64,14,8,2,0,3,-1,0,0,0,128-128-128,0-0-0,|0||128-128-128
10,74,natural mortality,606,409,48,28,8,2,0,3,-1,0,0,0,128-128-128,0-0-0,|0||128-128-128
10,75,total moths,161,398,64,14,8,2,0,3,-1,0,0,0,128-128-128,0-0-0,|0||128-128-128
10,76,total moths,188,207,64,14,8,2,0,3,-1,0,0,0,128-128-128,0-0-0,|0||128-128-128
1,77,1,40,1,0,0,0,0,64,0,-1--1--1,,1|(419,176)|
1,78,40,8,1,0,0,0,0,64,0,-1--1--1,,1|(348,159)|
1,79,76,9,1,0,0,0,0,64,0,-1--1--1,,1|(223,155)|
1,81,44,9,1,0,0,0,0,64,0,-1--1--1,,1|(127,177)|
1,82,76,40,1,0,0,0,0,64,0,-1--1--1,,1|(273,242)|
1,83,40,21,1,0,0,0,0,64,0,-1--1--1,,1|(343,275)|
1,84,44,20,1,0,0,0,0,64,0,-1--1--1,,1|(172,322)|
1,85,75,21,1,0,0,0,0,64,0,-1--1--1,,1|(217,355)|
1,86,45,21,1,0,0,0,0,64,0,-1--1--1,,1|(334,378)|
1,87,2,45,1,0,0,0,0,64,0,-1--1--1,,1|(417,400)|
1,88,75,45,1,0,0,0,0,64,0,-1--1--1,,1|(241,438)|
1,89,45,33,1,0,0,0,0,64,0,-1--1--1,,1|(342,478)|
1,90,3,45,1,0,0,0,0,64,0,-1--1--1,,1|(413,454)|
1,91,44,33,1,0,0,0,0,64,0,-1--1--1,,1|(79,444)|
1,92,73,33,1,0,0,0,0,64,0,-1--1--1,,1|(231,595)|
1,93,1,15,1,0,0,0,0,64,0,-1--1--1,,1|(592,143)|
1,94,41,15,1,0,0,0,0,64,0,-1--1--1,,1|(644,180)|
1,95,42,15,1,0,0,0,0,64,0,-1--1--1,,1|(720,184)|
1,96,43,42,1,0,0,0,0,64,0,-1--1--1,,1|(864,164)|
1,97,46,42,1,0,0,0,0,64,0,-1--1--1,,1|(845,335)|
1,98,46,47,1,0,0,0,0,64,0,-1--1--1,,1|(823,543)|
1,99,71,47,1,0,0,0,0,64,0,-1--1--1,,1|(835,635)|
1,100,47,39,1,0,0,0,0,64,0,-1--1--1,,1|(698,601)|
1,101,72,39,1,0,0,0,0,64,0,-1--1--1,,1|(636,597)|
1,102,3,39,1,0,0,0,0,64,0,-1--1--1,,1|(589,563)|
1,103,74,27,1,0,0,0,0,64,0,-1--1--1,,1|(670,405)|
1,104,42,27,1,0,0,0,0,64,0,-1--1--1,,1|(747,293)|
1,105,2,27,1,0,0,0,0,64,0,-1--1--1,,1|(595,358)|
1,106,2,40,1,0,0,0,0,64,0,-1--1--1,,1|(420,246)|
12,107,2949662,1363,271,308,207,3,188,0,0,1,0,0,0
Phenotype_Frequency
///---\\\
:GRAPH Phenotype_Frequency
:TITLE Phenotype Frequency
:SCALE
:VAR dark freq
:Y-MIN 0
:Y-MAX 1
:LINE-WIDTH 3
:VAR light freq
:Y-MIN 0
:Y-MAX 1
:LINE-WIDTH 3
:VAR pollution
:LINE-WIDTH 3
:L<%^E!@
1:Current.vdf
9:Current
22:$,Dollar,Dollars,$s
22:Hour,Hours
22:Month,Months
22:Person,People,Persons
22:Unit,Units
22:Week,Weeks
22:Year,Years
22:Day,Days
23:0
15:0,0,0,0,0,0
19:100,0
27:2,
34:0,
4:Time
5:repro rate
35:Date
36:YYYY-MM-DD
37:2000
38:1
39:1
40:0
41:0
24:0
25:200
26:200