TITLE kir
COMMENT
potassium inward rectifier after geukes foppen en siegenbeek
simply, with no time dynamics. (instantanous kinetics)
ENDCOMMENT

NEURON {
	SUFFIX kir
	USEION k READ ko, ek WRITE ik
	RANGE gbar, ik, qk, acck
	GLOBAL vs
}

UNITS {
	(molar) = (1/liter)
	(mV) =	(millivolt)
	(mA) =	(milliamp)
	(mM) =	(millimolar)
}

PARAMETER {
	gbar	= 0	(mho/cm2)
	vs	= 1	:voltage-dependency factor in boltzmann equation was 4
	vh	= 0 (mV)	:halfmaximal activation relative to ek was -10
}

ASSIGNED { 
	ik	(mA/cm2)
    acck (mA/cm2)
	v	(mV)
	ek	(mV)
	ko 	(mM)
}

BREAKPOINT {
	ik = gbar*pkir((v), ko, ek)*(v-ek)
    acck = acck + ik
}

INITIAL {
    acck = 0
}

FUNCTION pkir(vm,kout,erev) {
	pkir = 1 / ( sqrt(kout) * (1+exp((vm-vh-erev)/vs)) )
}