pulse_start = 240 pulse_duration = 160 print("pulse duration: {0} ms".format(pulse_duration * cell.dt)) pulse = np.zeros(n_tsteps) pulse[pulse_start:(pulse_start + pulse_duration)] = 1. # TO DETERMINE OR NOT, maybe just start from zmin = - max cortical thickness cortical_surface_height = 50 clamp = False # CLAMPING if clamp: if cell_id == 1: utils.clamp_ends(cell, 0, pulse_start + pulse_duration, -76.) # Parameters for the external field sigma = 0.3 # source_xs = np.array([-50, -50, -10, -10, 10, 10, 50, 50]) # source_ys = np.array([-50, 50, -10, 10, 10, -10, -50, 50]) # source_xs = np.array([-50, 0, 50, 0, 0]) # source_ys = np.array([0, 50, 0, -50, 0]) # source_geometry = np.array([0, 0, 1, 1, 1, 1, 0, 0]) # stim_amp = 1. # n_stim_amp = -stim_amp / 4 # source_geometry = np.array([0, 0, 0, 0, stim_amp]) # source_geometry = np.array([-stim_amp / 4, -stim_amp / 4, -stim_amp / 4, -stim_amp / 4, stim_amp]) # source_geometry = np.array([stim_amp, stim_amp, stim_amp, stim_amp, -stim_amp])
# cell.set_rotation(y=np.pi/2) # cell.set_rotation(z=np.pi/2) n_tsteps = int(cell.tstop / cell.dt + 1) # print("number of segments: ", cell.totnsegs) t = np.arange(n_tsteps) * cell.dt pulse_start = 240 pulse_duration = 160 # CLAMPING 1/2 if clamp: if cell_id == 0 or cell_id == 1: utils.clamp_ends(cell, pulse_start, pulse_start + pulse_duration) if RANK == 0: print("pulse duration: {0} ms".format(pulse_duration * cell.dt)) pulse = np.zeros(n_tsteps) pulse[pulse_start:(pulse_start + pulse_duration)] = 1. # TO DETERMINE OR NOT, maybe just start from zmin = - max cortical thickness cortical_surface_height = 50 # Parameters for the external field sigma = 0.3