//genesis
float PI = 3.14159
// =====================================================================
// PHYSIOLOGICAL PARAMETERS
// =====================================================================
// All units in SI and usually dimensions are given as well
// IONIC EQUILIBRIUM POTENTIALS
float ENA = 55.0e-3 // V
float ECL = -65.0e-3 // V
float EK = -90.0e-3 // V
float EREST = -70.0e-3 // V
float I_EREST = -70.0e-3 // V
// ACTIVE COMPARTMENT POTENTIALS
float EREST_ACT = -70.0e-3 // V
float ENA_ACT = 115.0e-3 + EREST_ACT // V
float EK_ACT = -12.0e-3 + EREST_ACT // V
float ELEAK_ACT = 10.6e-3 + EREST_ACT // V
// squid giant axon - Conti et al.(1975),J.Physiol(Lond),248,45-82
// gNa = 4pS
// squid giant axon - Conti and Neher(1980),Nature(Lond),285,140-143
// gK = 18pS
// mouse spinal neurone - Mathers and Barker(1982),Int.Rev.Neurobiol,23,1-34
// gCl = 18pS
// UNIT CONDUCTANCES
// (Siemens/channel)
float UNIT_GNA = 12.0e-12 // S
float UNIT_GCL = 18.0e-12 // S
float UNIT_GK = 4.0e-12 // S
float SUNIT_GNA = 8.0e-12 // S
float SUNIT_GK = 4.0e-12 // S
// squid giant axon - Conti et al.(1975),J.Physiol(Lond),248,45-82
// rhoNa = 330 channels/um^2
// squid giant axon - Conti and Neher(1980),Nature(Lond),285,140-143
// rhoK = 72 channels/um^2
// CHANNEL DENSITIES
// (channels/m^2)
float RHO_NA = 60.0e12
float RHO_CL = 30.0e12
float RHO_K = 30.0e12
float SRHO_NA = 330.0e12
float SRHO_K = 72.0e12
float IRHO_NA = 10.0e12
// SYNAPTIC AREA
// Haberly and Presto(1986),J.Comp.Neurol,248,464-474
// basal spine d = .40um
// distal spine d = .74um
// presynaptic to dendritic spines d = .61
// presynaptic to dendritic shafts d = .89
// presynaptic to initial segment d = .89
//
// using synaptic contact area A = pi*d^2/4
// (m^2/synapse)
float ASYN_LOCAL_NA = 0.12e-12
float ASYN_DISTAL_NA = 0.43e-12
float ASYN_CL = 0.62e-12
float ASYN_K = 0.43e-12
float IASYN_NA = 0.29e-12
// ACTIVE AREA
// fraction of somatic area containing active channels
float fAC = 0.1
// PEAK CONDUCTANCE
// (Siemens/synapse)
float I_GMAX_NA = UNIT_GNA * IRHO_NA * IASYN_NA
float LOCAL_GMAX_NA = UNIT_GNA * RHO_NA * ASYN_LOCAL_NA
float DISTAL_GMAX_NA = UNIT_GNA * RHO_NA * ASYN_DISTAL_NA
float GMAX_K = UNIT_GK * RHO_K * ASYN_K
float GMAX_CL = UNIT_GCL * RHO_CL * ASYN_CL
float SGMAX_NA = SUNIT_GNA * SRHO_NA // S/m^2
float SGMAX_K = SUNIT_GK * SRHO_K // S/m^2
// MEMBRANE PARAMETERS
float RM = 0.2 // ohm-m^2
float RA = 0.5 // ohm-m
float CM = 0.01 // F/m^2
// TYPICAL CELL DIMENSIONS
float SOMA_D = 20.0e-6 // m
float SOMA_L = 20.0e-6 // m
float SOMA_A = PI * SOMA_D * SOMA_L // m^2
float SOMA_XA = PI * SOMA_D * SOMA_D / 4 // m^2
float DEND_D = 1.5e-6 // m
float DEND_L = 100.0e-6 // m
float DEND_A = PI * DEND_D * DEND_L // m^2
float DEND_XA = PI * DEND_D * DEND_D / 4 // m^2
/*
** VARIABLES USED BY ACTIVE COMPONENTS
*/
int EXPONENTIAL = 1
int SIGMOID = 2
int LINOID = 3
|