-
Notifications
You must be signed in to change notification settings - Fork 1
/
IK1.em
executable file
·124 lines (102 loc) · 3.01 KB
/
IK1.em
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
@{
IK1_I = 24.29288398159507
IK1_Pspm = 0.6373407661728812
IK1_amp = {
"V" : 2.6,
"EMB" : 2.6,
"LAT" : 2.6,
"NEO" : 2.6,
"SAN" : 0.007079,
"A" : 2.6,
}
}
System System(/CELL/MEMBRANE/IK1)
{
StepperID ODE;
Variable Variable( I )
{
Value @IK1_I;
}
Variable Variable( POpen )
{
Value 4.26E-01;
}
Variable Variable( vPspm )
{
Value -3.73433916526e-07;
}
Variable Variable( Pspm )
{
Value @IK1_Pspm;
}
Variable Variable( fracO )
{
Value 0.414217508395;
}
Variable Variable( GX ){
Value @( Kir2_1[SimulationMode]);
}
Process IK1AssignmentProcess( I )
{
StepperID PSV;
Priority 20;
VariableReferenceList
[ Vm :..:Vm 0 ]
[ EK :..:EK 0 ]
[ Mg :../../CYTOPLASM:Mg 0 ]
[ SPM :../../CYTOPLASM:Spermine 0 ]
[ vPspm :.:vPspm 1 ]
[ Pspm :.:Pspm 0 ]
[ fracO :.:fracO 1 ]
[ pOpen :.:POpen 1 ]
[ GX :.:GX 0 ]
# [ GX :../../CYTOPLASM:Kir2_1 0 ]
[ Ko :/:K 0 ]
[ Cm :..:Cm 0 ]
[ I :.:I 1 ];
Phi 0.9;
amplitude @(IK1_amp[SimulationMode]);
constant 5.4;
power 0.6;
A -0.6;
B 3.1;
}
@{''' Vm.getValue() - EK.getValue() + 8000.0 * Mg.getMolarConc() の
計算手続きの関係でExpressionProcessモデルとの間に若干の誤差が生じる。
Priorityを上げても解消しないので、単純に計算順序の問題ではない？
'''}
Process ZeroVariableAsFluxProcess( Pspm )
{
Priority 15;
VariableReferenceList
[ Pspm :.:Pspm 1 ]
[ vPspm :.:vPspm 0 ];
}
@setCurrents( [ 'I' ], [ 'K', 'I' ] )
}
@{'''
Author Maria Takeuchi
Author Yasuhiro Naito
Version 0.2 2008-11-24 13:40:48 +0900
<IK1 name="IK1" initial_value="24.29288398159507" units="pA"
className="org.simBio.bio.kuzumoto_et_al_2007.current.IK1">
<KoDependencyIK1 name="permeabilityK" initial_value="2.5449343128168245" units="pA/mV"
className="org.simBio.bio.kuzumoto_et_al_2007.function.KoDependency_IK1">
<parameter name="amplitude" initial_value="2.6" />
<parameter name="constant" initial_value="5.4" units="mM" />
<parameter name="power" initial_value="0.6" units="dimension_less" />
<parameter name="A" initial_value="-0.6" units="dimension_less" />
<parameter name="B" initial_value="3.1" units="dimension_less" />
<link name="Ko" initial_value="../../../K" />
</KoDependencyIK1>
<variable name="Pspm" initial_value="0.6373407661728812" />
<parameter name="Phi" initial_value="0.9" />
<link name="Vm" initial_value="../Vm" units="mV" />
<link name="Cm" initial_value="../membrane capacitance" />
<link name="Mg" initial_value="../Mg" />
<link name="SPM" initial_value="../Spermine" />
<link name="reversalPotential" initial_value="../K reversal potential" units="mV" />
<link name="current" initial_value="../current" />
<link name="currentK" initial_value="../currentK" />
</IK1>
'''}