print "Loading Input/Output Routines for Experiment ......."
// Input/Output for Experiment:
proc WriteClampCurrents() { local TStepCount, VStepCount
sprint(iClampFileName,"ClampCurrents/%s.iCl",ExperimentName)
print " Writing ClampCurrents to ", iClampFileName //debug
wopen(iClampFileName) //m WriteOpen /ClampCurrents/ExperimentName.iCl
fprint(" ExperimentName: _%s_",ExperimentName)
fprint(" Number of recorded voltage steps:\t\t\t\t%g\n",NumVSteps)
fprint(" vclmp \t iclmp(time) \n \t")
for TStepCount=0,NumTSteps fprint("%g\t",MeasTimes[TStepCount])
fprint("\n")
for VStepCount=0,NumVSteps-1 { // writes: V I(V,t1) I(V,t2) ...
fprint("%g\t",vclmp[VStepCount])
for TStepCount=0,NumTSteps {
fprint("%g\t",iclmp[VStepCount][TStepCount])
}
fprint("\n")
} //m next v
wopen()
sprint(iClampFileName,"ClampCurrents/%s.vCl",ExperimentName)
print " Writing Clamp voltages to ", iClampFileName //debug
wopen(iClampFileName) //m WriteOpen /ClampCurrents/ExperimentName.iCl
fprint(" ExperimentName: _%s_",ExperimentName)
fprint(" Number of recorded voltage steps:\t\t\t\t%g\n",NumVSteps)
fprint(" vclmp \t iclmp(time) \n \t")
for TStepCount=0,NumTSteps fprint("%g\t",MeasTimes[TStepCount])
fprint("\n")
for VStepCount=0,NumVSteps-1 { // writes: V I(V,t1) I(V,t2) ...
fprint("%g\t",vclmp[VStepCount])
for TStepCount=0,NumTSteps {
fprint("%g\t",VClmp[VStepCount][TStepCount])
}
fprint("\n")
} //m next v
wopen()
sprint(iClampFileName,"ClampCurrents/%s.gCl",ExperimentName)
print " Writing Conductance estimate from ClampCurrents to ", iClampFileName //debug
wopen(iClampFileName) //m WriteOpen /ClampCurrents/ExperimentName.iCl
fprint("\t")
for VStepCount=0,NumVSteps-1 {
fprint("%g\t",vclmp[VStepCount])
}
fprint("\n")
for TStepCount=0,NumTSteps {
fprint("%g\t",MeasTimes[TStepCount])
for VStepCount=0,NumVSteps-1 {
if (vclmp[VStepCount]==Ek) fprint("0/t") else {
fprint("%g\t",iclmp[VStepCount][TStepCount]/(vclmp[VStepCount]-Ek))
}
}
fprint("\n")
}
wopen()
}
|