Back to index.

func_test.m
clear all
list_constants;
 
%% define frequency points
%N=10; fmin = 1; fmax = 100;
%frequency = linspace(fmin,fmax,N);
 
 
%% calling function
freq = frequency;     % freq = logspace(log10(fmin),log10(fmax),N); would be another interesting option
r_m = MirrorRad;
m = M4;
 
 
%Suspension Thermal Noise
susTN = total_sus_therm(freq);
susTN_strain = susTN / (arm_length/2);
 
susTN_x = x_therm(freq);
susTN_z = z_therm(freq);
susTN_v = v_therm(freq);
 
susTN_x_strain = susTN_x / (arm_length/2);
susTN_z_strain = susTN_z / (arm_length/2);
susTN_v_strain = susTN_v / (arm_length/2);
 
 
%Seismic (longitudinal transmissibility)
long_seis_ap = long_seis_approx(freq);
 
%Seismic (vertical transmissibility)
cross_coupling=1/1000; 
ver_seis_ap = vertical_seis_approx(freq);
 
%Seismic
BSC = BSC_req(freq);
 
long_trans_BSC_approx_strain = BSC .* long_seis_ap ./ (arm_length/2);
ver_trans_BSC_approx_strain = 0*cross_coupling .* BSC .* ver_seis_ap./ (arm_length/2);
total_seis_BSC_approx_strain = sqrt(long_trans_BSC_approx_strain.^2 + ver_trans_BSC_approx_strain.^2);
 
 
%Coating Thermal
%total_coating_thermal = total_coating_therm(freq)./ (arm_length/2);
total_coating_thermal = total_coating_therm_final(freq);%./ (arm_length/2);
 
 
%Clipping Loss
total_clipping_loss = total_clipping(r_m)
 
 
%Total
total_strain=total_seis_BSC_approx_strain + susTN_strain + total_coating_thermal;
 
 
%% Save/name
mass_string=num2str(m);
L4_string=num2str(F_l);
tot_len=num2str(total_length);
 
name_string=[mass_string,'_',L4_string,'_',tot_len,'_770MPa_06m'];
 
 
 
%% (1) aLIGO
% aLIGO.coat_aLIGO=total_coating_thermal;
% aLIGO.seis_aLIGO=total_seis_BSC_approx_strain;
% aLIGO.seis_BSC_aLIGO=BSC;
% aLIGO.sus_aLIGO=susTN_strain;
% aLIGO.sus_x_aLIGO=susTN_x_strain;
% aLIGO.sus_z_aLIGO=susTN_z_strain;
% aLIGO.sus_v_aLIGO=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/test/aLIGO.mat'],'aLIGO');
 
 
%% (1-1') aLIGO (1GPa)
% aLIGO_1GPa.coat_aLIGO_1GPa=total_coating_thermal;
% aLIGO_1GPa.seis_aLIGO_1GPa=total_seis_BSC_approx_strain;
% aLIGO_1GPa.sus_aLIGO_1GPa=susTN_strain;
% aLIGO_1GPa.sus_x_aLIGO_1GPa=susTN_x_strain;
% aLIGO_1GPa.sus_z_aLIGO_1GPa=susTN_z_strain;
% aLIGO_1GPa.sus_v_aLIGO_1GPa=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_1GPa.mat'],'aLIGO_1GPa');
% 
 
%% (1-2') aLIGO (1.5GPa)
% aLIGO_15GPa.coat_aLIGO_15GPa=total_coating_thermal;
% aLIGO_15GPa.seis_aLIGO_15GPa=total_seis_BSC_approx_strain;
% aLIGO_15GPa.sus_aLIGO_15GPa=susTN_strain;
% aLIGO_15GPa.sus_x_aLIGO_15GPa=susTN_x_strain;
% aLIGO_15GPa.sus_z_aLIGO_15GPa=susTN_z_strain;
% aLIGO_15GPa.sus_v_aLIGO_15GPa=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_15GPa.mat'],'aLIGO_15GPa');
 
 
%% (1-3') aLIGO (2GPa)
% aLIGO_2GPa.coat_aLIGO_2GPa=total_coating_thermal;
% aLIGO_2GPa.seis_aLIGO_2GPa=total_seis_BSC_approx_strain;
% aLIGO_2GPa.sus_aLIGO_2GPa=susTN_strain;
% aLIGO_2GPa.sus_x_aLIGO_2GPa=susTN_x_strain;
% aLIGO_2GPa.sus_z_aLIGO_2GPa=susTN_z_strain;
% aLIGO_2GPa.sus_v_aLIGO_2GPa=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_2GPa.mat'],'aLIGO_2GPa');
 
 
%% (1-4') aLIGO (2.5GPa)
% aLIGO_25GPa.coat_aLIGO_25GPa=total_coating_thermal;
% aLIGO_25GPa.seis_aLIGO_25GPa=total_seis_BSC_approx_strain;
% aLIGO_25GPa.sus_aLIGO_25GPa=susTN_strain;
% aLIGO_25GPa.sus_x_aLIGO_25GPa=susTN_x_strain;
% aLIGO_25GPa.sus_z_aLIGO_25GPa=susTN_z_strain;
% aLIGO_25GPa.sus_v_aLIGO_25GPa=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_25GPa.mat'],'aLIGO_25GPa');
 
 
%% (1-5') aLIGO (3GPa)
% aLIGO_3GPa.coat_aLIGO_3GPa=total_coating_thermal;
% aLIGO_3GPa.seis_aLIGO_3GPa=total_seis_BSC_approx_strain;
% aLIGO_3GPa.sus_aLIGO_3GPa=susTN_strain;
% aLIGO_3GPa.sus_x_aLIGO_3GPa=susTN_x_strain;
% aLIGO_3GPa.sus_z_aLIGO_3GPa=susTN_z_strain;
% aLIGO_3GPa.sus_v_aLIGO_3GPa=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_3GPa.mat'],'aLIGO_3GPa');
 
 
%% (1-6') aLIGO (3.5GPa)
% aLIGO_35GPa.coat_aLIGO_35GPa=total_coating_thermal;
% aLIGO_35GPa.seis_aLIGO_35GPa=total_seis_BSC_approx_strain;
% aLIGO_35GPa.sus_aLIGO_35GPa=susTN_strain;
% aLIGO_35GPa.sus_x_aLIGO_35GPa=susTN_x_strain;
% aLIGO_35GPa.sus_z_aLIGO_35GPa=susTN_z_strain;
% aLIGO_35GPa.sus_v_aLIGO_35GPa=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_35GPa.mat'],'aLIGO_35GPa');
 
 
 
%% (1-7') aLIGO (4GPa)
% aLIGO_4GPa.coat_aLIGO_4GPa=total_coating_thermal;
% aLIGO_4GPa.seis_aLIGO_4GPa=total_seis_BSC_approx_strain;
% aLIGO_4GPa.sus_aLIGO_4GPa=susTN_strain;
% aLIGO_4GPa.sus_x_aLIGO_4GPa=susTN_x_strain;
% aLIGO_4GPa.sus_z_aLIGO_4GPa=susTN_z_strain;
% aLIGO_4GPa.sus_v_aLIGO_4GPa=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_4GPa.mat'],'aLIGO_4GPa');
 
 
%% (1-8') aLIGO (4.5GPa)
% aLIGO_45GPa.coat_aLIGO_45GPa=total_coating_thermal;
% aLIGO_45GPa.seis_aLIGO_45GPa=total_seis_BSC_approx_strain;
% aLIGO_45GPa.sus_aLIGO_45GPa=susTN_strain;
% aLIGO_45GPa.sus_x_aLIGO_45GPa=susTN_x_strain;
% aLIGO_45GPa.sus_z_aLIGO_45GPa=susTN_z_strain;
% aLIGO_45GPa.sus_v_aLIGO_45GPa=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_45GPa.mat'],'aLIGO_45GPa');
% 
 
%% (1-9') aLIGO (5GPa)
% aLIGO_5GPa.coat_aLIGO_5GPa=total_coating_thermal;
% aLIGO_5GPa.seis_aLIGO_5GPa=total_seis_BSC_approx_strain;
% aLIGO_5GPa.sus_aLIGO_5GPa=susTN_strain;
% aLIGO_5GPa.sus_x_aLIGO_5GPa=susTN_x_strain;
% aLIGO_5GPa.sus_z_aLIGO_5GPa=susTN_z_strain;
% aLIGO_5GPa.sus_v_aLIGO_5GPa=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_5GPa.mat'],'aLIGO_5GPa');
 
 
%% (1-1) aLIGO + 2.14m
aLIGO_LT214.coat_aLIGO_LT214=total_coating_thermal;
aLIGO_LT214.seis_aLIGO_LT214=total_seis_BSC_approx_strain;
aLIGO_LT214.sus_aLIGO_LT214=susTN_strain;
save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_LT214.mat'],'aLIGO_LT214');
 
 
%% (2) aAH: aLIGO + 1.2GPa
% aLIGO12GPa.coat_aLIGO12GPa=total_coating_thermal;
% aLIGO12GPa.seis_aLIGO12GPa=total_seis_BSC_approx_strain;
% aLIGO12GPa.sus_aLIGO12GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/paper/aLIGO12GPa.mat'],'aLIGO12GPa');
 
 
%% (2-2) aAH: aLIGO + 1.2GPa + 1.1m + 2.14m
% aLIGO_opt_11_214_12GPa.coat_aLIGO_opt_11_214_12GPa=total_coating_thermal;
% aLIGO_opt_11_214_12GPa.seis_aLIGO_opt_11_214_12GPa=total_seis_BSC_approx_strain;
% aLIGO_opt_11_214_12GPa.sus_aLIGO_opt_11_214_12GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/mat/aLIGO_opt_11_214_12GPa.mat'],'aLIGO_opt_11_214_12GPa');
 
 
%% (2-1) aAH: aLIGO + 1.2GPa + 2.14m
% aLIGO_LT214_12GPa.coat_aLIGO_LT214_12GPa=total_coating_thermal;
% aLIGO_LT214_12GPa.seis_aLIGO_LT214_12GPa=total_seis_BSC_approx_strain;
% aLIGO_LT214_12GPa.sus_aLIGO_LT214_12GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_LT214_12GPa.mat'],'aLIGO_LT214_12GPa');
 
 
%% (3) aBL: 40kg + 770MPa + 1.1m
% aLIGO770MPa_11m.coat_aLIGO770MPa_11m=total_coating_thermal;
% aLIGO770MPa_11m.seis_aLIGO770MPa_11m=total_seis_BSC_approx_strain;
% aLIGO770MPa_11m.sus_aLIGO770MPa_11m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO770MPa_11m.mat'],'aLIGO770MPa_11m');
 
 
 
%% (4) aBH: 40kg + 1.2GPa + 1.1m
% aLIGO12GPa_11m.coat_aLIGO12GPa_11m=total_coating_thermal;
% aLIGO12GPa_11m.seis_aLIGO12GPa_11m=total_seis_BSC_approx_strain;
% aLIGO12GPa_11m.sus_aLIGO12GPa_11m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO12GPa_11m.mat'],'aLIGO12GPa_11m');
 
 
 
%% (1) 1AL: 80kg + 770MPa + 0.6m
% m80_770MPa_06m.coat_m80_770MPa_06m=total_coating_thermal;
% m80_770MPa_06m.seis_m80_770MPa_06m=total_seis_BSC_approx_strain;
% m80_770MPa_06m.sus_m80_770MPa_06m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/m80_770MPa_06m.mat'],'m80_770MPa_06m');
 
 
% m80_770MPa_0582m.coat_m80_770MPa_0582m=total_coating_thermal;
% m80_770MPa_0582m.seis_m80_770MPa_0582m=total_seis_BSC_approx_strain;
% m80_770MPa_0582m.sus_m80_770MPa_0582m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/m80_770MPa_0582m.mat'],'m80_770MPa_0582m');
 
 
%% (2) 1AH: 80kg + 1.2GPa + 0.6m
% m80_12GPa_06m.coat_m80_12GPa_06m=total_coating_thermal;
% m80_12GPa_06m.seis_m80_12GPa_06m=total_seis_BSC_approx_strain;
% m80_12GPa_06m.sus_m80_12GPa_06m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/m80_12GPa_06m.mat'],'m80_12GPa_06m');
 
 
 
 
%% (3) 1BL: 80kg + 770MPa + 1.1m
% m80_770MPa_11m.coat_m80_770MPa_11m=total_coating_thermal;
% m80_770MPa_11m.seis_m80_770MPa_11m=total_seis_BSC_approx_strain;
% m80_770MPa_11m.sus_m80_770MPa_11m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/m80_770MPa_11m.mat'],'m80_770MPa_11m');
 
 
%% (4) 1BH: 80kg + 1.2GPa + 1.1m
% m80_12GPa_11m.coat_m80_12GPa_11m=total_coating_thermal;
% m80_12GPa_11m.seis_m80_12GPa_11m=total_seis_BSC_approx_strain;
% m80_12GPa_11m.sus_m80_12GPa_11m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/paper/m80_12GPa_11m.mat'],'m80_12GPa_11m');
 
 
 
%% (1) 2AL: 160kg + 770MPa + 0.6m
% m160_770MPa_06m.coat_m160_770MPa_06m=total_coating_thermal;
% m160_770MPa_06m.seis_m160_770MPa_06m=total_seis_BSC_approx_strain;
% m160_770MPa_06m.sus_m160_770MPa_06m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/m160_770MPa_06m.mat'],'m160_770MPa_06m');
 
% m160_770MPa_0582m.coat_m160_770MPa_0582m=total_coating_thermal;
% m160_770MPa_0582m.seis_m160_770MPa_0582m=total_seis_BSC_approx_strain;
% m160_770MPa_0582m.sus_m160_770MPa_0582m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/m160_770MPa_0582m.mat'],'m160_770MPa_0582m');
 
 
%% (2) 2AH: 160kg + 1.2GPa + 0.6m
% m160_12GPa_06m.coat_m160_12GPa_06m=total_coating_thermal;
% m160_12GPa_06m.seis_m160_12GPa_06m=total_seis_BSC_approx_strain;
% m160_12GPa_06m.sus_m160_12GPa_06m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/m160_12GPa_06m.mat'],'m160_12GPa_06m');
 
 
 
 
%% (3) 2BL: 160kg + 770MPa + 1.1m
% m160_770MPa_11m.coat_m160_770MPa_11m=total_coating_thermal;
% m160_770MPa_11m.seis_m160_770MPa_11m=total_seis_BSC_approx_strain;
% m160_770MPa_11m.sus_m160_770MPa_11m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/m160_770MPa_11m.mat'],'m160_770MPa_11m');
 
 
%% (4) 2BH: 160kg + 1.2GPa + 1.1m
% m160_12GPa_11m.coat_m160_12GPa_11m=total_coating_thermal;
% m160_12GPa_11m.seis_m160_12GPa_11m=total_seis_BSC_approx_strain;
% m160_12GPa_11m.sus_m160_12GPa_11m=susTN_strain;
% save(['/Users/klee/Documents/Matlab/1718_Thesis/paper/m160_12GPa_11m.mat'],'m160_12GPa_11m');
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
%% aLIGO upgrade #1: 40/40/22/22, 0.6m, 1.6m, 1.54GPa, 5mm
% aLIGO_154GPa.coat_aLIGO_154GPa=total_coating_thermal;
% aLIGO_154GPa.seis_aLIGO_154GPa=total_seis_BSC_approx_strain;
% aLIGO_154GPa.sus_aLIGO_154GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/aLIGO_154GPa.mat'],'aLIGO_154GPa');
 
 
 
%% aLIGO upgrade #2: 40/opt, 0.6m, 2.14m, 1.54GPa, 5mm
% aLIGO_opt_06_214_154GPa.coat_aLIGO_opt_06_214_154GPa=total_coating_thermal;
% aLIGO_opt_06_214_154GPa.seis_aLIGO_opt_06_214_154GPa=total_seis_BSC_approx_strain;
% aLIGO_opt_06_214_154GPa.sus_aLIGO_opt_06_214_154GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/aLIGO_opt_06_214_154GPa.mat'],'aLIGO_opt_06_214_154GPa');
 
%% aLIGO upgrade #3: 40/opt, 1.1m, 2.14m, 1.54GPa, 5mm
% aLIGO_opt_11_214_154GPa.coat_aLIGO_opt_11_214_154GPa=total_coating_thermal;
% aLIGO_opt_11_214_154GPa.seis_aLIGO_opt_11_214_154GPa=total_seis_BSC_approx_strain;
% aLIGO_opt_11_214_154GPa.sus_aLIGO_opt_11_214_154GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/aLIGO_opt_11_214_154GPa.mat'],'aLIGO_opt_11_214_154GPa');
 
 
 
%% aLIGO upgrade #4: 80/opt, 1.1m, 2.14m, 1.54GPa, 5mm
% m80_11_214_154GPa.coat_80_11_214_154GPa=total_coating_thermal;
% m80_11_214_154GPa.seis_80_11_214_154GPa=total_seis_BSC_approx_strain;
% m80_11_214_154GPa.sus_80_11_214_154GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_11_214_154GPa.mat'],'m80_11_214_154GPa');
%save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_06_214_770MPa');
 
%% aLIGO upgrade #5: 160/opt, 1.1m, 2.14m, 1.54GPa, 5mm
% m160_11_214_154GPa.coat_160_11_214_154GPa=total_coating_thermal;
% m160_11_214_154GPa.seis_160_11_214_154GPa=total_seis_BSC_approx_strain;
% m160_11_214_154GPa.sus_160_11_214_154GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m160_11_214_154GPa.mat'],'m160_11_214_154GPa');
% %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_11_214_770MPa');
 
 
 
%% aLIGO upgrade #4-1: 80/opt, 0.6m, 1.6m, 1.54GPa, 5mm
% m80_06_16_154GPa.coat_80_06_16_154GPa=total_coating_thermal;
% m80_06_16_154GPa.seis_80_06_16_154GPa=total_seis_BSC_approx_strain;
% m80_06_16_154GPa.sus_80_06_16_154GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_06_16_154GPa.mat'],'m80_06_16_154GPa');
% %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_06_16_154GPa');
 
%% aLIGO upgrade #4-2: 80/opt, 1.1m, 1.6m, 1.54GPa, 5mm
% m80_11_16_154GPa.coat_80_11_16_154GPa=total_coating_thermal;
% m80_11_16_154GPa.seis_80_11_16_154GPa=total_seis_BSC_approx_strain;
% m80_11_16_154GPa.sus_80_11_16_154GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_11_16_154GPa.mat'],'m80_11_16_154GPa');
% %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_11_16_154GPa');
 
 
 
%% aLIGO upgrade #5-1: 80/opt, 0.6m, 2.14m, 1.54GPa, 5mm
% m80_06_214_154GPa.coat_80_06_214_154GPa=total_coating_thermal;
% m80_06_214_154GPa.seis_80_06_214_154GPa=total_seis_BSC_approx_strain;
% m80_06_214_154GPa.sus_80_06_214_154GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_06_214_154GPa.mat'],'m80_06_214_154GPa');
% %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_06_214_154GPa');
 
%% aLIGO upgrade #5-2: 80/opt, 1.1m, 2.14m, 1.54GPa, 5mm
% m80_11_214_154GPa.coat_80_11_214_154GPa=total_coating_thermal;
% m80_11_214_154GPa.seis_80_11_214_154GPa=total_seis_BSC_approx_strain;
% m80_11_214_154GPa.sus_80_11_214_154GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_11_214_154GPa.mat'],'m80_11_214_154GPa');
% %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_11_214_154GPa');
 
 
%% aLIGO upgrade 1.2GPa: 40/40/22/22, 0.6m, 1.6m, 1.2GPa, 5mm
% aLIGO_12GPa.coat_aLIGO_12GPa=total_coating_thermal;
% aLIGO_12GPa.seis_aLIGO_12GPa=total_seis_BSC_approx_strain;
% aLIGO_12GPa.sus_aLIGO_12GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/aLIGO_12GPa.mat'],'aLIGO_12GPa');
 
%% aLIGO upgrade 1.2GPa: 160/opt, 1.1m, 2.14m, 1.2GPa, 5mm
% m160_11_214_154GPa.coat_160_11_214_154GPa=total_coating_thermal;
% m160_11_214_154GPa.seis_160_11_214_154GPa=total_seis_BSC_approx_strain;
% m160_11_214_12GPa.sus_160_11_214_12GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m160_11_214_12GPa.mat'],'m160_11_214_12GPa');
% %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_11_214_770MPa');
 
%% aLIGO upgrade 1.2GPa: 80/opt, 1.1m, 2.14m, 1.2GPa, 5mm
% m80_11_214_154GPa.coat_80_11_214_154GPa=total_coating_thermal;
% m80_11_214_154GPa.seis_80_11_214_154GPa=total_seis_BSC_approx_strain;
% m80_11_214_12GPa.sus_80_11_214_12GPa=susTN_strain;
% save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/m80_11_214_12GPa.mat'],'m80_11_214_12GPa');
% %save(['/Users/klee/Documents/Matlab/150504_standard_figs/',name_string,'.mat'],'m80_11_214_154GPa');
 
 
%% aLIGO upgrade 1.2GPa: 40/opt, 1.1m, 2.14m, 1.2GPa, 5mm
% aLIGO_opt_11_214_154GPa.coat_aLIGO_opt_11_214_154GPa=total_coating_thermal;
% aLIGO_opt_11_214_154GPa.seis_aLIGO_opt_11_214_154GPa=total_seis_BSC_approx_strain;
%  aLIGO_opt_11_214_12GPa.sus_aLIGO_opt_11_214_12GPa=susTN_strain;
%  save(['/Users/klee/Documents/Matlab/2016LVC_Glasgow/aLIGO_opt_11_214_12GPa.mat'],'aLIGO_opt_11_214_12GPa');
 
 
 
%% (L4_test1) aLIGO 0.6, 2.14
% aLIGO_06_214.coat_aLIGO_06_214=total_coating_thermal;
% aLIGO_06_214.seis_aLIGO_06_214=total_seis_BSC_approx_strain;
% aLIGO_06_214.sus_aLIGO_06_214=susTN_strain;
% aLIGO_06_214.sus_x_aLIGO_06_214=susTN_x_strain;
% aLIGO_06_214.sus_z_aLIGO_06_214=susTN_z_strain;
% aLIGO_06_214.sus_v_aLIGO_06_214=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_06_214.mat'],'aLIGO_06_214');
 
 
%% (L4_test2) aLIGO 1.1, 2.14
% aLIGO_11_214.coat_aLIGO_11_214=total_coating_thermal;
% aLIGO_11_214.seis_aLIGO_11_214=total_seis_BSC_approx_strain;
% aLIGO_11_214.sus_aLIGO_11_214=susTN_strain;
% aLIGO_11_214.sus_x_aLIGO_11_214=susTN_x_strain;
% aLIGO_11_214.sus_z_aLIGO_11_214=susTN_z_strain;
% aLIGO_11_214.sus_v_aLIGO_11_214=susTN_v_strain;
% 
% save(['/Users/klee/Documents/Matlab/1718_Thesis/aLIGO_11_214.mat'],'aLIGO_11_214');
 
 
 
%% =========figure========
%{
%tic
figure(1)
title(['Frequency vs Strain for M4 = ',num2str(m),' kg'])
loglog(freq,total_seis_BSC_approx_strain,'c'); hold on
loglog(freq,susTN_strain,'k'); hold on
loglog(freq,total_coating_thermal,'r'); hold on
grid on
xlabel('Frequency (Hz)');
ylabel('Strain (1/rtHz)');
legend('seis-total','sus ther','coating')
axis([0 100 1e-30 1e-12])
 
figure(2)
title(['Frequency vs Strain for M4 = ',num2str(m),' kg'])
loglog(freq,total_seis_BSC_approx_strain,'b'); hold on
loglog(freq,susTN_strain,'g'); hold on
loglog(freq,total_coating_thermal,'r'); hold on
loglog(freq,total_strain,'k'); hold on
grid on
xlabel('Frequency (Hz)');
ylabel('Strain (1/rtHz)');
legend('seis-total','sus ther','coating','total')
axis([0 100 1e-30 1e-12])
 
%toc
%}