TITLE K/Cl cotransporter KCC2

NEURON {
	SUFFIX kcc2
	USEION k READ ek WRITE ik
	USEION cl READ ecl WRITE icl VALENCE -1
	RANGE ik, icl, g, acck, acccl
}

UNITS {
	(mV)	= (millivolt)
	(molar) = (1/liter)
	(mM)	= (millimolar)
	(um)	= (micron)
	(mA)	= (milliamp)
    FARADAY = (faraday) (coulomb)
    R = (k-mole) (joule/degC)
}

PARAMETER {
    : g = 0.00181 (mA/cm2) : based on 7e-5 (A/(V*cm2)) at 300K
    g = 0.0003 (mM/ms)
    surf = 615.7650056818251 (um2) : surface area of original model
    voli = 1436.8 (um3)            : volume of original model
    
}

ASSIGNED {
	ik		(mA/cm2)
	icl		(mA/cm2)
    acck    (mA/cm2)
    acccl   (mA/cm2)
  	ek      (mV)
  	ecl     (mV)
    celsius (degC)
    diam    (um)
}

INITIAL {
    acck = 0
    acccl = 0
}

BREAKPOINT {
	LOCAL cur 
	:cur = g * (FARADAY/(R * celsius)) * 0.001 * (ecl - ecl)
	cur = (g * FARADAY * voli/surf) * (FARADAY/(R * (celsius+273.15))) * (ek - ecl)
    : volume/area == diam/4.0 (um)
    : main Cl extruder in the cell
    ik = -(1e-7)*(cur)  
	icl = (1e-7)*(cur)
    acck = acck + ik
    acccl = acccl + icl
}