// LEAK-file, fit.par is loaded first to get appropriate time & voltage scales
// this file must be run first (and once per cell) to produce *.iL -leak current file for fit
strdef iClampFileName, loadProg
objref vC, sh, st
load_file("stdgui.hoc")
{load_proc("nrnmainmenu") nrnmainmenu()}
// VARIABLE CHECK: //
print " utilized Timesteps: "
for TStepCount=0, NumTSteps print TStepCount, ": " , MeasTimes[TStepCount], "ms"
print " utilized VoltageSteps: "
for VStepCount=0, NumVSteps-1 {
print VStepCount, ": " , vclmp[VStepCount]}
print "+++++ Experiment Name: ", ExperimentName, "++++++++"
print " "
print " Passive Parameters: (homogenous)"
print "Rm: ",rm
print "Ra: ",ra
print "c_m: ",c_m
print "Ek: ",Ek
print "v_init: ",v_init
// CELL (loaded by fit.par already):
forall {
// passive:
insert pas
Ra = ra
cm = c_m
g_pas = 1/rm
e_pas = Epas
// no active
}
// reducing cell model:
if (ReduceCell) {
xopen("ReduceModel.cell")
ReduceModel()
}
// VOLTAGE CLAMP PARAMETERS //
// all explicitly set by Experiment.par
stimSec.sec vC = new TwoElectrodeIdealVClamp() // perfect voltage clamp at soma
//vC.dur1 = MeasTStart + MaxMeasTime*2 // factor 2 is for safety reasons
// Amplitudes will be vector-played
//vC.rs = VClampSeriesR
stimSec.sec vC.xmeasure(stimPosV)
stimSec.sec vC.xinject(stimPosI)
xopen("Routines.leak") // Loading Clamp-Routines
xopen("InOut.leak") // Loading Input/Output-Routines
// Main Program:
CalculateLeak()
WriteLeakCurrents()
|