N_TEMPERATURES = 10
o = OmR()
o.set_params(m=6.75e-12,
f_mech=1e6,
q=1000.,
lambda_nm = 1064.,
losses=1.e-6,
transmission_input=300e-6,
length=0.00764,
i_incident=1e-3,
delta_hz=1.)
tau = 1./(2*pi*1.7e6)
o.tau_in = tau/0.4
o.tau_ex = tau/0.6
o.n_th = 0
o.freq = linspace(0.99e6, 1.01e6, 400)
o.delta_hz = 0
figure('SB asymmetry on resonance', figsize=(10,13))
subplot(211)
o.n_th = 0
o.cooperativity = 0.1
title('Cooperativity = 0.1, n_th = 0')
plot(o.freq, o.right.output.spectrum_ss('upper'), '-r', label='Upper sideband')
plot(o.freq, o.right.output.spectrum_ss('lower'), '-b', label='Lower sideband')
