// genesis
// allows one to document how many autapses each MLI receives
// rmaex  March 2017


int number_of_neurons
int msgcount, msgindex 
int i, j, sum, nautapses, nsynaps
int totalnautapses = 0
int nneurons_with_autapses = 0 
str list, name, target, sender 
str output 



  number_of_neurons = 800 // MLIs
  list = ""
  output = "Stc_stc_nautapses.dat"

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

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

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

           echo {i} {nautapses}
           if ({nautapses > 0}) 
                nneurons_with_autapses = {nneurons_with_autapses + 1}
           end // if
           totalnautapses = {totalnautapses + nautapses}

       list = (list) @ " " @ {i} @ " " @ {nautapses} @ " " @ {chr 10}
  end // for

  echo totalnautapses {totalnautapses}
  echo nnneurons_with_autapses {nneurons_with_autapses}
  list = (list) @ " " @ {nneurons_with_autapses} @ " " @ {chr 10}
  list = (list) @ " " @ {totalnautapses} @ " " @ {chr 10}
  echo {list}  > {output}