<PSICSRun timeStep="0.05ms" runTime="100ms" startPotential="-65mV"
	environment="environment"
	properties="membrane"
	access="recording"
	morphology = "cell"
	repeats="20"
	stochThreshold="0">

	<StructureDiscretization baseElementSize="25um"/>

<!--
    <ChannelStochThreshold channel="km" threshold="0"/>
    <ChannelStochThreshold channel="kna_nature" threshold="0"/>
	<ChannelStochThreshold channel="kvz_nature" threshold="0"/>
 -->

	<RunSet vary="stochThreshold" values="[0, 100000]"
		filepattern="out-$"/>

	<info></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="out-100000.txt" color="blue"/>
		<LineSet file="out-0.txt" color="red"/>


		 <View id="whole" xmin="-10." xmax="50." ymin="-100." ymax="80."/>
		  <View id="enlarged" xmin="-10." xmax="50." ymin="-70." ymax="-50."/>
		</LineGraph>
	</ViewConfig>


</PSICSRun>

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


</CellMorphology>

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


    <ChannelPopulation channel="leak-Na" density="1per_um2"/>
	<ChannelPopulation channel="leak-K" density="1per_um2"/>

	<ChannelPopulation channel="naz_nature" density="30per_um2"/>
    <ChannelPopulation channel="kvz_nature" density="18per_um2"/>


	<ChannelPopulation channel="km" density="0.1per_um2"/>



    <DensityAdjustment maintain="-65mV" vary="leak-Na, leak-K"/>

</CellProperties>

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

<Access id="recording">


	<ConductanceClamp at="p0" lineColor="red" hold="0.0nS" potential="-50mV">
	     <TimeSeries id="stim" file="stim1.txt" column="1" timeUnit="1ms" valueUnit="10pS"
	     	scaleFactor="1.0"/>
	</ConductanceClamp>


</Access>

<KSChannel id="leak-Na" gSingle="2pS" permeantIon="Na">
	<OpenState id="o1"/>
	<ClosedState id="c1"/>
	<FixedRateTransition from="o1" to="c1" forward="3.per_ms" reverse="7.per_ms"/>
</KSChannel>

<KSChannel id="leak-K" gSingle="2pS" permeantIon="K">
	<OpenState id="o1"/>
	<ClosedState id="c1"/>
	<FixedRateTransition from="o1" to="c1" forward="3.per_ms" reverse="7.per_ms"/>
</KSChannel>

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

	<About>
		Converted from naz_nature.mod (Zack Mainen, 1994)
	</About>

	<KSComplex id="m" instances="3">
		<ClosedState id="c"/>
		<OpenState id="o"/>
		<ExpLinearTransition from="c" to="o" rate="1.64per_ms" midpoint ="-35.mV" scale="9mV"
			baseTemperature="23Celsius" q10="2.3"/>
		<ExpLinearTransition from="o" to="c" rate="1.12per_ms" midpoint="-35.mV" scale="-9mV"
		  baseTemperature="23Celsius" q10="2.3"/>
	</KSComplex>


	<CodedTransitionFunction name="trap0" returnVariable="rate" type="double">
		<Argument name="v" type="double"/>
		<Argument name="th" type="double"/>
		<Argument name="a" type="double"/>
		<Argument name="q" type="double"/>
		<![CDATA[
			if (Math.abs(v - th) > 1.e-6) {
				rate =  a * (v - th) / (1 - Math.exp(-(v - th)/q));
			} else {
				rate = a * q;
			}
		]]>
	</CodedTransitionFunction>



	<KSComplex id="h">
		<ClosedState id="c"/>
		<OpenState id="o"/>


 	<TauInfCodedTransition from="c" to="o" tauvar="htau" infvar="hinf"
        	baseTemperature="23Celsius" q10="2.3">
 		<Constant name="thi1" value="-50"/>
 		<Constant name="thi2" value="-75"/>
 		<Constant name ="rd" value="0.024"/>
 		<Constant name="rg" value="0.0091"/>
 		<Constant name="qi" value="5"/>

 		<Constant name="thinf" value="-65"/>
 		<Constant name="qinf" value="6.2"/>
 	<![CDATA[

        double a = trap0(v, thi1, rd, qi);
        double b = trap0(-v, -thi2, rg, qi);
        htau = 1/(a+b);
        hinf = 1. / (1 + Math.exp((v - thinf)/qinf));
	]]>
 	</TauInfCodedTransition>


	</KSComplex>

</KSChannel>

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

	<About>
		Converted from kvz_nature.mod (Zach Mainen, 1995)
	</About>

	<KSComplex id="n" instances="1">
		<ClosedState id="c"/>
		<OpenState id="o"/>
		<ExpLinearTransition from="c" to="o" rate="0.18per_ms" midpoint ="25.mV" scale="9mV"
		baseTemperature="23Celsius" q10="2.3"/>
		 <ExpLinearTransition from="o" to="c" rate="0.018per_ms" midpoint ="25.mV" scale="-9mV"
		baseTemperature="23Celsius" q10="2.3"/>

	</KSComplex>


</KSChannel>

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

	<About>
		Converted from km.mod, slow non-inactivating potassium current
		based on I-M (Zach Mainen, 1995)
	</About>

	<KSComplex id="n" instances="1">
		<ClosedState id="c"/>
		<OpenState id="o"/>
		<ExpLinearTransition from="c" to="o" rate="0.09per_ms" midpoint ="-30.mV" scale="9mV"
		baseTemperature="23Celsius" q10="2.3"/>
		 <ExpLinearTransition from="o" to="c" rate="0.009per_ms" midpoint ="-30.mV" scale="-9mV"
		baseTemperature="23Celsius" q10="2.3"/>

	</KSComplex>


</KSChannel>