// genesis

int number_of_neurons, selected_neuron
int i, j, sum, nsynapses, all_synapses
int nafferents = 0
str list, name
str output 
int soma_synapses, dend1_synapses, dend2_synapses, dend3_synapses


  all_synapses = 0
  list = ""
  str output = "Stell_GABAA_nsynapses_soma-dend.dat"

  ce   /molecular_layer

       soma_synapses = 0
       foreach name ({el ./Stellate[]/s#[]/GABAA})
//         echo {name}
             ce {name}
             nsynapses = {getfield . nsynapses}
             echo  {name} {nsynapses}
             soma_synapses = {soma_synapses} + {nsynapses}
       end
       list = (list) @ "soma" @ " " @ {soma_synapses} @ " " @ {chr 10}
       all_synapses = {all_synapses + soma_synapses}

  echo total number of GABAA synapses on soma of stcs = {soma_synapses}

  ce   /molecular_layer

       dend1_synapses = 0
       foreach name ({el ./Stellate[]/dend1#[]/GABAA})
//         echo {name}
             ce {name}
             nsynapses = {getfield . nsynapses}
             echo  {name} {nsynapses}
             dend1_synapses = {dend1_synapses} + {nsynapses}
       end
       list = (list) @ "dend1" @ " " @ {dend1_synapses} @ " " @ {chr 10}
       all_synapses = {all_synapses + dend1_synapses}

  echo total number of GABAA synapses on dend1 of stcs = {dend1_synapses}


  ce   /molecular_layer

      dend2_synapses = 0
       foreach name ({el ./Stellate[]/dend2#[]/GABAA})
//         echo {name}
             ce {name}
             nsynapses = {getfield . nsynapses}
             echo  {name} {nsynapses}
             dend2_synapses = {dend2_synapses} + {nsynapses}
       end
       list = (list) @ "dend2" @ " " @ {dend2_synapses} @ " " @ {chr 10}
       all_synapses = {all_synapses + dend2_synapses}

  echo total number of GABAA synapses on dend2 of stcs = {dend2_synapses}


  ce   /molecular_layer

      dend3_synapses = 0
       foreach name ({el ./Stellate[]/dend3#[]/GABAA})
//         echo {name}
             ce {name}
             nsynapses = {getfield . nsynapses}
             echo  {name} {nsynapses}
             dend3_synapses = {dend3_synapses} + {nsynapses}
       end
       list = (list) @ "dend3" @ " " @ {dend3_synapses} @ " " @ {chr 10}
       all_synapses = {all_synapses + dend3_synapses}

  echo total number of GABAA synapses on dend3 of stcs = {dend3_synapses}


  echo    {all_synapses}  
  echo  soma {soma_synapses * 1.0 / all_synapses} %  
  echo dend1 {dend1_synapses * 1.0 / all_synapses} % 
  echo dend2 {dend2_synapses * 1.0 / all_synapses} % 
  echo dend3 {dend3_synapses * 1.0 / all_synapses} % 


  echo {list}  > {output} // {source_list}