#TF.save_eps_mu(savedir) #sys.exit() #------------------------------------------------------------------# #-------------------- Collector object settings -------------------# #------------------------------------------------------------------# loc1 = (120 * nm, 100 * nm) loc2 = (450 * nm, 300 * nm) loc3 = (90 * nm, 270 * nm) loc4 = (300 * nm, 480 * nm) loc5 = (250 * nm, 250 * nm) filename = "lamby{:05d}".format(int(round(wvlen / nm))) fap1 = collector.FieldAtPoint("fap1", savedir + "graph/{}".format(filename), TF, loc1, 'cupy') fap2 = collector.FieldAtPoint("fap2", savedir + "graph/{}".format(filename), TF, loc2, 'cupy') fap3 = collector.FieldAtPoint("fap3", savedir + "graph/{}".format(filename), TF, loc3, 'cupy') fap4 = collector.FieldAtPoint("fap4", savedir + "graph/{}".format(filename), TF, loc4, 'cupy') fap5 = collector.FieldAtPoint("fap5", savedir + "graph/{}".format(filename), TF, loc5, 'cupy') #------------------------------------------------------------------# #-------------------- Graphtool object settings -------------------# #------------------------------------------------------------------# # Set plotfield options plot_per = 100
########## Save eps, mu data. #TF.save_eps_mu(savedir) #sys.exit() #------------------------------------------------------------------# #-------------------- Collector object settings -------------------# #------------------------------------------------------------------# loc1 = (Lx / 2, 120 * nm, 40 * nm) loc2 = (Lx / 2, 200 * nm, 60 * nm) loc3 = (Lx / 2, 370 * nm, 100 * nm) loc4 = (Lx / 2, 500 * nm, 90 * nm) fap1 = collector.FieldAtPoint( "fap1", savedir + "graph/{:04d}".format(round(wvlen / nm)), TF, loc1, 'cupy') fap2 = collector.FieldAtPoint( "fap2", savedir + "graph/{:04d}".format(round(wvlen / nm)), TF, loc2, 'cupy') fap3 = collector.FieldAtPoint( "fap3", savedir + "graph/{:04d}".format(round(wvlen / nm)), TF, loc3, 'cupy') fap4 = collector.FieldAtPoint( "fap4", savedir + "graph/{:04d}".format(round(wvlen / nm)), TF, loc4, 'cupy') #------------------------------------------------------------------# #-------------------- Graphtool object settings -------------------# #------------------------------------------------------------------#
#center3 = ( Lx, 0) #center4 = ( Lx, Lz) #cylinder1 = structure.Cylinder(TF, 'x', radius, height, center1, 8.9, 1.) #cylinder2 = structure.Cylinder(TF, 'y', radius, height, center2, 8.9, 1.) #cylinder3 = structure.Cylinder(TF, 'y', radius, height, center3, 8.9, 1.) #cylinder4 = structure.Cylinder(TF, 'y', radius, height, center4, 8.9, 1.) ########## Save eps, mu data. TF.save_eps_mu(savedir) #------------------------------------------------------------------# #-------------------- Collector object settings -------------------# #------------------------------------------------------------------# loc = (Lx / 2, 200 * nm, 200 * nm) field_at_point = collector.FieldAtPoint("fap1", savedir + "graph/", TF, loc, 'cupy') #------------------------------------------------------------------# #-------------------- Graphtool object settings -------------------# #------------------------------------------------------------------# # Set plotfield options plot_per = 1000 TFgraphtool = plotfield.Graphtool(TF, 'TF', savedir) #------------------------------------------------------------------# #------------------------ Time loop begins ------------------------# #------------------------------------------------------------------# # Save what time the simulation begins. start_time = datetime.datetime.now()
#cylinder2 = structure.Cylinder(TF, 'y', radius, height, center2, 8.9, 1.) #cylinder3 = structure.Cylinder(TF, 'y', radius, height, center3, 8.9, 1.) #cylinder4 = structure.Cylinder(TF, 'y', radius, height, center4, 8.9, 1.) ########## Save eps, mu data. TF.save_eps_mu(savedir) #------------------------------------------------------------------# #-------------------- Collector object settings -------------------# #------------------------------------------------------------------# loc1 = (Lx / 2, 450 * nm, 450 * nm) loc2 = (Lx / 2, 70 * nm, 70 * nm) loc3 = (Lx / 2, 200 * nm, 170 * nm) loc4 = (Lx / 3, 380 * nm, 270 * nm) fap1 = collector.FieldAtPoint("fap1", savedir + "graph/", TF, loc1, 'cupy') fap2 = collector.FieldAtPoint("fap2", savedir + "graph/", TF, loc2, 'cupy') fap3 = collector.FieldAtPoint("fap3", savedir + "graph/", TF, loc3, 'cupy') fap4 = collector.FieldAtPoint("fap4", savedir + "graph/", TF, loc4, 'cupy') #------------------------------------------------------------------# #-------------------- Graphtool object settings -------------------# #------------------------------------------------------------------# # Set plotfield options plot_per = 100 TFgraphtool = plotfield.Graphtool(TF, 'TF', savedir) #------------------------------------------------------------------# #------------------------ Time loop begins ------------------------# #------------------------------------------------------------------#