import nsrunit;
unit conversion on;
unit s=1 second^1;
unit nM=1E-6 meter^(-3)*mole^1;
unit s=1 second^1;
unit nM=1E-6 meter^(-3)*mole^1;
unit per_s=1 second^(-1);
unit nM_per_s=1E-6 meter^(-3)*second^(-1)*mole^1;
unit s=1 second^1;
unit nM=1E-6 meter^(-3)*mole^1;
unit s=1 second^1;
unit nM=1E-6 meter^(-3)*mole^1;
unit per_s=1 second^(-1);
unit nM_per_s=1E-6 meter^(-3)*second^(-1)*mole^1;
unit per_nM_s=1E6 meter^3*second^(-1)*mole^(-1);

math main {
	realDomain t s;
	t.min=0;
	extern t.max;
	extern t.delta;
	real Ccn dimensionless;
	Ccn=50;
	real nuclearNFAT(t) nM;
	real NFATN_n(t) nM;
	when(t=t.min) NFATN_n=0.10902441;
	real NFATp_n(t) nM;
	when(t=t.min) NFATp_n=0.077917275;
	real stimEnd s;
	stimEnd=1800;
	real Ca(t) nM;
	real M nM;
	M=6000;
	real Ntot nM;
	Ntot=1000;
	real K_mN nM;
	K_mN=535;
//	Var below replaced by constant in model eqns to satisfy unit correction
//	real n dimensionless;
//	n=2.92;
	real K_dN nM;
	K_dN=1760;
	real act_N(t) dimensionless;
	real NFATN_c(t) nM;
	when(t=t.min) NFATN_c=0.0004806012;
	real NFATp_c(t) nM;
	when(t=t.min) NFATp_c=0.9957805651;
	real k_f1 per_nM_s;
	k_f1=7.68934e-6;
	real k_r1 per_s;
	k_r1=0.019256;
	real k_f2 per_s;
	k_f2=0.00144192;
	real k_f3 per_s;
	k_f3=0.000361944;
	real k_r3 per_nM_s;
	k_r3=4.70813e-5;
	real k_f4 per_s;
	k_f4=0.000444695;
	real J1(t) nM_per_s;
	real J2(t) nM_per_s;
	real J3(t) nM_per_s;
	real J4(t) nM_per_s;

	// <component name="environment">

	// <component name="geometry">

	// <component name="NFAT_Cycling_interface">
	nuclearNFAT=(NFATN_n+NFATp_n);

	// <component name="calcium">
	Ca=(if ((t>=(0 s)) and (t<=(180 s))) ((459 nM)-(576 nM))*t/(180 s)+(576 nM) else if ((t>(180 s)) and (t<=(360 s))) ((406 nM)-(459 nM))*(t-(180 s))/(180 s)+(459 nM) else if ((t>(360 s)) and (t<=(540 s))) ((377 nM)-(406 nM))*(t-(360 s))/(180 s)+(406 nM) else if ((t>(540 s)) and (t<=(720 s))) ((348 nM)-(377 nM))*(t-(540 s))/(180 s)+(377 nM) else if ((t>(720 s)) and (t<=(900 s))) ((336 nM)-(348 nM))*(t-(720 s))/(180 s)+(348 nM) else if ((t>(900 s)) and (t<=(1080 s))) ((334 nM)-(336 nM))*(t-(900 s))/(180 s)+(336 nM) else if ((t>(1080 s)) and (t<=stimEnd)) (334 nM) else (60 nM));

	// <component name="calcineurin">
	act_N=(Ca^2.92/(Ca^2.92+K_mN^2.92*(1+K_dN/M)));

	// <component name="NFAT_Cycling">
	J1=(k_f1*NFATp_c*Ntot*act_N-k_r1*NFATN_c*(1-act_N));
	J2=(k_f2*NFATN_c);
	J3=(k_f3*NFATN_n*(1-act_N)-k_r3*NFATp_n*Ntot*act_N);
	J4=(k_f4*NFATp_n);
	NFATp_c:t=(J4/Ccn-J1);
	NFATN_c:t=(J1-J2);
	NFATN_n:t=(J2*Ccn-J3);
	NFATp_n:t=(J3-J4);
}