// 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
// v1: s => Clb2 
// v2: Clb2 => p 
// v3: s => Cdc20 
// v4: Cdc20 => p 
// v5: s => Cdh1 
// v6: Cdh1 => p 
// v7: s => securinT securin 
// v8: securinT securin => p 
// v9: s => separaseT separase 
// v10: separaseT separase => p 
// v11: securin separase => securinseparase 
// v12: securinseparase => securin separase 
// v13: s => Net1 
// v14: Net1 => p 
// v15: Net1 => Net1Cdc14 
// v16: Net1Cdc14 => Net1 
// v17: Net1Cdc14 => p 
// v18: s => PoloT 
// v19: PoloT => p 
// v20: s => Polo 
// v21: Polo => p 
// v22: Polo => p 
// v23: s => Tem1 
// v24: Tem1 => p 
// v25: s => Cdc15 
// v26: Cdc15 => p 
// v27: s => MEN 
// v28: MEN => p 
// v29: MEN => p 
// v30: MEN => p 

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

  // variable definitions
  real compartment = 1 L;
  real PP2A(time);
  real Cdk(time);
  real Net1P(time);
  real Cdc14(time);
  real ksclb2 = .03;
  real kdclb2 = .03;
  real kadclb2 = .2;
  real kaadclb2 = 2;
  real kscdc20 = .015;
  real kdcdc20 = .05;
  real kadcdc20 = 2;
  real kdcdh = .01;
  real kadcdh = 1;
  real kapcdh = 1;
  real Jcdh = .0015;
  real Cdh1T = 1;
  real kssecurin = .03;
  real kdsecurin = .05;
  real kadsecurin = 2;
  real ksseparase = .001;
  real kdseparase = .004;
  real lasecurin = 500;
  real ldsecurin = 1;
  real kad = .1;
  real kd = .45;
  real kp = .4;
  real kap = 2;
  real Jnet = .2;
  real lanet = 500;
  real ldnet = 1;
  real Net1T = 1;
  real Cdc14T = .5;
  real kspolo = .01;
  real kdpolo = .01;
  real kadpolo = .25;
  real kapolo = 0;
  real kaapolo = .5;
  real kipolo = .1;
  real Jpolo = .25;
  real katem = 0;
  real kaatem = .5;
  real kitem = .1;
  real kaitem = 1;
  real Jtem1 = .005;
  real Tem1T = 1;
  real kacdc15 = .02;
  real kaacdc15 = .5;
  real kicdc15 = 0;
  real kaicdc15 = .2;
  real Jcdc15 = .2;
  real lamen = 10;
  real ldmen = .1;
  real Cdc15T = 1;
  real kpp = .1;
  real ki = 20;
  real PP2AT = 1;
  real Inh = 0;
  real s = 1 M;
  real Clb2(time) M;
  real p = 1 M;
  real Cdc20(time) M;
  real Cdh1(time) M;
  real securinT(time) M;
  real securin(time) M;
  real separaseT(time) M;
  real separase(time) M;
  real securinseparase(time) M;
  real Net1(time) M;
  real Net1Cdc14(time) M;
  real PoloT(time) M;
  real Polo(time) M;
  real Tem1(time) M;
  real Cdc15(time) M;
  real MEN(time) M;
  real v1(time) katal;
  real v2(time) katal;
  real v3(time) katal;
  real v4(time) katal;
  real v5(time) katal;
  real v6(time) katal;
  real v7(time) katal;
  real v8(time) katal;
  real v9(time) katal;
  real v10(time) katal;
  real v11(time) katal;
  real v12(time) katal;
  real v13(time) katal;
  real v14(time) katal;
  real v15(time) katal;
  real v16(time) katal;
  real v17(time) katal;
  real v18(time) katal;
  real v19(time) katal;
  real v20(time) katal;
  real v21(time) katal;
  real v22(time) katal;
  real v23(time) katal;
  real v24(time) katal;
  real v25(time) katal;
  real v26(time) katal;
  real v27(time) katal;
  real v28(time) katal;
  real v29(time) katal;
  real v30(time) katal;

  // equations
  PP2A = (1+kpp*ki*separase)/(1+ki*separase)*PP2AT;
  Cdk = Clb2/(1+Inh);
  Net1P = Net1T-Net1-Net1Cdc14;
  Cdc14 = Cdc14T-Net1Cdc14;
  when (time=time.min) Clb2 = .997;
  Clb2:time = (v1 + -1*v2)/compartment;
  when (time=time.min) Cdc20 = 0;
  Cdc20:time = (v3 + -1*v4)/compartment;
  when (time=time.min) Cdh1 = 3.8E-5;
  Cdh1:time = (v5 + -1*v6)/compartment;
  when (time=time.min) securinT = .6;
  securinT:time = (v7 + -1*v8)/compartment;
  when (time=time.min) securin = .352;
  securin:time = (v7 + -1*v8 + -1*v11 + v12)/compartment;
  when (time=time.min) separaseT = .25;
  separaseT:time = (v9 + -1*v10)/compartment;
  when (time=time.min) separase = .002;
  separase:time = (v9 + -1*v10 + -1*v11 + v12)/compartment;
  when (time=time.min) securinseparase = .248;
  securinseparase:time = (v11 + -1*v12)/compartment;
  when (time=time.min) Net1 = .098;
  Net1:time = (v13 + -1*v14 + -1*v15 + v16)/compartment;
  when (time=time.min) Net1Cdc14 = .485;
  Net1Cdc14:time = (v15 + -1*v16 + -1*v17)/compartment;
  when (time=time.min) PoloT = .99;
  PoloT:time = (v18 + -1*v19)/compartment;
  when (time=time.min) Polo = .945;
  Polo:time = (v20 + -1*v21 + -1*v22)/compartment;
  when (time=time.min) Tem1 = .00389;
  Tem1:time = (v23 + -1*v24)/compartment;
  when (time=time.min) Cdc15 = .0257;
  Cdc15:time = (v25 + -1*v26)/compartment;
  when (time=time.min) MEN = 8E-6;
  MEN:time = (v27 + -1*v28 + -1*v29 + -1*v30)/compartment;
  v1 = ksclb2;
  v2 = (kdclb2+kadclb2*Cdc20+kaadclb2*Cdh1)*Clb2;
  v3 = kscdc20;
  v4 = (kdcdc20+kadcdc20*Cdh1)*Cdc20;
  v5 = (kdcdh+kadcdh*Cdc14)*((Cdh1T-Cdh1)/(Jcdh+Cdh1T-Cdh1));
  v6 = kapcdh*Cdk*(Cdh1/(Jcdh+Cdh1));
  v7 = kssecurin;
  v8 = (kdsecurin+kadsecurin*Cdc20)*securinT;
  v9 = ksseparase;
  v10 = kdseparase*separaseT;
  v11 = lasecurin*securin*separase;
  v12 = (ldsecurin+kdseparase+kdsecurin+kadsecurin*Cdc20)*securinseparase;
  v13 = (kad*Cdc14+kd*PP2A)*(Net1P/(Jnet+Net1P));
  v14 = (kp*Cdk+kap*MEN)*(Net1/(Jnet+Net1+Net1Cdc14));
  v15 = lanet*Net1*Cdc14;
  v16 = ldnet*Net1Cdc14;
  v17 = (kp*Cdk+kap*MEN)*(Net1Cdc14/(Jnet+Net1+Net1Cdc14));
  v18 = kspolo;
  v19 = (kdpolo+kadpolo*Cdh1)*PoloT;
  v20 = (kapolo+kaapolo*Cdk)*((PoloT-Polo)/(Jpolo+PoloT-Polo));
  v21 = kipolo*(Polo/(Jpolo+Polo));
  v22 = (kdpolo+kadpolo*Cdh1)*Polo;
  v23 = (katem+kaatem*Polo)*((Tem1T-Tem1)/(Jtem1+Tem1T-Tem1));
  v24 = (kitem+kaitem*PP2A)*(Tem1/(Jtem1+Tem1));
  v25 = (kacdc15+kaacdc15*Cdc14)*((Cdc15T-Cdc15)/(Jcdc15+Cdc15T-Cdc15));
  v26 = (kicdc15+kaicdc15*Cdk)*(Cdc15/(Jcdc15+Cdc15));
  v27 = lamen*(Tem1-MEN)*(Cdc15-MEN);
  v28 = ldmen*MEN;
  v29 = (kitem+kaitem*PP2A)*(MEN/(Jtem1+Tem1));
  v30 = (kicdc15+kaicdc15*Cdk)*(MEN/(Jcdc15+Cdc15));

  // variable properties
  compartment.sbmlRole="compartment";
  PP2A.sbmlRole="parameter";
  Cdk.sbmlRole="parameter";
  Net1P.sbmlRole="parameter";
  Cdc14.sbmlRole="parameter";
  ksclb2.sbmlRole="parameter";
  kdclb2.sbmlRole="parameter";
  kadclb2.sbmlRole="parameter";
  kaadclb2.sbmlRole="parameter";
  kscdc20.sbmlRole="parameter";
  kdcdc20.sbmlRole="parameter";
  kadcdc20.sbmlRole="parameter";
  kdcdh.sbmlRole="parameter";
  kadcdh.sbmlRole="parameter";
  kapcdh.sbmlRole="parameter";
  Jcdh.sbmlRole="parameter";
  Cdh1T.sbmlRole="parameter";
  kssecurin.sbmlRole="parameter";
  kdsecurin.sbmlRole="parameter";
  kadsecurin.sbmlRole="parameter";
  ksseparase.sbmlRole="parameter";
  kdseparase.sbmlRole="parameter";
  lasecurin.sbmlRole="parameter";
  ldsecurin.sbmlRole="parameter";
  kad.sbmlRole="parameter";
  kd.sbmlRole="parameter";
  kp.sbmlRole="parameter";
  kap.sbmlRole="parameter";
  Jnet.sbmlRole="parameter";
  lanet.sbmlRole="parameter";
  ldnet.sbmlRole="parameter";
  Net1T.sbmlRole="parameter";
  Cdc14T.sbmlRole="parameter";
  kspolo.sbmlRole="parameter";
  kdpolo.sbmlRole="parameter";
  kadpolo.sbmlRole="parameter";
  kapolo.sbmlRole="parameter";
  kaapolo.sbmlRole="parameter";
  kipolo.sbmlRole="parameter";
  Jpolo.sbmlRole="parameter";
  katem.sbmlRole="parameter";
  kaatem.sbmlRole="parameter";
  kitem.sbmlRole="parameter";
  kaitem.sbmlRole="parameter";
  Jtem1.sbmlRole="parameter";
  Tem1T.sbmlRole="parameter";
  kacdc15.sbmlRole="parameter";
  kaacdc15.sbmlRole="parameter";
  kicdc15.sbmlRole="parameter";
  kaicdc15.sbmlRole="parameter";
  Jcdc15.sbmlRole="parameter";
  lamen.sbmlRole="parameter";
  ldmen.sbmlRole="parameter";
  Cdc15T.sbmlRole="parameter";
  kpp.sbmlRole="parameter";
  ki.sbmlRole="parameter";
  PP2AT.sbmlRole="parameter";
  Inh.sbmlRole="parameter";
  s.sbmlRole="species";
  s.sbmlCompartment="compartment";
  Clb2.sbmlRole="species";
  Clb2.sbmlCompartment="compartment";
  p.sbmlRole="species";
  p.sbmlCompartment="compartment";
  Cdc20.sbmlRole="species";
  Cdc20.sbmlCompartment="compartment";
  Cdh1.sbmlRole="species";
  Cdh1.sbmlCompartment="compartment";
  securinT.sbmlRole="species";
  securinT.sbmlCompartment="compartment";
  securin.sbmlRole="species";
  securin.sbmlCompartment="compartment";
  separaseT.sbmlRole="species";
  separaseT.sbmlCompartment="compartment";
  separase.sbmlRole="species";
  separase.sbmlCompartment="compartment";
  securinseparase.sbmlRole="species";
  securinseparase.sbmlCompartment="compartment";
  Net1.sbmlRole="species";
  Net1.sbmlCompartment="compartment";
  Net1Cdc14.sbmlRole="species";
  Net1Cdc14.sbmlCompartment="compartment";
  PoloT.sbmlRole="species";
  PoloT.sbmlCompartment="compartment";
  Polo.sbmlRole="species";
  Polo.sbmlCompartment="compartment";
  Tem1.sbmlRole="species";
  Tem1.sbmlCompartment="compartment";
  Cdc15.sbmlRole="species";
  Cdc15.sbmlCompartment="compartment";
  MEN.sbmlRole="species";
  MEN.sbmlCompartment="compartment";
  v1.sbmlRole="rate";
  v2.sbmlRole="rate";
  v3.sbmlRole="rate";
  v4.sbmlRole="rate";
  v5.sbmlRole="rate";
  v6.sbmlRole="rate";
  v7.sbmlRole="rate";
  v8.sbmlRole="rate";
  v9.sbmlRole="rate";
  v10.sbmlRole="rate";
  v11.sbmlRole="rate";
  v12.sbmlRole="rate";
  v13.sbmlRole="rate";
  v14.sbmlRole="rate";
  v15.sbmlRole="rate";
  v16.sbmlRole="rate";
  v17.sbmlRole="rate";
  v18.sbmlRole="rate";
  v19.sbmlRole="rate";
  v20.sbmlRole="rate";
  v21.sbmlRole="rate";
  v22.sbmlRole="rate";
  v23.sbmlRole="rate";
  v24.sbmlRole="rate";
  v25.sbmlRole="rate";
  v26.sbmlRole="rate";
  v27.sbmlRole="rate";
  v28.sbmlRole="rate";
  v29.sbmlRole="rate";
  v30.sbmlRole="rate";
}