phi (Phi-deltaPhi)/betaPhi "Dimensionless potential"
eps_delay "e0+(einf-e0)*exp( -exp( -zeta*(Phi-Phibar) ) )" "Delay function"
tact 50[ms]-Z/cL*50[ms] "Activation time"
betat ts*(1-tau0*(tact-t0)/(t1-t0)) "Time scaling factor"
lambda sqrt(solid.hmm1.fib1.IaCIe) "Stretch in fiber direction"
thetaSwitch if(lambda>1,1,0) "Activation parameter"
Ie Cm*betaPhi/betat*(c*phi*(phi-alpha)*(phi-1)+r*phi) "Ionic current (pure electric)"
Im Cm*betaPhi/betat*thetaSwitch*Gs*(lambda-1)*(phi-phis) "Ionic current (stretch induced)"
D_iso d_iso*(Cm*Chi_m) "Isotropic conductivity"
D_ani d_ani*(Cm*Chi_m) "Anisotropic conductivity, fiber direction"
DfibXX solid.hmm1.fib1.a0X*D_ani*solid.hmm1.fib1.a0X "Conductivity tensor, global coordinate system"
DfibXY solid.hmm1.fib1.a0X*D_ani*solid.hmm1.fib1.a0Y "Conductivity tensor, global coordinate system"
DfibXZ solid.hmm1.fib1.a0X*D_ani*solid.hmm1.fib1.a0Z "Conductivity tensor, global coordinate system"
DfibYY solid.hmm1.fib1.a0Y*D_ani*solid.hmm1.fib1.a0Y "Conductivity tensor, global coordinate system"
DfibYZ solid.hmm1.fib1.a0Y*D_ani*solid.hmm1.fib1.a0Z "Conductivity tensor, global coordinate system"
DfibZZ solid.hmm1.fib1.a0Z*D_ani*solid.hmm1.fib1.a0Z "Conductivity tensor, global coordinate system"