// 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 = 326;   i > 325; i = i - 1)
	   for (i = 450;   i > 0; i = i - 320)
	        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}

*/

//                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE itotal[8]
//                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE itotal[9]
	   end
//           addmsg /PID1 ^ SAVE  output
//           addmsg /PID2 ^ SAVE  output

/*
           for (i = 234;   i > 233; i = i - 1)
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve soma Vm}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
           end

// to this 238 is connected but 233 not, hence second-order connection
           for (i = 259;   i > 258; i = i - 1)
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve soma Vm}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
           end
*/

/*
           for (i = 226;   i > 225; i = i - 1)
//           for (i = 3;   i > 2; i = i - 1)
                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 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}
//                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE itotal[8]
//                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE itotal[9]
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve soma/Ca_pool Ca}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend3[3]/Ca_pool Ca}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
           end
           for (i = 126;   i > 125; i = i - 1)
//           for (i = 2;   i > 1; i = i - 1)
                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 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}
//                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE itotal[8]
//                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE itotal[9]
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve soma/Ca_pool Ca}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend3[3]/Ca_pool Ca}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
           end
           for (i = 26;   i > 25; i = i - 1)
//             for (i = 1;   i > 0; i = i - 1)
              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 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}
//                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE itotal[8]
//                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE itotal[9]
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve soma/Ca_pool Ca}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
                elem = {findsolvefield /molecular_layer/Stellate[{i-1}]/solve dend3[3]/Ca_pool Ca}
                addmsg /molecular_layer/Stellate[{i-1}]/solve ^ SAVE {elem}
           end
*/

   end

