// File for initialising the full human network model

load_file("nrngui.hoc")
load_file("Cell_positions.hoc")         // File for generating presynaptic network cell positions
load_file("create_Cells.hoc")           // File for generating objects for presynaptic cells
load_file("create_Connections.hoc")     // File for generating gap junctions and synapses
load_file("save_Matrices.hoc")          // File for saving matrices
load_file("save_responses.hoc")         // File for saving all cell responses
load_file("Parameters.hoc")
load_file("create_light_input.hoc")
load_file("create_ExtStim.hoc")

proc show_time() {
	xpanel("timestep")
	xlabel("Simulation Time")
	xlabel("")
	xvalue("Real Time", "realtime")
	xlabel("")
	xvalue("Current Time Step", "t")
	xpanel(0,0,0)
}

proc init() {
    finitialize()
    fcurrent()
    dt = Time_Step
}

proc start() {

    init()

    // Create the connections between cells
    generate_PR_Gap_Junctions()
    generate_AII_gap_junctions()
    generate_dendritic_fields()
    generate_PR2BC_Synapses()
    generate_BC2RGC_Synapse()
    generate_HZ_Synapses()
    generate_AII_Synapses()

    //Init_Stim()
    //save_matrices()
    create_Input()
    //create_Clamps()

    tstop = Simul_Time
    show_time()

}

start()
sav_matrix()
print_total_cells()
load_file("session1.ses")
record_response()
run()
save_response()