Python FieldAtPoint Beispiele

Beispiel #1

Datei: oblique_GtoC_SPSTD.py Projekt: steve1029/2nd_paper

#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

Beispiel #2

Datei: band_1D_free.py Projekt: steve1029/2nd_paper

########## 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 -------------------#
#------------------------------------------------------------------#

Beispiel #3

Datei: oblique_ExHyHz_yzplane_SHPF.py Projekt: steve1029/2nd_paper

#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()

Beispiel #4

#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 ------------------------#
#------------------------------------------------------------------#