// 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
// Reaction1: EmptySet <=> Perm 
// Reaction2: Perm <=> EmptySet 
// Reaction3: EmptySet <=> Timm 
// Reaction4: Timm <=> EmptySet 
// Reaction5: EmptySet <=> Clkm 
// Reaction6: Clkm <=> EmptySet 
// Reaction7: CCc <=> CCn 
// Reaction8: CCn <=> CCc 
// Reaction9: PTn <=> PTc 
// Reaction10: PTc <=> PTn 
// Reaction11: species_0000012 Clkc => CCc 
// Reaction12: Perc Timc => PTc 
// Reaction16: EmptySet <=> Timc 
// Reaction18: EmptySet <=> Clkc 
// Reaction19: EmptySet <=> Perc 
// Reaction20: Perc <=> EmptySet 
// Reaction21: PTc <=> EmptySet 
// Reaction23: PTn <=> EmptySet 
// Reaction24: CCc <=> EmptySet 
// Reaction25: Clkc <=> EmptySet 
// Reaction26: CCn <=> EmptySet 
// Reaction27: Timc <=> EmptySet 
// Reaction28: Perm <=> EmptySet 
// Reaction29: Perc <=> EmptySet 
// Reaction30: Timm <=> EmptySet 
// Reaction31: Timc <=> EmptySet 
// Reaction32: PTc <=> EmptySet 
// Reaction33: PTn <=> EmptySet 
// Reaction34: Clkm <=> EmptySet 
// Reaction35: Clkc <=> EmptySet 
// Reaction36: CCc <=> EmptySet 
// Reaction37: CCn <=> EmptySet 

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

  // variable definitions
  real Drosophilia = 1 L;
  real compartment_0000003 = 1 L;
  real compartment_0000002 = 1 L;
  real EmptySet = 0 nM;
  real CCc(time) nM;
  real CCn(time) nM;
  real Clkc(time) nM;
  real Clkm(time) nM;
  real Perc(time) nM;
  real Perm(time) nM;
  real PTc(time) nM;
  real PTn(time) nM;
  real Timc(time) nM;
  real Timm(time) nM;
  real species_0000012 nM;
  real species_0000013 nM;
  real Reaction1(time) nmol/hr;
  real a = 1;
  real A1 = .45;
  real B1 = 0;
  real c1 = 0;
  real r1 = 1.02;
  real s1 = 1.45;
  real r = 4;
  real Reaction2(time) nmol/hr;
  real D0 = .012;
  real Reaction3(time) nmol/hr;
  real A2 = .45;
  real B2 = 0;
  real c2 = 0;
  real r2 = 1.02;
  real s3 = 1.45;
  real Reaction4(time) nmol/hr;
  real Reaction5(time) nmol/hr;
  real A3 = .8;
  real B3 = .6;
  real c3 = 0;
  real r3 = .89;
  real s5 = 1.63;
  real Reaction6(time) nmol/hr;
  real Reaction7(time) nmol/hr;
  real k3 = 2;
  real T3 = 1.63;
  real Reaction8(time) nmol/hr;
  real k4 = 2;
  real T4 = .52;
  real Reaction9(time) nmol/hr;
  real k2 = 2;
  real T2 = .72;
  real Reaction10(time) nmol/hr;
  real k1 = 2;
  real T1 = 1.73;
  real Reaction11(time) nmol/hr;
  real v3 = 1.63;
  real parameter_0000073 = 1.63;
  real Reaction12(time) nmol/hr;
  real v1 = 1.45;
  real parameter_0000072 = 1.45;
  real Reaction16(time) nmol/hr;
  real s4 = .48;
  real Reaction18(time) nmol/hr;
  real s6 = .47;
  real Reaction19(time) nmol/hr;
  real s2 = .48;
  real Reaction20(time) nmol/hr;
  real Reaction21(time) nmol/hr;
  real Reaction23(time) nmol/hr;
  real Reaction24(time) nmol/hr;
  real Reaction25(time) nmol/hr;
  real Reaction26(time) nmol/hr;
  real Reaction27(time) nmol/hr;
  real Reaction28(time) nmol/hr;
  real D1 = .94;
  real L1 = .3;
  real Reaction29(time) nmol/hr;
  real D2 = .44;
  real L2 = .2;
  real Reaction30(time) nmol/hr;
  real D3 = .94;
  real L3 = .3;
  real Reaction31(time) nmol/hr;
  real D4 = .44;
  real L4 = .2;
  real Reaction32(time) nmol/hr;
  real D5 = .44;
  real L5 = .2;
  real Reaction33(time) nmol/hr;
  real D6 = .29;
  real L6 = .2;
  real Reaction34(time) nmol/hr;
  real D7 = .54;
  real L7 = .13;
  real Reaction35(time) nmol/hr;
  real D8 = .6;
  real L8 = .2;
  real Reaction36(time) nmol/hr;
  real D9 = .6;
  real L9 = .2;
  real Reaction37(time) nmol/hr;
  real D10 = .3;
  real L10 = .2;

  // equations
  when (time=time.min) CCc = .3/compartment_0000003;
  CCc:time = (-1*Reaction7 + Reaction8 + Reaction11 + -1*Reaction24 + -1*Reaction36)/compartment_0000003;
  when (time=time.min) CCn = .4/compartment_0000002;
  CCn:time = (Reaction7 + -1*Reaction8 + -1*Reaction26 + -1*Reaction37)/compartment_0000002;
  when (time=time.min) Clkc = .2/compartment_0000003;
  Clkc:time = (-1*Reaction11 + Reaction18 + -1*Reaction25 + -1*Reaction35)/compartment_0000003;
  when (time=time.min) Clkm = .1/compartment_0000003;
  Clkm:time = (Reaction5 + -1*Reaction6 + -1*Reaction34)/compartment_0000003;
  when (time=time.min) Perc = .6/compartment_0000003;
  Perc:time = (-1*Reaction12 + Reaction19 + -1*Reaction20 + -1*Reaction29)/compartment_0000003;
  when (time=time.min) Perm = .5/compartment_0000003;
  Perm:time = (Reaction1 + -1*Reaction2 + -1*Reaction28)/compartment_0000003;
  when (time=time.min) PTc = .9/compartment_0000003;
  PTc:time = (Reaction9 + -1*Reaction10 + Reaction12 + -1*Reaction21 + -1*Reaction32)/compartment_0000003;
  when (time=time.min) PTn = 1/compartment_0000002;
  PTn:time = (-1*Reaction9 + Reaction10 + -1*Reaction23 + -1*Reaction33)/compartment_0000002;
  when (time=time.min) Timc = .8/compartment_0000003;
  Timc:time = (-1*Reaction12 + Reaction16 + -1*Reaction27 + -1*Reaction31)/compartment_0000003;
  when (time=time.min) Timm = .7/compartment_0000003;
  Timm:time = (Reaction3 + -1*Reaction4 + -1*Reaction30)/compartment_0000003;
  when (time=time.min) species_0000012 = 1/compartment_0000003;
  when (time=time.min) species_0000013 = 1/compartment_0000003;
  Reaction1 = compartment_0000003*(c1+(B1+(CCn/A1)^a)*s1/(1+B1+(CCn/A1)^a+(PTn/r1)^r));
  Reaction2 = compartment_0000003*D0*Perm;
  Reaction3 = compartment_0000003*(c2+(B2+(CCn/A2)^a)*s3/(1+B2+(CCn/A2)^a+(PTn/r2)^r));
  Reaction4 = Drosophilia*D0*Timm;
  Reaction5 = compartment_0000003*(c3+(B3+(PTn/A3)^a)*s5/(1+B3+(PTn/A3)^a+(CCn/r3)^r));
  Reaction6 = Drosophilia*Clkm*D0;
  Reaction7 = compartment_0000003*(CCc*T3/(k3+CCc));
  Reaction8 = compartment_0000002*(CCn*T4/(k4+CCn));
  Reaction9 = compartment_0000002*(PTn*T2/(k2+PTn));
  Reaction10 = compartment_0000003*(PTc*T1/(k1+PTc));
  Reaction11 = compartment_0000003*(Clkc*v3*species_0000012-parameter_0000073*CCc);
  Reaction12 = compartment_0000003*(Perc*Timc*v1-parameter_0000072*PTc);
  Reaction16 = compartment_0000003*s4*Timm;
  Reaction18 = compartment_0000003*Clkm*s6;
  Reaction19 = compartment_0000003*s2*Perm;
  Reaction20 = Drosophilia*D0*Perc;
  Reaction21 = compartment_0000003*D0*PTc;
  Reaction23 = compartment_0000002*D0*PTn;
  Reaction24 = compartment_0000003*CCc*D0;
  Reaction25 = compartment_0000003*Clkc*D0;
  Reaction26 = compartment_0000002*CCn*D0;
  Reaction27 = compartment_0000003*D0*Timc;
  Reaction28 = compartment_0000003*(D1*Perm/(L1+Perm));
  Reaction29 = compartment_0000003*(D2*species_0000013*Perc/(L2+Perc));
  Reaction30 = compartment_0000003*(D3*Timm/(L3+Timm));
  Reaction31 = compartment_0000003*(D4*Timc/(L4+Timc));
  Reaction32 = compartment_0000003*(D5*PTc/(L5+PTc));
  Reaction33 = compartment_0000002*(D6*PTn/(L6+PTn));
  Reaction34 = compartment_0000003*(Clkm*D7/(Clkm+L7));
  Reaction35 = compartment_0000003*(Clkc*D8/(Clkc+L8));
  Reaction36 = compartment_0000003*(CCc*D9/(CCc+L9));
  Reaction37 = compartment_0000002*(CCn*D10/(CCn+L10));

  // variable properties
  Drosophilia.sbmlRole="compartment";
  compartment_0000003.sbmlRole="compartment";
  compartment_0000003.sbmlName="cytoplasm";
  compartment_0000002.sbmlRole="compartment";
  compartment_0000002.sbmlName="nucleus";
  EmptySet.sbmlRole="species";
  EmptySet.sbmlCompartment="Drosophilia";
  CCc.sbmlRole="species";
  CCc.sbmlCompartment="compartment_0000003";
  CCn.sbmlRole="species";
  CCn.sbmlCompartment="compartment_0000002";
  Clkc.sbmlRole="species";
  Clkc.sbmlCompartment="compartment_0000003";
  Clkm.sbmlRole="species";
  Clkm.sbmlCompartment="compartment_0000003";
  Perc.sbmlRole="species";
  Perc.sbmlCompartment="compartment_0000003";
  Perm.sbmlRole="species";
  Perm.sbmlCompartment="compartment_0000003";
  PTc.sbmlRole="species";
  PTc.sbmlCompartment="compartment_0000003";
  PTn.sbmlRole="species";
  PTn.sbmlCompartment="compartment_0000002";
  Timc.sbmlRole="species";
  Timc.sbmlCompartment="compartment_0000003";
  Timm.sbmlRole="species";
  Timm.sbmlCompartment="compartment_0000003";
  species_0000012.sbmlRole="species";
  species_0000012.sbmlCompartment="compartment_0000003";
  species_0000013.sbmlRole="species";
  species_0000013.sbmlCompartment="compartment_0000003";
  Reaction1.sbmlRole="rate";
  a.sbmlRole="parameter";
  A1.sbmlRole="parameter";
  B1.sbmlRole="parameter";
  c1.sbmlRole="parameter";
  c1.sbmlName="C1";
  r1.sbmlRole="parameter";
  r1.sbmlName="R1";
  s1.sbmlRole="parameter";
  s1.sbmlName="S1";
  r.sbmlRole="parameter";
  Reaction2.sbmlRole="rate";
  D0.sbmlRole="parameter";
  Reaction3.sbmlRole="rate";
  A2.sbmlRole="parameter";
  B2.sbmlRole="parameter";
  c2.sbmlRole="parameter";
  c2.sbmlName="C2";
  r2.sbmlRole="parameter";
  r2.sbmlName="R2";
  s3.sbmlRole="parameter";
  s3.sbmlName="S3";
  Reaction4.sbmlRole="rate";
  Reaction5.sbmlRole="rate";
  A3.sbmlRole="parameter";
  B3.sbmlRole="parameter";
  c3.sbmlRole="parameter";
  c3.sbmlName="C3";
  r3.sbmlRole="parameter";
  r3.sbmlName="R3";
  s5.sbmlRole="parameter";
  s5.sbmlName="S5";
  Reaction6.sbmlRole="rate";
  Reaction7.sbmlRole="rate";
  k3.sbmlRole="parameter";
  k3.sbmlName="K3";
  T3.sbmlRole="parameter";
  Reaction8.sbmlRole="rate";
  k4.sbmlRole="parameter";
  T4.sbmlRole="parameter";
  Reaction9.sbmlRole="rate";
  k2.sbmlRole="parameter";
  k2.sbmlName="K2";
  T2.sbmlRole="parameter";
  Reaction10.sbmlRole="rate";
  k1.sbmlRole="parameter";
  k1.sbmlName="K1";
  T1.sbmlRole="parameter";
  Reaction11.sbmlRole="rate";
  v3.sbmlRole="parameter";
  v3.sbmlName="V3";
  parameter_0000073.sbmlRole="parameter";
  parameter_0000073.sbmlName="V4";
  Reaction12.sbmlRole="rate";
  v1.sbmlRole="parameter";
  v1.sbmlName="V1";
  parameter_0000072.sbmlRole="parameter";
  parameter_0000072.sbmlName="V2";
  Reaction16.sbmlRole="rate";
  s4.sbmlRole="parameter";
  s4.sbmlName="S4";
  Reaction18.sbmlRole="rate";
  s6.sbmlRole="parameter";
  s6.sbmlName="S6";
  Reaction19.sbmlRole="rate";
  s2.sbmlRole="parameter";
  s2.sbmlName="S2";
  Reaction20.sbmlRole="rate";
  Reaction21.sbmlRole="rate";
  Reaction23.sbmlRole="rate";
  Reaction24.sbmlRole="rate";
  Reaction25.sbmlRole="rate";
  Reaction26.sbmlRole="rate";
  Reaction27.sbmlRole="rate";
  Reaction28.sbmlRole="rate";
  D1.sbmlRole="parameter";
  L1.sbmlRole="parameter";
  Reaction29.sbmlRole="rate";
  D2.sbmlRole="parameter";
  L2.sbmlRole="parameter";
  Reaction30.sbmlRole="rate";
  D3.sbmlRole="parameter";
  L3.sbmlRole="parameter";
  Reaction31.sbmlRole="rate";
  D4.sbmlRole="parameter";
  L4.sbmlRole="parameter";
  Reaction32.sbmlRole="rate";
  D5.sbmlRole="parameter";
  L5.sbmlRole="parameter";
  Reaction33.sbmlRole="rate";
  D6.sbmlRole="parameter";
  L6.sbmlRole="parameter";
  Reaction34.sbmlRole="rate";
  D7.sbmlRole="parameter";
  L7.sbmlRole="parameter";
  Reaction35.sbmlRole="rate";
  D8.sbmlRole="parameter";
  L8.sbmlRole="parameter";
  Reaction36.sbmlRole="rate";
  D9.sbmlRole="parameter";
  L9.sbmlRole="parameter";
  Reaction37.sbmlRole="rate";
  D10.sbmlRole="parameter";
  L10.sbmlRole="parameter";
}