Consumo=
	Tasa de agotamiento de OD*OD hipolimnético
	~	ppm
	~		|

Descomposición aerobia en función de OD(
	[(0,-0.001)-(10,0.1)],(0,0),(2,0.004),(3,0.007),(6,0.025),(8,0.05),(10,0.1))
	~	{/año}
	~	No he encontrado datos que relacionen la tasa de descomposición aerobia \
		con el OD. Los valores están inspirados en el modelo de eutrofización que \
		ehe encontrado en http://www.wsu.edu/ forda/lake1.html
	|

Mortalidad=
	Mortalidad en función de la biomasa(Biomasa de fitoplancton)
	~	{/año} [0,1,0.001]
	~	Este valor viene dado por la gráfica de la Mortalidad en función de la \
		biomasa
	|

Tasa de descomposición aerobia=
	Descomposición aerobia en función de OD(OD hipolimnético)
	~	mg cla/m3·año [0,1,0.0001]
	~	La tasa de descomposición aerobia estará en función de la cantidad de \
		oxígeno disponible
	|

Tasa de descomposición anaerobia=
	0.01-Tasa de descomposición aerobia
	~	mg cla/m3·año [0,1,0.0001]
	~	Asumiendo una tasa de descomposición constante, la tasa de descomposición \
		anaerobia va a aumentar lo que se reduzca la tasa de descomposición aerobia
	|

Mortalidad en función de la biomasa(
	[(0,0)-(80,1)],(0,0),(25.69,0.06),(52.6,0.2125),(75.841,0.7719))
	~	{/año}
	~	Al aumentar la densidad del fitoplancton el crecimiento se frena hasta \
		alcanzar una densidad máxima
	|

Aportes=
	Tasa producción de oxígeno
	~	ppm
	~		|

Tasa de agotamiento de OD=
	Agotamiento de OD en función de la biomasa muerta(Biomasa muerta)
	~	ppm/año
	~		|

Descomposición=
	(Tasa de descomposición aerobia+Tasa de descomposición anaerobia)*Biomasa muerta
	~	mg cla/m3
	~		|

Muertes=
	Biomasa de fitoplancton*Mortalidad
	~	mg cla/m3
	~	Representa el descenso de biomasa en mg cla/m3
	|

Biomasa muerta= INTEG (
	+Muertes-Descomposición,
		0)
	~	mg cla/m3
	~		|

Biomasa de fitoplancton= INTEG (
	+Reproducción-Muertes,
		4.2)
	~	mg cla/m3
	~	Empleo la concentración de clorofila A como indicador de la biomasa de fitoplancton.
		El valor inicial (4.2) corresponde al valor de un lago oligotrófico según \
		la clasificacióndel estado trófico de la OCDE
	|

Producción de oxígeno en función de la turbidez(
	[(0,0)-(20,6)],(0.122324,0.0307),(2.20183,0.789474),(3.73089,1.36842),(5.62691,1.97368\
		),(7.95107,2.63158),(10.0306,3.18421),(12.3547,3.76316),(14.8624,4.145),(17.0031,4.39\
		),(19.9388,4.39))
	~	ppm/año
	~	Estos datos habría que obtenerlos experimentalmente
	|

Tasa producción de oxígeno=
	Producción de oxígeno en función de la turbidez(Profundidad de Secchi)
	~	ppm/año
	~	Suponiendo que el único aporte de O2 al hipolimnion procede de la fotosíntesis.
		La tasa fotosintética está condicionada por la penetración de la luz y \
		esta por la turbidez que depende a su vez de la densidad del fitoplancton
	|

Agotamiento de OD en función de la biomasa muerta(
	[(0,0)-(10,1)],(0.1835,0.3145),(2.08,0.3465),(4.22,0.443),(5.535,0.5439),(6.85,0.7018\
		),(8.318,0.9912))
	~	ppm/año
	~	Estos datos habría que obtenerlos experimentalmente
	|

Profundidad de Secchi=
	Profundidad de Secchi en función de la biomasa de fitoplancton(Biomasa de fitoplancton\
		)
	~	Profundidad de Secchi
	~	Se considera que sólo depende de la densidad del fitoplancton
	|

Profundidad de Secchi en función de la biomasa de fitoplancton(
	[(0,0)-(200,20)],(0.611621,19.7368),(4.2,9.9),(7.33945,7.2807),(16.1,4.2),(28.1346,2.98246\
		),(42.6,2.45),(72.7829,1.66667),(111.927,1.14035),(199.388,0.175439))
	~	
	~	Gráfica construida a partir de los datos  de la clasificación preliminar \
		del estado trófico de los lagos de la OCDE
	|

Aportes de P=
	Aportes naturales+P procedente de contaminación
	~	mg/L·año
	~		|

Aportes naturales=
	0.00807
	~	mg/L·año [0,0.005,1e-005]
	~	Ligeramente superior a la Tasa de pérdida global de P para simular el \
		proceso de eutrofización natural
	|

OD hipolimnético
	= INTEG (Aportes-Consumo,10)
	~	mg/L [1,?]
	~	El valor inicial corresponde al de un lago oligotrófico
	~	:SUPPLEMENTARY 
	|

P procedente de contaminación=
	0
	~	mg/L [0,0.1,0.0001]
	~	Modificando este valor se simula la transformación del lago oligotrófico a \
		mesotrófico, eutrófico e incluso hipereutrófico
	|

Reproducción=
	Biomasa de fitoplancton*Tasa de división
	~	mg cla/m3
	~	Indica el incremento de biomasa de fitoplancton expresado en mg cla/m3
	|

Fósforo total= INTEG (
	Aportes de P-Pérdida neta de P,
		0.008)
	~	mg/L
	~	El valor inicial (0,008 ) corresponde al de un lago oligotrófico según la \
		clasificación de la OCDE (Organización de Cooperación y Desarrollo \
		Económicos)
	|

Pérdida neta de P=
	Tasa de pérdida global de P
	~	mg/L·año
	~		|

Tasa de pérdida global de P=
	0.008
	~	mg/L·año
	~		|

Tasa de división en función de PT(
	[(0,0)-(1,1)],(0,0.001),(0.0172,0.02),(0.0345,0.05),(0.0517,0.1),(0.069,0.135),(0.0862\
		,0.205),(0.0948012,0.232456),(0.103,0.26),(0.121,0.405),(0.138,0.495),(0.155,0.58),\
		(0.172,0.65),(0.19,0.7),(0.207,0.89),(0.224,0.755),(0.276,0.755),(0.5,0.755),(1,0.755\
		))
	~	{/año}
	~	Los datos de esta tabla los he extraido de un modelo de eutrofización que \
		he encontrado en la página web: Environmental Inquiry \
		(http://ei.cornell.edu/)
	|

Tasa de división=
	Tasa de división en función de PT(Fósforo total)
	~	{/año}
	~	Este valor viene dado por la gráfica de la Tasa de división en función de \
		PT
	|

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

FINAL TIME  = 10
	~	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.25
	~	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
*Gráfica
$192-192-192,0,Arial|12||0-0-0|0-0-0|0-0-255|-1--1--1|-1--1--1|96,96,100
12,1,0,309,283,288,262,3,188,0,0,1,0,0,0
Eutrofización_cultural
12,2,0,760,65,80,20,3,124,0,0,0,0,0,0
P procedente de contaminación,0,0.05,1e-006
12,3,0,777,179,141,36,8,132,0,0,-1,0,0,0
INSTRUCCIONES: Selecciona la herramienta "Automatically simulate on change" y utiliza el control deslizante para simular  el aporte extra de P
12,4,0,781,290,173,45,8,132,0,4,-1,0,0,0,-1--1--1,255-255-128,|12||0-0-0
Al aumentar la tasa de entrada de P al lago se aprecia como la biomasa del fitoplancton aumenta hasta alcanzar un máximo  y  el OD del hipolimnion  se reduce hasta agotarse en la situación de hipereutrofía  (PT>100 mg/m3)
12,5,0,249,587,227,27,8,132,0,4,-1,0,0,0,0-0-0,255-128-192,|0||0-0-0
Con las flechas que aparecen en la esquina inferior izquierda se puede cambiar de esta vista a la que muestra el modelo
\\\---/// Sketch information - do not modify anything except names
V300  Do not put anything below this section - it will be ignored
*Modelo
$192-192-192,0,Arial|12||64-64-128|0-128-0|0-128-0|-1--1--1|-1--1--1|96,96,5
10,1,Fósforo total,306,170,49,21,3,131,0,3,0,0,0,0,0-0-128,0-0-0,|12||0-0-128
10,2,Biomasa de fitoplancton,469,365,49,21,3,131,0,3,0,0,0,0,0-128-0,0-0-0,|12||0-128-0
10,3,OD hipolimnético,680,644,49,20,3,131,0,3,0,0,0,0,128-64-0,0-0-0,|12||128-64-0
12,4,48,62,166,10,8,0,3,0,0,-1,0,0,0
1,5,7,1,4,0,0,22,1,0,0,0-0-128,|12||0-0-0,1|(213,166)|
1,6,7,4,100,0,0,22,1,0,0,0-0-128,|12||0-0-0,1|(115,166)|
11,7,48,164,166,6,8,34,3,0,1,1,0,0,0,0-0-128,0-0-0,|12||0-0-0
10,8,Aportes de P,164,184,48,10,40,3,0,1,-1,0,0,0,0-0-128,0-0-0,|12||0-0-0
12,9,48,563,170,10,8,0,3,0,0,-1,0,0,0
1,10,12,9,4,0,0,22,1,0,0,0-0-128,|12||0-0-0,1|(506,170)|
1,11,12,1,100,0,0,22,1,0,0,0-0-128,|12||0-0-0,1|(401,170)|
11,12,48,454,170,6,8,34,3,0,1,1,0,0,0,0-64-128,0-0-0,|12||0-0-0
10,13,Pérdida neta de P,454,187,70,9,40,131,0,1,-1,0,0,0,0-64-128,0-0-0,|12||0-0-0
12,14,48,224,359,10,8,0,3,0,3,-1,0,0,0,0-128-0,0-0-0,|12||0-128-0
1,15,17,2,4,0,0,22,3,0,0,0-128-0,|12||64-64-128,1|(376,359)|
1,16,17,14,100,0,0,22,3,0,0,0-128-0,|12||64-64-128,1|(277,359)|
11,17,48,327,359,6,8,34,3,0,1,1,0,0,0,0-128-0,0-0-0,|12||0-0-0
10,18,Reproducción,327,377,51,10,40,3,0,3,-1,0,0,0,0-128-0,0-0-0,|12||0-128-0
12,19,48,436,644,10,8,0,3,0,3,-1,0,0,0,128-64-0,0-0-0,|12||128-64-0
1,20,22,3,4,0,0,22,1,0,0,128-64-0,|12||0-0-0,1|(587,644)|
1,21,22,19,100,0,0,22,1,0,0,128-64-0,|12||0-0-0,1|(489,644)|
11,22,48,538,644,6,8,34,3,0,1,1,0,0,0,128-64-0,0-0-0,|12||0-0-0
10,23,Aportes,538,662,34,10,40,3,0,3,-1,0,0,0,128-64-0,0-0-0,|12||128-64-0
12,24,48,903,645,10,8,0,3,0,3,-1,0,0,0,128-64-0,0-0-0,|12||128-64-0
1,25,27,24,4,0,0,22,1,0,0,128-64-0,|12||0-0-0,1|(855,645)|
1,26,27,3,100,0,0,22,1,0,0,128-64-0,|12||0-0-0,1|(767,645)|
11,27,48,811,645,6,8,34,3,0,1,1,0,0,0,128-64-0,0-0-0,|12||0-0-0
10,28,Consumo,811,664,37,11,40,131,0,3,-1,0,0,0,128-64-0,0-0-0,|12||128-64-0
10,29,Profundidad de Secchi,433,465,18,18,2,131,0,2,4,0,0,0,0-128-0,0-0-0,|12||0-128-0
10,30,Tasa de pérdida global de P,481,114,61,18,8,3,0,0,0,0,0,0
10,31,Aportes naturales,72,101,34,18,8,3,0,0,0,0,0,0
10,32,P procedente de contaminación,267,99,61,18,8,3,0,0,0,0,0,0
1,33,31,7,1,0,0,0,1,64,0,0-0-128,|12||0-0-0,1|(138,129)|
1,34,32,7,1,0,0,0,1,64,0,0-0-128,|12||0-0-0,1|(193,123)|
1,35,30,12,0,0,0,0,1,64,0,0-0-128,|12||0-0-0,1|(467,141)|
10,36,Tasa de división,240,287,20,20,2,3,0,3,4,0,0,0,0-128-0,0-0-0,|12||0-128-0
10,37,Mortalidad,588,264,39,10,8,3,0,3,0,0,0,0,0-128-0,0-0-0,|12||0-128-0
1,38,36,17,1,0,0,0,0,64,0,-1--1--1,,1|(266,327)|
10,39,Tasa producción de oxígeno,426,573,61,18,8,3,0,3,0,0,0,0,128-64-0,0-0-0,|12||128-64-0
1,40,29,39,1,0,0,0,0,64,0,-1--1--1,,1|(408,522)|
1,41,39,22,1,0,0,0,1,64,0,128-64-0,|12||0-0-0,1|(497,594)|
10,42,Tasa de agotamiento de OD,837,550,74,18,8,131,0,3,0,0,0,0,128-64-0,0-0-0,|12||128-64-0
1,43,42,27,1,0,0,0,1,64,0,128-64-0,|12||0-0-0,1|(837,596)|
1,44,1,36,1,0,0,0,1,64,0,0-0-128,|12||0-0-0,1|(242,255)|
10,45,Tasa de división en función de PT,85,287,71,18,8,131,0,2,0,0,0,0,0-0-128,0-0-0,|12||255-0-0
1,46,45,36,0,0,0,0,1,64,0,255-0-0,|12||0-0-0,1|(181,287)|
1,47,2,17,1,0,0,0,0,64,0,-1--1--1,,1|(384,325)|
10,48,Profundidad de Secchi en función de la biomasa de fitoplancton,229,465,95,27,8,131,0,2,0,0,0,0,0-0-0,0-0-0,|12||255-0-0
1,49,48,29,1,0,0,0,1,64,0,255-0-0,|12||0-0-0,1|(362,465)|
1,50,2,29,1,0,0,0,0,64,0,-1--1--1,,1|(469,417)|
10,51,Producción de oxígeno en función de la turbidez,215,573,84,27,8,3,0,2,0,0,0,0,0-0-0,0-0-0,|12||255-0-0
1,52,51,39,0,0,0,0,1,64,0,255-0-0,|12||0-0-0,1|(325,573)|
10,53,Agotamiento de OD en función de la biomasa muerta,627,549,84,27,8,3,0,2,0,0,0,0,0-0-0,0-0-0,|12||255-0-0
1,54,53,42,0,0,0,0,1,64,0,255-0-0,|12||0-0-0,1|(730,549)|
10,55,Biomasa muerta,729,364,40,20,3,3,0,2,0,0,0,0,0-128-0,0-0-0,|12||0-128-0
1,56,58,55,4,0,0,22,3,0,0,0-128-0,|12||64-64-128,1|(649,363)|
1,57,58,2,100,0,0,22,0,0,0,-1--1--1,,1|(557,363)|
11,58,412,603,363,6,8,34,3,0,1,1,0,0,0,0-128-0,0-0-0,|12||0-0-0
10,59,Muertes,603,381,30,10,40,3,0,3,-1,0,0,0,0-128-0,0-0-0,|12||0-128-0
12,60,48,936,365,10,8,0,3,0,3,-1,0,0,0,0-128-0,0-0-0,|12||0-128-0
1,61,63,60,4,0,0,22,0,0,0,-1--1--1,,1|(886,364)|
1,62,63,55,100,0,0,22,0,0,0,-1--1--1,,1|(801,364)|
11,63,48,840,364,6,8,34,3,0,1,1,0,0,0,0-128-0,0-0-0,|12||0-0-0
10,64,Descomposición,840,380,62,8,40,131,0,3,-1,0,0,0,0-128-0,0-0-0,|12||0-128-0
1,65,37,58,1,0,0,0,0,64,0,-1--1--1,,1|(628,302)|
10,66,Tasa de descomposición aerobia,943,466,61,27,8,3,0,3,0,0,0,0,0-128-0,0-0-0,|12||128-0-0
10,67,Tasa de descomposición anaerobia,790,467,61,27,8,3,0,2,0,0,0,0,0-128-0,0-0-0,|12||128-0-0
1,68,66,64,1,0,0,0,1,64,0,128-0-0,|12||0-0-0,1|(920,411)|
1,69,67,64,1,0,0,0,1,64,0,128-0-0,|12||0-0-0,1|(792,405)|
1,70,2,58,1,0,0,0,0,64,0,-1--1--1,,1|(540,326)|
1,71,55,42,1,0,0,0,0,64,0,-1--1--1,,1|(707,471)|
1,72,55,63,1,0,0,0,0,64,0,-1--1--1,,1|(789,330)|
1,73,2,37,1,0,0,0,0,64,0,-1--1--1,,1|(487,284)|
10,74,Mortalidad en función de la biomasa,744,264,50,27,8,3,0,2,0,0,0,0,0-0-0,0-0-0,|12||255-0-0
1,75,74,37,0,0,0,0,1,64,0,255-0-0,|12||0-0-0,1|(667,264)|
1,76,3,66,1,0,0,0,1,64,0,128-0-0,|12||0-0-0,1|(909,672)|
1,77,66,67,0,0,0,0,1,64,0,128-0-0,|12||0-0-0,1|(873,466)|
10,78,Descomposición aerobia en función de OD,1102,466,65,29,8,131,0,2,0,0,0,0,0-0-0,0-0-0,|12||255-0-0
1,79,78,66,0,0,0,0,1,64,0,255-0-0,|12||0-0-0,1|(1027,466)|
1,80,3,27,1,0,0,0,1,64,0,128-0-0,|12||0-0-0,1|(744,602)|
12,81,0,589,32,575,26,8,132,0,30,-1,0,0,0,0-128-0,128-192-0,|16|B|0-128-0
EUTROFIZACIÓN CULTURAL
\\\---/// Sketch information - do not modify anything except names
V300  Do not put anything below this section - it will be ignored
*
$192-192-192,0,Arial|12||0-0-0|0-0-0|0-0-255|-1--1--1|-1--1--1|96,96,100
12,1,0,349,312,311,285,3,188,0,0,1,0,0,0
Eutrofización_cultural_2
12,2,0,801,65,80,20,3,124,0,0,0,0,0,0
P procedente de contaminación,0,0.05,0.0001
///---\\\
:GRAPH Eutrofización_cultural
:TITLE Eutrofización cultural
:X-LABEL Tiempo (años)
:SCALE
:VAR Fósforo total
:SCALE
:VAR Biomasa de fitoplancton
:SCALE
:VAR OD hipolimnético

:GRAPH Eutrofización_cultural_2
:TITLE Eutrofización cultural 2
:X-LABEL Tiempo (años)
:SCALE
:VAR Fósforo total
:SCALE
:VAR Biomasa muerta
:SCALE
:VAR Profundidad de Secchi
:SCALE
:VAR Biomasa muerta
:SCALE
:VAR Tasa de descomposición aerobia
:SCALE
:VAR Tasa de descomposición anaerobia
:L<%^E!@
1:Current.vdf
9:Current
15:0,0,0,0,0,0
19:100,0
27:2,
34:0,
4:Time
5:Tasa de agotamiento de OD
24:0
25:10
26:10