Testing the c302 model based on hub and spoke structure as described in https://github.com/openworm/OpenWorm/issues/212

hub_spoke_testing.py

import os

if __name__ == '__main__':

    hub_cell = "RMGR"

    spoke_cells =["URXR","ASKR","ASJR","RMHR","AWBR","IL2R","ASHR","ADLR"]


# stimulating RMGR (hub cell) when paired with all spoke cells
    for cell in spoke_cells:
        cells_to_plot="["+hub_cell+","+cell+"]"
        cells_to_stimulate="["+hub_cell+"]"
        #print cells_to_plot
        #print cells_to_stimulate
        bashCommand = 'python c302.py '+cell+"h"+' parameters_A -cells '+cells_to_plot+' -cellstostimulate '+cells_to_stimulate+' -duration 500'
        os.system(bashCommand)

# stimulating each spoke cell when paired with RMGR
    for cell in spoke_cells:
        cells_to_plot="["+hub_cell+","+cell+"]"
        cells_to_stimulate="["+cell+"]"
        #print cells_to_plot
        #print cells_to_stimulate
        bashCommand = 'python c302.py '+cell+"s"+' parameters_A -cells '+cells_to_plot+' -cellstostimulate '+cells_to_stimulate+' -duration 500'
        os.system(bashCommand)