//bcell-const.g
/** Used with the terminal branches model **/
/* Units: Millimole, nanoamps, megaohms, microsiemens, milliseconds, 
   nanoFarads, Liters, centimeters */

float PI = 3.14159

float umole	=	1e-6	/* units are umoles and uM */
float mmole	=	1e-3	/* units are mmoles and mM */
int quant	=	1 /* integrate using quantity, not concentration */
int concen	=	0 /* integrate using concentration, not quantity */

/* RM*CM = ~ 35 msec => passive membrane time constant. */
float RI	=	100e-6	/* Mohm-cm, could be 200e-6 */
float RM	=	1e-2	/* Megaohms - cm^2; 5e-2 from Spruston et al, used with increased rhab SA */
float CM	=	1e3	/* nanoFarads/cm^2; 1e3 from Spruston et al. */

float gshunt	=	0.005
float Er	=	-60
float Vinit	=	-60

float gleak	=	250	/* produces ~0.029 uS in rhab (consistent with IKlight)*/
float gna	=	5e5
float g_ih	=	2000	/* produces 25nS in entire cell? */
float pca_p	=	0.1	/* consider 0.2 instead of 0.1*/
float pca_t	=	0.1
float gkca	=	64e3
float gka	=	100e3	/* uS/cm^2; consider 150e3 */
float g_gabaa	=	70e4	/* gamma*(mol/cm^2) */
float g_gabab	=	70e3	/* uS/cm^2 */

float somarad	=	10.4e-4 /* e-4 converts from microns to cm */
float somalen	=	24e-4	/* 24 microns */
int somacyls	=	24 
int somashells	=	2

float rhabrad	=	6.17e-4 
float rhabcorerad=  1.0e-4
float rhablen	=	12e-4
float rhabSA	=	0.000125  /*Assuming 5000 microvilli. 0.08 um radius X 10 um len or 0.08 diam x 5 um len */
float rhabxarea	=	3.14e-8  /* Core cross section assuming 2e-4 dia non-villi portion */
int rhabcyls	=	12
int rhabshells	=	2
/*float rhabvillen	=	5.17e-4  length of microvilli used prior to 08/16/02 */
float rhabvillen	=	10e-4	/* actual length, Eakin et al.*/
/*float rhabvilrad	=	0.1e-4 	 radius of microvil used prior to 08/16/02 */
float rhabvilrad	=	0.04e-4 /* actual radius of microvil, Eakin et al. */
int numvilli	=	5000 		/* number of microvilli */
float shellsize	=	1e-4

float neckrad	=	3e-4	/* 3 microns */
float necklen	=	1e-4	/* 1 micron neck width */
float axonlen	=	100e-4	/* 100 micron length */
/*float axonrada  =       1e-4    1 micron short axis  prior to 08/16/02*/
float axondiama	=	1.5e-4	/* 1.5 micron short axis diameter! */
float axondiamb	=	3e-4	/* 3 micron long axis  diameter!*/
int axonslice	=	4	/* four axon  voltage compartments */
int axoncyls	=	100

float syn_br_rad	=	0.63e-4		/*0.215e-4  for asym*/
float nosyn_br_rad	=	0.63e-4		/*0.93e-4 for asym */
float branchlen1	=	5e-4	/* branch between synapse & axon*/
float branchlen2	=	10e-4	/* terminal branch with synapse*/
int branchcyls		=	15

float Cacyt     =       0.11e-3      /* mM, for ca=0.0100 */
float CaER      =       0.020         /* mM for ca=0.0100 */
float bufcyt    =       0.1497       /* mM for ca=0.0100*/
float bufER     =       2.399           /* mM, for ca=0.0100 */

float Cadif     =       6.0e-9  /* cm^2/msec */
float bufcyttot =       153e-3          /* mM */
float bufERtot  =       12.0            /* mM, 2x previous value */
float buf_kf    =       1e2     /* per mM-msec, from Nowycky et al. and */
float buf_kb    =       0.5     /* per msec, from Blumenfeld et al. */

float init00	=	0.7434	/*  for ca=0.0110 */
float init10	=	0.1615	/*  for ca=0.0110 */
float init01	=	0.0781	/*  for ca=0.0110 */
float init000	=	0.325	/*  for ca=0.0110 */
float init100	=	0.0	/*  for ca=0.0110 */
float init010	=	0.435	/*  for ca=0.0110 */
float init001	=	0.103	/*  for ca=0.0110 */
float init101	=	0.0	/*  for ca=0.0110 */
float init011	=	0.137	/*  for ca=0.0110 */
float init110	=	0.0	/*  for ca=0.0110 */

float maxiicr   =       0.16             /* Units are /msec */
int iicrpower	=	3
float maxcicr	=	0.08
int cicrpower	=	1
float serca     =       0.00047           /* Units are mmole/msec */
float pumppower	=	2
float ERfactor	=	0.185	/*0.185 ratio of ER volume to cyt vol */

float ip3dif    =       2.83e-9 /* cm^2/msec - Allbritton et al. 1992  */
float ip3init   =       0.0e-3  /* */
float ip3degrad =       0.7e-3 /* IP3 degradation - per msec - Allbritton et al. 1992*/

float Rhodtot	=	10000	/* uM*/
float krhoF	=	1e-6	/* per msec */
int duration	=	3000	/* msec*/
float intensity	=	10.0	/* isomerizations per msec */
			/* 0.01 = ND3, 0.1 = ND2, 1.0 = ND1, 10.0 = ND0 */
int lightdelay	=	10	/*msec*/
float deplete_power	=	0.16	/* decrease in mrho effectiveness = 1/ time^deplete_factor */


float RKArrtot	=	30	/* uM */
float Krkf	=	0.5e-3	/* per msec - uM */
float Krkb	=	5.0e-3	/* per msec */
float Krkcat	=	5e-3	/* per msec */

float Gtot	=	1000	/* uM */
float Kgf	=	6e-3	/* per msec - uM */
float Kgb	=	10e-3	/* per msec */
float Kgcat	=	1920e-3	/* per msec */
float Khyd	=	0.0057e-3  /* per msec */

float Kplcf	=	100e-3 /* per uM per msec */
float Kplcb	=	0.2e-3 /* per msec */
float plctot	=	100	/* uM */
float piptot	=	160	/* uM */
float Kpif	=	0.83e-3	/* per uM per msec */
float Kpib	=	0.1e-3	/* per msec */
float Kpicat1   =       10.0e-3   
float Kpicat2	=	6.0e-3  /* per msec */
float Kgap	=	3e-3	/* per msec */

float Gprot_syn	=	100
float plc_syn	=	10
float gabab_kf 	=	0.06	/* modified to produce EC50 of 1.5 uM */
float gabab_kb 	=	0.5
float recept_tot 	=	1.0	/* unitless or uM */
float gabab_gf	=	2.0
float gabab_gb	=	0.5		/* g5 modified to slow down peak curren
t, g4 not significant */
float gabab_gcat =	0.5