-
Notifications
You must be signed in to change notification settings - Fork 0
/
reduced.pDE
executable file
·368 lines (296 loc) · 16 KB
/
reduced.pDE
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
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
# 06Apr10
# Copied from 1spk-th-102606b-n420-Passive-v7c.pDE
# Modified to have what 1spk-th-10606b-n420+Fix-h+Naf+KDR-3iALpaper-v7.pDE ha
#
#================================================================================
# No sensitivity analysis
Fit-SA 0 0
Tmp-dir MyTemp
#Fit 102606b 311 200 10 3 0.85 1.0 randE /Volumes/d1/Users/ximing/Projects/ParSims/Runs_pDE/thin-1stspike/main-n400-Vbtim-simL-v8d-reduced.hoc
#Fit 102606b 199 110 25 3 0.7 1.0 randE /Volumes/d1/Users/ximingProjects/ParSims/Runs_pDE/thin-1stspike/main-n400-Vbtim-simL-v8d-reduced.hoc
#Fit 102606b 713 120 50 3 0.85 1.0 randE /Volumes/d1/Users/ximing/Projects/ParSims/Runs_pDE/thin-1stspike/main-n400-Vbtim-simL-v8d-reduced.hoc
Fit 102606b 713 100 80 3 0.85 1.0 randE /Volumes/d1/Users/ximing/Projects/ParSims/Runs_pDE/thin-1stspike/main-n400-Vbtim-simL-v8d-reduced.hoc
#================================================================================
# 10oct08 4:40pm
# Set up traces for 102606b
# 10oct08 4pm
# Stage 0; select traces; time v=0 post spike for +1500 800 400 200
# Stich 800 with 3 4 5 spikes; select v of end trace, stitch at same volt after spike
# Determine bl_avg with 13 traces; BL_avg=-74.3052
# Determine BBF; need closer fit first.
# Run Passive fit with 9 traces
# Use lambda_f_d=0.1
Pop-init uniform noset
# Determined BBF=1.6 from run so far. Resume above
#================================================================================
VERBATIM
//cell_file_name="n400td-pc2-mM-axn-reduced.hoc"
cell_file_name="n400td-pc2-sp-axn.hoc"
obfunc dfc_set_stims_etc(){
print "check pDE dfc_set_stims++++++++++++++++++++"
/*
if(vdef("aX_hill_nseg"))\
execute( "if(axon_hill.nseg<aX_hill_nseg) axon_hill.nseg=aX_hill_nseg" )
if(vdef("aX_iseg_nseg")) \
execute( "for i=0,aX_iseg_n-1 {if(axon_iseg[i].nseg<aX_iseg_nseg) axon_iseg[i].nseg= aX_iseg_nseg }" )
//if(vdef("aX_prox_nseg")) \
// execute( "for i=0,aX_prox_n-1 {if(axon_prox[i].nseg<aX_prox_nseg) axon_prox[i].nseg= aX_prox_nseg }" )
if( vdef("aX_mnode_nseg" )) \
execute( "forsec \"mnode\" if(nseg<aX_mnode_nseg) nseg= aX_mnode_nseg" )
*/
// execute( "axon_prox.nseg = aX_prox_nseg" )
// execute( "for i=0,aX_iseg_n-1 {axon_iseg[i].nseg= aX_iseg_nseg }" )
stim.amp=$1
stim.dur=$2
stim.del=$3
print "check stim.del +++++++++++++++++++++++++++"
dfc_options_all( dfc_eD )
return stim
}
obfunc dfc_set_tgt(){
print "dfc_set_tgt()+++++++++++++++++++++++++++"
dfc_tgt.stim_amp=$1
dfc_tgt.stim_dur=$2
dfc_tgt.stim_del=$3
if( numarg()>3 ) dfc_tgt.cut_at2($4,$5)
return dfc_tgt
}
proc dfc_sim_BBF(){
dfc_simL_i.o(0).dc_shift(stim.del+BBFdelay,stim.del+stim.dur+BBFdelay,-BBF*stim.amp)
}
//================================================================================
// New stuff for pkbeg posfix
obfunc dfc_seg_look_options(){ local i localobj dI
dI = new eMatch() if( numarg()>0)dI =$o1
dI.seg_option = 1101
dI.seg_dt = 0.09
dI.seg_slp_dt = 0.09
dI.seg_slp_th = 10
dI.seg_v_th = 0
dI.seg_tiny_last = 0
if( strcmp(dfc_i.id,"102606b-2ms+1500pA-pkbeg") == 0 ){ dI.seg_dt=dI.seg_slp_dt= 0.02 }
return dI
}
obfunc dfc_options_all(){ local i localobj eD
if( numarg()==0 ) eD = dfc_seg_look_options()
if( numarg()> 0 ) eD = dfc_seg_look_options($o1)
eD.match_option = 5.2
eD.alpha = eD.beta = eD.gamma = eD.delta = eD.eta = eD.mF1 = eD.mF2 = eD.pkP = eD.pkPs = 0
eD.ilog = 1
return eD
}
// rough match to train;
// necessary to avoid traces that don't return to baseline
// also check spikes also appear in axon
func dfc_trim(){ local i, t0_tail, f1, f2, r1, r2, res localobj eD
f1=$1 f2=$2
eD = dfc_options_all(dfc_eD)
eD.eta = 1
eD.mF2 = 200
eD.pkP = 1500*5
eD.pkPs = 1
eD.pkPsA = 1
t0_tail = stim.del + stim.dur // end of stimulus
eD.tmin = t0_tail // compare tails after end stimulus: axon vs tgt
r2 = eD.m0( dfc_simL_i.o(1), dfc_tgt_i )
eD.tmin = 0
printf( "\n\t >>dfc_trim: DfpwPkdVa_aXe: " )
r1 = eD.DfpwPkdVa_aXe(dfc_simL_i,dfc_tgt_i,-1,-1) // 800=8.90, 400=15.2, 200=30.82, 1500=6.7
res = r1*f1 + r2*f2
printf( "\n\t dfc_trim t0_tail %g r1,2 %g %g f1,f2 %g %g ri*fi %g %g\n\t res %g<<\n",\
t0_tail, r1,r2, f1,f2, r1*f1, r2*f2, res )
return res
}
//================================================================================
// match first spike (up to tmax)
// gamma and delta to match dV and ddV; largely independent of trace length/stim.amp
// eta scaled by trace length
//================================================================================
func dfc_pkbeg(){ local i, tmax, dV_f, eta_f, r1 localobj eD
tmax = $1
dV_f = 45*2
eta_f = 300 + 100/(tmax-stim.del-1.1) // A + B/(time from stim.del to base of spike) A=for spike; B for before spike
eD = dfc_eD // eD = dfc_options_all(dfc_eD)
// eD.mF2 = 200*10
eD.mF2 = 6000 * stim.amp // peak time error proportional to stimulation amplitude
eD.pkP = 1500*5
eD.pkPs = 1
eD.phi = 600/116/10 // range_dV/range_V old=0.4 * 10
eD.phi_1 = 0.5+0.15
eD.phi_2 = 0.45-0.02
eD.chi = 1000/116/100 // range_ddV/range_V old=0.02 * 8
eD.chi_1 = 0.5+0.15
eD.chi_2 = 0.15
if( numarg()>1) eD.chi_2=$2 // adhoc inserted to override
eD.phi*=dV_f eD.chi*=dV_f
eD.eta = eta_f
// eD.mF2=0 eD.pkP=1000*eta_f // for initial pop for Naf & KDR parms
eD.tmax = tmax
printf( "\n\t >>dfc_pkbeg: DfpwPkdVa: " )
r1 = eD.DfpwPkdVa(dfc_simL_i.o(0),dfc_tgt_i) // 800=8.90, 400=15.2, 200=30.82, 1500=6.7
printf( "\n\t dfc_pkbeg: r1 %g \n\t<<\n", r1 )
return r1
}
ENDVERBATIM
#================================================================================
# Fit-case name template func scale hoc_var hoc_before hoc_after tgt_files tgt_before tgt_after
# Active traces 20kHz
# For Naf & KDR initial parm search
#Fit-case 102606b-2ms+1500pA eMatch dfc_eD.DfpwPkdVa(dfc_simL_i.o(0),dfc_tgt_i) 5 soma dfc_set_stims_etc(1500/1000,2,5);dfc_eD.eta=1;dfc_eD.pkP=1500*5;dfc_eD.pkPs=1 dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF/2ms+1500pA.htf blMax=4 dfc_set_tgt(1500/1000,2,21,136,16)
#Fit-case 102606b-100ms+800pA eMatch dfc_eD.DfpwPkdVa(dfc_simL_i.o(0),dfc_tgt_i) 5 soma dfc_set_stims_etc(800/1000,100,5);dfc_eD.eta=1;dfc_eD.pkP=1500*5;dfc_eD.pkPs=1 dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF/100ms+800pA.htf blMax=4 dfc_set_tgt(800/1000,100,21,200,16)
# End for Naf & KDR initial parm search
# Note dfc_eD.seg_dt=dfc_eD.seg_slp_dt=0.02 cause end_stim almost on peak
# 2:active set ================================================================================
#Fit-case 102606b-2ms+1500pA-pkbeg eMatch dfc_pkbeg(7.36) 1 soma dfc_set_stims_etc(1500/1000,2,5);dfc_eD.seg_dt=dfc_eD.seg_slp_dt=0.02;dfc_eD.pkPs=3 dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF_20kHz/2ms+1500pA.htf blMax=4 dfc_set_tgt(1500/1000,2,21,7.36+1+16,16)
# later runs use the following instead
#modified
#Fit-case 2ms+1500pA_072006b_002_tr2_chopped.dat eMatch dfc_pkbeg(7.36,0.05) 1 soma dfc_set_stims_etc(1500/1000,2,5);dfc_eD.pkPs=3; dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/072006b/ACSF_0kHz/2ms+1500pA.htf blMax=4 dfc_set_tgt(1500/1000,2,21,7.36+1+16,16)
Fit-case 500ms+200pA_short.dat eMatch dfc_pkbeg(7.36,0.05) 1 soma dfc_set_stims_etc(200/1000,500,21.07);dfc_eD.pkPs=3; dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF_20kHz/2ms+1500pA.htf blMax=4 dfc_set_tgt(1500/1000,2,21,7.36+1+16,16)
# Note about TRIMMED file. Each trace was cut at a diffent time from same trace 2 of 102606b-100ms+800pA
# So, all trimmed traces should be compared to trace 2:102606b-100ms+800pA
# Trimmed #1; 26.39
#Fit-case 102606b-100ms+800pA-tr-pkbeg eMatch (dfc_trim(1,1)+dfc_pkbeg(10.03)) 1 soma;axon_inode[4] dfc_set_stims_etc(800/1000,26.39-5,5) dfc_sim_BBF() 1:cNeuro/Data/Grover-2006/102606b/ACSF_20kHz/100ms+800pA-trimmed.htf blMax=4 dfc_set_tgt(800/1000,26.39-5,5)
# Trimmed #2; 37.49
#Fit-case 102606b-100ms+800pA-tr-pkbeg eMatch (dfc_trim(1,1)+dfc_pkbeg(10.03)) 1 soma;axon_inode[4] dfc_set_stims_etc(800/1000,37.49-5,5) dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF_20kHz/100ms+800pA-trimmed.htf blMax=4 dfc_set_tgt(800/1000,37.49-5,5)
#Fit-case 102606b-100ms+400pA-pkbeg eMatch dfc_pkbeg(16.56) 1 soma dfc_set_stims_etc(400/1000,100,5) dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF_20kHz/100ms+400pA.htf blMax=4 dfc_set_tgt(400/1000,100,21,16.56+1+16,16)
#Fit-case 102606b-100ms+200pA-pkbeg eMatch dfc_pkbeg(34.01) 1 soma dfc_set_stims_etc(200/1000,100,5) dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF_20kHz/100ms+200pA.htf blMax=4 dfc_set_tgt(200/1000,100,21,34.01+1+16,16)
#--------------------------------------------------------------------------------
# Passive traces (all 10kHz)
# +750 has no spike; we could use m0; but we want to penalize unwanted spikes
#Fit-case 102606b-2ms+750pA eMatch dfc_eD.DfpwPkdVa(dfc_simL_i.o(0),dfc_tgt_i) 6.4*dfc_PassWgt*100/tstop soma dfc_set_stims_etc(750/1000,2,5);dfc_eD.eta=1 dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF/2ms+750pA.htf blMax=4 dfc_set_tgt(750/1000,2,21,136,16)
# m0
#Fit-case 102606b-2ms+375pA eMatch dfc_eD.m0(dfc_simL_i.o(0),dfc_tgt_i) 4*dfc_PassWgt*100/tstop soma dfc_set_stims_etc(375/1000,2,5) dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF/2ms+375pA.htf blMax=4 dfc_set_tgt(375/1000,2,21,136,16)
# m0
#Fit-case 102606b-100ms+100pA eMatch dfc_eD.m0(dfc_simL_i.o(0),dfc_tgt_i) 7.9*dfc_PassWgt*100/tstop soma dfc_set_stims_etc(100/1000,100,5) dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF_20kHz/100ms+100pA.htf blMax=4 dfc_set_tgt(100/1000,100,21,300,16)
#====== 20ms
# m0
#Fit-case 102606b-20ms+50pA eMatch dfc_eD.m0(dfc_simL_i.o(0),dfc_tgt_i) 2*dfc_PassWgt*100/tstop soma dfc_set_stims_etc(50/1000,20,5) dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF/20ms+50pA.htf blMax=4 dfc_set_tgt(50/1000,20,21,120,16)
# m0
#Fit-case 102606b-20ms-50pA eMatch dfc_eD.m0(dfc_simL_i.o(0),dfc_tgt_i) 2*dfc_PassWgt*100/tstop soma dfc_options_all();dfc_set_stims_etc(-50/1000,20,5) dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF/20ms-50pA.htf blMax=4 dfc_set_tgt(-50/1000,20,21,120,16)
# m0
#Fit-case 102606b-20ms-200pA eMatch dfc_eD.m0(dfc_simL_i.o(0),dfc_tgt_i) 8.3*dfc_PassWgt*100/tstop soma dfc_options_all();dfc_set_stims_etc(-200/1000,20,5) dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF/20ms-200pA.htf blMax=4 dfc_set_tgt(-200/1000,20,21,120,16)
#====== 500ms
# m0
#Fit-case 102606b-500ms+50pA eMatch dfc_eD.m0(dfc_simL_i.o(0),dfc_tgt_i) 4.1*dfc_PassWgt*100/tstop soma dfc_set_stims_etc(50/1000,500,21) dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF/500ms+50pA.htf blMax=19 dfc_set_tgt(50/1000,500,21,800,0)
# m0
#Fit-case 102606b-500ms-50pA eMatch dfc_eD.m0(dfc_simL_i.o(0),dfc_tgt_i) 3.9*dfc_PassWgt*100/tstop soma dfc_set_stims_etc(-50/1000,500,21) dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF/500ms-50pA.htf blMax=19 dfc_set_tgt(-50/1000,500,21,800,0)
# m0
#Fit-case 102606b-500ms-200pA eMatch dfc_eD.m0(dfc_simL_i.o(0),dfc_tgt_i) 12.7*dfc_PassWgt*100/tstop soma dfc_set_stims_etc(-200/1000,500,21) dfc_sim_BBF() 2:cNeuro/Data/Grover-2006/102606b/ACSF/500ms-200pA.htf blMax=19 dfc_set_tgt(-200/1000,500,21,800,0)
# --------------------------------------------------------------------------------
#
# Fit-parm name set fit-lims
Set-parm dfc_PassWgt 1 0:1
Set-parm RmCm_spine_gbar 0 0:1
Set-parm Spine_number 30000 0:40000
Set-parm BL_avg -74.3052 -80:-60
Set-parm BBF 1.60 0:20
Set-parm BBFdelay 0 0:0.1
Set-parm lambda_f_d 0.1 0.1:2
#Set-parm G_ena 65.4 60:70
#Set-parm G_ek -89.5 -95:-85
Set-parm G_eh -30 -40:-10
Set-parm celsius 35 22:37
#Fit-parm soma_diam_f 1 0.1:3
Set-parm G_Rm 25.9668 10:400
Set-parm G_cm 0.570817 0.25:2.5
Set-parm G_Ra 40.2087 10:150
#Set-parm G_e_pas -77.7604 -100:-65
#---------------------------------------------------------------------------
#Set-parm soma_cm 1.0 0:3
#Set-parm soma_Rm 40 1:100
#Set-parm soma_Ra 100 1:500
#Set-parm soma_l 17 4:30
#Set-parm soma_diam 13.5 2:30
#Set-parm aX_hill_nseg 5 1:100
#Set-parm aX_hill_l 48.7698 1:50
#Set-parm aX_iseg_nseg 2 1:100
Set-parm aX_iseg_n 2 11:50
#Set-parm aX_iseg_l 58.30 1:100
#Set-parm aX_iseg_d 0.5 0.2:1.0
#Set-parm aX_prox_nseg 11 1:100
#Set-parm aX_prox_n 3 1:50
#Fit-parm aX_prox_l 150 10:500
#Fit-parm aX_prox_d 0.5 0.2:1.0
#Set-parm aX_prox_l 150 10:500
#Set-parm aX_prox_d 0.5 0.2:1.0
#Set-parm aX_inode_n 2 1:50
#Set-parm aX_inode_d 0.4 0.1:2
#Set-parm aX_inode_l 1 0.1:5
#Set-parm aX_mnode_nseg 3 1:100
#Set-parm aX_mnode_n 3 1:100
#Set-parm aX_mnode_d 0.4 1:100
#Set-parm aX_mnode_l 100 1:1000
#Set-parm aX_cm 1.0 0:3
#Set-parm aX_Rm 40 1:100
#Set-parm aX_Ra 100 1:500
#Set-parm aXm_cm 0.04 0:1
#Set-parm aXm_Rm 100 1:200
#--------------------------------------------------------------------------------------
Fit-parm Gbar_Naf 268.743 10:300
Fit-parm s0_Gbar_Naf 377.784 50:1000
Fit-parm aXh_Gbar_Naf 5130.56 50:2000
Fit-parm aXi0_Gbar_Naf 3930.61 100:7000
Fit-parm aXi1_Gbar_Naf 508.49 100:8000
Set-parm aXr_Gbar_Naf 1000 800:2000
Set-parm aXm_Gbar_Naf 5 1:100
Set-parm gates_n_Naf_i0 3 1:3
Fit-parm vhalf_n_Naf_i0 -30 -45:-29
Fit-parm slope_n_Naf_i0 -5 -10:-2
Fit-parm tauA_n_Naf_i0 0.1 0.0001:0.9
Set-parm tauG_n_Naf_i0 0.5 0.03:0.97
Fit-parm vhalf_h_Naf_i0 -64 -75:-59
Fit-parm slope_h_Naf_i0 8 5:11
Fit-parm tauA_h_Naf_i0 10 5:50
Set-parm tauG_h_Naf_i0 0.5 0.2:0.8
Set-parm gates_n_Naf_i1 3 1:3
Fit-parm vhalf_n_Naf_i1 -38.143 -40:-35
Fit-parm slope_n_Naf_i1 -5 -10:-3
Fit-parm tauA_n_Naf_i1 0.2 0.0001:1
Set-parm tauG_n_Naf_i1 0.5 0.03:0.97
Fit-parm vhalf_h_Naf_i1 -67 -75:-62
Fit-parm slope_h_Naf_i1 8 3:14
Fit-parm tauA_h_Naf_i1 20 5:50
Set_parm tauG_h_Naf_i1 0.5 0.2:0.8
#---------------------------------------------------------------------------------
Set-parm Gbar_KDR 40.4103 10:300
Fit-parm s0_Gbar_KDR 350.697 100:800
Fit-parm aXh_Gbar_KDR 4016.54 100:2000
Fit-parm aXi0_Gbar_KDR 2163.57 500:8000
Fit-parm aXi1_Gbar_KDR 3059.39 500:7000
Set-parm aXm_Gbar_KDR 5 0.1:100
Set-parm aXr_Gbar_KDR 800 10:800
Set-parm Gbar_KDR_dslope 0 0:0.008
Set-parm gates_n_KDR_i1 4 1:4
Fit-parm vhalf_n_KDR_i1 2.98 -15:8
Fit-parm slope_n_KDR_i1 -14.72 -15:-10
Fit-parm tauA_n_KDR_i1 20 1:40
Set-parm tauG_n_KDR_i1 0.5 0.03:0.97
#---------------------------------------------------------------------------------
Fit-parm Gbar_KA 7.44833 0.1:15
Set-parm vhalf_n_KA_i1 -31 -35:5
Set-parm slope_n_KA_i1 -10 -40:-5
Fit-parm tauA_n_KA_i1 2 0.1:10
Set-parm vhalfD_h_KA_i1 -35 -60:-20
Set-parm slope_h_KA_i1 8 3:20
Fit-parm tauA_h_KA_i1 45.64 8:50
Set-parm Gbar_KA_dslope 0.008 0.0001:0.03
#---------------------------------------------------------------------------------
Set-parm Gbar_h 0.103157 0.01:4
Set-parm Gbar_h_dslope 0.008 0.0001:0.03
Set-parm gates_n_h_i0 1 1:3
Set-parm vhalf_n_h_i0 -92.07 -110:-70
Set-parm slope_n_h_i0 10 4:12
Set-parm tauA_n_h_i0 100 10:100
Set-parm tau0_n_h_i0 0.001 0.001:0.1
Set-parm tauF_n_h_i0 0 -3.5:2
Set-parm tauG_n_h_i0 0.5 0.05:0.95
Set-parm tauDv_n_h_i0 0 -55:55
#---------------------------------------------------------------------------------
Set-parm Gbar_KS1 11.82 0:150
Set-parm Gbar_KS2 82 0:150
Set-parm tauA_n_KS1 47 10:80
Set-parm tauA_n_KS2 115.43 80:200
Set-parm vhalf_n_KS1 -30 -40:-25
Set-parm vhalf_n_KS2 -30 -40:-25
Set-parm slope_n_KS1 -10 -12:-5
Set-parm slope_n_KS2 -10 -12:-5