/* Model version description

*/

// Quick simulation
cvode.active(1)
cvode.use_daspk(1)

// Temperature

celsius = 25

// Soma

Granule[0].soma.gkbar_GRC_KIR = 0.002
Granule[0].soma.gl_GRC_LKG1 = 0.0000916
Granule[0].soma.gkbar_GRC_KA = 0.04
Granule[0].soma.gkbar_GRC_KM = 0.006
Granule[0].soma.gkbar_GRC_KCA = 0
Granule[0].soma.gnabar_GRC_NA = 0

// Dendrites
for i = 0,3 {Granule[0].dend_1[i].gkbar_GRC_KCA(0.5) = 0.0014}
for i = 0,3 {Granule[0].dend_2[i].gkbar_GRC_KCA(0.5) = 0.0014}
for i = 0,3 {Granule[0].dend_3[i].gkbar_GRC_KCA(0.5) = 0.0014}
for i = 0,3 {Granule[0].dend_4[i].gkbar_GRC_KCA(0.5) = 0.0014}

// Hillock

for i = 0,4 {Granule[0].hillock[i].gl_GRC_LKG1(0.5) = 0} 
for i = 0,4 {Granule[0].hillock[i].gnabar_GRC_NA(0.5) = 0.13}
for i = 0,4 {Granule[0].hillock[i].gkbar_GRC_KV(0.5) = 3}
for i = 0,4 {Granule[0].hillock[i].Aepsilon_GRC_NA(0.5) = 0}
for i = 0,4 {Granule[0].hillock[i].Aalpha_n_GRC_KV(0.5) = -0.00222}
for i = 0,4 {Granule[0].hillock[i].Abeta_n_GRC_KV(0.5) = 0.0285}
for i = 0,4 {Granule[0].hillock[i].Aalfa_GRC_NA(0.5) = 8.85}
for i = 0,4 {Granule[0].hillock[i].Valfa_GRC_NA(0.5) = 10}
for i = 0,4 {Granule[0].hillock[i].Abeta_GRC_NA(0.5) = 0.0318}
for i = 0,4 {Granule[0].hillock[i].Vbeta_GRC_NA(0.5) = 10}
for i = 0,4 {Granule[0].hillock[i].ACon_GRC_NA(0.5) = 0.002}
for i = 0,4 {Granule[0].hillock[i].ACoff_GRC_NA(0.5) = 0.24}
for i = 0,4 {Granule[0].hillock[i].AOon_GRC_NA(0.5) = 1.05}
for i = 0,4 {Granule[0].hillock[i].AOoff_GRC_NA(0.5) = 0.004}
for i = 0,4 {Granule[0].hillock[i].Vshift_GRC_NA(0.5) = -22}
for i = 0,4 {
    Granule[0].hillock[i].use_threshold_GRC_NA(0.5) = 0
    Granule[0].hillock[i].V_threshold_GRC_NA(0.5) = -72
}
for i = 0,4 {Granule[0].hillock[i].Ra = 500}


// ++++++++ Ascending axon ++++++++


for i = 0,2 {Granule[0].axon[i].gkbar_GRC_KV(0.5) = 10}
for i = 0,2 {Granule[0].axon[i].Aalpha_n_GRC_KV(0.5) = -0.00222}
for i = 0,2 {Granule[0].axon[i].Abeta_n_GRC_KV(0.5) = 0.0285}
for i = 0,2 {Granule[0].axon[i].gnabar_GRC_NA(0.5) = 0.75}
for i = 0,2 {Granule[0].axon[i].Aalfa_GRC_NA(0.5) = 8.85}
for i = 0,2 {Granule[0].axon[i].Valfa_GRC_NA(0.5) = 10}
for i = 0,2 {Granule[0].axon[i].Abeta_GRC_NA(0.5) = 0.0318}
for i = 0,2 {Granule[0].axon[i].Vbeta_GRC_NA(0.5) = 10}
for i = 0,2 {Granule[0].axon[i].ACon_GRC_NA(0.5) = 0.002}
for i = 0,2 {Granule[0].axon[i].Aepsilon_GRC_NA(0.5) = 0}
for i = 0,2 {Granule[0].axon[i].ACoff_GRC_NA(0.5) = 0.24}
for i = 0,2 {Granule[0].axon[i].AOon_GRC_NA(0.5) = 1.05}
for i = 0,2 {Granule[0].axon[i].AOoff_GRC_NA(0.5) = 0.004}
for i = 0,2 {Granule[0].axon[i].Vshift_GRC_NA(0.5) = -22}
for i = 0,2 {
    Granule[0].axon[i].use_threshold_GRC_NA(0.5) = 0
    Granule[0].axon[i].V_threshold_GRC_NA(0.5) = -72
}

for i = 0,2 {Granule[0].axon[i].Ra = 500}
for i = 0,2 {Granule[0].axon[i].gl_GRC_LKG3(0.5) = 0.001}
for i = 0,2 {Granule[0].axon[i].el_GRC_LKG3(0.5) = -74.5}

// FHF WT


for i = 3,55 {Granule[0].axon[i].Aalfa_GRC_NA(0.5) = 8.85}
for i = 3,55 {Granule[0].axon[i].Valfa_GRC_NA(0.5) = 10}
for i = 3,55 {Granule[0].axon[i].Abeta_GRC_NA(0.5) = 0.0318}
for i = 3,55 {Granule[0].axon[i].Vbeta_GRC_NA(0.5) = 10}
for i = 3,55 {Granule[0].axon[i].ACon_GRC_NA(0.5) = 0.002}
for i = 3,55 {Granule[0].axon[i].Aepsilon_GRC_NA(0.5) = 0}
for i = 3,55 {Granule[0].axon[i].AOon_GRC_NA(0.5) = 1.05}
for i = 3,55 {Granule[0].axon[i].ACoff_GRC_NA(0.5) = 0.24}
for i = 3,55 {Granule[0].axon[i].AOoff_GRC_NA(0.5) = 0.004}
for i = 3,55 {Granule[0].axon[i].Vshift_GRC_NA(0.5) = -22}
for i = 3,55 {Granule[0].axon[i].gnabar_GRC_NA(0.5) = 0.0135}
for i = 3,55 {
    Granule[0].axon[i].use_threshold_GRC_NA(0.5) = 0
    Granule[0].axon[i].V_threshold_GRC_NA(0.5) = -60
}

for i = 3,55 {Granule[0].axon[i].gkbar_GRC_KV(0.5) = 0.0185}
for i = 3,55 {Granule[0].axon[i].Abeta_n_GRC_KV(0.5) = 0.0285}
for i = 3,55 {Granule[0].axon[i].Aalpha_n_GRC_KV(0.5) = -0.00222}
for i = 3,55 {Granule[0].axon[i].Ra = 100}
for i = 3,55 {Granule[0].axon[i].gl_GRC_LKG3(0.5) = 0.000005}
for i = 3,55 {Granule[0].axon[i].el_GRC_LKG3(0.5) = -74.5}

// ++++++  PFs ++++++++

for i = 0,499 {Granule[0].pf[i].L = 10}
for i = 0,499 {Granule[0].pf[i].nseg = 5}

// FHF WT


for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].Aepsilon_GRC_NA(s/10) =  0}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].Aalfa_GRC_NA(s/10) =  8.85}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].Valfa_GRC_NA(s/10) =  10}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].Abeta_GRC_NA(s/10) =  0.0318}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].Vbeta_GRC_NA(s/10) =  10}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].ACon_GRC_NA(s/10) =  0.002}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].AOon_GRC_NA(s/10) =  1.05}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].ACoff_GRC_NA(s/10) = 0.24}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].AOoff_GRC_NA(s/10) = 0.004}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].Vshift_GRC_NA(s/10) = -22}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].gnabar_GRC_NA(s/10) = 0.0135}}
for i = 0,499 {
    for (s=1;s<10;s=s+2) {
        Granule[0].pf[i].use_threshold_GRC_NA(s/10) = 0
        Granule[0].pf[i].V_threshold_GRC_NA(s/10) = -60
        }
 }

for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].gkbar_GRC_KV(s/10) = 0.0185}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].Abeta_n_GRC_KV(s/10) = 0.0285}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].Aalpha_n_GRC_KV(s/10) = -0.00222}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].Ra = 100}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].gl_GRC_LKG3(s/10) = 0.000005}}
for i = 0,499 {for (s=1;s<10;s=s+2) {Granule[0].pf[i].el_GRC_LKG3(s/10) = -74.5}}