zebrafish_cb1r.ode_sub_am251.set

## Set file for zebrafish_cb1r.ode on Fri Apr 9 2021
12   Number of equations and auxiliaries
54   Number of parameters
# Numerical stuff
1    nout
40    nullcline mesh
8   Qual RK
70000  total
0.02  DeltaT
0  T0
10000  Transient
500001  Bound
1e-12  DtMin
1  DtMax
0.001  Tolerance
0.001  Abs. Tolerance
0  Max Delay
100   Eigenvector iterates
0.001  Eigenvector tolerance
0.001  Newton tolerance
0  Poincare plane
1e-05  Boundary value tolerance
1e-05  Boundary value epsilon
20   Boundary value iterates
0   Poincare None
1   Poincare variable
1   Poincare sign
0   Stop on Section
0   Delay flag
46870.7  Current time
46870.7  Last Time
1   MyStart
1   INFLAG
# Delays
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
0.0
# Bndry conds
0
0
0
0
0
0
0
0
0
0
0
0
# Old ICs
-33.269  MV
0.0048411  MW
3.3  MCA
0.1  IGSYN
-29  EV
0.007  EW
4.7  ECA
0.01  ES
-29  I2V
0.007  I2W
4.7  I2CA
0.01  I2S
# Ending  ICs
-42.13341190895153  MV
0.001711400771223228  MW
1.43224952950728  MCA
1.716601931422579  IGSYN
-33.0107098529168  EV
0.004974115881352404  EW
6.440678748456065  ECA
0.01272894404551673  ES
-41.69695790687987  I2V
0.001798206102459405  I2W
18.04025272833771  I2CA
0.00531446653509741  I2S
# Parameters
-1.2  vf1
18  vf2
12  vf3
17  vf4
4  gca
120  vca
2  gl
8  gk
-60  vl
-84  vk
0.23  phi
1  cm
1  kca
0.25  gkca
0.005  eps
0.19  mu
10  ca0
31  mbiapp
20000  it0
1  fre
2  dur
0  thetas
4  ss
0.5  gsynmm
0.029  sm2
-50  mvsyn
20  igsynmax
10  mca0
10000  tauagm
4  mgca
0.3  mgkca
0.9  mkca
0  cb1rim
0  cb1rem
0.15  gsynem
0.5  gsynim
43.9  ebiapp
60  eiapp
30  evsyn
15  ealpha
0.3  ebeta
20  ecm
0.5  egsyni
0.005  eeps
36  ibiapp
-50  ivsyn
3  ialpha
0.3  ibeta
20  icm
36  i2biapp
8.5  i2alpha
0.046  i2beta
0  cb1rei
0.7  gsynei
# Graphics
0.7071097251469876  rm
-0.4999999999826661  rm
0.4999958366025516  rm
0.7071038372138505  rm
0.5000041633974484  rm
-0.4999999999826661  rm
0  rm
0.7071038372138505  rm
0.7071097251469876  rm
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
0    
1    
1    
1    
0    
-1000   
1000   
0    
0   3DFlag
1   Timeflag
0   Colorflag
0   Type
1  color scale
0   minscale
70000   xmax
0   xmin
39.8256950378418   ymax
-53.27024459838867   ymin
12   zmax
-12   zmin
35000   
2.857142857142857e-05   
-6.722274780273438   
0.02148321406728305   
0   
0.08333333333333333   
45   Theta
45   Phi
0    xshft
0    yshft
0    zshft
0   xlo
-53.27024459838867   ylo
10000   
-40   
70000   xhi
39.8256950378418   yhi
70000   
30   
# Transpose variables etc
MV
2   n columns
1   n rows
1   row skip
1   col skip
1   row 0
# Coupling stuff for H funs
0
0
0
0
0
0
0
0
0
0
0
0
# Array plot stuff

1   NCols
0   Row 1
50   NRows
8   RowSkip
0  Zmin
1  Zmax
# Torus information 
0    Torus flag 1=ON
6.283185307179586  Torus period
# Range information
vf1
-1   eq-range stab col
0   shoot flag 1=on
10   eq-range steps
0  eq_range low
1  eq_range high
MV
MV
20   Range steps
0   Cycle color 1=on
1   Reset data 1=on
1   Use old I.C.s 1=yes
0  Par1 low
0  Par2 low
1  Par1 high
1  Par2 high
vf1
0   BVP side
0   color cycle flag 1=on
10   BVP range steps
0  BVP range low
1  BVP range high
RHS etc ...
dMV/dT=(MBIAPP-E2MSYN-I2MSYN-MGCA*MINF(MV)*(MV-VCA)-GK*MW*(MV-VK)-GL*(MV-VL)-MGKCA*(MCA/(MCA+CA0))*(MV-VK)-GSYNMM*SM2*(MV-MVSYN))/CM
dMW/dT=PHI*(WINF(MV)-MW)/TAUW(MV)
dMCA/dT=EPS*(-MU*GCA*MINF(MV)*(MV-VCA)-MKCA*MCA)
dIGSYN/dT=(IGSYNMAX/(MCA+MCA0)-IGSYN)/TAUAGM
dEV/dT=(EBIAPP+EIAPP*(HEAV(T-IT0)*HEAV(SIN(PI*2*FRE/1000*T)-SIN(5*PI/2-PI*DUR*FRE/1000)))-GCA*MINF(EV)*(EV-VCA)-GK*EW*(EV-VK)-GL*(EV-VL)-GKCA*(ECA/(ECA+CA0))*(EV-VK))/ECM
dEW/dT=PHI*(WINF(EV)-EW)/TAUW(EV)
dECA/dT=EPS*(-MU*GCA*MINF(EV)*(EV-VCA)-KCA*ECA)
dES/dT=EALPHA*(1-ES)*ESINF(EV)-EBETA*ES
dI2V/dT=(I2BIAPP-E2ISYN2-GCA*MINF(I2V)*(I2V-VCA)-GK*I2W*(I2V-VK)-GL*(I2V-VL)-GKCA*(I2CA/(I2CA+CA0))*(I2V-VK))/ICM
dI2W/dT=PHI*(WINF(I2V)-I2W)/TAUW(I2V)
dI2CA/dT=EPS*(-MU*GCA*MINF(I2V)*(I2V-VCA)-KCA*I2CA)
dI2S/dT=I2ALPHA*(1-I2S)*ISINF(I2V)-I2BETA*I2S

where ...
E2MSYN =  GSYNEM*IGSYN*(1+CB1REM)*ES*(MV-EVSYN)
   
I2MSYN =  GSYNIM*IGSYN*(1-CB1RIM)*I2S*(MV-IVSYN)
   
E2ISYN2 =  GSYNEI*IGSYN*ES*(1+CB1REI)*(I2V-EVSYN)
   

User-defined functions:
MINF(V) = 0.5*(1+TANH((V-VF1)/VF2))
  
WINF(V) = 0.5*(1+TANH((V-VF3)/VF4))
  
TAUW(V) = 1/COSH((V-VF3)/(2*VF4))
  
ESINF(V) = 1./(1.+EXP(-(V+THETAS)/SS))
  
ISINF(V) = 1./(1.+EXP(-(V+THETAS)/SS))