close all
clear all
clc
%%
freq = logspace(-2,2,100);

susTypes = {'BHSS';'OPOS';'OFIS';'HAUX';'HTTS';'HATS';'HATS'};
isoStage = {'M1';  'M1';  'M1';  'M1';  'M1';  'M1';  'M3'};
osemType = { 'A';   'A';   'A';   'A';   'B';   'B';  'A'};
drvrType = { '-v1 (100 Ohm)';'-v2 (1.2 kOhm)';'-v2 (1.2 kOhm)';'-v2 (1.2 kOhm)';'-v2 (1.2 kOhm)';'-v1 (100 Ohm)';'-v2 (1.2 kOhm)'};
magSize  = {'2.0x6.0'; '2.0x6.0'; '3.0x6.0'; '1.9x3.2'; '2x3'; '4.0x4.0';'2.0x0.5'};

%%
for iSusStage = 1:length(susTypes)
    disp([susTypes{iSusStage} isoStage{iSusStage}])
    [~,~,osemModel,coilDriverModel,~] = make_OSEM_filter_model(freq,susTypes{iSusStage},isoStage{iSusStage},'lp');
    output{1,1} = 'susType';
    output{1,2} = 'isoStage';
    output{1,3} = 'Osem Type';
    output{1,4} = 'Magnet Size';
    output{1,5} = 'Force Co [N/A]';
    output{1,6} = 'Driver Type';
    output{1,7} = 'Driver TC [mA/V]';
    output{1,8} = 'Total Strength [N/V]';
    output{iSusStage+1,1} = susTypes{iSusStage};
    output{iSusStage+1,2} = isoStage{iSusStage};
    output{iSusStage+1,3} = osemType{iSusStage};
    output{iSusStage+1,4} = magSize{iSusStage};
    output{iSusStage+1,5} = osemModel.NperA;
    output{iSusStage+1,6} = drvrType{iSusStage};
    output{iSusStage+1,7} = 1e3*coilDriverModel.dcTransConductance;
    output{iSusStage+1,8} = coilDriverModel.dcTransConductance * osemModel.NperA;
end

%%
output