for ii, condition in enumerate(conditions):
photoreceptor = FlyFactory.DrosophilaR16(shift="none")
DepolarisePhotoreceptor.WithLight(photoreceptor, V=V)
if condition == 'PIP2':
ShiftConductances.WithLight(photoreceptor,
change_LIC_to_keep_depolarisation=
change_LIC_to_keep_depolarisation)
elif condition == 'Serotonin':
ShiftConductances.WithSerotonin(photoreceptor,
change_LIC_to_keep_depolarisation=
change_LIC_to_keep_depolarisation)
elif condition == 'Shab-50':
Experiment.modify_conductance(photoreceptor,
"Shab",
.5,
change_LIC_to_keep_depolarisation=
change_LIC_to_keep_depolarisation)
elif condition == 'Shab+50':
Experiment.modify_conductance(photoreceptor,
"Shab",
1.5,
change_LIC_to_keep_depolarisation=
change_LIC_to_keep_depolarisation)
else:
print("Warning, option not recognised")
Vr_new_r[ii, i] = photoreceptor.body.V_m
Z_cut = photoreceptor.body.impedance(
f_from_medium
for i,V in enumerate(Vr_continuous):
DepolarisePhotoreceptor.WithLight(HH,V)
GBWP_continuous_[i] = GBWP(HH.body.impedance, f_min = f_medium)
Experiment.freeze_conductances(HH)
GBWP_RC_continuous_[i] = GBWP(HH.body.impedance, f_min = f_medium)
Experiment.unfreeze_conductances(HH)
for ii in range(3):
HH_shifted = copy.deepcopy(HH)
if ii==0:
ShiftConductances.WithLight(HH_shifted, change_LIC_to_keep_depolarisation = change_LIC_to_keep_depolarisation)
elif ii==2:
ShiftConductances.WithSerotonin(HH_shifted, change_LIC_to_keep_depolarisation = change_LIC_to_keep_depolarisation)
elif ii==1:
Experiment.modify_conductance(HH_shifted, "Shab", .5, change_LIC_to_keep_depolarisation = change_LIC_to_keep_depolarisation)
else:
raise Exception
GBWP_shift_continuous_[ii,i] = GBWP(HH_shifted.body.impedance, f_min = f_medium)
Vr_shift_continuous[ii,i] = HH_shifted.body.V_m
for ii in range(3):
ax_GBWP[ii].plot(Vr_continuous,GBWP_continuous_/1e3,'k',zorder=0,alpha=0.2)
ax_GBWP[ii].plot(Vr_shift_continuous[ii],GBWP_shift_continuous_[ii,:]/1e3,'k',zorder=0)
ax_GBWP[ii].plot(Vr_continuous,GBWP_RC_continuous_/1e3,'k--',zorder=0,alpha=0.2)
GBWP_ = np.zeros_like(Vr)
GBWP_RC_ = np.zeros_like(Vr)
GBWP_shift_ = np.zeros((3,len(Vr)))
