// genesis



str elem
int i, j

str ascii_label = {label}

str Purkinje_cells_ascii_filename    =  {filename} @ "PCs_"    @ {ascii_label} @ ".ascii"
str Stellate_cells_ascii_filename =  {filename} @ "stcs_" @ {ascii_label} @ ".ascii" 


/*
storing leak in itotal[0]
storing InNa in itotal[1]
storing KDr in itotal[2]
storing KA in itotal[3]
storing CaHVA in itotal[4]
storing H in itotal[5]
storing Moczyd_KC in itotal[6]
storing mf_AMPA in itotal[7]
storing pf_AMPA in itotal[8]
storing GABAA in itotal[9]
storing Ca_pool in itotal[10]
storing Ca_nernst in itotal[11]
storing spike in itotal[12]
*/


// stellate cells

   if ({{number_stellate_cells} > 0})
	   echo   creating ascii-file for stellate cells

           create asc_file /output/plot_out_stcs
           setclock 8 4e-5
           useclock ^ 8
           enable /output
           enable ^

	   setfield ^ filename {Stellate_cells_ascii_filename} \
        	    initialize 1 leave_open 1 flush 1
	   echo Output to {Stellate_cells_ascii_filename}

//           for (i = {number_stellate_cells};   i > 0 ; i = i - 10)
//        for (i = 447;   i > 446; i = i - 1)
	   for (i = 447;   i > 0; i = i - 225)
	        elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve soma Vm}
	        addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}


/*
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend1  Vm}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend2 Vm}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend3 Vm}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}

                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve soma/GABAA Gk}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve soma/GABAA Ik}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}

                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend1/GABAA Gk}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend1/GABAA Ik}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}

                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend1/pf_AMPA Gk}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend1/pf_AMPA Ik}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend2/pf_AMPA Gk}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend2/pf_AMPA Ik}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend3/pf_AMPA Gk}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend3/pf_AMPA Ik}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}

*/

	   end
           addmsg /PID1 ^ SAVE  output
           addmsg /PID2 ^ SAVE  output

   end