updated inhibitor initial conditions

mjowen · Jan 17, 2024 · df76634 · df76634
1 parent a08ff7f
commit df76634
Show file tree

Hide file tree

Showing 5 changed files with 5 additions and 5 deletions.
diff --git a/chatterjee.py b/chatterjee.py
@@ -96,7 +96,7 @@ def setIC(ICvector = np.array([10e-12, 1.4e-6, 2e-8, 1e-8, 1e-8/100, 7e-10, 9e-8
     IC = np.zeros(81);
     IC[[0,13,20,1,3,14,10,7,50,28,25,80]] = np.append(ICvector,eps0); # np.append adds the initial value for eps onto the input. Indices correspond to factors in initial condition vector
     if includeExtras:
-        IC[[56,58,47]] = [4e-5, 9.5e-7, 2.1e-6];
+        IC[[56,58,47]] = [4e-5, 9.75e-7, 1.7e-6]; #a1AT, a2AP, C1-inh
     return IC
 
 def getRates(includeExtras = False):

diff --git a/lakshmanan.py b/lakshmanan.py
@@ -57,7 +57,7 @@ def setIC(ICvector = np.array([10e-12, 1.4e-6, 2e-8, 1e-8, 1e-8/100, 7e-10, 9e-8
     IC = np.zeros(42);
     IC[[0,13,20,1,3,14,10,7,34,28,25]] = ICvector; # Indices correspond to factors in initial condition vector
     if includeExtras:
-        IC[39] = 1.7e-6;
+        IC[39] = 1.7e-6; #C1-inh
     return IC
 
 def getRates(includeExtras = False):

diff --git a/panteleev.py b/panteleev.py
@@ -55,7 +55,7 @@ def getRates(includeExtras = False):
     n = np.array([260, 750, 16000, 2700]);
     h = np.array([0.44*1e9/60, 6*1e9/60, 8.2e-06*1e9/60, 0.00015*1e9/60, 4e-05*1e9/60, 1.36e-05*1e9/60, 0.0012*1e9/60, 2.2e-05*1e9/60, 0.00041*1e9/60, 0.0001*1e9/60, 3e-06*1e9/60, 0.00037*1e9/60, 6.3e-05*1e9/60, 0.35/60, 7.7*1e9/60, 1.9e-05*1e9/60, 2.6e-05*1e9/60, 6e-06*1e9/60, 0.0054*1e9/60, 0.00014*1e9/60, 6e-06*1e9/60, 6e-06*1e9/60, 7e-07*1e9/60, 0.00039*1e9/60]);
     if includeExtras:
-        i = np.array([3250/1e9, 40000/1e9, 950/1e9, 0/1e9, 70/1e9, 2100/1e9, 0/1e9]);
+        i = np.array([3000/1e9, 40000/1e9, 975/1e9, 0/1e9, 79/1e9, 1700/1e9, 0/1e9]); #a2M, a1AT, a2AP, hep, PCI, C1-inh, PS
     else:
         i = np.array([0, 0, 0, 0, 0, 0, 0]);
     p = np.array([7.5e-5/1e9]);

diff --git a/tyurin.py b/tyurin.py
@@ -48,7 +48,7 @@ def setIC(ICvector = np.array([10e-12, 1.4e-6, 2e-8, 1e-8, 1e-8/100, 7e-10, 9e-8
     IC = np.zeros(33);
     IC[[0,18,8,4,3,15,13,6,11,25,23]] = ICvector; # Indices correspond to factors in initial condition vector
     if includeExtras:
-        IC[[31, 28, 30, 29, 27, 26]] = [7e-8, 4e-5, 3.25e-6, 9.5e-7, 4.6e-10, 2.1e-6];
+        IC[[31, 28, 30, 29, 27, 26]] = [7.9e-8, 4e-5, 3e-6, 9.75e-7, 4.6e-10, 1.7e-6]; #PCI, a1AT, a2M, a2AP, PAI1, C1inh
     return IC
 
 def getRates(includeExtras = False):

diff --git a/zhu.py b/zhu.py
@@ -52,7 +52,7 @@ def setIC(ICvector = np.array([10e-12, 1.4e-6, 2e-8, 1e-8, 1e-8/100, 7e-10, 9e-8
     IC = np.zeros(37);
     IC[[0,17,7,3,2,14,12,5,10,24,22]] = ICvector; # Indices correspond to factors in initial condition vector
     if includeExtras:
-        IC[[26, 28, 27, 25, 34]] = [4e-5, 3.25e-6, 9.5e-7, 4.6e-10, 2.1e-6];
+        IC[[26, 28, 27, 25, 34]] = [4e-5, 3e-6, 9.75e-7, 4.6e-10, 1.7e-6]; #a1AT, a2M, a2AP, PAI1, C1-inh
     return IC
 
 def getRates(includeExtras = False):