// genesis
// allows one to document the number of electrical synapses to selected neuron types
// rmaex  March 2017

int number_of_neurons, selected_neuron
int msgcountin, msgindex 
int i, j, sum, nsynapses
int nmsgs = 0
str list, name
str output 



  sum = 0
  number_of_neurons = 800 // MLIs
  list = ""
  str output = "Stell_el_nsynapses_torus.dat"

  for (i = 1; {i <= number_of_neurons}; i = i + 1)

       list = (list) @ {i} @ " "
       ce /molecular_layer/Stellate[{i-1}]
       nsynapses = 0
       foreach name ({el ./d#[]})
             ce {name}
             nmsgs = {getmsg . -out  -count}
             for (msgindex = 0; {msgindex < nmsgs}; msgindex = msgindex + 1)
                if ({! {strcmp {getmsg . -out -type {msgindex}} "PLUS"}})
                  nsynapses = {nsynapses} + 1
                end // if
             end // for
       echo  Stellate cell  {i}   {nsynapses}
       end

       list = (list) @ {nsynapses} @ " " @ {chr 10}
       sum = {sum} + {nsynapses}
  end
  echo total number of electrical synapses on stcs = {sum}
  echo {list}  > {output} // {source_list}