import nsrunit;
// Warning: unit conversion turned off due to unit errors in 1 equation(s)
unit conversion off;
unit hour=3600 second^1;
unit per_hour=2.7777778E-4 second^(-1);
unit mole_per_litre=1E3 meter^(-3)*mole^1;
unit per_mole_per_litre=.001 meter^3*mole^(-1);
unit mole_per_litre_per_hour=.27777778 meter^(-3)*second^(-1)*mole^1;
unit per_mole_per_litre_per_hour=2.7777778E-7 meter^3*second^(-1)*mole^(-1);

math main {
	realDomain time hour;
	time.min=0;
	extern time.max;
	extern time.delta;
	real CycE(time) mole_per_litre;
	when(time=time.min) CycE=0;
	real v91 per_hour;
	v91=0.8;
	real HBEGF(time) mole_per_litre;
	when(time=time.min) HBEGF=0;
	real k91 mole_per_litre;
	k91=40;
	real v92 per_hour;
	v92=0.16;
	real CycD(time) mole_per_litre;
	when(time=time.min) CycD=0;
	real k92 mole_per_litre;
	k92=10;
	real d11 per_hour;
	d11=800;
	real CycECKI(time) mole_per_litre;
	when(time=time.min) CycECKI=0;
	real d9 per_hour;
	d9=0.0002;
	real deg(time) mole_per_litre;
	when(time=time.min) deg=0;
	real v11 per_mole_per_litre_per_hour;
	v11=800;
	real CKI(time) mole_per_litre;
	when(time=time.min) CKI=0;
	real v8 per_hour;
	v8=0.8;
	real IEG(time) mole_per_litre;
	when(time=time.min) IEG=0;
	real k8 mole_per_litre;
	k8=0.005;
	real d10 per_hour;
	d10=8.1;
	real d8 per_hour;
	d8=0.002;
	real v10 per_mole_per_litre_per_hour;
	v10=0.8;
	real CycDCKI(time) mole_per_litre;
	when(time=time.min) CycDCKI=0;
	real v7 per_hour;
	v7=1.6;
	real PIC(time) mole_per_litre;
	when(time=time.min) PIC=0;
	real d7 per_hour;
	d7=0.3;
	real HGF(time) mole_per_litre;
	when(time=time.min) HGF=0;
	real v5 per_mole_per_litre;
	v5=8.1;
	real k5 mole_per_litre;
	k5=0.125;
	real km5 mole_per_litre;
	km5=0.18;
	real PAI(time) mole_per_litre;
	when(time=time.min) PAI=0;
	real f5 mole_per_litre_per_hour;
	f5=5;
	real d5 per_hour;
	d5=0.008;
	real PGE(time) mole_per_litre;
	when(time=time.min) PGE=0;
	real v2 per_hour;
	v2=0.8;
	real d2 per_hour;
	d2=0.01;
	real v1 mole_per_litre_per_hour;
	v1=2;
	real km1 mole_per_litre;
	km1=1;
	real D mole_per_litre;
	real k1 mole_per_litre;
	k1=0.5;
	real d1 per_hour;
	d1=0.7;
	real I0 mole_per_litre;
	I0=4;
	real d0 per_hour;
	d0=0.05;
	real t hour;
	t=0;
	real v12 mole_per_litre_per_hour;
	v12=0.56;
	real d12 per_hour;
	d12=0.025;
	real v6 mole_per_litre_per_hour;
	v6=24;
	real k61 mole_per_litre;
	k61=0.9;
	real k62 mole_per_litre;
	k62=0.0005;
	real k63 mole_per_litre;
	k63=1;
	real d6 per_hour;
	d6=0.1;
	real v3 mole_per_litre_per_hour;
	v3=6;
	real k3 mole_per_litre;
	k3=0.5;
	real d3 per_hour;
	d3=10;
	real v4 mole_per_litre_per_hour;
	v4=60;
	real k4 mole_per_litre;
	k4=0.5;
	real d4 per_hour;
	d4=10;

	// <component name="environment">

	// <component name="CycE">
	CycE:time=((1 mole_per_litre)*v91*HBEGF^2/(k91^2+HBEGF^2)+(1 mole_per_litre)*v92*CycD^4/(k92^4+CycD^4)+d11*CycECKI-(d9*deg^2+v11*CKI)*CycE);

	// <component name="CycD">
	CycD:time=(v8*IEG*HBEGF/(k8+IEG)+d10*CycDCKI-(d8*(deg/(1 mole_per_litre))^1.5+v10*CKI)*CycD);

	// <component name="CycECKI">
	CycECKI:time=(v11*CKI*CycE-d11*CycECKI);

	// <component name="CycDCKI">
	CycDCKI:time=(v10*CKI*CycD-d10*CycDCKI);

	// <component name="CKI">
	CKI:time=(v7*PIC+d10*CycDCKI+d11*CycECKI-(d7+v10*CycD+v11*CycE)*CKI);

	// <component name="HGF">
	HGF:time=(v5*IEG/(1+(PAI/km5)^4)*((1 mole_per_litre_per_hour)+f5*HGF/(k5+HGF))-d5*HGF);

	// <component name="PGE">
	PGE:time=(v2*PIC-d2*PGE);

	// <component name="PIC">
	PIC:time=(v1/(1+(PGE/km1)^2)*(D/(k1+D))-d1*PIC);

	// <component name="D">
	D=(I0*exp((-1)*d0*t));

	// <component name="deg">
	deg:time=(v12*sqrt(CycE/(1 mole_per_litre))-d12*deg);

	// <component name="HBEGF">
	HBEGF:time=(v6*(D^4/(k61^4+D^4))*(IEG/(k62+IEG))*(HGF/(k63+HGF))-d6*HBEGF);

	// <component name="IEG">
	IEG:time=(v3*PIC^3/(k3^3+PIC^3)-d3*IEG);

	// <component name="PAI">
	PAI:time=(v4*PIC^3/(k4^3+PIC^3)-d4*PAI);

	// <component name="model_parameters">
}