clearvars; close all;
addpath(genpath('Classes'));

G1 = Grid(128,0.075);

%nominal waist = 1 cm with R=-400 - the beam is going toward the waist 
% E_input = E_Field(G1,'w',0.01025,'R',-400);

%10% mismatching waist
%E_input = E_Field(G1,'w',0.0068,'R',-400);

%4% mismatching waist
E_input = E_Field(G1,'w',0.008,'R',-400);

%9% mismatching waist
%E_input = E_Field(G1,'w',0.013,'R',-400);

%ROC of the IM and EM mirrors = 400 m 
IM = Interface(G1,'RoC',400,'CA',0.10,'T',0.0014);
EM = Interface(G1,'RoC',400,'CA',0.10,'T',1e-5);

% Misalignment of the end mirror
EM = Add_tilt(EM,5E-6,'y');

% 300 m cavity with mirrors IM and EM 
C1 = Cavity1(IM,EM,300,E_input);
C1.Laser_start_on_input = false ;
C1 = Cavity_resonance_phase(C1);
C1 = Calculate_fields(C1); 
Display_results(C1);

% remove the following command to have a colored colormap
colormap gray;