Equations for capump.mod
INITIAL {
ica = 0
ica_pmp = 0
ica_pmp_last = 0
SOLVE pmp STEADYSTATE sparse
}
BREAKPOINT {
SOLVE pmp METHOD sparse
ica_pmp_last = ica_pmp
ica = ica_pmp
}
KINETIC pmp {
~ cabulk <-> cam (width/tau, width/tau)
~ cam + pump <-> capump ((1e7)*k1, (1e10)*k2)
~ capump <-> cao + pump ((1e10)*k3, (1e10)*k4)
ica_pmp = (1e-7)*2*FARADAY*(f_flux - b_flux)
: ica_pmp_last vs ica_pmp needed because
: of STEADYSTATE calculation
~ cam << (-(ica - ica_pmp_last)/(2*FARADAY)*(1e7))
CONSERVE pump + capump = (1e13)*pump0
COMPARTMENT width {cam} : volume has dimensions of um
COMPARTMENT (1e13) {pump capump} : area is dimensionless
COMPARTMENT 1(um) {cabulk}
COMPARTMENT (1e3)*1(um) {cao}
cai = (0.001)*cam
}