<PSICSRun   timeStep="0.02ms" runTime="1000ms" startPotential="-65mV"
	morphology="cell"
	environment="environment"
	properties="membrane"
	access="recording"
	stochThreshold="1000000"
	oneByOne="true"
	quickChannels="true">

<RunSet vary="p0:r" values="[0.89206, 1.99471, 2.8209, 3.98942]um" filepattern="radius-$"/>

	<info>Chow and White (1996), Biophys J. Showing spontaneous spiking with varying numbers of channels</info>

	<ViewConfig>
		<LineGraph width="500" height="400">
		<XAxis min="0" max="250" label="time / ms"/>
		<YAxis min="-80" max="60" label="potential / mV"/>


		<LineSet file="radius-0.89206.txt" color="orange"/>
		<LineSet file="radius-1.99471.txt" color="blue"/>
		<LineSet file="radius-2.8209.txt" color="green"/>
		<LineSet file="radius-3.98942.txt" color="yellow"/>
	<!--
		<LineSet file="radius-5.64189.txt" color="red"/>
	-->

		<View id="whole" xmin="-10." xmax="10000." ymin="-100." ymax="80."/>
		<View id="start"   xmin="0." xmax="500." ymin="-100." ymax="80."/>
		 <View id="end"   xmin="500." xmax="600." ymin="-100." ymax="80."/>

		</LineGraph>
	</ViewConfig>
</PSICSRun>

<CellMorphology id="cell">
	<Point id="p0" x="0" y="0" z="0" r="5.4"/>
</CellMorphology>

<CellProperties id="membrane"
      cytoplasmResistivity="100ohm_cm"
      membraneCapacitance="1uF_per_cm2">

	<ChannelPopulation channel="NaHH" density="60per_um2"/>
    <ChannelPopulation channel="KHH" density="18per_um2"/>
    <ChannelPopulation channel="Leak" density="300per_um2"/>

</CellProperties>

<CellEnvironment id="environment">
	<Ion id="Na" name="Sodium" reversalPotential="50mV"/>
	<Ion id="K" name="Potassium" reversalPotential="-77mV"/>
	<Ion id="LEAK" name="Non-specific leak" reversalPotential="-54.4mV"/>
</CellEnvironment>

<Access id="recording">

    <VoltageRecorder at="p0" lineColor="blue"/>

</Access>

<KSChannel id="NaHH" permeantIon="Na" gSingle="20pS">

	<KSComplex id="m" instances="3">
		<ClosedState id="c"/>
		<OpenState id="o"/>
		<ExpLinearTransition from="c" to="o" rate="1per_ms" midpoint ="-40mV" scale="10mV"/>
		<ExpTransition from="o" to = "c" rate="4per_ms" midpoint="-65mV" scale="-18mV"/>
	</KSComplex>

	<KSComplex id="h">
		<ClosedState id="c"/>
		<OpenState id="o"/>
		<ExpTransition from="c" to="o" rate="0.07per_ms" midpoint="-65mV" scale="-20mV"/>
		<SigmoidTransition from="o" to="c" rate="1per_ms" midpoint="-35mV" scale="10mV"/>
	</KSComplex>

</KSChannel>

<KSChannel id="KHH" permeantIon="K" gSingle="20pS">

	<KSComplex id="n" instances="4">
		<ClosedState id="c"/>
		<OpenState id="o"/>
		<ExpLinearTransition from="c" to="o" rate="0.1per_ms" midpoint ="-55mV" scale="10mV"/>
		<ExpTransition from="o" to = "c" rate="0.125per_ms" midpoint="-65mV" scale="-80mV"/>
	</KSComplex>


</KSChannel>

<KSChannel id="Leak" gSingle="0.01pS" permeantIon="LEAK">
	<OpenState id="o1"/>
</KSChannel>