/*
 * Limit cycle models for circadian rhythms based on transcriptional
 * regulation in Drosophila and Neurospora
 * 
 * Model Status
 * 
 * This is an import for the model "leloup_1999_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
 * 
 * The original paper reference is cited below:
 * 
 * Limit Cycle Models for Circadian Rhythms Based on Transcriptional
 * Regulation in Drosophila and Neurospora, Jean-Christophe Leloup,
 * Didier Gonze, and Albert Goldbeter, 1999, Journal of Biological
 * Rhythms, 14, 433-448. PubMed ID: 10643740
 */

import nsrunit;
unit conversion on;
// unit nanomolar predefined
unit hour=3600 second^1;
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:  vd, KdP, V3, V4, k3, k4, K3, K4, kd, Y1, Z2,
	//  C, Y2_init
	realDomain time hour;
	time.min=0;
	extern time.max;
	extern time.delta;
	extern real vd flux;
	extern real KdP nanomolar;
	extern real V3 flux;
	extern real V4 flux;
	extern real k3 second_order_rate_constant;
	extern real k4 first_order_rate_constant;
	extern real K3 nanomolar;
	extern real K4 nanomolar;
	extern real kd first_order_rate_constant;
	extern real Y1 nanomolar;
	extern real Z2 nanomolar;
	extern real C nanomolar;
	extern real Y2_init nanomolar;
	real Y2(time) nanomolar;
	when(time=time.min) Y2=Y2_init;

	// <component name="interface">
	Y2:time=(V3*(Y1/(K3+Y1))+k4*C-(V4*(Y2/(K4+Y2))+k3*Y2*Z2+vd*(Y2/(KdP+Y2))+kd*Y2));
}