- Author:
- Megan Guidry <mgui210@aucklanduni.ac.nz>
- Date:
- 2017-10-17 13:34:15+13:00
- Desc:
- First commit. Uploading models and protocols
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/4ca/rawfile/9680cc4249799c3501eec685cbf076d80f3b8ece/src/Model_Versions/main_WL_fixedCai.py
import sys
import FixedCai_parameters
#preset_afterloads = [0.135, 0.1911, 0.2576, 0.3359, 0.4312, 0.6, 0.9]
preset_CaiSL = [2.3]
preset_afterloads = [0.135, 0.1911, 0.2576, 0.3359, 0.4312, 0.6, 0.9]
print range(len(preset_afterloads))
CaiSL_count = 0
while CaiSL_count < len(preset_CaiSL):
Afterload_count = 0
value_CaiSL = preset_CaiSL[CaiSL_count]
if value_CaiSL == 1.955:
parameters = [0.089002126, 0.283359682, 39.423358, 2.84495395, 1.6649028, 37.6686509, 9.995708, 49.5660573, 100.759014, 50, 70.3621179]
elif value_CaiSL == 2.024:
parameters = [0.08968451, 0.27926453, 37.66516463, 2.80900835, 1.67887931, 36.03943529, 5.42060899, 45.08829658, 82.71763373, 50, 72.90385783]
elif value_CaiSL == 2.093:
parameters = [0.08981397, 0.23605426, 37.19706502, 2.36169272, 1.76341758, 31.40873567, 9.99030452, 40.1918327, 74.94501223, 50 ,79.1466021]
elif value_CaiSL == 2.162:
parameters = [0.08982691, 0.23170048, 35.4328648, 2.32764163, 1.78087606, 29.92922172, 9.99824712, 36.10844273, 70.54464675, 50, 87.58238235]
elif value_CaiSL == 2.231:
parameters = [0.089699541, 0.225542405, 34.1045897, 2.28321127, 1.80277288, 28.3291654, 0.016845165, 32.1908114, 66.8007906, 50, 97.3577797]
elif value_CaiSL == 2.3:
parameters = [0.0894001, 0.21621596, 33.1456627, 2.21842617, 1.83302661, 26.491962, 9.97086786, 28.7470311, 63.8698117, 50, 108.854814]
value_tstart = parameters[3]
value_Tau1 = parameters[5]
value_Tau2 = parameters[7]
value_Tau3 = parameters[10]
value_Ca_dia = parameters[0]
value_a1 = parameters[1]
value_b = parameters[2]
value_k1 = parameters[4]
value_c = parameters[8]
CaiSL_count += 1
while Afterload_count < len(preset_afterloads):
value_afterload = preset_afterloads[Afterload_count]
user_number = Afterload_count + 1
execfile(r'Workloop_protocol_fixedCai.py')
Afterload_count += 1
print value_CaiSL
print value_afterload
#from mlabwrap import mlab
#mlab.plotting_figures('wrkloop.csv')
sys.exit()
#Grab the csv workloops generated above (they will be in the same folder
#as main.py and Workloop_protocol_with_passive_cycling2.py) and plot
