import nsrunit;
unit conversion on;
unit s=1 second^1;
unit uM=1E-3 meter^(-3)*mole^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);
unit per_uM_per_s=1E3 meter^3*second^(-1)*mole^(-1);
unit um2_per_s=1E-12 meter^2*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:  L, Gd
	realDomain t s;
	t.min=0;
	extern t.max;
	extern t.delta;
	extern real L uM;
	extern real Gd per_um2;
	real J2(t) per_um2_per_s;
	real J3(t) per_um2_per_s;
	real J6(t) per_um2_per_s;
	real J1(t) per_um2_per_s;
	real kf1 per_uM_per_s;
	kf1=0.0003;
	real kr1 per_s;
	real Kd1 uM;
	Kd1=3.00E-5;
	real kf2 um2_per_s;
	kf2=2.75E-4;
	real kr2 per_s;
	real Kd2 per_um2;
	Kd2=27500;
	real kf3 um2_per_s;
	kf3=1.00;
	real kr3 per_s;
	kr3=0.00100;
	real J4(t) per_um2_per_s;
	real kf4 per_uM_per_s;
	kf4=0.3;
	real kr4 per_s;
	real Kd4 uM;
	Kd4=3.00E-5;
	real J5(t) per_um2_per_s;
	real kf5 per_s;
	kf5=0.0004;
	real kf6 per_s;
	kf6=1.00;
	real R(t) per_um2;
	when(t=t.min) R=13.9;
	real Rl(t) per_um2;
	when(t=t.min) Rl=0.00;
	real Rg(t) per_um2;
	when(t=t.min) Rg=5.06;
	real Rlg(t) per_um2;
	when(t=t.min) Rlg=0.00;
	real Rlgp(t) per_um2;
	when(t=t.min) Rlgp=0.00;

	// <component name="GPCR_Cycle">
	kr1=(kf1*Kd1);
	J1=(kf1*R*L-kr1*Rl);
	kr2=(kf2*Kd2);
	J2=(kf2*R*Gd-kr2*Rg);
	R:t=((-1)*J1-J2);
	J3=(kf3*Rl*Gd-kr3*Rlg);
	Rl:t=(J1+J6-J3);
	kr4=(kf4*Kd4);
	J4=(kf4*L*Rg-kr4*Rlg);
	Rg:t=(J2-J4);
	J5=(kf5*Rlg);
	Rlgp:t=J5;
	J6=(kf6*Rlg);
	Rlg:t=(J3-J5+(J4-J6));
}