import nsrunit;
unit conversion on;
unit s=1 second^1;
unit um2=1E-12 meter^2;
unit per_um2=1E12 meter^(-2);
unit per_um2_per_s=1E12 meter^(-2)*second^(-1);
unit um2_per_s=1E-12 meter^2*second^(-1);
unit per_s=1 second^(-1);
unit uM=1E-3 meter^(-3)*mole^1;
unit per_uM_per_s=1E3 meter^3*second^(-1)*mole^(-1);

math main {
	//Warning:  the following variables were set 'extern' or given
	//  an initial value of '0' because the model would otherwise be
	//  underdetermined:  Gt, Ca
	realDomain t s;
	t.min=0;
	extern t.max;
	extern t.delta;
	real Pc(t) per_um2;
	when(t=t.min) Pc=9.09;
	real Pcg(t) per_um2;
	when(t=t.min) Pcg=0.00;
	real J2(t) per_um2_per_s;
	real J3(t) per_um2_per_s;
	extern real Gt per_um2;
	extern real Ca uM;
	real J1(t) per_um2_per_s;
	real J4(t) per_um2_per_s;
	real J5(t) per_um2_per_s;
	real J6(t) per_um2_per_s;
	real kf1 per_uM_per_s;
	kf1=0.0167;
	real kr1 per_s;
	kr1=0.0167;
	real kf2 um2_per_s;
	kf2=0.00420;
	real kr2 per_s;
	kr2=1.00;
	real kf3 um2_per_s;
	kf3=0.0420;
	real kr3 per_s;
	kr3=1.00;
	real kf4 per_uM_per_s;
	kf4=0.0334;
	real kr4 per_s;
	real Kd4 uM;
	Kd4=0.1;
	real kf5 per_s;
	kf5=6.00;
	real kf6 per_s;
	kf6=6.00;
	real P(t) per_um2;
	when(t=t.min) P=90.9;
	real Pg(t) per_um2;
	when(t=t.min) Pg=0.00;

	// <component name="PLC_Cycle">
	J2=(kf2*P*Gt-kr2*Pg);
	J1=(kf1*P*Ca-kr1*Pc);
	J3=(kf3*Pc*Gt-kr3*Pcg);
	kr4=(kf4*Kd4);
	J4=(kf4*Pg*Ca-kr4*Pcg);
	J5=(kf5*Pcg);
	J6=(kf6*Pg);
	P:t=(J6-(J2+J1));
	Pg:t=(J2-(J4+J6));
	Pc:t=(J1+J5-J3);
	Pcg:t=(J3+J4-J5);
}