<PSICSRun   lib="." timeStep="0.01ms" runTime="40ms" startPotential="-80mV"
	morphology="cell"
	environment="environment"
	properties="membrane"
	access="recording-iv"
	stochThreshold="0">

	<StructureDiscretization baseElementSize="3um"/>

<RunSet vary="napop:density" values="[16]per_um2"
                filepattern="out-$.txt">
 		<CommandSet vary="vstep:to" values="[-60, -40, -30, -20, -10, 0, 10, 20, 30, 40]mV"/>
	</RunSet>



<About>
  	I-V curves for the same model as is used for mean variance analysis.
 </About>



	<ViewConfig>
		<LineGraph width="500" height="400">
		<XAxis min="0" max="250" label="time / ms"/>
		<YAxis min="-80" max="60" label="current / nA"/>
		<LineSet file="out-16-f.txt" color="green"/>
		 <View id="16" xmin="0." xmax="40." ymin="-500." ymax="20."/>
		</LineGraph>
	</ViewConfig>
</PSICSRun>

<CellMorphology id="cell">

    <Point id="p0" x="0" y="0" z="0" r="1."/>

</CellMorphology>

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

	<ChannelPopulation id="napop" channel="na1" density="12per_um2"/>


</CellProperties>

<CellEnvironment id="environment">
	<Ion id="Na" name="Sodium" reversalPotential="80mV"/>
	<Ion id="K" name="Potassium" reversalPotential="-80mV"/>
</CellEnvironment>

<Access id="recording-iv" saveInterval="0.01ms">

	<VoltageClamp at="p0" lineColor="red" hold="-80.0mV">
		<VoltagePulse id="vstep" start="5ms" to="0mV"
			      duration="20ms"/>

	</VoltageClamp>

</Access>

<KSChannel id="na1" gSingle="30pS" permeantIon="Na">
	<ClosedState id="c1"/>
	<OpenState id="o1"/>
	<ClosedState id="c2"/>

	<VHalfTransition from="c1" to="o1" vHalf = "-35mV" z="4.5e" gamma="0.8" tau="0.55ms" tauMin="0.1ms"/>

	<VHalfTransition from="o1" to="c2" vHalf = "-70mV" z="1.1e" gamma="0.90" tau="12.0ms" tauMin="1.5ms"/>

</KSChannel>