// genesis

int number_of_neurons, selected_neuron, number_of_afferents
int msgcount, msgindex 
int i, j, sum, nsynapses, nsynaps
int totalnsynapses = 0
int nafferents 
int totalnafferents = 0
float meantotalnsynapses,meantotalnafferents
str list, name, target, sender 
str output 




  number_of_neurons = 800 // 84 // 50 // PCs
  number_of_afferents = 800 // stcs
  list = ""
  output = "Stc_stc_newtable.dat"
  population_sum = 0

//  for (i = 1; {i <= number_of_neurons}; i = i + 1)
//       list = (list) @ {i} @ " "
//  end

  list = (list) @ " " @ {chr 10}

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

       nsynapses = 0
       nafferents = 0
       for (j = 1; {j <= number_of_afferents}; j = j + 1)
            ce /molecular_layer/Stellate[{j-1}]/soma/spike
            msgcount = {getmsg . -outgoing -count}  
            echo {i} {j} {msgcount}
           sum = 0

           foreach name ({el /molecular_layer/Stellate[{i-1}]/##[][TYPE=compartment]/GABAA})
                   nsynaps = {getsyncount /molecular_layer/Stellate[{j-1}]/soma/spike  {name}}
                   sum = {sum} + {nsynaps}
           end

//           list = (list) @ {sum} @ " "
           echo {i} {j} {sum}
           if ({sum > 0}) 
                nafferents = {nafferents + 1}
           end
           totalnsynapses = {totalnsynapses + sum}
           nsynapses = {nsynapses + sum}

       end

       list = (list) @ " " @ {nsynapses} @ " " @ {nafferents} @ " " @ {chr 10}
       totalnafferents = {totalnafferents + nafferents}
  end
  echo totalnafferents {totalnafferents}
  meantotalnafferents = {totalnafferents * 1.0 / number_of_neurons}
  echo mean totalnafferents {meantotalnafferents}
  list = (list) @ " " @ {meantotalnafferents} @ " " @ {chr 10}
  echo totalnsynapses {totalnsynapses}
  meantotalnsynapses = {totalnsynapses * 1.0 / number_of_neurons}
  echo mean totalnsynapses {meantotalnsynapses}
  list = (list) @ " " @ {meantotalnsynapses} @ " " @ {chr 10}
  echo {list}  > {output}