diff --git a/src/new/thermal/examples/FTE_aLIGO_thermal_lenses.m b/src/new/thermal/examples/FTE_aLIGO_thermal_lenses.m
index 22e4a26e83c214bbed80503914c6d4c5610ceda9..cc62bc3f07a8c72e3cc208ea2ed3644438e51316 100644
--- a/src/new/thermal/examples/FTE_aLIGO_thermal_lenses.m
+++ b/src/new/thermal/examples/FTE_aLIGO_thermal_lenses.m
@@ -20,17 +20,9 @@ clear all;
% Choose coating absorption
option = 'coating';
-
-% Filenames for map files
-file1 = 'etm_thermal_lens.txt';
-file2 = 'etm_thermal_res_lens.txt';
-file3 = 'itm_thermal_lens.txt';
-file4 = 'itm_thermal_res_lens.txt';
-filea1 = 'etm_aperture.txt';
-filea2 = 'itm_aperture.txt';
% Number
-num_terms = 100;
+num_terms = 50;
% Beam parameters (aLIGO parameters)
% Two beam parameters at the two cavity mirrors
@@ -55,7 +47,7 @@ Y = 7.3e10;
pois = 0.17;
a = 0.168;
h = 0.2;
-dn_dT = -1.1e-5;
+dn_dT = 1.1e-5;
% For 1 ppm absorption
Babs = 1e-6/(P*h);
@@ -63,7 +55,7 @@ Babs = 1e-6/(P*h);
mirror = FT_init_thermal_mirror(dn_dT,Cabs,Babs,K,E,alpha,Y,pois,a,h);
% Radial coordinate
-r = linspace(-a,a,128);
+r = linspace(-a,a,256);
% Thermal lens: etm
[Z_etm,c_etm,p_etm] = FT_thermal_lens(beam_etm,mirror,r,option,num_terms);
@@ -76,12 +68,6 @@ Zr_etm = Z_etm - Zp_etm;
% Creating mirror maps representing overall and residual lenses: etm
[Zmap_etm,rmap_etm,Rc_etm] = FT_phase_map_from_thermal_lens(beam_etm,mirror,r,option,num_terms);
-
-% Creating aperture map and writing maps to file: etm
-Amap_etm = FT_create_aperture_map(Zmap_etm);
-FT_write_surface_map(file1,Zmap_etm);
-FT_write_surface_map(file2,rmap_etm);
-FT_write_surface_map(filea1,Amap_etm);
% Thermal lens: itm
[Z_itm,c_itm,p_itm] = FT_thermal_lens(beam_itm,mirror,r,option,num_terms);
@@ -95,33 +81,27 @@ Zr_itm = Z_itm - Zp_itm;
% Creating mirror maps representing overall and residual lenses: itm
[Zmap_itm,rmap_itm,Rc_itm] = FT_phase_map_from_thermal_lens(beam_itm,mirror,r,option,num_terms);
-% Creating aperture maps and writing maps to file: itm
-Amap_itm = FT_create_aperture_map(Zmap_itm);
-FT_write_surface_map(file3,Zmap_itm);
-FT_write_surface_map(file4,rmap_itm);
-FT_write_surface_map(filea2,Amap_itm);
-
% Plotting results
figure
% Plotting etm thermal lens and approximate parabaloid
hold on;
- plot(r,Z_etm,'r','LineWidth',3)
- plot(r,Zp_etm,'b--','LineWidth',3)
+ plot(1e2*r,1e9*Z_etm,'r','LineWidth',3)
+ plot(1e2*r,1e9*Zp_etm,'b--','LineWidth',3)
grid on;
- xlabel('r [m] ')
- ylabel('Thermal lens [m] ')
+ xlabel('r [cm] ')
+ ylabel('Thermal lens [nm] ')
title('Thermal distortion: ETM')
legend('Thermal lens','Approximate paraboloid')
- set(gca,'XLim',[-a a])
+ set(gca,'XLim',[-a a]*1e2)
figure
% Plotting etm residual lens
- plot(r,Zr_etm,'r','LineWidth',3)
+ plot(1e2*r,1e9*Zr_etm,'r','LineWidth',3)
grid on;
- xlabel('r [m] ')
- ylabel('Thermal lens [m] ')
+ xlabel('r [cm] ')
+ ylabel('Thermal lens [nm] ')
title('Residual thermal lens: ETM')
- set(gca,'XLim',[-a a])
+ set(gca,'XLim',[-a a]*1e2)
% Plotting thermal distortion as map: etm
FT_plot_mirror_map(Zmap_etm,0)
@@ -132,23 +112,23 @@ FT_plot_mirror_map(rmap_etm,0)
figure
% Plotting itm thermal distortion and approximate parabaloid
hold on;
- plot(r,Z_itm,'r','LineWidth',3)
- plot(r,Zp_itm,'b--','LineWidth',3)
+ plot(1e2*r,1e9*Z_itm,'r','LineWidth',3)
+ plot(1e2*r,1e9*Zp_itm,'b--','LineWidth',3)
grid on;
- xlabel('r [m] ')
- ylabel('Thermal lens [m] ')
+ xlabel('r [cm] ')
+ ylabel('Thermal lens [nm] ')
title('Thermal lens: ITM')
legend('Thermal lens','Approximate paraboloid')
- set(gca,'XLim',[-a a])
+ set(gca,'XLim',[-a a]*1e2)
figure
% Plotting itm residual distortion
- plot(r,Zr_itm,'r','LineWidth',3)
+ plot(1e2*r,1e9*Zr_itm,'r','LineWidth',3)
grid on;
- xlabel('r [m] ')
- ylabel('Thermal lens [m] ')
+ xlabel('r [cm] ')
+ ylabel('Thermal lens [nm] ')
title('Residual thermal lens: ITM')
- set(gca,'XLim',[-a a])
+ set(gca,'XLim',[-a a]*1e2)
% Plotting thermal lens as map: itm
FT_plot_mirror_map(Zmap_itm,0)
diff --git a/src/new/thermal/examples/FTE_aLIGO_thermal_mirrors.m b/src/new/thermal/examples/FTE_aLIGO_thermal_mirrors.m
index ca64a9424756aea413b114d6c21c197e83731404..1dafc4ba3181294b580f89ffab526ca87f7c9537 100644
--- a/src/new/thermal/examples/FTE_aLIGO_thermal_mirrors.m
+++ b/src/new/thermal/examples/FTE_aLIGO_thermal_mirrors.m
@@ -21,14 +21,6 @@ clear all;
% Choose coating absorption
option = 'coating';
-% Filenames for map files
-file1 = 'etm_thermal_dis_map.txt';
-file2 = 'etm_thermal_res_map.txt';
-file3 = 'itm_thermal_dis_map.txt';
-file4 = 'itm_thermal_res_map.txt';
-filea1 = 'etm_aperture.txt';
-filea2 = 'itm_aperture.txt';
-
% Number
num_terms = 100;
@@ -78,12 +70,6 @@ Ur_etm = U_etm - Up_etm;
[Umap_etm,rmap_etm,Rc_etm] = FT_mirror_map_from_thermal_distortion(beam_etm,mirror,r,option,num_terms);
Umap_etm = FT_remove_offset_from_mirror_map(Umap_etm,a);
rmap_etm = FT_remove_offset_from_mirror_map(rmap_etm,a);
-
-% Creating aperture map and writing maps to file: etm
-Amap_etm = FT_create_aperture_map(Umap_etm);
-FT_write_surface_map(file1,Umap_etm);
-FT_write_surface_map(file2,rmap_etm);
-FT_write_surface_map(filea1,Amap_etm);
% Thermal distortion: itm
[U_itm,c_itm,p_itm] = FT_thermal_distortion(beam_itm,mirror,r,option,num_terms);
@@ -98,34 +84,28 @@ Ur_itm = U_itm - Up_itm;
[Umap_itm,rmap_itm,Rc_itm] = FT_mirror_map_from_thermal_distortion(beam_itm,mirror,r,option,num_terms);
Umap_itm = FT_remove_offset_from_mirror_map(Umap_itm,a);
rmap_itm = FT_remove_offset_from_mirror_map(rmap_itm,a);
-
-% Creating aperture maps and writing maps to file: itm
-Amap_itm = FT_create_aperture_map(Umap_itm);
-FT_write_surface_map(file3,Umap_itm);
-FT_write_surface_map(file4,rmap_itm);
-FT_write_surface_map(filea2,Amap_itm);
% Plotting results
figure
% Plotting etm thermal distortion and approximate parabaloid
hold on;
- plot(r,U_etm,'r','LineWidth',3)
- plot(r,Up_etm,'b--','LineWidth',3)
+ plot(1e2*r,1e9*U_etm,'r','LineWidth',3)
+ plot(1e2*r,1e9*Up_etm,'b--','LineWidth',3)
grid on;
- xlabel('r [m] ')
- ylabel('Thermal distortion [m] ')
+ xlabel('r [cm] ')
+ ylabel('Thermal distortion [nm] ')
title('Thermal distortion: ETM')
legend('Thermal distortion','Approximate paraboloid')
- set(gca,'XLim',[-a a])
+ set(gca,'XLim',[-a a]*1e2)
figure
% Plotting etm residual distortion
- plot(r,Ur_etm,'r','LineWidth',3)
+ plot(1e2*r,1e9*Ur_etm,'r','LineWidth',3)
grid on;
- xlabel('r [m] ')
- ylabel('Thermal distortion [m] ')
+ xlabel('r [cm] ')
+ ylabel('Thermal distortion [nm] ')
title('Residual thermal distortion: ETM')
- set(gca,'XLim',[-a a])
+ set(gca,'XLim',[-a a]*1e2)
% Plotting thermal distortion as map: etm
FT_plot_mirror_map(Umap_etm,0)
@@ -136,23 +116,23 @@ FT_plot_mirror_map(rmap_etm,0)
figure
% Plotting itm thermal distortion and approximate parabaloid
hold on;
- plot(r,U_itm,'r','LineWidth',3)
- plot(r,Up_itm,'b--','LineWidth',3)
+ plot(1e2*r,1e9*U_itm,'r','LineWidth',3)
+ plot(1e2*r,1e9*Up_itm,'b--','LineWidth',3)
grid on;
- xlabel('r [m] ')
- ylabel('Thermal distortion [m] ')
+ xlabel('r [cm] ')
+ ylabel('Thermal distortion [nm] ')
title('Thermal distortion: ITM')
legend('Thermal distortion','Approximate paraboloid')
- set(gca,'XLim',[-a a])
+ set(gca,'XLim',[-a a]*1e2)
figure
% Plotting itm residual distortion
- plot(r,Ur_itm,'r','LineWidth',3)
+ plot(1e2*r,1e9*Ur_itm,'r','LineWidth',3)
grid on;
- xlabel('r [m] ')
- ylabel('Thermal distortion [m] ')
+ xlabel('r [cm] ')
+ ylabel('Thermal distortion [nm] ')
title('Residual thermal distortion: ITM')
- set(gca,'XLim',[-a a])
+ set(gca,'XLim',[-a a]*1e2)
% Plotting thermal distortion as map: itm
FT_plot_mirror_map(Umap_itm,0)