% set operating point to RSE
disp('Option 1 - low power - define interferometer');
opt1 = adVirgo1_LP;
pos = zeros(opt1.Ndrive, 1);
nSRM = getDriveIndex(opt1, 'SRM');
pos(nSRM) =  1.064e-6 / 4;

% frequency bins
f = logspace(-1, 3.7, 1000);

disp('Option 1 - low power - run simulation');
[fDC1lp, sigDC1lp, sigAC1lp, mMech1lp, noiseAC1lp] = tickle(opt1, pos, f);
disp('Option 1 - low power - saving');
save advirgo1_LP.mat f fDC1lp sigDC1lp sigAC1lp mMech1lp noiseAC1lp
