// 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
// votc: 148445cp 148445ornmat => 148445citmat 148445pi 
// votcR: 148445cpR 148445ornmat => 148445citmatR 148445pi 
// vcitmi: 148445citmat => 244266citims 
// vcitmiR: 148445citmatR => 244266citimsR 
// vcitic: 244266citims => 40395.6cit 
// vciticR: 244266citimsR => 40395.6citR 
// vass: 40395.6asp 40395.6atp 40395.6cit => 40395.6amp 40395.6as 40395.6ppi 
// vassR: 40395.6asp 40395.6atp 40395.6citR => 40395.6amp 40395.6asR 40395.6ppi 
// vasl: 40395.6as => 40395.6arg 40395.6fum 
// vaslR: 40395.6asR => 40395.6argR 40395.6fum 
// va: 40395.6arg => 40395.6orn 40395.6ure 
// vaR: 40395.6argR => 40395.6orn 40395.6ureR 
// vaQ: 40395.6argQ => 40395.6orn 40395.6ure 
// vaE: 509.861argE => 509.861ornE 509.861ureE 
// vaRE: 509.861argRE => 509.861ornE 509.861ureRE 
// vaQE: 509.861argQE => 509.861ornE 509.861ureE 
// vornci: 40395.6orn => 244266ornims 
// vornim: 244266ornims => 148445ornmat 
// vcpin: 509.861cppool => 509.861cpeE 
// vcpci: 40395.6cpe => 244266cpims 
// vcpim: 244266cpims => 148445cp 
// vcpciR: 40395.6cpeR => 244266cpimsR 
// vcpimR: 244266cpimsR => 148445cpR 
// vatpin: 40395.6atppool => 40395.6atp 
// vaspin: 40395.6asppool => 40395.6asp 
// vppisink: 40395.6ppi => XX 
// vfumsink: 40395.6fum => XX 
// vpisink: 148445pi => XX 
// vampsink: 40395.6amp => XX 
// vusink: 40395.6ure => XX 
// vusinkR: 40395.6ureR => XX 
// vcitCE: 40395.6cit => 509.861citE 
// vcitRCE: 40395.6citR => 509.861citRE 
// vatpCE: 40395.6atp => 509.861atpE 
// vaspCE: 40395.6asp => 509.861aspE 
// vampCE: 40395.6amp => 509.861ampE 
// vppiCE: 40395.6ppi => 509.861ppiE 
// vasCE: 40395.6as => 509.861asE 
// vasRCE: 40395.6asR => 509.861asRE 
// vfumCE: 40395.6fum => 509.861fumE 
// vargQCE: 40395.6argQ => 509.861argQE 
// vargCE: 40395.6arg => 509.861argE 
// vargRCE: 40395.6argR => 509.861argRE 
// vureCE: 40395.6ure => 509.861ureE 
// vureRCE: 40395.6ureR => 509.861ureRE 
// vornCE: 40395.6orn => 509.861ornE 
// vcpeCE: 40395.6cpe => 509.861cpeE 
// vcpeRCE: 40395.6cpeR => 509.861cpeRE 

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

  // variable definitions
  real compartment = 1 L;
  extern real denom1;
  real vel1FOR(time);
  real vel1REV(time);
  real vel1FORRc(time);
  real vel1REVRc(time);
  real vel20(time);
  real vel21(time);
  real vel20Rc(time);
  real vel21Rc(time);
  real vel19(time);
  real vel18(time);
  real vel18Rc(time);
  real vel19Rc(time);
  real denom2(time);
  real vel2FOR(time);
  real vel2REV(time);
  real vel2FORRc(time);
  real vel2REVRa(time);
  real denom3(time);
  real vel3FOR(time);
  real vel3REV(time);
  real vel3FORRa(time);
  real vel3REVRa(time);
  real denom4(time);
  real vel4FORa(time);
  real vel4FORRa(time);
  real vel4FORQ(time);
  real denom4E(time);
  real vel4FORE(time);
  real vel4FORRaE(time);
  real vel4FORQE(time);
  real vel14(time);
  real vel15(time);
  real vel22(time);
  real vel23(time);
  real vel8(time);
  real vel108(time);
  real vel109(time);
  real vel208(time);
  real vel209(time);
  real vel108R(time);
  real vel109R(time);
  real vel208R(time);
  real vel209R(time);
  real vel5(time);
  real vel9(time);
  real vel6(time);
  real vel7(time);
  real vel11(time);
  real vel10(time);
  real vel12(time);
  real vel12Ru(time);
  real vcitCE(time);
  real vcitCEREV(time);
  real vcitRCE(time);
  real vcitRCEREV(time);
  real vatpCE(time);
  real vatpCEREV(time);
  real vaspCE(time);
  real vaspCEREV(time);
  real vampCE(time);
  real vampCEREV(time);
  real vppiCE(time);
  real vppiCEREV(time);
  real vasCE(time);
  real vasCEREV(time);
  real vasRCE(time);
  real vasRCEREV(time);
  real vfumCE(time);
  real vfumCEREV(time);
  real vargQCE(time);
  real vargQCEREV(time);
  real vargCE(time);
  real vargCEREV(time);
  real vargRCE(time);
  real vargRCEREV(time);
  real vureCE(time);
  real vureCEREV(time);
  real vureRCE(time);
  real vureRCEREV(time);
  real vornCE(time);
  real vornCEREV(time);
  real vcpeCE(time);
  real vcpeCEREV(time);
  real vcpeRCE(time);
  real vcpeRCEREV(time);
  real sampleV = .002;
  real cellV = 3.868E-5;
  real mitoV = 1.08304E-5;
  real nucV = 3.0944E-6;
  real cytV = 2.47552E-5;
  real imsV = 4.09389E-6;
  real matV = 6.73651E-6;
  real eV = .00196132;
  real k1e1 = 1.7E7;
  real k2e1 = 63;
  real k3e1 = 2.1E6;
  real k4e1 = 1E3;
  real k5e1 = 3E3;
  real k6e1 = 9E4;
  real k7e1 = 2600;
  real k8e1 = 5E5;
  real k1e2 = 2.4E5;
  real k2e2 = 2.3;
  real k3e2 = 3.5E5;
  real k4e2 = 10;
  real k5e2 = 4.8E5;
  real k6e2 = 10;
  real k7e2 = 20;
  real k8e2 = 8.9E5;
  real k9e2 = 50;
  real k10e2 = 6.4E5;
  real k11e2 = 50;
  real k12e2 = 1.7E5;
  real k1e3 = 2.7E6;
  real k2e3 = 70;
  real k3e3 = 75;
  real k4e3 = 1.5E6;
  real k5e3 = 1100;
  real k6e3 = 7E5;
  real k1e4 = 1E7;
  real k2e4 = 54000;
  real k3e4 = 5300;
  real k4e4 = 3E4;
  real k5e4 = 1E7;
  real steadyStateUreaProdn = 7.80047E-12;
  real e10 = 1.49288E-5;
  real e20 = 6.25E-6;
  real e30 = 3.4375E-6;
  real e40 = 1.20984E-5;
  real e40E = 2.28177E-8;
  real katp = 1.57552E-4;
  real kasp = 3.15104E-4;
  real kcp = 2.65143E-7;
  real ppifinal = 1E-5;
  real fumfinal = 1E-5;
  real ampfinal = 5E-5;
  real pifinal = 5E-5;
  real urefinal = 1E-6;
  real kppi = .0315104;
  real kfum = .0315104;
  real kamp = .00630208;
  real kpi = .0231588;
  real kure = .315104;
  real Kiarg = 7E-4;
  real kornci = 1E4;
  real kornic = 60468.6;
  real kcitci = 1E4;
  real kcitic = 60468.6;
  real kornim = 110;
  real kornmi = 66.8489;
  real kcitim = 100;
  real kcitmi = 60.7717;
  real kcpei = 1E4;
  real kcpie = 60468.6;
  real kcpmi = .5;
  real kcpim = .822751;
  real kcitCE = 1E7;
  real kcitRCE = 1E7;
  real katpCE = 1E7;
  real kaspCE = 1E7;
  real kampCE = 1E7;
  real kppiCE = 1E7;
  real kasCE = 1E7;
  real kasRCE = 1E7;
  real kfumCE = 1E7;
  real kargQCE = 1E7;
  real kargCE = 1E7;
  real kargRCE = 1E7;
  real kureCE = 1E7;
  real kureRCE = 1E7;
  real kornCE = 1E7;
  real kcpeCE = 1E7;
  real kcpeRCE = 1E7;
  real kcitCEREV = 126217;
  real kcitRCEREV = 126217;
  real katpCEREV = 126217;
  real kaspCEREV = 126217;
  real kampCEREV = 126217;
  real kppiCEREV = 126217;
  real kasCEREV = 126217;
  real kasRCEREV = 126217;
  real kfumCEREV = 126217;
  real kargQCEREV = 126217;
  real kargCEREV = 126217;
  real kargRCEREV = 126217;
  real kureCEREV = 126217;
  real kureRCEREV = 126217;
  real kornCEREV = 126217;
  real kcpeCEREV = 126217;
  real kcpeRCEREV = 126217;
  real otcAdj = 1;
  real assAdj = 2;
  real aslAdj = 1;
  real argfac = 1;
  real fm = 1;
  real fmE = 1;
  real cp(time) M;
  real ornmat(time) M;
  real citmat(time) M;
  real pi(time) M;
  real cpR(time) M;
  real citmatR(time) M;
  real citims(time) M;
  real citimsR(time) M;
  real cit(time) M;
  real citR(time) M;
  real asp(time) M;
  real atp(time) M;
  real amp(time) M;
  real as(time) M;
  real ppi(time) M;
  real asR(time) M;
  real arg(time) M;
  real fum(time) M;
  real argR(time) M;
  real orn(time) M;
  real ure(time) M;
  real ureR(time) M;
  real argQ(time) M;
  real argE(time) M;
  real ornE(time) M;
  real ureE(time) M;
  real argRE(time) M;
  real ureRE(time) M;
  real argQE(time) M;
  real ornims(time) M;
  real cppool = .015 M;
  real cpeE(time) M;
  real cpe(time) M;
  real cpims(time) M;
  real cpeR(time) M;
  real cpimsR(time) M;
  real atppool = .002 M;
  real asppool = .001 M;
  real XX = 0 M;
  real citE(time) M;
  real citRE(time) M;
  real atpE(time) M;
  real aspE(time) M;
  real ampE(time) M;
  real ppiE(time) M;
  real asE(time) M;
  real asRE(time) M;
  real fumE(time) M;
  real cpeRE(time) M;
  real votc(time) katal;
  real votcR(time) katal;
  real vcitmi(time) katal;
  real vcitmiR(time) katal;
  real vcitic(time) katal;
  real vciticR(time) katal;
  real vass(time) katal;
  real vassR(time) katal;
  real vasl(time) katal;
  real vaslR(time) katal;
  real va(time) katal;
  real vaR(time) katal;
  real vaQ(time) katal;
  real vaE(time) katal;
  real vaRE(time) katal;
  real vaQE(time) katal;
  real vornci(time) katal;
  real vornim(time) katal;
  real vcpin(time) katal;
  real vcpci(time) katal;
  real vcpim(time) katal;
  real vcpciR(time) katal;
  real vcpimR(time) katal;
  real vatpin(time) katal;
  real vaspin(time) katal;
  real vppisink(time) katal;
  real vfumsink(time) katal;
  real vpisink(time) katal;
  real vampsink(time) katal;
  real vusink(time) katal;
  real vusinkR(time) katal;

  // equations
  vel1FOR = otcAdj*k1e1*k3e1*k5e1*k7e1*cp*ornmat*e10/denom1;
  vel1REV = otcAdj*k2e1*k4e1*k6e1*k8e1*citmat*3.141592653589793*e10/denom1;
  vel1FORRc = otcAdj*k1e1*k3e1*k5e1*k7e1*cpR*ornmat*e10/denom1;
  vel1REVRc = otcAdj*k2e1*k4e1*k6e1*k8e1*citmatR*3.141592653589793*e10/denom1;
  vel20 = kcitim*citims;
  vel21 = kcitmi*citmat;
  vel20Rc = kcitim*citimsR;
  vel21Rc = kcitmi*citmatR;
  vel19 = kcitic*citims;
  vel18 = kcitci*cit;
  vel18Rc = kcitci*citR;
  vel19Rc = kcitic*citimsR;
  denom2 = k2e2*k4e2*k9e2*k11e2*(k6e2+k7e2)+k1e2*k4e2*k6e2*k8e2*k11e2*(cit+citR)*ppi+k1e2*k4e2*k9e2*k11e2*(k6e2+k7e2)*(cit+citR)+k2e2*k5e2*k7e2*k9e2*k12e2*asp*(as+asR)+k2e2*k5e2*k7e2*k9e2*k11e2*asp+k1e2*k3e2*k6e2*k8e2*k11e2*(cit+citR)*atp*ppi+k1e2*k3e2*k9e2*k11e2*(k6e2+k7e2)*(cit+citR)*atp+k1e2*k4e2*k6e2*k8e2*k10e2*(cit+citR)*ppi*amp+k1e2*k5e2*k7e2*k9e2*k11e2*(cit+citR)*asp+k3e2*k5e2*k7e2*k9e2*k12e2*atp*asp*(as+asR)+k3e2*k5e2*k7e2*k9e2*k11e2*atp*asp+k2e2*k5e2*k7e2*k10e2*k12e2*asp*amp*(as+asR)+k1e2*k3e2*k5e2*(k7e2*k9e2+k7e2*k11e2+k9e2*k11e2)*(cit+citR)*atp*asp+k1e2*k3e2*k5e2*k8e2*k11e2*(cit+citR)*atp*asp*ppi+k2e2*k4e2*k6e2*k8e2*k11e2*ppi+k1e2*k3e2*k5e2*k7e2*k10e2*(cit+citR)*atp*asp*amp+k2e2*k4e2*k9e2*k12e2*(k6e2+k7e2)*(as+asR)+k1e2*k3e2*k6e2*k8e2*k10e2*(cit+citR)*atp*ppi*amp+k2e2*k4e2*k6e2*k8e2*k10e2*ppi*amp+k3e2*k5e2*k7e2*k10e2*k12e2*atp*asp*amp*(as+asR)+k2e2*k4e2*k6e2*k8e2*k12e2*ppi*(as+asR)+k3e2*k6e2*k8e2*k10e2*k12e2*atp*ppi*amp*(as+asR)+k2e2*k4e2*k10e2*k12e2*(k6e2+k7e2)*amp*(as+asR)+k2e2*k5e2*k8e2*k10e2*k12e2*asp*ppi*amp*(as+asR)+k8e2*k10e2*k12e2*(k2e2*k4e2+k2e2*k6e2+k4e2*k6e2)*ppi*amp*(as+asR)+k1e2*k3e2*k5e2*k8e2*k10e2*(cit+citR)*atp*asp*ppi*amp+k3e2*k5e2*k8e2*k10e2*k12e2*atp*asp*ppi*amp*(as+asR);
  vel2FOR = k1e2*k3e2*k5e2*k7e2*k9e2*k11e2*cit*atp*asp*assAdj*e20*(1+(arg+argR+argQ)/Kiarg)^(-1)/denom2;
  vel2REV = k2e2*k4e2*k6e2*k8e2*k10e2*k12e2*ppi*amp*as*assAdj*e20*(1+(arg+argR+argQ)/Kiarg)^(-1)/denom2;
  vel2FORRc = k1e2*k3e2*k5e2*k7e2*k9e2*k11e2*citR*atp*asp*assAdj*e20*(1+(arg+argR+argQ)/Kiarg)^(-1)/denom2;
  vel2REVRa = k2e2*k4e2*k6e2*k8e2*k10e2*k12e2*ppi*amp*asR*assAdj*e20*(1+(arg+argR+argQ)/Kiarg)^(-1)/denom2;
  denom3 = k5e3*(k2e3+k3e3)+(as+asR)*k1e3*(k3e3+k5e3)+fum*k2e3*k4e3+(arg+argR+fm*argQ)*k6e3*(k2e3+k3e3)+fum*(arg+argR+fm*argQ)*k4e3*k6e3+(as+asR+fm*argQ)*fum*k1e3*k4e3;
  vel3FOR = aslAdj*k1e3*k3e3*k5e3*as*e30/denom3;
  vel3REV = aslAdj*k2e3*k4e3*k6e3*fum*arg*e30/denom3;
  vel3FORRa = aslAdj*k1e3*k3e3*k5e3*asR*e30/denom3;
  vel3REVRa = aslAdj*k2e3*k4e3*k6e3*fum*argR*e30/denom3;
  denom4 = k4e4*(k2e4+k3e4)+k5e4*(k2e4+k3e4)*orn+k1e4*(k3e4+k4e4)*(arg+argR+fm*argQ);
  vel4FORa = argfac*k1e4*k3e4*k4e4*arg*e40/denom4;
  vel4FORRa = argfac*k1e4*k3e4*k4e4*argR*e40/denom4;
  vel4FORQ = argfac*fm*k1e4*k3e4*k4e4*argQ*e40/denom4;
  denom4E = k4e4*(k2e4+k3e4)+k5e4*(k2e4+k3e4)*ornE+k1e4*(k3e4+k4e4)*(argE+argRE+fmE*argQE);
  vel4FORE = argfac*k1e4*k3e4*k4e4*argE*e40E/denom4E;
  vel4FORRaE = argfac*k1e4*k3e4*k4e4*argRE*e40E/denom4E;
  vel4FORQE = argfac*fm*k1e4*k3e4*k4e4*argQE*e40E/denom4E;
  vel14 = kornci*orn;
  vel15 = kornic*ornims;
  vel22 = kornim*ornims;
  vel23 = kornmi*ornmat;
  vel8 = kcp*cppool;
  vel108 = kcpei*cpe;
  vel109 = kcpie*cpims;
  vel208 = kcpim*cpims;
  vel209 = kcpmi*cp;
  vel108R = kcpei*cpeR;
  vel109R = kcpie*cpimsR;
  vel208R = kcpim*cpimsR;
  vel209R = kcpmi*cpR;
  vel5 = katp*atppool;
  vel9 = kasp*asppool;
  vel6 = kppi*ppi;
  vel7 = kfum*fum;
  vel11 = kpi*3.141592653589793;
  vel10 = kamp*amp;
  vel12 = kure*ure;
  vel12Ru = kure*ureR;
  vcitCE = kcitCE*cit;
  vcitCEREV = kcitCEREV*citE;
  vcitRCE = kcitRCE*citR;
  vcitRCEREV = kcitRCEREV*citRE;
  vatpCE = katpCE*atp;
  vatpCEREV = katpCEREV*atpE;
  vaspCE = kaspCE*asp;
  vaspCEREV = kaspCEREV*aspE;
  vampCE = kampCE*amp;
  vampCEREV = kampCEREV*ampE;
  vppiCE = kppiCE*ppi;
  vppiCEREV = kppiCEREV*ppiE;
  vasCE = kasCE*as;
  vasCEREV = kasCEREV*asE;
  vasRCE = kasRCE*asR;
  vasRCEREV = kasRCEREV*asRE;
  vfumCE = kfumCE*fum;
  vfumCEREV = kfumCEREV*fumE;
  vargQCE = kargQCE*argQ;
  vargQCEREV = kargQCEREV*argQE;
  vargCE = kargCE*arg;
  vargCEREV = kargCEREV*argE;
  vargRCE = kargRCE*argR;
  vargRCEREV = kargRCEREV*argRE;
  vureCE = kureCE*ure;
  vureCEREV = kureCEREV*ureE;
  vureRCE = kureRCE*ureR;
  vureRCEREV = kureRCEREV*ureRE;
  vornCE = kornCE*orn;
  vornCEREV = kornCEREV*ornE;
  vcpeCE = kcpeCE*cpe;
  vcpeCEREV = kcpeCEREV*cpeE;
  vcpeRCE = kcpeRCE*cpeR;
  vcpeRCEREV = kcpeRCEREV*cpeRE;
  when (time=time.min) cp = 1E-5;
  cp:time = (-148445*votc + 148445*vcpim)/compartment;
  when (time=time.min) ornmat = 1E-6;
  ornmat:time = (-148445*votc + -148445*votcR + 148445*vornim)/compartment;
  when (time=time.min) citmat = 1E-6;
  citmat:time = (148445*votc + -148445*vcitmi)/compartment;
  when (time=time.min) pi = 1E-4;
  pi:time = (148445*votc + 148445*votcR + -148445*vpisink)/compartment;
  when (time=time.min) cpR = 0;
  cpR:time = (-148445*votcR + 148445*vcpimR)/compartment;
  when (time=time.min) citmatR = 0;
  citmatR:time = (148445*votcR + -148445*vcitmiR)/compartment;
  when (time=time.min) citims = 1E-6;
  citims:time = (244266*vcitmi + -244266*vcitic)/compartment;
  when (time=time.min) citimsR = 0;
  citimsR:time = (244266*vcitmiR + -244266*vciticR)/compartment;
  when (time=time.min) cit = 1E-6;
  cit:time = (40395.6*vcitic + -40395.6*vass + -40395.6*vcitCE)/compartment;
  when (time=time.min) citR = 0;
  citR:time = (40395.6*vciticR + -40395.6*vassR + -40395.6*vcitRCE)/compartment;
  when (time=time.min) asp = 1E-4;
  asp:time = (-40395.6*vass + -40395.6*vassR + 40395.6*vaspin + -40395.6*vaspCE)/compartment;
  when (time=time.min) atp = .0015;
  atp:time = (-40395.6*vass + -40395.6*vassR + 40395.6*vatpin + -40395.6*vatpCE)/compartment;
  when (time=time.min) amp = 1E-6;
  amp:time = (40395.6*vass + 40395.6*vassR + -40395.6*vampsink + -40395.6*vampCE)/compartment;
  when (time=time.min) as = 1E-6;
  as:time = (40395.6*vass + -40395.6*vasl + -40395.6*vasCE)/compartment;
  when (time=time.min) ppi = 1E-6;
  ppi:time = (40395.6*vass + 40395.6*vassR + -40395.6*vppisink + -40395.6*vppiCE)/compartment;
  when (time=time.min) asR = 0;
  asR:time = (40395.6*vassR + -40395.6*vaslR + -40395.6*vasRCE)/compartment;
  when (time=time.min) arg = 1E-6;
  arg:time = (40395.6*vasl + -40395.6*va + -40395.6*vargCE)/compartment;
  when (time=time.min) fum = 1E-6;
  fum:time = (40395.6*vasl + 40395.6*vaslR + -40395.6*vfumsink + -40395.6*vfumCE)/compartment;
  when (time=time.min) argR = 0;
  argR:time = (40395.6*vaslR + -40395.6*vaR + -40395.6*vargRCE)/compartment;
  when (time=time.min) orn = 2E-4;
  orn:time = (40395.6*va + 40395.6*vaR + 40395.6*vaQ + -40395.6*vornci + -40395.6*vornCE)/compartment;
  when (time=time.min) ure = 1E-6;
  ure:time = (40395.6*va + 40395.6*vaQ + -40395.6*vusink + -40395.6*vureCE)/compartment;
  when (time=time.min) ureR = 0;
  ureR:time = (40395.6*vaR + -40395.6*vusinkR + -40395.6*vureRCE)/compartment;
  when (time=time.min) argQ = .001;
  argQ:time = (-40395.6*vaQ + -40395.6*vargQCE)/compartment;
  when (time=time.min) argE = 1E-6;
  argE:time = (-509.861*vaE + 509.861*vargCE)/compartment;
  when (time=time.min) ornE = 2E-4;
  ornE:time = (509.861*vaE + 509.861*vaRE + 509.861*vaQE + 509.861*vornCE)/compartment;
  when (time=time.min) ureE = 1E-6;
  ureE:time = (509.861*vaE + 509.861*vaQE + 509.861*vureCE)/compartment;
  when (time=time.min) argRE = 0;
  argRE:time = (-509.861*vaRE + 509.861*vargRCE)/compartment;
  when (time=time.min) ureRE = 0;
  ureRE:time = (509.861*vaRE + 509.861*vureRCE)/compartment;
  when (time=time.min) argQE = .001;
  argQE:time = (-509.861*vaQE + 509.861*vargQCE)/compartment;
  when (time=time.min) ornims = 4.5E-4;
  ornims:time = (244266*vornci + -244266*vornim)/compartment;
  when (time=time.min) cpeE = 5E-4;
  cpeE:time = (509.861*vcpin + 509.861*vcpeCE)/compartment;
  when (time=time.min) cpe = 1E-5;
  cpe:time = (-40395.6*vcpci + -40395.6*vcpeCE)/compartment;
  when (time=time.min) cpims = 1E-5;
  cpims:time = (244266*vcpci + -244266*vcpim)/compartment;
  when (time=time.min) cpeR = 0;
  cpeR:time = (-40395.6*vcpciR + -40395.6*vcpeRCE)/compartment;
  when (time=time.min) cpimsR = 0;
  cpimsR:time = (244266*vcpciR + -244266*vcpimR)/compartment;
  when (time=time.min) citE = 1E-6;
  citE:time = (509.861*vcitCE)/compartment;
  when (time=time.min) citRE = 0;
  citRE:time = (509.861*vcitRCE)/compartment;
  when (time=time.min) atpE = .0015;
  atpE:time = (509.861*vatpCE)/compartment;
  when (time=time.min) aspE = 1E-4;
  aspE:time = (509.861*vaspCE)/compartment;
  when (time=time.min) ampE = 1E-6;
  ampE:time = (509.861*vampCE)/compartment;
  when (time=time.min) ppiE = 1E-6;
  ppiE:time = (509.861*vppiCE)/compartment;
  when (time=time.min) asE = 1E-6;
  asE:time = (509.861*vasCE)/compartment;
  when (time=time.min) asRE = 0;
  asRE:time = (509.861*vasRCE)/compartment;
  when (time=time.min) fumE = 1E-6;
  fumE:time = (509.861*vfumCE)/compartment;
  when (time=time.min) cpeRE = 2.26667E-6;
  cpeRE:time = (509.861*vcpeRCE)/compartment;
  votc = vel1FOR*matV-vel1REV*matV;
  votcR = vel1FORRc*matV-vel1REVRc*matV;
  vcitmi = vel21*matV-vel20*imsV;
  vcitmiR = vel21Rc*matV-vel20Rc*imsV;
  vcitic = vel19*imsV-vel18*cytV;
  vciticR = vel19Rc*imsV-vel18Rc*cytV;
  vass = vel2FOR*cytV-vel2REV*cytV;
  vassR = vel2FORRc*cytV-vel2REVRa*cytV;
  vasl = vel3FOR*cytV-vel3REV*cytV;
  vaslR = vel3FORRa*cytV-vel3REVRa*cytV;
  va = vel4FORa*cytV;
  vaR = vel4FORRa*cytV;
  vaQ = vel4FORQ*cytV;
  vaE = vel4FORE*eV;
  vaRE = vel4FORRaE*eV;
  vaQE = vel4FORQE*eV;
  vornci = vel14*cytV-vel15*imsV;
  vornim = vel22*imsV-vel23*matV;
  vcpin = vel8*eV;
  vcpci = vel108*cytV-vel109*imsV;
  vcpim = vel208*imsV-vel209*matV;
  vcpciR = vel108R*cytV-vel109R*imsV;
  vcpimR = vel208R*imsV-vel209R*matV;
  vatpin = vel5*cytV;
  vaspin = vel9*cytV;
  vppisink = vel6*cytV;
  vfumsink = vel7*cytV;
  vpisink = vel11*matV;
  vampsink = vel10*cytV;
  vusink = vel12*cytV;
  vusinkR = 0*vel12Ru*cytV;

  // variable properties
  compartment.sbmlRole="compartment";
  denom1.sbmlRole="parameter";
  vel1FOR.sbmlRole="parameter";
  vel1REV.sbmlRole="parameter";
  vel1FORRc.sbmlRole="parameter";
  vel1REVRc.sbmlRole="parameter";
  vel20.sbmlRole="parameter";
  vel21.sbmlRole="parameter";
  vel20Rc.sbmlRole="parameter";
  vel21Rc.sbmlRole="parameter";
  vel19.sbmlRole="parameter";
  vel18.sbmlRole="parameter";
  vel18Rc.sbmlRole="parameter";
  vel19Rc.sbmlRole="parameter";
  denom2.sbmlRole="parameter";
  vel2FOR.sbmlRole="parameter";
  vel2REV.sbmlRole="parameter";
  vel2FORRc.sbmlRole="parameter";
  vel2REVRa.sbmlRole="parameter";
  denom3.sbmlRole="parameter";
  vel3FOR.sbmlRole="parameter";
  vel3REV.sbmlRole="parameter";
  vel3FORRa.sbmlRole="parameter";
  vel3REVRa.sbmlRole="parameter";
  denom4.sbmlRole="parameter";
  vel4FORa.sbmlRole="parameter";
  vel4FORRa.sbmlRole="parameter";
  vel4FORQ.sbmlRole="parameter";
  denom4E.sbmlRole="parameter";
  vel4FORE.sbmlRole="parameter";
  vel4FORRaE.sbmlRole="parameter";
  vel4FORQE.sbmlRole="parameter";
  vel14.sbmlRole="parameter";
  vel15.sbmlRole="parameter";
  vel22.sbmlRole="parameter";
  vel23.sbmlRole="parameter";
  vel8.sbmlRole="parameter";
  vel108.sbmlRole="parameter";
  vel109.sbmlRole="parameter";
  vel208.sbmlRole="parameter";
  vel209.sbmlRole="parameter";
  vel108R.sbmlRole="parameter";
  vel109R.sbmlRole="parameter";
  vel208R.sbmlRole="parameter";
  vel209R.sbmlRole="parameter";
  vel5.sbmlRole="parameter";
  vel9.sbmlRole="parameter";
  vel6.sbmlRole="parameter";
  vel7.sbmlRole="parameter";
  vel11.sbmlRole="parameter";
  vel10.sbmlRole="parameter";
  vel12.sbmlRole="parameter";
  vel12Ru.sbmlRole="parameter";
  vcitCE.sbmlRole="parameter";
  vcitCEREV.sbmlRole="parameter";
  vcitRCE.sbmlRole="parameter";
  vcitRCEREV.sbmlRole="parameter";
  vatpCE.sbmlRole="parameter";
  vatpCEREV.sbmlRole="parameter";
  vaspCE.sbmlRole="parameter";
  vaspCEREV.sbmlRole="parameter";
  vampCE.sbmlRole="parameter";
  vampCEREV.sbmlRole="parameter";
  vppiCE.sbmlRole="parameter";
  vppiCEREV.sbmlRole="parameter";
  vasCE.sbmlRole="parameter";
  vasCEREV.sbmlRole="parameter";
  vasRCE.sbmlRole="parameter";
  vasRCEREV.sbmlRole="parameter";
  vfumCE.sbmlRole="parameter";
  vfumCEREV.sbmlRole="parameter";
  vargQCE.sbmlRole="parameter";
  vargQCEREV.sbmlRole="parameter";
  vargCE.sbmlRole="parameter";
  vargCEREV.sbmlRole="parameter";
  vargRCE.sbmlRole="parameter";
  vargRCEREV.sbmlRole="parameter";
  vureCE.sbmlRole="parameter";
  vureCEREV.sbmlRole="parameter";
  vureRCE.sbmlRole="parameter";
  vureRCEREV.sbmlRole="parameter";
  vornCE.sbmlRole="parameter";
  vornCEREV.sbmlRole="parameter";
  vcpeCE.sbmlRole="parameter";
  vcpeCEREV.sbmlRole="parameter";
  vcpeRCE.sbmlRole="parameter";
  vcpeRCEREV.sbmlRole="parameter";
  sampleV.sbmlRole="parameter";
  cellV.sbmlRole="parameter";
  mitoV.sbmlRole="parameter";
  nucV.sbmlRole="parameter";
  cytV.sbmlRole="parameter";
  imsV.sbmlRole="parameter";
  matV.sbmlRole="parameter";
  eV.sbmlRole="parameter";
  k1e1.sbmlRole="parameter";
  k2e1.sbmlRole="parameter";
  k3e1.sbmlRole="parameter";
  k4e1.sbmlRole="parameter";
  k5e1.sbmlRole="parameter";
  k6e1.sbmlRole="parameter";
  k7e1.sbmlRole="parameter";
  k8e1.sbmlRole="parameter";
  k1e2.sbmlRole="parameter";
  k2e2.sbmlRole="parameter";
  k3e2.sbmlRole="parameter";
  k4e2.sbmlRole="parameter";
  k5e2.sbmlRole="parameter";
  k6e2.sbmlRole="parameter";
  k7e2.sbmlRole="parameter";
  k8e2.sbmlRole="parameter";
  k9e2.sbmlRole="parameter";
  k10e2.sbmlRole="parameter";
  k11e2.sbmlRole="parameter";
  k12e2.sbmlRole="parameter";
  k1e3.sbmlRole="parameter";
  k2e3.sbmlRole="parameter";
  k3e3.sbmlRole="parameter";
  k4e3.sbmlRole="parameter";
  k5e3.sbmlRole="parameter";
  k6e3.sbmlRole="parameter";
  k1e4.sbmlRole="parameter";
  k2e4.sbmlRole="parameter";
  k3e4.sbmlRole="parameter";
  k4e4.sbmlRole="parameter";
  k5e4.sbmlRole="parameter";
  steadyStateUreaProdn.sbmlRole="parameter";
  e10.sbmlRole="parameter";
  e20.sbmlRole="parameter";
  e30.sbmlRole="parameter";
  e40.sbmlRole="parameter";
  e40E.sbmlRole="parameter";
  katp.sbmlRole="parameter";
  kasp.sbmlRole="parameter";
  kcp.sbmlRole="parameter";
  ppifinal.sbmlRole="parameter";
  fumfinal.sbmlRole="parameter";
  ampfinal.sbmlRole="parameter";
  pifinal.sbmlRole="parameter";
  urefinal.sbmlRole="parameter";
  kppi.sbmlRole="parameter";
  kfum.sbmlRole="parameter";
  kamp.sbmlRole="parameter";
  kpi.sbmlRole="parameter";
  kure.sbmlRole="parameter";
  Kiarg.sbmlRole="parameter";
  kornci.sbmlRole="parameter";
  kornic.sbmlRole="parameter";
  kcitci.sbmlRole="parameter";
  kcitic.sbmlRole="parameter";
  kornim.sbmlRole="parameter";
  kornmi.sbmlRole="parameter";
  kcitim.sbmlRole="parameter";
  kcitmi.sbmlRole="parameter";
  kcpei.sbmlRole="parameter";
  kcpie.sbmlRole="parameter";
  kcpmi.sbmlRole="parameter";
  kcpim.sbmlRole="parameter";
  kcitCE.sbmlRole="parameter";
  kcitRCE.sbmlRole="parameter";
  katpCE.sbmlRole="parameter";
  kaspCE.sbmlRole="parameter";
  kampCE.sbmlRole="parameter";
  kppiCE.sbmlRole="parameter";
  kasCE.sbmlRole="parameter";
  kasRCE.sbmlRole="parameter";
  kfumCE.sbmlRole="parameter";
  kargQCE.sbmlRole="parameter";
  kargCE.sbmlRole="parameter";
  kargRCE.sbmlRole="parameter";
  kureCE.sbmlRole="parameter";
  kureRCE.sbmlRole="parameter";
  kornCE.sbmlRole="parameter";
  kcpeCE.sbmlRole="parameter";
  kcpeRCE.sbmlRole="parameter";
  kcitCEREV.sbmlRole="parameter";
  kcitRCEREV.sbmlRole="parameter";
  katpCEREV.sbmlRole="parameter";
  kaspCEREV.sbmlRole="parameter";
  kampCEREV.sbmlRole="parameter";
  kppiCEREV.sbmlRole="parameter";
  kasCEREV.sbmlRole="parameter";
  kasRCEREV.sbmlRole="parameter";
  kfumCEREV.sbmlRole="parameter";
  kargQCEREV.sbmlRole="parameter";
  kargCEREV.sbmlRole="parameter";
  kargRCEREV.sbmlRole="parameter";
  kureCEREV.sbmlRole="parameter";
  kureRCEREV.sbmlRole="parameter";
  kornCEREV.sbmlRole="parameter";
  kcpeCEREV.sbmlRole="parameter";
  kcpeRCEREV.sbmlRole="parameter";
  otcAdj.sbmlRole="parameter";
  assAdj.sbmlRole="parameter";
  aslAdj.sbmlRole="parameter";
  argfac.sbmlRole="parameter";
  fm.sbmlRole="parameter";
  fmE.sbmlRole="parameter";
  cp.sbmlRole="species";
  cp.sbmlCompartment="compartment";
  ornmat.sbmlRole="species";
  ornmat.sbmlCompartment="compartment";
  citmat.sbmlRole="species";
  citmat.sbmlCompartment="compartment";
  pi.sbmlRole="species";
  pi.sbmlCompartment="compartment";
  cpR.sbmlRole="species";
  cpR.sbmlCompartment="compartment";
  citmatR.sbmlRole="species";
  citmatR.sbmlCompartment="compartment";
  citims.sbmlRole="species";
  citims.sbmlCompartment="compartment";
  citimsR.sbmlRole="species";
  citimsR.sbmlCompartment="compartment";
  cit.sbmlRole="species";
  cit.sbmlCompartment="compartment";
  citR.sbmlRole="species";
  citR.sbmlCompartment="compartment";
  asp.sbmlRole="species";
  asp.sbmlCompartment="compartment";
  atp.sbmlRole="species";
  atp.sbmlCompartment="compartment";
  amp.sbmlRole="species";
  amp.sbmlCompartment="compartment";
  as.sbmlRole="species";
  as.sbmlCompartment="compartment";
  ppi.sbmlRole="species";
  ppi.sbmlCompartment="compartment";
  asR.sbmlRole="species";
  asR.sbmlCompartment="compartment";
  arg.sbmlRole="species";
  arg.sbmlCompartment="compartment";
  fum.sbmlRole="species";
  fum.sbmlCompartment="compartment";
  argR.sbmlRole="species";
  argR.sbmlCompartment="compartment";
  orn.sbmlRole="species";
  orn.sbmlCompartment="compartment";
  ure.sbmlRole="species";
  ure.sbmlCompartment="compartment";
  ureR.sbmlRole="species";
  ureR.sbmlCompartment="compartment";
  argQ.sbmlRole="species";
  argQ.sbmlCompartment="compartment";
  argE.sbmlRole="species";
  argE.sbmlCompartment="compartment";
  ornE.sbmlRole="species";
  ornE.sbmlCompartment="compartment";
  ureE.sbmlRole="species";
  ureE.sbmlCompartment="compartment";
  argRE.sbmlRole="species";
  argRE.sbmlCompartment="compartment";
  ureRE.sbmlRole="species";
  ureRE.sbmlCompartment="compartment";
  argQE.sbmlRole="species";
  argQE.sbmlCompartment="compartment";
  ornims.sbmlRole="species";
  ornims.sbmlCompartment="compartment";
  cppool.sbmlRole="species";
  cppool.sbmlCompartment="compartment";
  cpeE.sbmlRole="species";
  cpeE.sbmlCompartment="compartment";
  cpe.sbmlRole="species";
  cpe.sbmlCompartment="compartment";
  cpims.sbmlRole="species";
  cpims.sbmlCompartment="compartment";
  cpeR.sbmlRole="species";
  cpeR.sbmlCompartment="compartment";
  cpimsR.sbmlRole="species";
  cpimsR.sbmlCompartment="compartment";
  atppool.sbmlRole="species";
  atppool.sbmlCompartment="compartment";
  asppool.sbmlRole="species";
  asppool.sbmlCompartment="compartment";
  XX.sbmlRole="species";
  XX.sbmlCompartment="compartment";
  citE.sbmlRole="species";
  citE.sbmlCompartment="compartment";
  citRE.sbmlRole="species";
  citRE.sbmlCompartment="compartment";
  atpE.sbmlRole="species";
  atpE.sbmlCompartment="compartment";
  aspE.sbmlRole="species";
  aspE.sbmlCompartment="compartment";
  ampE.sbmlRole="species";
  ampE.sbmlCompartment="compartment";
  ppiE.sbmlRole="species";
  ppiE.sbmlCompartment="compartment";
  asE.sbmlRole="species";
  asE.sbmlCompartment="compartment";
  asRE.sbmlRole="species";
  asRE.sbmlCompartment="compartment";
  fumE.sbmlRole="species";
  fumE.sbmlCompartment="compartment";
  cpeRE.sbmlRole="species";
  cpeRE.sbmlCompartment="compartment";
  votc.sbmlRole="rate";
  votcR.sbmlRole="rate";
  vcitmi.sbmlRole="rate";
  vcitmiR.sbmlRole="rate";
  vcitic.sbmlRole="rate";
  vciticR.sbmlRole="rate";
  vass.sbmlRole="rate";
  vassR.sbmlRole="rate";
  vasl.sbmlRole="rate";
  vaslR.sbmlRole="rate";
  va.sbmlRole="rate";
  vaR.sbmlRole="rate";
  vaQ.sbmlRole="rate";
  vaE.sbmlRole="rate";
  vaRE.sbmlRole="rate";
  vaQE.sbmlRole="rate";
  vornci.sbmlRole="rate";
  vornim.sbmlRole="rate";
  vcpin.sbmlRole="rate";
  vcpci.sbmlRole="rate";
  vcpim.sbmlRole="rate";
  vcpciR.sbmlRole="rate";
  vcpimR.sbmlRole="rate";
  vatpin.sbmlRole="rate";
  vaspin.sbmlRole="rate";
  vppisink.sbmlRole="rate";
  vfumsink.sbmlRole="rate";
  vpisink.sbmlRole="rate";
  vampsink.sbmlRole="rate";
  vusink.sbmlRole="rate";
  vusinkR.sbmlRole="rate";
}