import nsrunit;
unit conversion on;
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 {
	//Warning:  the following variables were set 'extern' or given
	//  an initial value of '0' because the model would otherwise be
	//  underdetermined:  Ccn, act_N, Ntot
	realDomain t s;
	t.min=0;
	extern t.max;
	extern t.delta;
	real NFATN_c(t) nM;
	when(t=t.min) NFATN_c=0.0017047398;
	real NFATN_n(t) nM;
	when(t=t.min) NFATN_n=0.50910553;
	real NFATp_c(t) nM;
	when(t=t.min) NFATp_c=0.9825855442;
	real NFATp_n(t) nM;
	when(t=t.min) NFATp_n=0.27638027;
	extern real Ccn dimensionless;
	extern real act_N dimensionless;
	extern real Ntot nM;
	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="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);
}