// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;
unit item=scalar;
unit time = 3600 second;
unit substance = mole;
unit volume = litre;
unit area = metre^2;
unit length = metre;

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

// SBML reactions
// Photoreceptor_activation: Pfr <=> Pr 
// Photoreceptor_inactivation: Pr <=> Pfr 
// Transducer_activation: Xi <=> Xa 
// Transducer_inactivation: Xa <=> Xi 
// preS_formation: prepreS <=> preS 
// S_generation: preS <=> S 
// Glucose_sensor_inactivation: Ya Gluc <=> Yi 
// S_formation: <=> S 
// V_formation: <=> V 
// S_degradation: S 
// V_degradation: V 
// Photoreceptor_decay: Pr <=> Pi 

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

  // variable definitions
  real compartment = 1 volume;
  private real Pfr.amt(time) substance;
  real Pfr(time) substance/volume;
  real Pfr.init substance;
  private real Pr.amt(time) substance;
  real Pr(time) substance/volume;
  real Pr.init substance;
  private real Xi.amt(time) substance;
  real Xi(time) substance/volume;
  real Xi.init substance;
  private real Xa.amt(time) substance;
  real Xa(time) substance/volume;
  real Xa.init substance;
  private real prepreS.amt(time) substance;
  real prepreS(time) substance/volume;
  real prepreS.init substance;
  private real preS.amt(time) substance;
  real preS(time) substance/volume;
  real preS.init substance;
  private real Ya.amt(time) substance;
  real Ya(time) substance/volume;
  real Ya.init substance;
  private real S.amt(time) substance;
  real S(time) substance/volume;
  real S.init substance;
  private real Gluc.amt(time) substance;
  real Gluc substance/volume;
  real Gluc.init substance;
  private real Yi.amt(time) substance;
  real Yi(time) substance/volume;
  real Yi.init substance;
  private real V.amt(time) substance;
  real V(time) substance/volume;
  real V.init substance;
  private real Pi.amt(time) substance;
  real Pi(time) substance/volume;
  real Pi.init substance;
  real Photoreceptor_activation.rate(time) substance/time;
  real IfrSfrPfr = .1;
  real Photoreceptor_inactivation.rate(time) substance/time;
  real IrSrPr = 0;
  real Transducer_activation.rate(time) substance/time;
  real kia = .1;
  real Transducer_inactivation.rate(time) substance/time;
  real kai = .8;
  real preS_formation.rate(time) substance/time;
  real kx = .2;
  real S_generation.rate(time) substance/time;
  real ky = 1;
  real Glucose_sensor_inactivation.rate(time) substance/time;
  real kG = .1;
  real S_formation.rate(time) substance/time;
  real alpha1 = 30;
  real V_formation.rate(time) substance/time;
  real alpha2 = 50;
  real S_degradation.rate(time) substance/time;
  real V_degradation.rate(time) substance/time;
  real Photoreceptor_decay.rate(time) substance/time;
  real kd = .1;

  // equations
  when (time=time.min) Pfr.amt = Pfr.init;
  Pfr.amt:time = -1*Photoreceptor_activation.rate + Photoreceptor_inactivation.rate;
  Pfr = Pfr.amt/compartment;
  Pfr.init = 10;
  when (time=time.min) Pr.amt = Pr.init;
  Pr.amt:time = Photoreceptor_activation.rate + -1*Photoreceptor_inactivation.rate + -1*Photoreceptor_decay.rate;
  Pr = Pr.amt/compartment;
  Pr.init = 0;
  when (time=time.min) Xi.amt = Xi.init;
  Xi.amt:time = -1*Transducer_activation.rate + Transducer_inactivation.rate;
  Xi = Xi.amt/compartment;
  Xi.init = 6;
  when (time=time.min) Xa.amt = Xa.init;
  Xa.amt:time = Transducer_activation.rate + -1*Transducer_inactivation.rate;
  Xa = Xa.amt/compartment;
  Xa.init = 0;
  when (time=time.min) prepreS.amt = prepreS.init;
  prepreS.amt:time = -1*preS_formation.rate;
  prepreS = prepreS.amt/compartment;
  prepreS.init = 200;
  when (time=time.min) preS.amt = preS.init;
  preS.amt:time = preS_formation.rate + -1*S_generation.rate;
  preS = preS.amt/compartment;
  preS.init = 0;
  when (time=time.min) Ya.amt = Ya.init;
  Ya.amt:time = -1*Glucose_sensor_inactivation.rate;
  Ya = Ya.amt/compartment;
  Ya.init = .9;
  when (time=time.min) S.amt = S.init;
  S.amt:time = S_generation.rate + S_formation.rate + -1*S_degradation.rate;
  S = S.amt/compartment;
  S.init = 0;
  Gluc.amt = Gluc*compartment;
  Gluc = Gluc.init/compartment;
  Gluc.init = 0;
  when (time=time.min) Yi.amt = Yi.init;
  Yi.amt:time = Glucose_sensor_inactivation.rate;
  Yi = Yi.amt/compartment;
  Yi.init = 0;
  when (time=time.min) V.amt = V.init;
  V.amt:time = V_formation.rate + -1*V_degradation.rate;
  V = V.amt/compartment;
  V.init = 30;
  when (time=time.min) Pi.amt = Pi.init;
  Pi.amt:time = Photoreceptor_decay.rate;
  Pi = Pi.amt/compartment;
  Pi.init = 0;
  Photoreceptor_activation.rate = compartment*Pfr*IfrSfrPfr;
  Photoreceptor_inactivation.rate = IrSrPr*Pr*compartment;
  Transducer_activation.rate = Xi*kia*Pr*compartment;
  Transducer_inactivation.rate = kai*Xa*compartment;
  preS_formation.rate = prepreS*kx*Xa*compartment;
  S_generation.rate = preS*ky*Ya*compartment;
  Glucose_sensor_inactivation.rate = kG*Ya*Gluc*compartment;
  S_formation.rate = compartment*(alpha1/(1+V^3));
  V_formation.rate = compartment*(alpha2/(1+S^3));
  S_degradation.rate = S*compartment;
  V_degradation.rate = compartment*V;
  Photoreceptor_decay.rate = compartment*kd*Pr;

  // variable properties
  compartment.sbmlRole="compartment";
  Pfr.amt.sbmlRole="speciesAmount";
  Pfr.sbmlRole="speciesConcentration";
  Pfr.sbmlCompartment="compartment";
  Pfr.init.sbmlRole="speciesInitialAmount";
  Pr.amt.sbmlRole="speciesAmount";
  Pr.sbmlRole="speciesConcentration";
  Pr.sbmlCompartment="compartment";
  Pr.init.sbmlRole="speciesInitialAmount";
  Xi.amt.sbmlRole="speciesAmount";
  Xi.sbmlRole="speciesConcentration";
  Xi.sbmlCompartment="compartment";
  Xi.init.sbmlRole="speciesInitialAmount";
  Xa.amt.sbmlRole="speciesAmount";
  Xa.sbmlRole="speciesConcentration";
  Xa.sbmlCompartment="compartment";
  Xa.init.sbmlRole="speciesInitialAmount";
  prepreS.amt.sbmlRole="speciesAmount";
  prepreS.sbmlRole="speciesConcentration";
  prepreS.sbmlCompartment="compartment";
  prepreS.init.sbmlRole="speciesInitialAmount";
  preS.amt.sbmlRole="speciesAmount";
  preS.sbmlRole="speciesConcentration";
  preS.sbmlCompartment="compartment";
  preS.init.sbmlRole="speciesInitialAmount";
  Ya.amt.sbmlRole="speciesAmount";
  Ya.sbmlRole="speciesConcentration";
  Ya.sbmlCompartment="compartment";
  Ya.init.sbmlRole="speciesInitialAmount";
  S.amt.sbmlRole="speciesAmount";
  S.sbmlRole="speciesConcentration";
  S.sbmlCompartment="compartment";
  S.init.sbmlRole="speciesInitialAmount";
  Gluc.amt.sbmlRole="speciesAmount";
  Gluc.sbmlRole="speciesConcentration";
  Gluc.sbmlCompartment="compartment";
  Gluc.init.sbmlRole="speciesInitialAmount";
  Yi.amt.sbmlRole="speciesAmount";
  Yi.sbmlRole="speciesConcentration";
  Yi.sbmlCompartment="compartment";
  Yi.init.sbmlRole="speciesInitialAmount";
  V.amt.sbmlRole="speciesAmount";
  V.sbmlRole="speciesConcentration";
  V.sbmlCompartment="compartment";
  V.init.sbmlRole="speciesInitialAmount";
  Pi.amt.sbmlRole="speciesAmount";
  Pi.sbmlRole="speciesConcentration";
  Pi.sbmlCompartment="compartment";
  Pi.init.sbmlRole="speciesInitialAmount";
  Photoreceptor_activation.rate.sbmlRole="rate";
  IfrSfrPfr.sbmlRole="parameter";
  Photoreceptor_inactivation.rate.sbmlRole="rate";
  IrSrPr.sbmlRole="parameter";
  Transducer_activation.rate.sbmlRole="rate";
  kia.sbmlRole="parameter";
  Transducer_inactivation.rate.sbmlRole="rate";
  kai.sbmlRole="parameter";
  preS_formation.rate.sbmlRole="rate";
  kx.sbmlRole="parameter";
  S_generation.rate.sbmlRole="rate";
  ky.sbmlRole="parameter";
  Glucose_sensor_inactivation.rate.sbmlRole="rate";
  kG.sbmlRole="parameter";
  S_formation.rate.sbmlRole="rate";
  alpha1.sbmlRole="parameter";
  V_formation.rate.sbmlRole="rate";
  alpha2.sbmlRole="parameter";
  S_degradation.rate.sbmlRole="rate";
  V_degradation.rate.sbmlRole="rate";
  Photoreceptor_decay.rate.sbmlRole="rate";
  kd.sbmlRole="parameter";
}