// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// rPer: EmptySet <=> Per 
// rdClk: EmptySet <=> dClk 
// rPerD: Per <=> EmptySet 
// rdClkD: dClk <=> EmptySet 

math main {
  realDomain time second;
  time.min=0;
  extern time.max;
  extern time.delta;

  // variable definitions
  real CELL = 1E-15 L;
  real dClkF_tau1(time) nM;
  real dClkF_tau2(time) nM;
  real tau1 = 10 hr;
  real tau2 = 10 hr;
  real Vsp = .5 nmol/L/hr;
  real K1 = .3 nM;
  real Vsc = .25 nmol/L/hr;
  real K2 = .1 nM;
  real kdc = .5 1/hr;
  real kdp = .5 1/hr;
  real EmptySet = 0 nM;
  real Per(time) nM;
  real dClk(time) nM;
  real dClkF(time) nM;
  real rPer(time) nmol/hr;
  real rdClk(time) nmol/hr;
  real rPerD(time) nmol/hr;
  real rdClkD(time) nmol/hr;

  // equations
  dClkF_tau1 = if ((dClk(if (time<(time.min+tau1)) time.min else time-(tau1))-Per(if (time<(time.min+tau1)) time.min else time-(tau1)))<0) 0 else dClk(if (time<(time.min+tau1)) time.min else time-(tau1))-Per(if (time<(time.min+tau1)) time.min else time-(tau1));
  dClkF_tau2 = if ((dClk(if (time<(time.min+tau2)) time.min else time-(tau2))-Per(if (time<(time.min+tau2)) time.min else time-(tau2)))<0) 0 else dClk(if (time<(time.min+tau2)) time.min else time-(tau2))-Per(if (time<(time.min+tau2)) time.min else time-(tau2));
  when (time=time.min) Per = 5E-16;
  Per:time = (rPer + -1*rPerD)/CELL;
  when (time=time.min) dClk = 1E-16;
  dClk:time = (rdClk + -1*rdClkD)/CELL;
  dClkF = if ((dClk-Per)<0) 0 else dClk-Per;
  rPer = Vsp*(dClkF_tau1/(K1+dClkF_tau1))*CELL;
  rdClk = CELL*Vsc*(K2/(K2+dClkF_tau2));
  rPerD = kdc*Per*CELL;
  rdClkD = kdp*dClk*CELL;

  // variable properties
  CELL.sbmlRole="compartment";
  dClkF_tau1.sbmlRole="parameter";
  dClkF_tau2.sbmlRole="parameter";
  tau1.sbmlRole="parameter";
  tau2.sbmlRole="parameter";
  Vsp.sbmlRole="parameter";
  K1.sbmlRole="parameter";
  Vsc.sbmlRole="parameter";
  K2.sbmlRole="parameter";
  kdc.sbmlRole="parameter";
  kdp.sbmlRole="parameter";
  EmptySet.sbmlRole="species";
  EmptySet.sbmlCompartment="CELL";
  Per.sbmlRole="species";
  Per.sbmlCompartment="CELL";
  dClk.sbmlRole="species";
  dClk.sbmlCompartment="CELL";
  dClkF.sbmlRole="species";
  dClkF.sbmlCompartment="CELL";
  rPer.sbmlRole="rate";
  rdClk.sbmlRole="rate";
  rPerD.sbmlRole="rate";
  rdClkD.sbmlRole="rate";
}