current_dir = '.' 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')