################################################################
# LCD_phased_array.lsf
#
# Used with LCD_phased_array.fsp
# This script sets up the spatially-varying LC orientation along
# the y-direction in order to steer a beam.
#
# Copyright 2016 Lumerical Solutions Inc.
################################################################

# delete existing attribute
select('LC attribute');
delete;

# define x/y/z
x = 0;
y = linspace(-80e-6,80e-6,100);
z = 0;

X = meshgrid3dx(x,y,z);
Y = meshgrid3dy(x,y,z);
Z = meshgrid3dz(x,y,z);
n = matrix(length(x),length(y),length(z),3);

# define the orientation function
n(1:length(x),1:length(y),1:length(z),1) = cos((Y+80e-6)/160e-6*pi/2);
n(1:length(x),1:length(y),1:length(z),2) = sin((Y+80e-6)/160e-6*pi/2);
n(1:length(x),1:length(y),1:length(z),3) = Z;

# add LC import grid attribute
LC=rectilineardataset("LC",x,y,z);
LC.addattribute("u",n);
addgridattribute("lc orientation",LC);
setnamed("LC attribute","ny",100); # set resolution