################################################################
# MODE_simulation_width.lsf
#
# Run with GPW.lms
# Calculates the effective index and loss of the GPW mode as a
# function of arm width, and plots the data.
#
# Copyright 2016 Lumerical Solutions Inc.
################################################################

# Calculate neff and loss as a function of arm width
npts=11;
antenna_width=linspace(50e-9,150e-9,npts);
neff=matrix(npts);

for (i=1:npts) {
  switchtolayout;
  structures;
  select("GPW");
  set("height",100e-9);
  set("gap",30e-9);
  set("antenna_width",antenna_width(i));

  ?"Calculation " + num2str(i) + " of " + num2str(npts);
  findmodes;

  m="mode1";
  neff(i)=getdata(m,"neff");
}
plot(antenna_width*1e9,real(neff),"Antenna width (nm)","real(neff)");
plot(antenna_width*1e9,imag(neff),"Antenna width (nm)","imag(neff)");