/*
 * A Model for Circadian Rhythms in Drosophila Incorporating the
 * Formation of a Complex between the PER and TIM Proteins
 * 
 * Model Status
 * 
 * This is an import for the model "leloup_1998a_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.
 */

import nsrunit;
// Warning: unit conversion turned off due to unit errors in 1 equation(s)
unit conversion off;
unit hour=3600 second^1;
unit per_hour=2.7777778E-4 second^(-1);
unit per_litre=1E3 meter^(-3);
// unit nanomolar predefined
unit nanomolar_hour=2.7777778E-10 meter^(-3)*second^(-1)*mole^1;
unit per_nanomolar_hour=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:  v_s, v_m, K_I, K_m, C_N, n, k_d, M_init
	realDomain time hour;
	time.min=0;
	extern time.max;
	extern time.delta;
	extern real v_s nanomolar_hour;
	extern real v_m nanomolar_hour;
	extern real K_I nanomolar;
	extern real K_m nanomolar;
	extern real C_N nanomolar;
	extern real n dimensionless;
	extern real k_d per_hour;
	extern real M_init nanomolar;
	real M(time) nanomolar;
	when(time=time.min) M=M_init;

	// <component name="interface">
	M:time=(v_s*(K_I^n/(K_I^n+C_N^n))-v_m*(M/(K_m+M))-k_d*M);
}