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 per_s=1 second^(-1);

math main {
	//Warning:  the following variables were set 'extern' or given
	//  an initial value of '0' because the model would otherwise be
	//  underdetermined:  J_gain_Gd, J_gain_Gt
	realDomain t s;
	t.min=0;
	extern t.max;
	extern t.delta;
	real Gd(t) per_um2;
	when(t=t.min) Gd=10000;
	extern real J_gain_Gd per_um2_per_s;
	real Gt(t) per_um2;
	when(t=t.min) Gt=0.00;
	extern real J_gain_Gt per_um2_per_s;
	real J1(t) per_um2_per_s;
	real kf1 per_s;
	kf1=0.150;

	// <component name="Galpha">
	Gd:t=(J1+J_gain_Gd);
	Gt:t=(J_gain_Gt-J1);
	J1=(kf1*Gt);
}