################################################################################
# graphene_solition_mode_calculation.lsf
#
# This script file calculates modes of a slab wavegide which embeds 
# a graphene sheet.
# After finding the modes, a TE mode profile is output to 
# a file named mode_profie.txt.
# 
# Copyright 2013 Lumerical Solutions Inc
################################################################################

findmodes;

#mode1 is a TE mode
Ex=pinch(getdata("MODE::data::mode1","Ex"));
Ey=pinch(getdata("MODE::data::mode1","Ey"));
Ez=pinch(getdata("MODE::data::mode1","Ez"));

y=pinch(getdata("MODE::data::mode1","y"));
f=pinch(getdata("MODE::data::mode1","f"));

plot(y,real(Ex),"y","Re(Ex)");

if (fileexists("mode_profile.txt")){
rm("mode_profile.txt");
}

for(i=1:length(y)){
  write("mode_profile.txt",num2str(y(i))+", "+num2str(real(Ex(i))));
}