- Author:
- WeiweiAi <wai484@aucklanduni.ac.nz>
- Date:
- 2022-07-12 17:11:03+12:00
- Desc:
- Change the stimulation protocol to double pulse with conditioning; add steady state equations
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/8af/rawfile/113b9d7fbcebee0427016d8f7d9b5b181656ab2e/Components/buildsrc/ON_S2.txt
def model ON_S2 as
def import using "../cellLib/Components/units.cellml" for
unit ms using unit ms;
unit per_ms using unit per_ms;
enddef;
def import using "../cellLib/Components/MarkovS.cellml" for
comp S2 using comp S2;
enddef;
def comp ON_S2 as
var k_L: per_ms {pub: in};
var k_Lminus: per_ms {pub: in};
var k_C: per_ms {pub: in};
var k_Cminus: per_ms {pub: in};
var f: dimensionless {pub: in};
var k_i1: per_ms {pub: out, priv: out};
var k_i2: per_ms {pub: out, priv: out};
var k_1i: per_ms {pub: out, priv: out};
var k_2i: per_ms {pub: out, priv: out};
var norm: dimensionless {pub: in, priv: out};
var S_1: dimensionless {pub: in, priv: out};
var S_2: dimensionless {pub: in, priv: out};
var S_i: dimensionless {pub: out, priv: in};
var dS_i: per_ms {pub: out, priv: in};
var Si_init: dimensionless {pub: in, priv: out};
var N: dimensionless {pub: in};
var t: ms {pub: in, priv: out};
k_i1 = k_Cminus*N*f;
k_i2 = k_Lminus*pow(f, N);
k_1i = k_C/f;
k_2i = k_L/pow(f, N);
enddef;
def group as encapsulation for
comp ON_S2 incl
comp S2;
endcomp;
enddef;
def map between S2 and ON_S2 for
vars S_i and S_i;
vars dS_i and dS_i;
vars t and t;
vars Si_init and Si_init;
vars k_i1 and k_i1;
vars k_i2 and k_i2;
vars k_1i and k_1i;
vars k_2i and k_2i;
vars norm and norm;
vars S_1 and S_1;
vars S_2 and S_2;
enddef;
enddef;