z
Sys.HStrain = [1 1 1]*220; % MHz
Sys.lwEndor = 0.2; % MHz

% EPR experiment settings
Exp.mwFreq = 9.5; % GHz
Exp.Range = [250 400]; % mT
Exp.Harmonic = 0;

[B,spec] = pepper(Sys,Exp);

% ENDOR experiment settings
ExpEndor.mwFreq = Exp.mwFreq;
ExpEndor.ExciteWidth = 100;
ExpEndor.Range = [8 20]; % MHz

Fields = [295 310 323 332 339]; % mT

for iField = 1:numel(Fields)   % loop over all field values
  ExpEndor.Field = Fields(iField);
  [freq,spectra(iField,:)] = salt(Sys,ExpEndor);
end

% Plot EPR spectrum
