dipole-on-surf/add_plasmon_monitors.lsf

groupscope("::model");
select("FDTD");
min_mesh_step=get("min mesh step");
z_mon = get("z");
z_span = get("z span");
select("sub-Au");
sub_Au_z_max = get("z max");
#remove old monitors
groupscope("::model::mon");
selectall;
delete;
#########################################
mkm = 1e-6;
z_arr = [5, 20];
for(x=1:10) {
    for (media=1:2){
        if (media == 1){shift_sign = 1;} else {shift_sign = -1;}
        for(i=1:length(z_arr)){
            addpower;
            zshift_all = z_arr(i)*1e-9;
            set("name","mon_x"+num2str(x)+"mkm_media"+num2str(media)+"_zshift"
            +num2str(z_arr(i))+"nm");
            set("monitor type",1);  # 1 = point, 2 = linear x, 3 = linear y, 4 = linear z, 5 = 2D x-normal, 6 = 2D y-normal, 7 = 2D z-normal, 8 = 3D
            set("x",x*mkm);
            set("y",0);
            set("z",sub_Au_z_max + shift_sign*zshift_all);
        }
    }
}
groupscope("::model");
for(i=1:2){
    select("monitor_"+num2str(i));
    delete;
    addpower;
    set("name","monitor_"+num2str(i));
    set("monitor type",4);  # 1 = point, 2 = linear x, 3 = linear y, 4 = linear z, 5 = 2D x-normal, 6 = 2D y-normal, 7 = 2D z-normal, 8 = 3D
    set("x",0.15*mkm+0.650*(i-1)*mkm);
    set("y",0);
    set("z", z_mon);
    set("z span", z_span);
}