TITLE KATP channel based on model from Knowlton et al 2018
 
                 
NEURON {
	SUFFIX katp
	USEION k WRITE ik
  USEION adp READ adpi VALENCE 0
	RANGE  km, gkbar, ik, sinf, n, gk
}

INDEPENDENT {v FROM -100 TO 50 WITH 50 (mV)}

UNITS {
	(molar) = (1/liter)
	(mM) =  (millimolar)
	(um) =  (micron)
	(mA) =  (milliamp)
	FARADAY = (faraday) (coulomb) 
	PI = (pi) (1)
}

PARAMETER {
         dt (ms)
         adpi (mM)
         km = 0.0077 (mM)
         n = 1.3 (1)
         gkbar = 20e-6 (S/cm2)
         ek  =  -90.0 (mV)
         }

ASSIGNED {
  ik  (mA/cm2)
  gk  (mS/cm2)
  sinf
}


BREAKPOINT {
   if(adpi > 1e-6){
     sinf = 1.0/(1.0 +  pow(km/adpi,n))
   } else {
     sinf = 0
   }
   gk = gkbar*sinf
   ik = gk*(v-ek)
}