Example_Pcirc.m

clear all; close all; clear classes
addpath(genpath('Classes'));

disp('---------------------------------------------------------------------------')
disp('                  OSCAR V3.0                                      ')
disp('  ')


% Define the grid for the simulation: 128 X 128, 30 cm X 30 cm
G1 = Grid(128,0.3);

% Define the incoming beam on the input mirror surface (beam radius 2 cm,
% wavefront curvature 2500 m, propagating toward the waist)
E_input = E_Field(G1,0.02,-2500);

% Define the 2 mirrors, RofC = 2500m, 10 cm in diameter, transmission 2%,
% no loss

IM = Interface(G1,2500,0.10,0.02,0);
EM = Interface(G1,2500,0.10,0.02,0);

% Use the 2 previous Interfaces and the input beam to defing a cavity 1000
% meter long
C1 = Cavity1(IM,EM,1000,E_input);

% Calculate the resonance length
 C1 = Cavity_resonance_phase(C1);

 % Display information about the cavity
Get_info(C1);