/*
 * A Detailed Computational Model for the Mammalian Circadian Clock
 * 
 * Model Status
 * 
 * This is an import for the model "leloup_goldbeter_2003_1.1model.cellml".
 * It is used to define a reoccuring structure within the model
 * and may be used for multiple equations in the model. This file
 * is known to run only in OpenCell in conjunction with the base
 * file.
 * 
 * Model Structure
 * 
 * Toward a detailed computational model for the mammalian circadian
 * clock, Jean-Christophe Leloup and Albert Goldbeter, 2003,PNAS,
 * 100, 7051-7056. PubMed ID: 12775757
 */

import nsrunit;
unit conversion on;
unit hour=3600 second^1;
// unit nanomolar predefined
unit flux=2.7777778E-10 meter^(-3)*second^(-1)*mole^1;
unit first_order_rate_constant=2.7777778E-4 second^(-1);
unit second_order_rate_constant=2.7777778E2 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:  Kp, Kdp, kc, ke, kb, ka, kdn, V1, V2, YP,
	//  E, A, C, Y_init
	realDomain time hour;
	time.min=0;
	extern time.max;
	extern time.delta;
	extern real Kp nanomolar;
	extern real Kdp nanomolar;
	extern real kc second_order_rate_constant;
	extern real ke first_order_rate_constant;
	extern real kb first_order_rate_constant;
	extern real ka first_order_rate_constant;
	extern real kdn first_order_rate_constant;
	extern real V1 flux;
	extern real V2 flux;
	extern real YP nanomolar;
	extern real E nanomolar;
	extern real A nanomolar;
	extern real C nanomolar;
	extern real Y_init nanomolar;
	real Y(time) nanomolar;
	when(time=time.min) Y=Y_init;

	// <component name="interface">
	Y:time=(V2*(YP/(Kdp+YP))+kb*A+ke*E-(V1*(Y/(Kp+Y))+ka*Y+kc*C*Y+kdn*Y));
}