{"id":326278,"date":"2023-12-05T00:51:29","date_gmt":"2023-12-05T00:51:29","guid":{"rendered":"\/forum\/forums\/topic\/troubleshooting-a-couple-of-errors-in-lumerical-scripting\/"},"modified":"2023-12-05T00:51:29","modified_gmt":"2023-12-05T00:51:29","slug":"troubleshooting-a-couple-of-errors-in-lumerical-scripting","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/troubleshooting-a-couple-of-errors-in-lumerical-scripting\/","title":{"rendered":"Troubleshooting a couple of errors in Lumerical scripting"},"content":{"rendered":"

I am attempting to recreate the Lumerical RCWA portion of the small metalens tutorial<\/a> using only Lumerical scripting rather than starting from the pre-built .fsp example file. When scripting Lumerical, I am having trouble with the command:
addrcwa;

Yet, this command appears in the online knowledge base here<\/a>. 
If I choose to use the command “addfdtd;” instead of “addrcwa;” the line executes, but I get an error further down the script when using the command:
addsweepresult(“thickness_sweep_script”, result_1);
<\/p>\n

This seems strange to me because the code worked last week without errors when using the “addfdtd;” command.
A couple of notes:
1) I initially built this script with FDTD in mind but decided to try switching to RCWA since that is what appears in the online tutorial linked above. I realize that I may need to change some parts of my code accordingly. However, “addrcwa;” fails even when run as a lone command in the command window.
2) The screenshots above show the problematic lines being run alone in the command window. The errors are identical when running the full code. Only the line number in the error message changes<\/p>\n

I have included my code below. (I am still trying to work out how best to add sweeps and collect the results to get the phase response info, but that’s another topic.) Thanks in advance for any advice!<\/p>\n

SCRIPT: <\/p>\n

deleteall;<\/div>\n
if(exist(“cylinder_height_sweep”)) {deletesweep(“cylinder_height_sweep”);}<\/div>\n
if(exist(“cylinder_radius_sweep”)) {deletesweep(“cylinder_radius_sweep”);}<\/div>\n
### VARIABLES ###<\/div>\n
# Constants and Units<\/div>\n
nm = 1e-9;<\/div>\n
um = 1e-6;<\/div>\n
# Nanopillar<\/div>\n
height_np_min = 600*nm;<\/div>\n
height_np_max = 1.2*um;<\/div>\n
radius_np_min = 75*nm;<\/div>\n
radius_np_max = 125*nm;<\/div>\n
material_np = “Si (Silicon) – Palik”;<\/div>\n
# Substrate<\/div>\n
period = 400*nm;<\/div>\n
thick_sub = 100*nm;<\/div>\n
material_sub = “SiO2 (Glass) – Palik”;<\/div>\n
# FDTD<\/div>\n
lda_min = 400*nm;<\/div>\n
lda_max = 800*nm;<\/div>\n
airgap = lda_max\/2;<\/div>\n
#height_cell = height_np_max+thick_sub+2*airgap;<\/div>\n
zmin_cell = -1*(thick_sub + airgap);<\/div>\n
zmax_cell = height_np_max + airgap;<\/div>\n
 <\/div>\n
# Cylindrical Nanopillar<\/div>\n
addcircle;<\/div>\n
set(“x”,0);<\/div>\n
set(“y”,0);<\/div>\n
set(“radius”,radius_np_min);<\/div>\n
set(“z min”,0);<\/div>\n
set(“z max”,height_np_min);<\/div>\n
set(“material”,material_np);<\/div>\n
 <\/div>\n
# Substrate<\/div>\n
addrect;<\/div>\n
set(“x”,0);<\/div>\n
set(“x span”,period);<\/div>\n
set(“y”,0);<\/div>\n
set(“y span”,period);<\/div>\n
set(“z min”,-1*thick_sub);<\/div>\n
set(“z max”,0);<\/div>\n
set(“material”,material_sub);<\/div>\n
 <\/div>\n
# FDTD Region<\/div>\n
#addrcwa;<\/div>\n
addfdtd;<\/div>\n
set(“dimension”,2);  #  1 = 2D, 2 = 3D<\/div>\n
set(“x”,0);<\/div>\n
set(“x span”,period);<\/div>\n
set(“y”,0);<\/div>\n
set(“y span”,period);<\/div>\n
set(“z min”,zmin_cell);<\/div>\n
set(“z max”,zmax_cell);<\/div>\n
set(“x min bc”,”Periodic”);<\/div>\n
set(“y min bc”,”Periodic”);<\/div>\n
set(“pml profile”,3); # 1. standard, 2. stabilized<\/div>\n
                     # 3. steep angle, 4. custom<\/div>\n
set(“allow symmetry on all boundaries”,true);<\/div>\n
set(“set simulation bandwidth”,true);<\/div>\n
set(“simulation wavelength min”,lda_min);<\/div>\n
set(“simulation wavelength max”,lda_max);<\/div>\n
set(“always use complex fields”,true);<\/div>\n
 <\/div>\n
### Mesh Override Region<\/div>\n
addmesh;<\/div>\n
set(“dx”,10*nm);<\/div>\n
set(“dy”,10*nm);<\/div>\n
set(“dz”,10*nm);<\/div>\n
set(“based on a structure”,true);<\/div>\n
set(“structure”,”circle”);<\/div>\n
 <\/div>\n
### Source<\/div>\n
addplane;<\/div>\n
set(“injection axis”,”z”);<\/div>\n
set(“x”,0);<\/div>\n
set(“x span”,period);<\/div>\n
set(“y”,0);<\/div>\n
set(“y span”,period);<\/div>\n
set(“z”,zmin_cell);<\/div>\n
# If Frequency Domain<\/div>\n
set(“wavelength start”,lda_min);<\/div>\n
set(“wavelength stop”,lda_max);<\/div>\n
# Else if Time Domain<\/div>\n
set(“set time domain”,true);<\/div>\n
 <\/div>\n
### Monitor<\/div>\n
addpower;<\/div>\n
set(“name”,”Transmission”);<\/div>\n
set(“monitor type”,”2D Z-normal”);<\/div>\n
set(“x”,0);<\/div>\n
set(“x span”,period);<\/div>\n
set(“y”,0);<\/div>\n
set(“y span”,period);<\/div>\n
set(“z”,zmax_cell-airgap\/2);<\/div>\n
 <\/div>\n
### Parameter Sweep<\/div>\n
# create parameter sweep “height”<\/div>\n
addsweep;<\/div>\n
setsweep(“sweep”, “name”, “cylinder_height_sweep”);<\/div>\n
setsweep(“cylinder_height_sweep”, “type”, “Ranges”);<\/div>\n
setsweep(“cylinder_height_sweep”, “number of points”, 10);<\/div>\n
# define the parameter height<\/div>\n
para_h = struct;<\/div>\n
para_h.Name = “height”;<\/div>\n
para_h.Parameter = “::model::height”;<\/div>\n
para_h.Type = “Length”;<\/div>\n
para_h.Start = height_np_min;<\/div>\n
para_h.Stop = height_np_max;<\/div>\n
para_h.Units = “microns”;<\/div>\n
addsweepparameter(“cylinder_height_sweep”, para_h);<\/div>\n
# define results<\/div>\n
result_1 = struct;<\/div>\n
result_1.Name = “S”;<\/div>\n
result_1.Result = “S”;<\/div>\n
result_2 = struct;<\/div>\n
result_2.Name = “T”;<\/div>\n
result_2.Result = “T”;<\/div>\n
# add the results R & T to the sweep<\/div>\n
addsweepresult(“thickness_sweep_script”, result_1);<\/div>\n
addsweepresult(“thickness_sweep_script”, result_2);<\/div>\n
 <\/div>\n
# create parameter sweep “radius”<\/div>\n
addsweep;<\/div>\n
setsweep(“sweep”,”name”,”cylinder_radius_sweep”);<\/div>\n
# define the parameter radius<\/div>\n
para_r = struct;<\/div>\n
para_r.Name = “radius”;<\/div>\n
para_r.Parameter = “::model::radius”;<\/div>\n
para_r.Type = “Length”;<\/div>\n
para_r.Start = radius_np_min;<\/div>\n
para_r.Stop = radius_np_max;<\/div>\n
para_r.Units = “microns”;<\/div>\n
# add the parameters to the sweep<\/div>\n
addsweepparameter(“cylinder_height_sweep”, para_r);<\/div>\n
# define results<\/div>\n
result_1 = struct;<\/div>\n
result_1.Name = “R”;<\/div>\n
result_1.Result = “::model::R::T”;<\/div>\n
 <\/div>\n
result_2 = struct;<\/div>\n
result_2.Name = “T”;<\/div>\n
result_2.Result = “::model::T::T”;<\/div>\n
 <\/div>\n
# add the results R & T to the sweep<\/div>\n
addsweepresult(“thickness_sweep_script”, result_1);<\/div>\n
addsweepresult(“thickness_sweep_script”, result_2);<\/div>\n

 <\/p>\n","protected":false},"template":"","class_list":["post-326278","topic","type-topic","status-closed","hentry"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t