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