/*
 * 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;
// Warning: unit conversion turned off due to unit errors in 1 equation(s)
unit conversion off;
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:  vs, vm, kdm, B, Km, c, A, M_init
	realDomain time hour;
	time.min=0;
	extern time.max;
	extern time.delta;
	extern real vs flux;
	extern real vm flux;
	extern real kdm first_order_rate_constant;
	extern real B nanomolar;
	extern real Km nanomolar;
	extern real c dimensionless;
	extern real A nanomolar;
	extern real M_init nanomolar;
	real M(time) nanomolar;
	when(time=time.min) M=M_init;

	// <component name="interface">
	M:time=(vs*(A^c/(B^c+A^c))-(vm*(M/(Km+M))+kdm*M));
}