In [1]:
# -*- coding: utf-8 -*-

import meep as mp
import math
import numpy as np
import numpy.matlib
import matplotlib.pyplot as plt

resolution = 50                 # pixels/um

dpml = 1.0                      # PML thickness
sz = 10 + 2*dpml
cell_size = mp.Vector3(z=sz)
pml_layers = [mp.PML(dpml)]

wvl_min = 0.4                   # min wavelength
wvl_max = 0.8                   # max wavelength
fmin = 1/wvl_max                # min frequency
fmax = 1/wvl_min                # max frequency
fcen = 0.5*(fmin+fmax)          # center frequency
df = fmax-fmin                  # frequency width
nfreq = 50                      # number of frequency bins

def planar_reflectance(theta):    
    # rotation angle (in degrees) of source: CCW around Y axis, 0 degrees along +Z axis
    theta_r = math.radians(theta)

    # plane of incidence is XZ; rotate counter clockwise (CCW) about y-axis
    k = mp.Vector3(z=fmin).rotate(mp.Vector3(y=1), theta_r)
    
    # if normal incidence, force number of dimensions to be 1
    if theta_r == 0:
        dimensions = 1
    else:
        dimensions = 3
    
    sources = [mp.Source(mp.GaussianSource(fcen,fwidth=df), component=mp.Ex, center=mp.Vector3(z=-0.5*sz+dpml))]

    sim = mp.Simulation(cell_size=cell_size,
                        boundary_layers=pml_layers,
                        sources=sources,
                        k_point=k,
                        dimensions=dimensions,
                        resolution=resolution)

    refl_fr = mp.FluxRegion(center=mp.Vector3(z=-0.25*sz))
    refl = sim.add_flux(fcen, df, nfreq, refl_fr)
    
    sim.run(until_after_sources=mp.stop_when_fields_decayed(50, mp.Ex, mp.Vector3(z=-0.5*sz+dpml), 1e-9))

    empty_flux = mp.get_fluxes(refl)
    empty_data = sim.get_flux_data(refl)

    sim.reset_meep()

    # add a block with n=3.5 for the air-dielectric interface
    geometry = [mp.Block(mp.Vector3(mp.inf,mp.inf,0.5*sz), center=mp.Vector3(z=0.25*sz), material=mp.Medium(index=3.5))]

    sim = mp.Simulation(cell_size=cell_size,
                        geometry=geometry,
                        boundary_layers=pml_layers,
                        sources=sources,
                        k_point=k,
                        dimensions=dimensions,
                        resolution=resolution)

    refl = sim.add_flux(fcen, df, nfreq, refl_fr)
    sim.load_minus_flux_data(refl, empty_data)

    sim.run(until_after_sources=mp.stop_when_fields_decayed(50, mp.Ex, mp.Vector3(z=-0.5*sz+dpml), 1e-9))

    refl_flux = mp.get_fluxes(refl)
    freqs = mp.get_flux_freqs(refl)

    wvls = np.empty(nfreq)
    theta_out = np.empty(nfreq)
    R = np.empty(nfreq)
    for i in range(nfreq):
        wvls[i] = 1/freqs[i]
        theta_out[i] = math.degrees(math.asin(k.x/freqs[i]))
        R[i] = -refl_flux[i]/empty_flux[i]
        print("refl:, {}, {}, {}, {}".format(k.x,wvls[i],theta_out[i],R[i]))
        
    return k.x*np.ones(nfreq), wvls, theta_out, R 
In [2]:
theta_in = np.arange(0,85,5)
wvl = np.empty(nfreq)
kxs = np.empty((nfreq,theta_in.size))
thetas = np.empty((nfreq,theta_in.size))
Rmeep = np.empty((nfreq,theta_in.size))

for j in range(theta_in.size):
  kxs[:,j], wvl, thetas[:,j], Rmeep[:,j] = planar_reflectance(theta_in[j])

# create a 2d matrix for the wavelength by repeating the column vector for each angle
wvls = np.transpose(np.matlib.repmat(wvl,theta_in.size,1))
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.25332329653323415 / 0.25332329653323415 = 1.0
field decay(t = 100.01): 6.806395978139867e-16 / 0.25332329653323415 = 2.6868417043700194e-15
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.25332329652480207 / 0.25332329652480207 = 1.0
field decay(t = 100.01): 1.9736380723733672e-11 / 0.25332329652480207 = 7.790985272371642e-11
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 0.0, 0.8, 0.0, 0.29476330384323207
refl:, 0.0, 0.784, 0.0, 0.29416578611798405
refl:, 0.0, 0.7686274509803922, 0.0, 0.2935560666132851
refl:, 0.0, 0.7538461538461539, 0.0, 0.29293290546540196
refl:, 0.0, 0.739622641509434, 0.0, 0.2922951429790378
refl:, 0.0, 0.7259259259259259, 0.0, 0.29164324229177097
refl:, 0.0, 0.7127272727272727, 0.0, 0.29097731021514417
refl:, 0.0, 0.7, 0.0, 0.29029575759541565
refl:, 0.0, 0.6877192982456141, 0.0, 0.2895975033066778
refl:, 0.0, 0.6758620689655173, 0.0, 0.2888834125810835
refl:, 0.0, 0.664406779661017, 0.0, 0.2881539488987047
refl:, 0.0, 0.6533333333333333, 0.0, 0.2874077841401458
refl:, 0.0, 0.6426229508196721, 0.0, 0.2866443691325259
refl:, 0.0, 0.632258064516129, 0.0, 0.2858651486082286
refl:, 0.0, 0.6222222222222222, 0.0, 0.28507064287177303
refl:, 0.0, 0.6124999999999999, 0.0, 0.2842593661161744
refl:, 0.0, 0.6030769230769231, 0.0, 0.28343094611850356
refl:, 0.0, 0.593939393939394, 0.0, 0.28258688014584116
refl:, 0.0, 0.5850746268656717, 0.0, 0.2817272035905085
refl:, 0.0, 0.5764705882352942, 0.0, 0.2808500622387931
refl:, 0.0, 0.5681159420289855, 0.0, 0.27995515394321197
refl:, 0.0, 0.56, 0.0, 0.27904384994339343
refl:, 0.0, 0.552112676056338, 0.0, 0.27811574176551984
refl:, 0.0, 0.5444444444444444, 0.0, 0.2771688207349352
refl:, 0.0, 0.536986301369863, 0.0, 0.2762029589907331
refl:, 0.0, 0.5297297297297298, 0.0, 0.2752195125027274
refl:, 0.0, 0.5226666666666667, 0.0, 0.2742178237240486
refl:, 0.0, 0.5157894736842105, 0.0, 0.27319583895650085
refl:, 0.0, 0.509090909090909, 0.0, 0.27215367552975717
refl:, 0.0, 0.5025641025641026, 0.0, 0.27109276268116617
refl:, 0.0, 0.4962025316455696, 0.0, 0.2700121717612113
refl:, 0.0, 0.49, 0.0, 0.26890977948431494
refl:, 0.0, 0.4839506172839506, 0.0, 0.26778602305927346
refl:, 0.0, 0.47804878048780486, 0.0, 0.26664234270279846
refl:, 0.0, 0.47228915662650606, 0.0, 0.2654773959072893
refl:, 0.0, 0.4666666666666667, 0.0, 0.2642890605183494
refl:, 0.0, 0.4611764705882353, 0.0, 0.26307815170848525
refl:, 0.0, 0.4558139534883721, 0.0, 0.2618459695575972
refl:, 0.0, 0.4505747126436782, 0.0, 0.2605907641320623
refl:, 0.0, 0.4454545454545454, 0.0, 0.25931062280334477
refl:, 0.0, 0.44044943820224725, 0.0, 0.25800677969709146
refl:, 0.0, 0.43555555555555553, 0.0, 0.2566804027177118
refl:, 0.0, 0.4307692307692308, 0.0, 0.25532959438320507
refl:, 0.0, 0.4260869565217391, 0.0, 0.2539529088676918
refl:, 0.0, 0.421505376344086, 0.0, 0.2525520862693137
refl:, 0.0, 0.41702127659574467, 0.0, 0.2511282531123407
refl:, 0.0, 0.4126315789473684, 0.0, 0.24967934554669954
refl:, 0.0, 0.4083333333333333, 0.0, 0.24820404444956207
refl:, 0.0, 0.4041237113402062, 0.0, 0.24670386250022297
refl:, 0.0, 0.4, 0.0, 0.24517862401765392
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.25242167342001043 / 0.25242167342001043 = 1.0
field decay(t = 100.01): 1.8867425254108467e-14 / 0.25242167342001043 = 7.474566267815881e-14
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.2524216734361254 / 0.2524216734361254 = 1.0
field decay(t = 100.01): 2.0310851043564926e-11 / 0.2524216734361254 = 8.046397429777174e-11
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 0.1089446784345727, 0.8, 5.0, 0.2934227118859898
refl:, 0.1089446784345727, 0.784, 4.899752997934953, 0.29287850682357097
refl:, 0.1089446784345727, 0.7686274509803922, 4.803451415694315, 0.2923194417293526
refl:, 0.1089446784345727, 0.7538461538461539, 4.710866569098618, 0.29174404751628646
refl:, 0.1089446784345727, 0.739622641509434, 4.621787131270349, 0.2911512998642181
refl:, 0.1089446784345727, 0.7259259259259259, 4.536017514803643, 0.2905417123133207
refl:, 0.1089446784345727, 0.7127272727272727, 4.45337643175598, 0.2899156301931104
refl:, 0.1089446784345727, 0.7, 4.373695609047488, 0.2892720597047358
refl:, 0.1089446784345727, 0.6877192982456141, 4.296818640028281, 0.288609702000128
refl:, 0.1089446784345727, 0.6758620689655173, 4.222599955651541, 0.28792940285545326
refl:, 0.1089446784345727, 0.664406779661017, 4.150903900954991, 0.28723211458816683
refl:, 0.1089446784345727, 0.6533333333333333, 4.081603904475513, 0.28651676808837245
refl:, 0.1089446784345727, 0.6426229508196721, 4.014581729857545, 0.2857826486835516
refl:, 0.1089446784345727, 0.632258064516129, 3.9497268003123533, 0.28503114405854746
refl:, 0.1089446784345727, 0.6222222222222222, 3.886935587780479, 0.2842631868294254
refl:, 0.1089446784345727, 0.6124999999999999, 3.826111059675504, 0.28347749048119336
refl:, 0.1089446784345727, 0.6030769230769231, 3.7671621769699377, 0.28267341154462916
refl:, 0.1089446784345727, 0.593939393939394, 3.7100034381453693, 0.28185242101019087
refl:, 0.1089446784345727, 0.5850746268656717, 3.6545544641874574, 0.2810149315091566
refl:, 0.1089446784345727, 0.5764705882352942, 3.6007396203768605, 0.2801592191834353
refl:, 0.1089446784345727, 0.5681159420289855, 3.5484876711229183, 0.2792846926505535
refl:, 0.1089446784345727, 0.56, 3.4977314645183832, 0.2783927047249737
refl:, 0.1089446784345727, 0.552112676056338, 3.4484076436699387, 0.2774832537476156
refl:, 0.1089446784345727, 0.5444444444444444, 3.4004563821883713, 0.27655440028581263
refl:, 0.1089446784345727, 0.536986301369863, 3.353821141510494, 0.2756056895508914
refl:, 0.1089446784345727, 0.5297297297297298, 3.3084484479779683, 0.2746385302061814
refl:, 0.1089446784345727, 0.5226666666666667, 3.2642876878206, 0.2736526713566295
refl:, 0.1089446784345727, 0.5157894736842105, 3.221290918387684, 0.2726460490596169
refl:, 0.1089446784345727, 0.509090909090909, 3.1794126941439105, 0.27161844709937716
refl:, 0.1089446784345727, 0.5025641025641026, 3.1386099060992034, 0.27057141563726533
refl:, 0.1089446784345727, 0.4962025316455696, 3.098841633477275, 0.2695044090077289
refl:, 0.1089446784345727, 0.49, 3.060069006547709, 0.2684152052151545
refl:, 0.1089446784345727, 0.4839506172839506, 3.0222550796530636, 0.26730394106050426
refl:, 0.1089446784345727, 0.47804878048780486, 2.985364713557375, 0.2661722224160734
refl:, 0.1089446784345727, 0.47228915662650606, 2.9493644663270095, 0.2650190380989821
refl:, 0.1089446784345727, 0.4666666666666667, 2.914222492030252, 0.26384212802287793
refl:, 0.1089446784345727, 0.4611764705882353, 2.879908446609464, 0.26264201676437865
refl:, 0.1089446784345727, 0.4558139534883721, 2.84639340033999, 0.2614202058647917
refl:, 0.1089446784345727, 0.4505747126436782, 2.8136497563440477, 0.26017524659219626
refl:, 0.1089446784345727, 0.4454545454545454, 2.7816511746763815, 0.258905029027019
refl:, 0.1089446784345727, 0.44044943820224725, 2.750372501542037, 0.2576105253535685
refl:, 0.1089446784345727, 0.43555555555555553, 2.719789703245832, 0.256293125487968
refl:, 0.1089446784345727, 0.4307692307692308, 2.689879804508449, 0.2549512279458179
refl:, 0.1089446784345727, 0.4260869565217391, 2.660620830815881, 0.2535831260087174
refl:, 0.1089446784345727, 0.421505376344086, 2.631991754497782, 0.25219031084070415
refl:, 0.1089446784345727, 0.41702127659574467, 2.603972444256211, 0.25077428012312447
refl:, 0.1089446784345727, 0.4126315789473684, 2.576543617889841, 0.24933314662117714
refl:, 0.1089446784345727, 0.4083333333333333, 2.549686797979975, 0.2478653741944476
refl:, 0.1089446784345727, 0.4041237113402062, 2.5233842703240543, 0.24637235281997996
refl:, 0.1089446784345727, 0.4, 2.4976190449198983, 0.2448541987123429
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.24974540035413878 / 0.24974540035413878 = 1.0
field decay(t = 100.01): 6.006906595531682e-14 / 0.24974540035413878 = 2.40521210281106e-13
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.24974540044608928 / 0.24974540044608928 = 1.0
field decay(t = 100.01): 2.1519308385946042e-11 / 0.24974540044608928 = 8.61649838095467e-11
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 0.2170602220836629, 0.8, 10.0, 0.2893781746253207
refl:, 0.2170602220836629, 0.784, 9.798006528153513, 0.288998732433702
refl:, 0.2170602220836629, 0.7686274509803922, 9.604050171837292, 0.2885938951899122
refl:, 0.2170602220836629, 0.7538461538461539, 9.417658416993296, 0.28816419107846264
refl:, 0.2170602220836629, 0.739622641509434, 9.238395240840497, 0.28770900366948954
refl:, 0.2170602220836629, 0.7259259259259259, 9.06585764090149, 0.28722834374816003
refl:, 0.2170602220836629, 0.7127272727272727, 8.899672554443574, 0.2867235913729682
refl:, 0.2170602220836629, 0.7, 8.739494117841588, 0.28619525258422535
refl:, 0.2170602220836629, 0.6877192982456141, 8.585001222725978, 0.28564218452120266
refl:, 0.2170602220836629, 0.6758620689655173, 8.435895331947279, 0.2850647891546287
refl:, 0.2170602220836629, 0.664406779661017, 8.291898523577625, 0.2844650049217599
refl:, 0.2170602220836629, 0.6533333333333333, 8.152751735551202, 0.2838430637467787
refl:, 0.2170602220836629, 0.6426229508196721, 8.018213187256704, 0.28319797552226944
refl:, 0.2170602220836629, 0.632258064516129, 7.88805695754783, 0.2825306580041294
refl:, 0.2170602220836629, 0.6222222222222222, 7.762071701325296, 0.2818429693994039
refl:, 0.2170602220836629, 0.6124999999999999, 7.640059489140416, 0.2811345740291295
refl:, 0.2170602220836629, 0.6030769230769231, 7.521834756238996, 0.28040433725297315
refl:, 0.2170602220836629, 0.593939393939394, 7.407223349155971, 0.27965322172446067
refl:, 0.2170602220836629, 0.5850746268656717, 7.296061659428915, 0.2788825653708164
refl:, 0.2170602220836629, 0.5764705882352942, 7.188195835257705, 0.27809144849085615
refl:, 0.2170602220836629, 0.5681159420289855, 7.083481063027822, 0.2772785981174172
refl:, 0.2170602220836629, 0.56, 6.981780911561047, 0.2764449243213644
refl:, 0.2170602220836629, 0.552112676056338, 6.882966732780441, 0.2755914805667849
refl:, 0.2170602220836629, 0.5444444444444444, 6.786917113194022, 0.27471697321444166
refl:, 0.2170602220836629, 0.536986301369863, 6.693517371228444, 0.27382008373306266
refl:, 0.2170602220836629, 0.5297297297297298, 6.602659095992853, 0.2729019522035499
refl:, 0.2170602220836629, 0.5226666666666667, 6.514239723534344, 0.27196350015165877
refl:, 0.2170602220836629, 0.5157894736842105, 6.428162147069652, 0.27100303851505764
refl:, 0.2170602220836629, 0.509090909090909, 6.3443343580501015, 0.27001939086671356
refl:, 0.2170602220836629, 0.5025641025641026, 6.262669115245526, 0.26901406912439724
refl:, 0.2170602220836629, 0.4962025316455696, 6.1830836393232005, 0.2679876733928022
refl:, 0.2170602220836629, 0.49, 6.105499330654852, 0.2669381026643071
refl:, 0.2170602220836629, 0.4839506172839506, 6.029841508312739, 0.265864537999668
refl:, 0.2170602220836629, 0.47804878048780486, 5.956039168418179, 0.26476870360049726
refl:, 0.2170602220836629, 0.47228915662650606, 5.884024760185935, 0.26365066570667506
refl:, 0.2170602220836629, 0.4666666666666667, 5.813733978168244, 0.26250810411848724
refl:, 0.2170602220836629, 0.4611764705882353, 5.745105569345401, 0.26134058873292715
refl:, 0.2170602220836629, 0.4558139534883721, 5.678081153837623, 0.26014984735445634
refl:, 0.2170602220836629, 0.4505747126436782, 5.612605058127373, 0.25893548643505515
refl:, 0.2170602220836629, 0.4454545454545454, 5.5486241597837695, 0.2576951510120501
refl:, 0.2170602220836629, 0.44044943820224725, 5.486087742772728, 0.2564288095400968
refl:, 0.2170602220836629, 0.43555555555555553, 5.424947362519055, 0.2551382775804875
refl:, 0.2170602220836629, 0.4307692307692308, 5.365156719961124, 0.25382296989852066
refl:, 0.2170602220836629, 0.4260869565217391, 5.306671543905553, 0.25248068938520546
refl:, 0.2170602220836629, 0.421505376344086, 5.249449481049872, 0.2511120329000034
refl:, 0.2170602220836629, 0.41702127659574467, 5.19344999309547, 0.2497191778959394
refl:, 0.2170602220836629, 0.4126315789473684, 5.138634260422533, 0.2483012320567707
refl:, 0.2170602220836629, 0.4083333333333333, 5.08496509184317, 0.2468560224651321
refl:, 0.2170602220836629, 0.4041237113402062, 5.032406839989342, 0.24538435149961998
refl:, 0.2170602220836629, 0.4, 4.980925321928872, 0.2438873823711191
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.24537918139687384 / 0.24537918139687384 = 1.0
field decay(t = 100.01): 1.2755953357084837e-13 / 0.24537918139687384 = 5.198466016745522e-13
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.24537918161848085 / 0.24537918161848085 = 1.0
field decay(t = 100.01): 2.2796495621609012e-11 / 0.24537918161848085 = 9.290313657111033e-11
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 0.3235238063781509, 0.8, 14.999999999999998, 0.28256373042287475
refl:, 0.3235238063781509, 0.784, 14.693171512000124, 0.28247093802487827
refl:, 0.3235238063781509, 0.7686274509803922, 14.398780921441814, 0.28233355058613946
refl:, 0.3235238063781509, 0.7538461538461539, 14.116078899389818, 0.2821545309608155
refl:, 0.3235238063781509, 0.739622641509434, 13.844375746673084, 0.2819358770069733
refl:, 0.3235238063781509, 0.7259259259259259, 13.583035518835887, 0.2816776835689441
refl:, 0.3235238063781509, 0.7127272727272727, 13.331470838933798, 0.28138133257604653
refl:, 0.3235238063781509, 0.7, 13.089138305160036, 0.281049646363028
refl:, 0.3235238063781509, 0.6877192982456141, 12.855534414514684, 0.28068351331308133
refl:, 0.3235238063781509, 0.6758620689655173, 12.630191935533823, 0.28028281405469785
refl:, 0.3235238063781509, 0.664406779661017, 12.412676672931594, 0.2798495314094937
refl:, 0.3235238063781509, 0.6533333333333333, 12.202584575236058, 0.27938611173524625
refl:, 0.3235238063781509, 0.6426229508196721, 11.999539143408532, 0.2788927304404999
refl:, 0.3235238063781509, 0.632258064516129, 11.803189104258525, 0.27836932361278227
refl:, 0.3235238063781509, 0.6222222222222222, 11.613206317390016, 0.27781781071798484
refl:, 0.3235238063781509, 0.6124999999999999, 11.429283888592414, 0.27723982281143944
refl:, 0.3235238063781509, 0.6030769230769231, 11.25113446614539, 0.27663483656101373
refl:, 0.3235238063781509, 0.593939393939394, 11.078488699542484, 0.2760025164435567
refl:, 0.3235238063781509, 0.5850746268656717, 10.91109384273803, 0.27534443882671783
refl:, 0.3235238063781509, 0.5764705882352942, 10.748712486253877, 0.2746616477877146
refl:, 0.3235238063781509, 0.5681159420289855, 10.591121404404543, 0.2739530013321228
refl:, 0.3235238063781509, 0.56, 10.438110505558328, 0.27321798287921667
refl:, 0.3235238063781509, 0.552112676056338, 10.289481874787974, 0.2724582796019277
refl:, 0.3235238063781509, 0.5444444444444444, 10.145048899510067, 0.27167450266463894
refl:, 0.3235238063781509, 0.536986301369863, 10.004635469795673, 0.2708649570485523
refl:, 0.3235238063781509, 0.5297297297297298, 9.868075245978739, 0.27002944005333146
refl:, 0.3235238063781509, 0.5226666666666667, 9.735210987013529, 0.2691698987799361
refl:, 0.3235238063781509, 0.5157894736842105, 9.60589393375409, 0.26828628327229254
refl:, 0.3235238063781509, 0.509090909090909, 9.479983241961918, 0.2673766390145325
refl:, 0.3235238063781509, 0.5025641025641026, 9.3573454604044, 0.266441343276572
refl:, 0.3235238063781509, 0.4962025316455696, 9.237854049896542, 0.2654822435375998
refl:, 0.3235238063781509, 0.49, 9.121388939570695, 0.26449858196569564
refl:, 0.3235238063781509, 0.4839506172839506, 9.00783611704105, 0.2634884799544352
refl:, 0.3235238063781509, 0.47804878048780486, 8.897087249467761, 0.26245270185974406
refl:, 0.3235238063781509, 0.47228915662650606, 8.789039332825531, 0.26139274110511584
refl:, 0.3235238063781509, 0.4666666666666667, 8.683594366947904, 0.26030735036150837
refl:, 0.3235238063781509, 0.4611764705882353, 8.5806590541555, 0.2591947457863499
refl:, 0.3235238063781509, 0.4558139534883721, 8.480144519487721, 0.25805589551088276
refl:, 0.3235238063781509, 0.4505747126436782, 8.381966050745921, 0.256892047194308
refl:, 0.3235238063781509, 0.4454545454545454, 8.286042856724524, 0.2557015480946698
refl:, 0.3235238063781509, 0.44044943820224725, 8.192297842157469, 0.2544827165503608
refl:, 0.3235238063781509, 0.43555555555555553, 8.100657398042404, 0.25323697315876975
refl:, 0.3235238063781509, 0.4307692307692308, 8.011051206126568, 0.25196549708688976
refl:, 0.3235238063781509, 0.4260869565217391, 7.923412056447124, 0.2506663282693837
refl:, 0.3235238063781509, 0.421505376344086, 7.837675676917115, 0.2493383700330171
refl:, 0.3235238063781509, 0.41702127659574467, 7.753780574036344, 0.24798388324916928
refl:, 0.3235238063781509, 0.4126315789473684, 7.671667883886533, 0.24660386238549814
refl:, 0.3235238063781509, 0.4083333333333333, 7.591281232641979, 0.2451961050843521
refl:, 0.3235238063781509, 0.4041237113402062, 7.512566605892175, 0.2437600099079589
refl:, 0.3235238063781509, 0.4, 7.435472226131853, 0.24229747168145319
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.23946072431081186 / 0.23946072431081186 = 1.0
field decay(t = 100.01): 2.143367304660103e-13 / 0.23946072431081186 = 8.950809410724429e-13
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.23946072470582727 / 0.23946072470582727 = 1.0
field decay(t = 100.01): 2.2924410648414684e-11 / 0.23946072470582727 = 9.573348897434795e-11
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 0.4275251791570859, 0.8, 20.0, 0.27287026974445927
refl:, 0.4275251791570859, 0.784, 19.583468236198428, 0.2732024101218738
refl:, 0.4275251791570859, 0.7686274509803922, 19.184283570762837, 0.2734622444200681
refl:, 0.4275251791570859, 0.7538461538461539, 18.80136282780684, 0.27365329059595683
refl:, 0.4275251791570859, 0.739622641509434, 18.433712921276232, 0.27378076185161915
refl:, 0.4275251791570859, 0.7259259259259259, 18.080421505473552, 0.27384851240423785
refl:, 0.4275251791570859, 0.7127272727272727, 17.74064878668832, 0.2738584846379507
refl:, 0.4275251791570859, 0.7, 17.413620328364388, 0.27381376974644933
refl:, 0.4275251791570859, 0.6877192982456141, 17.098620709730763, 0.27371850538265774
refl:, 0.4275251791570859, 0.6758620689655173, 16.794987920274295, 0.2735750966910372
refl:, 0.4275251791570859, 0.664406779661017, 16.50210839086207, 0.2733847949229058
refl:, 0.4275251791570859, 0.6533333333333333, 16.21941257752279, 0.2731503266875855
refl:, 0.4275251791570859, 0.6426229508196721, 15.946371026493122, 0.2728749017452174
refl:, 0.4275251791570859, 0.632258064516129, 15.682490859619808, 0.2725596258354739
refl:, 0.4275251791570859, 0.6222222222222222, 15.427312627971409, 0.2722049153191827
refl:, 0.4275251791570859, 0.6124999999999999, 15.18040748886694, 0.2718130795196188
refl:, 0.4275251791570859, 0.6030769230769231, 14.941374667722823, 0.2713863599041548
refl:, 0.4275251791570859, 0.593939393939394, 14.709839171355561, 0.2709245158308668
refl:, 0.4275251791570859, 0.5850746268656717, 14.485449723819883, 0.2704275378957415
refl:, 0.4275251791570859, 0.5764705882352942, 14.267876899642369, 0.26989771202120866
refl:, 0.4275251791570859, 0.5681159420289855, 14.056811432538773, 0.26933628796739983
refl:, 0.4275251791570859, 0.56, 13.851962680467771, 0.26874214027712073
refl:, 0.4275251791570859, 0.552112676056338, 13.653057230248155, 0.2681155761222093
refl:, 0.4275251791570859, 0.5444444444444444, 13.459837627011849, 0.2674588673074579
refl:, 0.4275251791570859, 0.536986301369863, 13.272061215531336, 0.2667723151469476
refl:, 0.4275251791570859, 0.5297297297297298, 13.089499081989581, 0.26605454446111854
refl:, 0.4275251791570859, 0.5226666666666667, 12.911935086088024, 0.2653063511402134
refl:, 0.4275251791570859, 0.5157894736842105, 12.73916497454351, 0.2645297343354968
refl:, 0.4275251791570859, 0.509090909090909, 12.570995568032437, 0.2637242645202348
refl:, 0.4275251791570859, 0.5025641025641026, 12.407244014521105, 0.2628885752055156
refl:, 0.4275251791570859, 0.4962025316455696, 12.247737102692609, 0.26202377900308255
refl:, 0.4275251791570859, 0.49, 12.092310629857906, 0.2611315306153532
refl:, 0.4275251791570859, 0.4839506172839506, 11.94080881933422, 0.2602107728685866
refl:, 0.4275251791570859, 0.47804878048780486, 11.793083782798943, 0.25926014291141586
refl:, 0.4275251791570859, 0.47228915662650606, 11.648995023590732, 0.25828106391345235
refl:, 0.4275251791570859, 0.4666666666666667, 11.508408977339421, 0.25727477704208235
refl:, 0.4275251791570859, 0.4611764705882353, 11.37119858666983, 0.2562395153461572
refl:, 0.4275251791570859, 0.4558139534883721, 11.23724290704701, 0.2551740956437644
refl:, 0.4275251791570859, 0.4505747126436782, 11.106426741117268, 0.2540803028372474
refl:, 0.4275251791570859, 0.4454545454545454, 10.978640299154753, 0.25295873317085843
refl:, 0.4275251791570859, 0.44044943820224725, 10.853778883451254, 0.2518070732573948
refl:, 0.4275251791570859, 0.43555555555555553, 10.731742594690267, 0.2506247305646683
refl:, 0.4275251791570859, 0.4307692307692308, 10.61243605852877, 0.249413825082604
refl:, 0.4275251791570859, 0.4260869565217391, 10.495768170772996, 0.2481743775271508
refl:, 0.4275251791570859, 0.421505376344086, 10.381651859681224, 0.24690419443990608
refl:, 0.4275251791570859, 0.41702127659574467, 10.270003864057921, 0.2456037441531416
refl:, 0.4275251791570859, 0.4126315789473684, 10.160744525922071, 0.24427554305470142
refl:, 0.4275251791570859, 0.4083333333333333, 10.053797596639106, 0.2429193465329389
refl:, 0.4275251791570859, 0.4041237113402062, 9.949090055502005, 0.24153322642643202
refl:, 0.4275251791570859, 0.4, 9.846551939834079, 0.2401180039660742
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.23217567515545764 / 0.23217567515545764 = 1.0
field decay(t = 100.01): 3.1895846890533654e-13 / 0.23217567515545764 = 1.3737807317315728e-12
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.23217567572617967 / 0.23217567572617967 = 1.0
field decay(t = 100.01): 2.511386584143695e-11 / 0.23217567572617967 = 1.0816751480484724e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 0.5282728271758743, 0.8, 25.0, 0.2601481976279196
refl:, 0.5282728271758743, 0.784, 24.46679999189225, 0.2610696372423531
refl:, 0.5282728271758743, 0.7686274509803922, 23.95662859236144, 0.2618770325163737
refl:, 0.5282728271758743, 0.7538461538461539, 23.467975961566616, 0.2625788276586943
refl:, 0.5282728271758743, 0.739622641509434, 22.99946566139384, 0.2631817618626755
refl:, 0.5282728271758743, 0.7259259259259259, 22.54983979533518, 0.2636926798419164
refl:, 0.5282728271758743, 0.7127272727272727, 22.117946146436154, 0.2641189902441729
refl:, 0.5282728271758743, 0.7, 21.70272699899909, 0.2644658127354963
refl:, 0.5282728271758743, 0.6877192982456141, 21.303209386217617, 0.26473707694833154
refl:, 0.5282728271758743, 0.6758620689655173, 20.91849655102253, 0.26493869102120726
refl:, 0.5282728271758743, 0.664406779661017, 20.54776044367494, 0.26507564513451415
refl:, 0.5282728271758743, 0.6533333333333333, 20.19023510896709, 0.2651500052068409
refl:, 0.5282728271758743, 0.6426229508196721, 19.84521083974841, 0.26516501686911953
refl:, 0.5282728271758743, 0.632258064516129, 19.51202899301187, 0.26512536746580095
refl:, 0.5282728271758743, 0.6222222222222222, 19.190077380825585, 0.26503310672375896
refl:, 0.5282728271758743, 0.6124999999999999, 18.878786161658894, 0.26488908132632305
refl:, 0.5282728271758743, 0.6030769230769231, 18.577624168665142, 0.26469630209606193
refl:, 0.5282728271758743, 0.593939393939394, 18.28609562066815, 0.26445767784614754
refl:, 0.5282728271758743, 0.5850746268656717, 18.003737169291743, 0.26417370471931184
refl:, 0.5282728271758743, 0.5764705882352942, 17.730115242139682, 0.26384505077481407
refl:, 0.5282728271758743, 0.5681159420289855, 17.46482364739322, 0.26347427561024556
refl:, 0.5282728271758743, 0.56, 17.207481409818527, 0.26306325708398504
refl:, 0.5282728271758743, 0.552112676056338, 16.957730812108316, 0.262611771262809
refl:, 0.5282728271758743, 0.5444444444444444, 16.715235618835717, 0.2621203557572436
refl:, 0.5282728271758743, 0.536986301369863, 16.479679463167955, 0.2615915488530281
refl:, 0.5282728271758743, 0.5297297297297298, 16.250764378950514, 0.26102646172552363
refl:, 0.5282728271758743, 0.5226666666666667, 16.028209462892256, 0.26042407309811766
refl:, 0.5282728271758743, 0.5157894736842105, 15.811749653412038, 0.25978529806673434
refl:, 0.5282728271758743, 0.509090909090909, 15.60113461429131, 0.2591127130188696
refl:, 0.5282728271758743, 0.5025641025641026, 15.396127712651408, 0.25840617419006323
refl:, 0.5282728271758743, 0.4962025316455696, 15.196505081970068, 0.2576643632602327
refl:, 0.5282728271758743, 0.49, 15.002054761894165, 0.2568889467859571
refl:, 0.5282728271758743, 0.4839506172839506, 14.81257590751694, 0.2560818324629057
refl:, 0.5282728271758743, 0.47804878048780486, 14.627878061586326, 0.2552417285733339
refl:, 0.5282728271758743, 0.47228915662650606, 14.44778048381156, 0.25436772779199046
refl:, 0.5282728271758743, 0.4666666666666667, 14.272111532051735, 0.2534618054882236
refl:, 0.5282728271758743, 0.4611764705882353, 14.100708090713411, 0.2525248442883521
refl:, 0.5282728271758743, 0.4558139534883721, 13.933415042164041, 0.25155500702726324
refl:, 0.5282728271758743, 0.4505747126436782, 13.770084777392734, 0.25055182001671567
refl:, 0.5282728271758743, 0.4454545454545454, 13.610576742526133, 0.24951717022490025
refl:, 0.5282728271758743, 0.44044943820224725, 13.454757018141454, 0.24845117408859577
refl:, 0.5282728271758743, 0.43555555555555553, 13.30249792861589, 0.24735169966979012
refl:, 0.5282728271758743, 0.4307692307692308, 13.153677679016674, 0.2462187296769865
refl:, 0.5282728271758743, 0.4260869565217391, 13.008180017272103, 0.24505436171970374
refl:, 0.5282728271758743, 0.421505376344086, 12.865893919575512, 0.24385825418409707
refl:, 0.5282728271758743, 0.41702127659574467, 12.726713297162874, 0.2426284381212972
refl:, 0.5282728271758743, 0.4126315789473684, 12.590536722774571, 0.24136614460408662
refl:, 0.5282728271758743, 0.4083333333333333, 12.457267175263851, 0.24007405600657153
refl:, 0.5282728271758743, 0.4041237113402062, 12.326811800951434, 0.238751417728497
refl:, 0.5282728271758743, 0.4, 12.199081690448809, 0.2373967326334257
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.22375100718842 / 0.22375100718842 = 1.0
field decay(t = 100.01): 4.4054832270861915e-13 / 0.22375100718842 = 1.9689221883038714e-12
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.22375100783149995 / 0.22375100783149995 = 1.0
field decay(t = 100.01): 2.604590852088861e-11 / 0.22375100783149995 = 1.1640577074183724e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 0.6249999999999999, 0.8, 29.999999999999993, 0.2442050337697129
refl:, 0.6249999999999999, 0.784, 29.34058157502373, 0.24592026654412946
refl:, 0.6249999999999999, 0.7686274509803922, 28.711017527148794, 0.2474593095717218
refl:, 0.6249999999999999, 0.7538461538461539, 28.10922128260952, 0.24883736533494075
refl:, 0.6249999999999999, 0.739622641509434, 27.53330580109674, 0.25006920079221073
refl:, 0.6249999999999999, 0.7259259259259259, 26.98155921981659, 0.25116546178265714
refl:, 0.6249999999999999, 0.7127272727272727, 26.452424118557172, 0.2521360799590939
refl:, 0.6249999999999999, 0.7, 25.944479772370002, 0.2529918131418104
refl:, 0.6249999999999999, 0.6877192982456141, 25.45642688540588, 0.2537410511682661
refl:, 0.6249999999999999, 0.6758620689655173, 24.98707439783817, 0.2543901348468849
refl:, 0.6249999999999999, 0.664406779661017, 24.53532803474182, 0.2549460427718285
refl:, 0.6249999999999999, 0.6533333333333333, 24.10018032643974, 0.25541571682136277
refl:, 0.6249999999999999, 0.6426229508196721, 23.680701877992792, 0.25580383669879514
refl:, 0.6249999999999999, 0.632258064516129, 23.276033704033082, 0.2561140687592272
refl:, 0.6249999999999999, 0.6222222222222222, 22.885380476158563, 0.2563515425795029
refl:, 0.6249999999999999, 0.6124999999999999, 22.508004555237044, 0.2565208309064006
refl:, 0.6249999999999999, 0.6030769230769231, 22.14322070144805, 0.2566238743420462
refl:, 0.6249999999999999, 0.593939393939394, 21.79039137167391, 0.25666334314703
refl:, 0.6249999999999999, 0.5850746268656717, 21.448922527676533, 0.256643740595326
refl:, 0.6249999999999999, 0.5764705882352942, 21.118259889941786, 0.25656730578199227
refl:, 0.6249999999999999, 0.5681159420289855, 20.797885581592283, 0.25643465337781024
refl:, 0.6249999999999999, 0.56, 20.487315114722662, 0.256248873599916
refl:, 0.6249999999999999, 0.552112676056338, 20.186094678183196, 0.25601329709669296
refl:, 0.6249999999999999, 0.5444444444444444, 19.89379869145801, 0.25572829145150733
refl:, 0.6249999999999999, 0.536986301369863, 19.610027594036424, 0.255394616552182
refl:, 0.6249999999999999, 0.5297297297297298, 19.33440584370919, 0.2550154350698738
refl:, 0.6249999999999999, 0.5226666666666667, 19.066580100655866, 0.25459252830565576
refl:, 0.6249999999999999, 0.5157894736842105, 18.806217577124535, 0.25412542518907427
refl:, 0.6249999999999999, 0.509090909090909, 18.553004535020655, 0.2536152551480315
refl:, 0.6249999999999999, 0.5025641025641026, 18.306644915884704, 0.2530647648651237
refl:, 0.6249999999999999, 0.4962025316455696, 18.066859089603533, 0.25247451149799194
refl:, 0.6249999999999999, 0.49, 17.833382709813034, 0.251843602956665
refl:, 0.6249999999999999, 0.4839506172839506, 17.605965665348197, 0.25117351956545597
refl:, 0.6249999999999999, 0.47804878048780486, 17.38437111831227, 0.2504665488589744
refl:, 0.6249999999999999, 0.47228915662650606, 17.168374620396104, 0.2497220477831976
refl:, 0.6249999999999999, 0.4666666666666667, 16.957763300004142, 0.24893902021832587
refl:, 0.6249999999999999, 0.4611764705882353, 16.752335113553887, 0.24811948817045582
refl:, 0.6249999999999999, 0.4558139534883721, 16.551898155026578, 0.24726486187036864
refl:, 0.6249999999999999, 0.4505747126436782, 16.35627001847215, 0.24637335126615545
refl:, 0.6249999999999999, 0.4454545454545454, 16.165277208722518, 0.24544445758208924
refl:, 0.6249999999999999, 0.44044943820224725, 15.978754596053776, 0.24448044797494067
refl:, 0.6249999999999999, 0.43555555555555553, 15.796544910968182, 0.24348141406546195
refl:, 0.6249999999999999, 0.4307692307692308, 15.618498275648548, 0.24244508607372425
refl:, 0.6249999999999999, 0.4260869565217391, 15.44447176897603, 0.24137187009594935
refl:, 0.6249999999999999, 0.421505376344086, 15.274329022304034, 0.2402639528246127
refl:, 0.6249999999999999, 0.41702127659574467, 15.107939843449019, 0.23912065251888115
refl:, 0.6249999999999999, 0.4126315789473684, 14.94517986659885, 0.23794038016947253
refl:, 0.6249999999999999, 0.4083333333333333, 14.78593022605342, 0.23672496742401775
refl:, 0.6249999999999999, 0.4041237113402062, 14.630077251904048, 0.23547697979370008
refl:, 0.6249999999999999, 0.4, 14.477512185929921, 0.23419574038931623
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.2144470639002981 / 0.2144470639002981 = 1.0
field decay(t = 100.01): 5.749510870152026e-13 / 0.2144470639002981 = 2.6810863089387505e-12
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.21444706445458073 / 0.21444706445458073 = 1.0
field decay(t = 100.01): 2.737725032991902e-11 / 0.21444706445458073 = 1.276643744205575e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 0.7169705454388076, 0.8, 35.0, 0.22482605012853205
refl:, 0.7169705454388076, 0.784, 34.201491482224874, 0.2275879297659293
refl:, 0.7169705454388076, 0.7686274509803922, 33.4413597291606, 0.2300861771421993
refl:, 0.7169705454388076, 0.7538461538461539, 32.71669424141867, 0.23234663839124975
refl:, 0.7169705454388076, 0.739622641509434, 32.02489215108529, 0.2343911194485937
refl:, 0.7169705454388076, 0.7259259259259259, 31.363616172201723, 0.23623793864435802
refl:, 0.7169705454388076, 0.7127272727272727, 30.73075960358516, 0.23790584943527182
refl:, 0.7169705454388076, 0.7, 30.12441698853809, 0.2394088473062314
refl:, 0.7169705454388076, 0.6877192982456141, 29.542859352190995, 0.24075882698160853
refl:, 0.7169705454388076, 0.6758620689655173, 28.984513173223203, 0.24196961848162934
refl:, 0.7169705454388076, 0.664406779661017, 28.447942424882076, 0.24305087670480616
refl:, 0.7169705454388076, 0.6533333333333333, 27.931833156170086, 0.24400999332366932
refl:, 0.7169705454388076, 0.6426229508196721, 27.43498018882287, 0.24485687989558083
refl:, 0.7169705454388076, 0.632258064516129, 26.95627558715663, 0.24559913215899493
refl:, 0.7169705454388076, 0.6222222222222222, 26.49469862174186, 0.246241502598612
refl:, 0.7169705454388076, 0.6124999999999999, 26.049306998353728, 0.24679042717466868
refl:, 0.7169705454388076, 0.6030769230769231, 25.619229163859046, 0.2472523465378225
refl:, 0.7169705454388076, 0.593939393939394, 25.20365753294407, 0.24763123330505407
refl:, 0.7169705454388076, 0.5850746268656717, 24.80184250561243, 0.24793067603660995
refl:, 0.7169705454388076, 0.5764705882352942, 24.41308716651858, 0.24815554480891827
refl:, 0.7169705454388076, 0.5681159420289855, 24.036742574466203, 0.24831017863638788
refl:, 0.7169705454388076, 0.56, 23.672203564580393, 0.24839671574872801
refl:, 0.7169705454388076, 0.552112676056338, 23.318904997368072, 0.24841768322180216
refl:, 0.7169705454388076, 0.5444444444444444, 22.976318398592397, 0.24837749388988659
refl:, 0.7169705454388076, 0.536986301369863, 22.643948941980693, 0.24827859346807957
refl:, 0.7169705454388076, 0.5297297297297298, 22.321332733561132, 0.24812142120897318
refl:, 0.7169705454388076, 0.5226666666666667, 22.008034362119382, 0.24790902685773963
refl:, 0.7169705454388076, 0.5157894736842105, 21.703644685073712, 0.24764492885992392
refl:, 0.7169705454388076, 0.509090909090909, 21.407778823139967, 0.24732917319523862
refl:, 0.7169705454388076, 0.5025641025641026, 21.120074340620857, 0.24696243326222161
refl:, 0.7169705454388076, 0.4962025316455696, 20.840189591108683, 0.24654821236606225
refl:, 0.7169705454388076, 0.49, 20.56780221092012, 0.246087935342834
refl:, 0.7169705454388076, 0.4839506172839506, 20.302607744753697, 0.2455807741452668
refl:, 0.7169705454388076, 0.47804878048780486, 20.044318389931625, 0.24502832199568722
refl:, 0.7169705454388076, 0.47228915662650606, 19.79266184720379, 0.24443328421700983
refl:, 0.7169705454388076, 0.4666666666666667, 19.54738026749198, 0.24379556106465222
refl:, 0.7169705454388076, 0.4611764705882353, 19.308229285168313, 0.24311440973651832
refl:, 0.7169705454388076, 0.4558139534883721, 19.07497712952083, 0.2423917032614313
refl:, 0.7169705454388076, 0.4505747126436782, 18.847403806983923, 0.24162921270758458
refl:, 0.7169705454388076, 0.4454545454545454, 18.62530034751981, 0.24082577322829152
refl:, 0.7169705454388076, 0.44044943820224725, 18.408468109247096, 0.23998068018388435
refl:, 0.7169705454388076, 0.43555555555555553, 18.196718136035763, 0.23909605307915388
refl:, 0.7169705454388076, 0.4307692307692308, 17.98987056333767, 0.23817255818097147
refl:, 0.7169705454388076, 0.4260869565217391, 17.787754068005896, 0.23720797538105634
refl:, 0.7169705454388076, 0.421505376344086, 17.590205358285754, 0.236202366418755
refl:, 0.7169705454388076, 0.41702127659574467, 17.39706870053973, 0.23515820136339463
refl:, 0.7169705454388076, 0.4126315789473684, 17.20819547960629, 0.2340750107196232
refl:, 0.7169705454388076, 0.4083333333333333, 17.02344378999232, 0.2329510888033509
refl:, 0.7169705454388076, 0.4041237113402062, 16.842678055366356, 0.23178864386805378
refl:, 0.7169705454388076, 0.4, 16.665768674058118, 0.2305906700422077
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.20454876789491885 / 0.20454876789491885 = 1.0
field decay(t = 100.01): 7.178635175771892e-13 / 0.20454876789491885 = 3.509498125874663e-12
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.20454876813558592 / 0.20454876813558592 = 1.0
field decay(t = 100.01): 2.8174572450318613e-11 / 0.20454876813558592 = 1.3774012284270024e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 0.8034845121081741, 0.8, 39.99999999999999, 0.20179462162966672
refl:, 0.8034845121081741, 0.784, 39.04509528455022, 0.2059250719323383
refl:, 0.8034845121081741, 0.7686274509803922, 38.139646206191365, 0.20966397204966036
refl:, 0.8034845121081741, 0.7538461538461539, 37.27949726994566, 0.21305427397378518
refl:, 0.8034845121081741, 0.739622641509434, 36.46098976062103, 0.21613269723029083
refl:, 0.8034845121081741, 0.7259259259259259, 35.680884053707004, 0.21892741270212848
refl:, 0.8034845121081741, 0.7127272727272727, 34.93629688867566, 0.22146683922801672
refl:, 0.8034845121081741, 0.7, 34.22465020933728, 0.22377279980470224
refl:, 0.8034845121081741, 0.6877192982456141, 33.54362905484555, 0.22586624586940138
refl:, 0.8034845121081741, 0.6758620689655173, 32.89114661043291, 0.22776517254801126
refl:, 0.8034845121081741, 0.664406779661017, 32.265314978971055, 0.2294835019692795
refl:, 0.8034845121081741, 0.6533333333333333, 31.664420565727855, 0.2310373437439053
refl:, 0.8034845121081741, 0.6426229508196721, 31.086903214629274, 0.23243828166083527
refl:, 0.8034845121081741, 0.632258064516129, 30.531338419096144, 0.23369542692345546
refl:, 0.8034845121081741, 0.6222222222222222, 29.996422070856152, 0.2348212532271984
refl:, 0.8034845121081741, 0.6124999999999999, 29.480957317802616, 0.23582404434823864
refl:, 0.8034845121081741, 0.6030769230769231, 28.983843185365757, 0.23670970324155213
refl:, 0.8034845121081741, 0.593939393939394, 28.504064681021994, 0.2374874618810027
refl:, 0.8034845121081741, 0.5850746268656717, 28.0406841528977, 0.23816414677923553
refl:, 0.8034845121081741, 0.5764705882352942, 27.592833714170844, 0.2387438174110003
refl:, 0.8034845121081741, 0.5681159420289855, 27.159708577552564, 0.239232508084756
refl:, 0.8034845121081741, 0.56, 26.740561170354134, 0.23963612123153596
refl:, 0.8034845121081741, 0.552112676056338, 26.33469592188663, 0.23995828312103898
refl:, 0.8034845121081741, 0.5444444444444444, 25.94146463225064, 0.24020224568501714
refl:, 0.8034845121081741, 0.536986301369863, 25.560262345758233, 0.24037244343095185
refl:, 0.8034845121081741, 0.5297297297297298, 25.190523663916437, 0.2404730149043355
refl:, 0.8034845121081741, 0.5226666666666667, 24.831719442577725, 0.2405057199519815
refl:, 0.8034845121081741, 0.5157894736842105, 24.483353825914165, 0.24047273277761147
refl:, 0.8034845121081741, 0.509090909090909, 24.144961576600384, 0.2403784518025041
refl:, 0.8034845121081741, 0.5025641025641026, 23.816105667237974, 0.24022489861203003
refl:, 0.8034845121081741, 0.4962025316455696, 23.496375102814216, 0.24001212440412875
refl:, 0.8034845121081741, 0.49, 23.185382948015043, 0.2397434129524054
refl:, 0.8034845121081741, 0.4839506172839506, 22.882764536632674, 0.2394219863711472
refl:, 0.8034845121081741, 0.47804878048780486, 22.58817584322306, 0.23904730465890137
refl:, 0.8034845121081741, 0.47228915662650606, 22.301291999661412, 0.23862033792003395
refl:, 0.8034845121081741, 0.4666666666666667, 22.021805941382443, 0.23814466238694304
refl:, 0.8034845121081741, 0.4611764705882353, 21.749427169932794, 0.23762099874005177
refl:, 0.8034845121081741, 0.4558139534883721, 21.483880620051718, 0.23704845905896132
refl:, 0.8034845121081741, 0.4505747126436782, 21.224905620871482, 0.2364290006894749
refl:, 0.8034845121081741, 0.4454545454545454, 20.972254942022712, 0.23576489053070165
refl:, 0.8034845121081741, 0.44044943820224725, 20.725693916468796, 0.23505540516189377
refl:, 0.8034845121081741, 0.43555555555555553, 20.48499963280006, 0.23429986728741037
refl:, 0.8034845121081741, 0.4307692307692308, 20.249960190511292, 0.23350021386178307
refl:, 0.8034845121081741, 0.4260869565217391, 20.020374012481035, 0.23265749904183425
refl:, 0.8034845121081741, 0.421505376344086, 19.79604920948226, 0.23176986721335244
refl:, 0.8034845121081741, 0.41702127659574467, 19.576802992091316, 0.23083693098169072
refl:, 0.8034845121081741, 0.4126315789473684, 19.362461125837058, 0.22986107450928206
refl:, 0.8034845121081741, 0.4083333333333333, 19.15285742585141, 0.2288423958877941
refl:, 0.8034845121081741, 0.4041237113402062, 18.9478332876545, 0.22777913143869027
refl:, 0.8034845121081741, 0.4, 18.747237251037504, 0.22667336344348024
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.1943564114319847 / 0.1943564114319847 = 1.0
field decay(t = 100.01): 8.613369505850488e-13 / 0.1943564114319847 = 4.431739319731548e-12
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.1943564112167818 / 0.1943564112167818 = 1.0
field decay(t = 100.01): 2.929219089217826e-11 / 0.1943564112167818 = 1.5071378766870859e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 0.8838834764831843, 0.8, 44.99999999999999, 0.1749506546010766
refl:, 0.8838834764831843, 0.784, 43.865246854678944, 0.18084614262894314
refl:, 0.8838834764831843, 0.7686274509803922, 42.79498686880298, 0.18616490888473136
refl:, 0.8838834764831843, 0.7538461538461539, 41.78306957363623, 0.1909834896778329
refl:, 0.8838834764831843, 0.739622641509434, 40.82419653278361, 0.1953532231998348
refl:, 0.8838834764831843, 0.7259259259259259, 39.913765803587594, 0.19932045676242707
refl:, 0.8838834764831843, 0.7127272727272727, 39.047751315763335, 0.20293092176670602
refl:, 0.8838834764831843, 0.7, 38.2226079814105, 0.20621465898807328
refl:, 0.8838834764831843, 0.6877192982456141, 37.435196089128915, 0.20920558346668153
refl:, 0.8838834764831843, 0.6758620689655173, 36.682720369951866, 0.21192936258982145
refl:, 0.8838834764831843, 0.664406779661017, 35.962680378531154, 0.2144091185242741
refl:, 0.8838834764831843, 0.6533333333333333, 35.272829708778914, 0.21666666414632094
refl:, 0.8838834764831843, 0.6426229508196721, 34.61114218453038, 0.2187182335956619
refl:, 0.8838834764831843, 0.632258064516129, 33.975783613579615, 0.22058228827610923
refl:, 0.8838834764831843, 0.6222222222222222, 33.36508802077569, 0.22227211810191635
refl:, 0.8838834764831843, 0.6124999999999999, 32.77753751829758, 0.22379917822267162
refl:, 0.8838834764831843, 0.6030769230769231, 32.21174515294742, 0.225178014297656
refl:, 0.8838834764831843, 0.593939393939394, 31.666440208036935, 0.22641701701757366
refl:, 0.8838834764831843, 0.5850746268656717, 31.14045554291693, 0.22752429105750807
refl:, 0.8838834764831843, 0.5764705882352942, 30.6327166347442, 0.22851103000945724
refl:, 0.8838834764831843, 0.5681159420289855, 30.142232050688204, 0.2293833309292736
refl:, 0.8838834764831843, 0.56, 29.66808512880701, 0.23014665030198314
refl:, 0.8838834764831843, 0.552112676056338, 29.20942668547734, 0.2308091072410071
refl:, 0.8838834764831843, 0.5444444444444444, 28.765468598924116, 0.23137604508998047
refl:, 0.8838834764831843, 0.536986301369863, 28.33547814384704, 0.23185120425684555
refl:, 0.8838834764831843, 0.5297297297297298, 27.91877297273401, 0.23223967284122596
refl:, 0.8838834764831843, 0.5226666666666667, 27.514716656212954, 0.23254649427327514
refl:, 0.8838834764831843, 0.5157894736842105, 27.12271470851516, 0.23277497605931805
refl:, 0.8838834764831843, 0.509090909090909, 26.742211035417114, 0.2329274877105852
refl:, 0.8838834764831843, 0.5025641025641026, 26.372684751370738, 0.23300817835325607
refl:, 0.8838834764831843, 0.4962025316455696, 26.013647320299032, 0.23302096475520206
refl:, 0.8838834764831843, 0.49, 25.66463998102006, 0.23296661309174865
refl:, 0.8838834764831843, 0.4839506172839506, 25.32523142370283, 0.2328473819749458
refl:, 0.8838834764831843, 0.47804878048780486, 24.99501568834094, 0.23266777547883646
refl:, 0.8838834764831843, 0.47228915662650606, 24.673610260104642, 0.2324286704369094
refl:, 0.8838834764831843, 0.4666666666666667, 24.360654339720863, 0.23213016162250638
refl:, 0.8838834764831843, 0.4611764705882353, 24.055807269832574, 0.23177609164136181
refl:, 0.8838834764831843, 0.4558139534883721, 23.75874710068368, 0.2313686559958002
refl:, 0.8838834764831843, 0.4505747126436782, 23.46916928052956, 0.23090689658521463
refl:, 0.8838834764831843, 0.4454545454545454, 23.18678545794031, 0.23039248392782377
refl:, 0.8838834764831843, 0.44044943820224725, 22.911322384688706, 0.22982847953509086
refl:, 0.8838834764831843, 0.43555555555555553, 22.64252090923418, 0.22921468760037644
refl:, 0.8838834764831843, 0.4307692307692308, 22.380135051959574, 0.2285503145110825
refl:, 0.8838834764831843, 0.4260869565217391, 22.123931154313652, 0.22783728199610298
refl:, 0.8838834764831843, 0.421505376344086, 21.873687094881706, 0.22707701516938034
refl:, 0.8838834764831843, 0.41702127659574467, 21.629191566166806, 0.226267933168201
refl:, 0.8838834764831843, 0.4126315789473684, 21.390243406530644, 0.22540929123705422
refl:, 0.8838834764831843, 0.4083333333333333, 21.156650982328358, 0.22450324859888127
refl:, 0.8838834764831843, 0.4041237113402062, 20.928231615787418, 0.2235503997974889
refl:, 0.8838834764831843, 0.4, 20.704811054635428, 0.22254911459777174
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.18417632037527257 / 0.18417632037527257 = 1.0
field decay(t = 100.01): 1.0178427595748204e-12 / 0.18417632037527257 = 5.5264583280895835e-12
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.18417631975521828 / 0.18417631975521828 = 1.0
field decay(t = 100.01): 2.8713012611847576e-11 / 0.18417631975521828 = 1.5589958931750262e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 0.9575555538987225, 0.8, 50.0, 0.14429067301668
refl:, 0.9575555538987225, 0.784, 48.6530933185261, 0.15241249489296654
refl:, 0.9575555538987225, 0.7686274509803922, 47.39207738436693, 0.15970237410745958
refl:, 0.9575555538987225, 0.7538461538461539, 46.207396358038395, 0.16628917742963503
refl:, 0.9575555538987225, 0.739622641509434, 45.091066324150106, 0.1722363131532337
refl:, 0.9575555538987225, 0.7259259259259259, 44.036334278275795, 0.1776288135399773
refl:, 0.9575555538987225, 0.7127272727272727, 43.037427456538545, 0.1825236470768417
refl:, 0.9575555538987225, 0.7, 42.089365210987516, 0.1869718164862598
refl:, 0.9575555538987225, 0.6877192982456141, 41.18781524255363, 0.19102552097133876
refl:, 0.9575555538987225, 0.6758620689655173, 40.32898196453324, 0.19471651743120957
refl:, 0.9575555538987225, 0.664406779661017, 39.50951857909891, 0.1980846376246996
refl:, 0.9575555538987225, 0.6533333333333333, 38.726456948394606, 0.20115681830261414
refl:, 0.9575555538987225, 0.6426229508196721, 37.97715101972185, 0.20395988101926588
refl:, 0.9575555538987225, 0.632258064516129, 37.25923071465613, 0.20651807398114963
refl:, 0.9575555538987225, 0.6222222222222222, 36.57056399547784, 0.20884957043569513
refl:, 0.9575555538987225, 0.6124999999999999, 35.909225393217994, 0.21097550520798009
refl:, 0.9575555538987225, 0.6030769230769231, 35.273469693564564, 0.21290924565401623
refl:, 0.9575555538987225, 0.593939393939394, 34.6617097783444, 0.214665416754943
refl:, 0.9575555538987225, 0.5850746268656717, 34.07249784378941, 0.2162593463941121
refl:, 0.9575555538987225, 0.5764705882352942, 33.504509384470694, 0.21769885837479563
refl:, 0.9575555538987225, 0.5681159420289855, 32.95652945897567, 0.2189956192414337
refl:, 0.9575555538987225, 0.56, 32.42744085087325, 0.22016058268759847
refl:, 0.9575555538987225, 0.552112676056338, 31.91621381391884, 0.22119898455404
refl:, 0.9575555538987225, 0.5444444444444444, 31.42189714930561, 0.22211935917204448
refl:, 0.9575555538987225, 0.536986301369863, 30.943610409083956, 0.2229298821236265
refl:, 0.9575555538987225, 0.5297297297297298, 30.480537056602333, 0.22363493657078803
refl:, 0.9575555538987225, 0.5226666666666667, 30.031918444163345, 0.2242400133418431
refl:, 0.9575555538987225, 0.5157894736842105, 29.597048491686717, 0.22475115535497284
refl:, 0.9575555538987225, 0.509090909090909, 29.17526896927061, 0.2251729335780506
refl:, 0.9575555538987225, 0.5025641025641026, 28.76596530209651, 0.22550858491249842
refl:, 0.9575555538987225, 0.4962025316455696, 28.36856282886028, 0.22576196217784258
refl:, 0.9575555538987225, 0.49, 27.982523455399697, 0.22593807040959396
refl:, 0.9575555538987225, 0.4839506172839506, 27.60734265386835, 0.22603917633476883
refl:, 0.9575555538987225, 0.47804878048780486, 27.24254676502359, 0.22606684493527063
refl:, 0.9575555538987225, 0.47228915662650606, 26.88769056722596, 0.22602594884033364
refl:, 0.9575555538987225, 0.4666666666666667, 26.54235508080712, 0.22591913511748662
refl:, 0.9575555538987225, 0.4611764705882353, 26.206145580726076, 0.2257461894573964
refl:, 0.9575555538987225, 0.4558139534883721, 25.878689794039474, 0.2255106337747593
refl:, 0.9575555538987225, 0.4505747126436782, 25.55963626177365, 0.22521611643853037
refl:, 0.9575555538987225, 0.4454545454545454, 25.248652847395146, 0.2248621048576318
refl:, 0.9575555538987225, 0.44044943820224725, 24.945425376307522, 0.2244497571914091
refl:, 0.9575555538987225, 0.43555555555555553, 24.64965639271632, 0.22398274531115286
refl:, 0.9575555538987225, 0.4307692307692308, 24.361064021851533, 0.22346171326862457
refl:, 0.9575555538987225, 0.4260869565217391, 24.079380926958596, 0.22288585592315466
refl:, 0.9575555538987225, 0.421505376344086, 23.804353351700417, 0.22225690920257898
refl:, 0.9575555538987225, 0.41702127659574467, 23.535740239681235, 0.22157652716566317
refl:, 0.9575555538987225, 0.4126315789473684, 23.27331242373393, 0.22084351574148542
refl:, 0.9575555538987225, 0.4083333333333333, 23.016851878423925, 0.220057130518136
refl:, 0.9575555538987225, 0.4041237113402062, 22.76615102993319, 0.21921919579654955
refl:, 0.9575555538987225, 0.4, 22.521012118111, 0.21832995493951468
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.17431196343224772 / 0.17431196343224772 = 1.0
field decay(t = 100.01): 1.226599619152291e-12 / 0.17431196343224772 = 7.0368068547919865e-12
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.17431196264184795 / 0.17431196264184795 = 1.0
field decay(t = 100.01): 2.9177363356244976e-11 / 0.17431196264184795 = 1.6738589201817764e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 1.0239400553612397, 0.8, 54.99999999999999, 0.11015750327021501
refl:, 1.0239400553612397, 0.784, 53.39534223243568, 0.12100519786278613
refl:, 1.0239400553612397, 0.7686274509803922, 51.90867481651075, 0.13067356321855755
refl:, 1.0239400553612397, 0.7538461538461539, 50.524208813748004, 0.1393743175816963
refl:, 1.0239400553612397, 0.739622641509434, 49.22931269112476, 0.14718461760941978
refl:, 1.0239400553612397, 0.7259259259259259, 48.013685667151286, 0.15425442212255983
refl:, 1.0239400553612397, 0.7127272727272727, 46.86879216202182, 0.16064116852684127
refl:, 1.0239400553612397, 0.7, 45.787462696319906, 0.1664466525245801
refl:, 1.0239400553612397, 0.6877192982456141, 44.7636047370949, 0.17171830451504444
refl:, 1.0239400553612397, 0.6758620689655173, 43.79198840588911, 0.17651889499651638
refl:, 1.0239400553612397, 0.664406779661017, 42.868084512177944, 0.1808955749197593
refl:, 1.0239400553612397, 0.6533333333333333, 41.98794000564178, 0.18488590864155194
refl:, 1.0239400553612397, 0.6426229508196721, 41.148080731197126, 0.1885334768453189
refl:, 1.0239400553612397, 0.632258064516129, 40.345434464399176, 0.19186399028554912
refl:, 1.0239400553612397, 0.6222222222222222, 39.57726925315509, 0.1949105801621511
refl:, 1.0239400553612397, 0.6124999999999999, 38.841143478370526, 0.19769436000574783
refl:, 1.0239400553612397, 0.6030769230769231, 38.13486500383326, 0.20023789123980396
refl:, 1.0239400553612397, 0.593939393939394, 37.45645745905701, 0.2025624060128163
refl:, 1.0239400553612397, 0.5850746268656717, 36.80413218014341, 0.204680856177325
refl:, 1.0239400553612397, 0.5764705882352942, 36.176264682924305, 0.2066128346051995
refl:, 1.0239400553612397, 0.5681159420289855, 35.57137479946729, 0.20837010840175924
refl:, 1.0239400553612397, 0.56, 34.98810980028866, 0.2099625172331732
refl:, 1.0239400553612397, 0.552112676056338, 34.42522996870841, 0.21140487263605512
refl:, 1.0239400553612397, 0.5444444444444444, 33.881596203503065, 0.2127045522560479
refl:, 1.0239400553612397, 0.536986301369863, 33.35615931039313, 0.21386953098052403
refl:, 1.0239400553612397, 0.5297297297297298, 32.847950708397136, 0.21491013773280657
refl:, 1.0239400553612397, 0.5226666666666667, 32.3560743283589, 0.21583212552616557
refl:, 1.0239400553612397, 0.5157894736842105, 31.87969952142079, 0.2166414381263308
refl:, 1.0239400553612397, 0.509090909090909, 31.41805482739263, 0.21734485598596862
refl:, 1.0239400553612397, 0.5025641025641026, 30.97042247873359, 0.2179480717640861
refl:, 1.0239400553612397, 0.4962025316455696, 30.536133536635912, 0.21845549976557374
refl:, 1.0239400553612397, 0.49, 30.114563572550164, 0.21887075023784372
refl:, 1.0239400553612397, 0.4839506172839506, 29.70512882224086, 0.21919941164026552
refl:, 1.0239400553612397, 0.47804878048780486, 29.307282750744406, 0.21944546389616373
refl:, 1.0239400553612397, 0.47228915662650606, 28.920512975908146, 0.21961022544356806
refl:, 1.0239400553612397, 0.4666666666666667, 28.544338505905493, 0.2196984674330614
refl:, 1.0239400553612397, 0.4611764705882353, 28.178307252549438, 0.21971428386147993
refl:, 1.0239400553612397, 0.4558139534883721, 27.821993787605123, 0.21965764862748693
refl:, 1.0239400553612397, 0.4505747126436782, 27.474997313821564, 0.21953193965864698
refl:, 1.0239400553612397, 0.4454545454545454, 27.13693982621647, 0.21934178652625613
refl:, 1.0239400553612397, 0.44044943820224725, 26.80746444237855, 0.21908700093353123
refl:, 1.0239400553612397, 0.43555555555555553, 26.486233883298294, 0.21876853077786598
refl:, 1.0239400553612397, 0.4307692307692308, 26.172929088582443, 0.2183905958409472
refl:, 1.0239400553612397, 0.4260869565217391, 25.867247951913267, 0.21795488284012238
refl:, 1.0239400553612397, 0.421505376344086, 25.56890416433822, 0.21746063743615796
refl:, 1.0239400553612397, 0.41702127659574467, 25.277626154460066, 0.21690850348145124
refl:, 1.0239400553612397, 0.4126315789473684, 24.993156115881433, 0.21630059246127992
refl:, 1.0239400553612397, 0.4083333333333333, 24.715249113370163, 0.21563737630418256
refl:, 1.0239400553612397, 0.4041237113402062, 24.443672260178424, 0.2149184652144797
refl:, 1.0239400553612397, 0.4, 24.178203959791162, 0.2141434469115452
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.16505586379118362 / 0.16505586379118362 = 1.0
field decay(t = 100.01): 1.6796657955389861e-12 / 0.16505586379118362 = 1.0176347310289891e-11
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.16505586311279807 / 0.16505586311279807 = 1.0
field decay(t = 100.01): 3.3223715038456985e-11 / 0.16505586311279807 = 2.0128769988468766e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 1.0825317547305482, 0.8, 59.99999999999999, 0.07360773441443474
refl:, 1.0825317547305482, 0.784, 58.07108488083593, 0.08756183349395384
refl:, 1.0825317547305482, 0.7686274509803922, 56.311309260669624, 0.09993186958234106
refl:, 1.0825317547305482, 0.7538461538461539, 54.69254540505164, 0.11102792920544496
refl:, 1.0825317547305482, 0.739622641509434, 53.19365179333026, 0.12095671757064261
refl:, 1.0825317547305482, 0.7259259259259259, 51.798237196004614, 0.12989111490169047
refl:, 1.0825317547305482, 0.7127272727272727, 50.49327371888257, 0.13795385843502186
refl:, 1.0825317547305482, 0.7, 49.268194122622916, 0.14525656081749277
refl:, 1.0825317547305482, 0.6877192982456141, 48.11428082873218, 0.15187584382721003
refl:, 1.0825317547305482, 0.6758620689655173, 47.024238635081495, 0.15790192109849083
refl:, 1.0825317547305482, 0.664406779661017, 45.99188753936469, 0.16337692994168027
refl:, 1.0825317547305482, 0.6533333333333333, 45.01193662107726, 0.16837576220748254
refl:, 1.0825317547305482, 0.6426229508196721, 44.07981414374115, 0.1729308563822558
refl:, 1.0825317547305482, 0.632258064516129, 43.191537588393835, 0.17709787538099556
refl:, 1.0825317547305482, 0.6222222222222222, 42.34361264736431, 0.18090702518494378
refl:, 1.0825317547305482, 0.6124999999999999, 41.53295361356457, 0.18439374819992094
refl:, 1.0825317547305482, 0.6030769230769231, 40.756819839298586, 0.1875864032458676
refl:, 1.0825317547305482, 0.593939393939394, 40.01276444425686, 0.1905057507798468
refl:, 1.0825317547305482, 0.5850746268656717, 39.298592486049785, 0.19318030905511865
refl:, 1.0825317547305482, 0.5764705882352942, 38.61232652957572, 0.1956239605857035
refl:, 1.0825317547305482, 0.5681159420289855, 37.9521780657052, 0.19785784772625742
refl:, 1.0825317547305482, 0.56, 37.31652360111115, 0.19989839276598212
refl:, 1.0825317547305482, 0.552112676056338, 36.70388451304266, 0.2017547526768604
refl:, 1.0825317547305482, 0.5444444444444444, 36.112909964600526, 0.20344360330889855
refl:, 1.0825317547305482, 0.536986301369863, 35.54236232752769, 0.20497480625666326
refl:, 1.0825317547305482, 0.5297297297297298, 34.991104674473426, 0.20635699484209216
refl:, 1.0825317547305482, 0.5226666666666667, 34.4580899908132, 0.20760151265764157
refl:, 1.0825317547305482, 0.5157894736842105, 33.94235182430718, 0.20871528450622134
refl:, 1.0825317547305482, 0.509090909090909, 33.442996144127704, 0.20970524738496094
refl:, 1.0825317547305482, 0.5025641025641026, 32.95919422270091, 0.2105783514069682
refl:, 1.0825317547305482, 0.4962025316455696, 32.4901763870572, 0.2113414215404359
refl:, 1.0825317547305482, 0.49, 32.035226512949734, 0.21200027362523924
refl:, 1.0825317547305482, 0.4839506172839506, 31.593677156370052, 0.21255807961996584
refl:, 1.0825317547305482, 0.47804878048780486, 31.164905234389202, 0.21302071809138926
refl:, 1.0825317547305482, 0.47228915662650606, 30.748328181342448, 0.21339356617448996
refl:, 1.0825317547305482, 0.4666666666666667, 30.343400517915825, 0.2136782430141135
refl:, 1.0825317547305482, 0.4611764705882353, 29.949610780196405, 0.21387969509149066
refl:, 1.0825317547305482, 0.4558139534883721, 29.566478763613983, 0.2140025648813409
refl:, 1.0825317547305482, 0.4505747126436782, 29.193553043244272, 0.21404664550288702
refl:, 1.0825317547305482, 0.4454545454545454, 28.830408737409662, 0.2140158741099825
refl:, 1.0825317547305482, 0.44044943820224725, 28.47664548610104, 0.2139163515040305
refl:, 1.0825317547305482, 0.43555555555555553, 28.131885619609587, 0.213747755347228
refl:, 1.0825317547305482, 0.4307692307692308, 27.79577249602797, 0.2135103221579689
refl:, 1.0825317547305482, 0.4260869565217391, 27.46796898905743, 0.21320790783224955
refl:, 1.0825317547305482, 0.421505376344086, 27.14815610992485, 0.2128440342182022
refl:, 1.0825317547305482, 0.41702127659574467, 26.836031749237822, 0.21242171734687695
refl:, 1.0825317547305482, 0.4126315789473684, 26.531309526343414, 0.21193780892543815
refl:, 1.0825317547305482, 0.4083333333333333, 26.23371773525153, 0.21138950860192895
refl:, 1.0825317547305482, 0.4041237113402062, 25.94299837747481, 0.21078013578978685
refl:, 1.0825317547305482, 0.4, 25.65890627325528, 0.2101171783272214
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.15668265201002152 / 0.15668265201002152 = 1.0
field decay(t = 100.01): 2.2446823295144107e-12 / 0.15668265201002152 = 1.4326297779098349e-11
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.1566826516265948 / 0.1566826516265948 = 1.0
field decay(t = 100.01): 3.452976998428892e-11 / 0.1566826516265948 = 2.2038030136597423e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 1.1328847337958123, 0.8, 64.99999999999999, 0.037170886043112084
refl:, 1.1328847337958123, 0.784, 62.645633431180784, 0.05362783985573467
refl:, 1.1328847337958123, 0.7686274509803922, 60.547811567584965, 0.06875709351287543
refl:, 1.1328847337958123, 0.7538461538461539, 58.65172469781183, 0.082354616478198
refl:, 1.1328847337958123, 0.739622641509434, 56.919780417243786, 0.09462755396270182
refl:, 1.1328847337958123, 0.7259259259259259, 55.32480276439999, 0.10557645638030469
refl:, 1.1328847337958123, 0.7127272727272727, 53.84635215217889, 0.11544950273548941
refl:, 1.1328847337958123, 0.7, 52.468573108244414, 0.1243168343725702
refl:, 1.1328847337958123, 0.6877192982456141, 51.17885777022349, 0.13235703881433808
refl:, 1.1328847337958123, 0.6758620689655173, 49.96697533403755, 0.139626360552005
refl:, 1.1328847337958123, 0.664406779661017, 48.82448244349733, 0.14625471025031767
refl:, 1.1328847337958123, 0.6533333333333333, 47.744310665402864, 0.1522788790416092
refl:, 1.1328847337958123, 0.6426229508196721, 46.720469815981424, 0.1577852862490647
refl:, 1.1328847337958123, 0.632258064516129, 45.7478295128327, 0.16280555175366354
refl:, 1.1328847337958123, 0.6222222222222222, 44.82195500406973, 0.167398621072885
refl:, 1.1328847337958123, 0.6124999999999999, 43.93898156096055, 0.17159777718924688
refl:, 1.1328847337958123, 0.6030769230769231, 43.09551684594588, 0.17544391308281032
refl:, 1.1328847337958123, 0.593939393939394, 42.288563952640665, 0.1789701508099767
refl:, 1.1328847337958123, 0.5850746268656717, 41.51545997434592, 0.18219866439447
refl:, 1.1328847337958123, 0.5764705882352942, 40.77382641078957, 0.185160771243366
refl:, 1.1328847337958123, 0.5681159420289855, 40.06152872079658, 0.18786997681737858
refl:, 1.1328847337958123, 0.56, 39.37664302675111, 0.19035022781311492
refl:, 1.1328847337958123, 0.552112676056338, 38.717428473384835, 0.19261993394590182
refl:, 1.1328847337958123, 0.5444444444444444, 38.08230410219336, 0.19469110743273724
refl:, 1.1328847337958123, 0.536986301369863, 37.4698293655767, 0.1965812755931606
refl:, 1.1328847337958123, 0.5297297297297298, 36.87868759977799, 0.19830004161335185
refl:, 1.1328847337958123, 0.5226666666666667, 36.307671922074825, 0.19985824848269654
refl:, 1.1328847337958123, 0.5157894736842105, 35.75567312877807, 0.20126711563403166
refl:, 1.1328847337958123, 0.509090909090909, 35.22166925577522, 0.20253568251243403
refl:, 1.1328847337958123, 0.5025641025641026, 34.70471652928195, 0.2036726801968963
refl:, 1.1328847337958123, 0.4962025316455696, 34.203941485939055, 0.20468320295652365
refl:, 1.1328847337958123, 0.49, 33.7185340819115, 0.20557517176695547
refl:, 1.1328847337958123, 0.4839506172839506, 33.247741642788455, 0.20635542059630924
refl:, 1.1328847337958123, 0.47804878048780486, 32.790863531764046, 0.2070277254966668
refl:, 1.1328847337958123, 0.47228915662650606, 32.347246434236666, 0.20759964094140226
refl:, 1.1328847337958123, 0.4666666666666667, 31.91628017368787, 0.20807496149136276
refl:, 1.1328847337958123, 0.4611764705882353, 31.497393987322308, 0.20845469225613245
refl:, 1.1328847337958123, 0.4558139534883721, 31.090053201105746, 0.20874562394291427
refl:, 1.1328847337958123, 0.4505747126436782, 30.693756253024553, 0.20895386380475825
refl:, 1.1328847337958123, 0.4454545454545454, 30.308032020994055, 0.20908023176898172
refl:, 1.1328847337958123, 0.44044943820224725, 29.9324374181665, 0.209125738596599
refl:, 1.1328847337958123, 0.43555555555555553, 29.5665552236732, 0.20909456295065154
refl:, 1.1328847337958123, 0.4307692307692308, 29.209992121269647, 0.20899049887482632
refl:, 1.1328847337958123, 0.4260869565217391, 28.86237692208816, 0.2088190433084917
refl:, 1.1328847337958123, 0.421505376344086, 28.523358950864424, 0.2085833129971406
refl:, 1.1328847337958123, 0.41702127659574467, 28.192606577687982, 0.2082742041204916
refl:, 1.1328847337958123, 0.4126315789473684, 27.86980587961545, 0.2078930131964556
refl:, 1.1328847337958123, 0.4083333333333333, 27.554659418441407, 0.20745662776433935
refl:, 1.1328847337958123, 0.4041237113402062, 27.246885122601544, 0.20697561960434535
refl:, 1.1328847337958123, 0.4, 26.946215262627685, 0.20643519600186538
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.14944315649486126 / 0.14944315649486126 = 1.0
field decay(t = 100.01): 4.987937476864823e-12 / 0.14944315649486126 = 3.3376820952228335e-11
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.14944315644072245 / 0.14944315644072245 = 1.0
field decay(t = 100.01): 3.807851695141726e-11 / 0.14944315644072245 = 2.548026812222836e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 1.1746157759823854, 0.8, 70.0, 0.008715934286147569
refl:, 1.1746157759823854, 0.784, 67.05783140972835, 0.024210799209089452
refl:, 1.1746157759823854, 0.7686274509803922, 64.53417775149677, 0.04029873111181387
refl:, 1.1746157759823854, 0.7538461538461539, 62.31059707934406, 0.0555678814867085
refl:, 1.1746157759823854, 0.739622641509434, 60.31629899672389, 0.06970125775811391
refl:, 1.1746157759823854, 0.7259259259259259, 58.50481213939229, 0.08255524971343331
refl:, 1.1746157759823854, 0.7127272727272727, 56.843594941553754, 0.09413263635827869
refl:, 1.1746157759823854, 0.7, 55.30875739675561, 0.10463698983404623
refl:, 1.1746157759823854, 0.6877192982456141, 53.88211612099151, 0.11407493400247079
refl:, 1.1746157759823854, 0.6758620689655173, 52.54943143507829, 0.1226472707657511
refl:, 1.1746157759823854, 0.664406779661017, 51.29929235382265, 0.13039930819744874
refl:, 1.1746157759823854, 0.6533333333333333, 50.122379445158494, 0.13744242808130103
refl:, 1.1746157759823854, 0.6426229508196721, 49.010959419193355, 0.1438468426182844
refl:, 1.1746157759823854, 0.632258064516129, 47.958527855818005, 0.14968921112177416
refl:, 1.1746157759823854, 0.6222222222222222, 46.9595500055931, 0.1550186900782557
refl:, 1.1746157759823854, 0.6124999999999999, 46.00926848902442, 0.15990623687364083
refl:, 1.1746157759823854, 0.6030769230769231, 45.103557823882895, 0.1643748532843572
refl:, 1.1746157759823854, 0.593939393939394, 44.23881248012816, 0.16847752780425043
refl:, 1.1746157759823854, 0.5850746268656717, 43.41185942214757, 0.1722356163059057
refl:, 1.1746157759823854, 0.5764705882352942, 42.619888854495166, 0.17567897390885354
refl:, 1.1746157759823854, 0.5681159420289855, 41.860398715356425, 0.17883996497200733
refl:, 1.1746157759823854, 0.56, 41.131149701242414, 0.18173480153232927
refl:, 1.1746157759823854, 0.552112676056338, 40.430128463310105, 0.18439121641848574
refl:, 1.1746157759823854, 0.5444444444444444, 39.75551721886283, 0.18682378756973045
refl:, 1.1746157759823854, 0.536986301369863, 39.10566845306894, 0.18904674104030247
refl:, 1.1746157759823854, 0.5297297297297298, 38.479083699202135, 0.19107686788902303
refl:, 1.1746157759823854, 0.5226666666666667, 37.87439561623927, 0.19292625454975332
refl:, 1.1746157759823854, 0.5157894736842105, 37.29035275442276, 0.19460884663783543
refl:, 1.1746157759823854, 0.509090909090909, 36.72580652885898, 0.1961345571897698
refl:, 1.1746157759823854, 0.5025641025641026, 36.179700019842116, 0.1975122933988443
refl:, 1.1746157759823854, 0.4962025316455696, 35.651058294458764, 0.19875051751344916
refl:, 1.1746157759823854, 0.49, 35.138980002927255, 0.19985541870413245
refl:, 1.1746157759823854, 0.4839506172839506, 34.64263004924646, 0.20083789883449643
refl:, 1.1746157759823854, 0.47804878048780486, 34.161233172132036, 0.2017049645889793
refl:, 1.1746157759823854, 0.47228915662650606, 33.69406830116886, 0.20245751259203246
refl:, 1.1746157759823854, 0.4666666666666667, 33.24046357629389, 0.20310354230391978
refl:, 1.1746157759823854, 0.4611764705882353, 32.79979193741594, 0.2036476720462249
refl:, 1.1746157759823854, 0.4558139534883721, 32.37146720614281, 0.20409307625590972
refl:, 1.1746157759823854, 0.4505747126436782, 31.954940593960732, 0.20445096223605416
refl:, 1.1746157759823854, 0.4454545454545454, 31.549697581364857, 0.20471845410715214
refl:, 1.1746157759823854, 0.44044943820224725, 31.155255120816925, 0.20489444063169535
refl:, 1.1746157759823854, 0.43555555555555553, 30.771159123350074, 0.2049889592193817
refl:, 1.1746157759823854, 0.4307692307692308, 30.396982194426784, 0.20500542736591415
refl:, 1.1746157759823854, 0.4260869565217391, 30.032321589495837, 0.20495768982428397
refl:, 1.1746157759823854, 0.421505376344086, 29.67679736376329, 0.20483056853588974
refl:, 1.1746157759823854, 0.41702127659574467, 29.33005069412437, 0.20461216911883964
refl:, 1.1746157759823854, 0.4126315789473684, 28.991742354111988, 0.2043374679527374
refl:, 1.1746157759823854, 0.4083333333333333, 28.66155132518973, 0.20400220026620766
refl:, 1.1746157759823854, 0.4041237113402062, 28.339173529827534, 0.2036399273981529
refl:, 1.1746157759823854, 0.4, 28.024320673604695, 0.20318419566870832
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.14355855111904026 / 0.14355855111904026 = 1.0
field decay(t = 100.01): 9.667874474337293e-12 / 0.14355855111904026 = 6.734446954901761e-11
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.1435585513271557 / 0.1435585513271557 = 1.0
field decay(t = 100.01): 5.780725980413138e-11 / 0.1435585513271557 = 4.0267374719040156e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 1.2074072828613354, 0.8, 75.0, 0.002279010064127535
refl:, 1.2074072828613354, 0.784, 71.19254305121969, 0.004052585814133056
refl:, 1.2074072828613354, 0.7686274509803922, 68.13230108184366, 0.017669454881364767
refl:, 1.2074072828613354, 0.7538461538461539, 65.5329128763873, 0.03336169435169541
refl:, 1.2074072828613354, 0.739622641509434, 63.25596068464814, 0.048654509225169895
refl:, 1.2074072828613354, 0.7259259259259259, 61.22160180379117, 0.06275216167717797
refl:, 1.2074072828613354, 0.7127272727272727, 59.378629248099124, 0.07576309243511313
refl:, 1.2074072828613354, 0.7, 57.691773022139905, 0.08758135787027703
refl:, 1.2074072828613354, 0.6877192982456141, 56.13545775731428, 0.098248929360899
refl:, 1.2074072828613354, 0.6758620689655173, 54.69040263962552, 0.10794664183105998
refl:, 1.2074072828613354, 0.664406779661017, 53.3416232015308, 0.11670121483117171
refl:, 1.2074072828613354, 0.6533333333333333, 52.0771858356392, 0.12465585169755833
refl:, 1.2074072828613354, 0.6426229508196721, 50.88739414232009, 0.13187389976516842
refl:, 1.2074072828613354, 0.632258064516129, 49.76423655717566, 0.13846172295217293
refl:, 1.2074072828613354, 0.6222222222222222, 48.70099913863975, 0.14446173011718028
refl:, 1.2074072828613354, 0.6124999999999999, 47.69198665804526, 0.14995965765204525
refl:, 1.2074072828613354, 0.6030769230769231, 46.73231696062462, 0.1549825384470223
refl:, 1.2074072828613354, 0.593939393939394, 45.817766249324485, 0.15958402411708034
refl:, 1.2074072828613354, 0.5850746268656717, 44.94465059915575, 0.16380217089082738
refl:, 1.2074072828613354, 0.5764705882352942, 44.109733785683, 0.1676654538139623
refl:, 1.2074072828613354, 0.5681159420289855, 43.31015457776125, 0.17121662693717593
refl:, 1.2074072828613354, 0.56, 42.543368664412704, 0.1744741445938872
refl:, 1.2074072828613354, 0.552112676056338, 41.80710174663229, 0.17746630440999525
refl:, 1.2074072828613354, 0.5444444444444444, 41.099311260706344, 0.1802095438402015
refl:, 1.2074072828613354, 0.536986301369863, 40.418154855029364, 0.18272173169827843
refl:, 1.2074072828613354, 0.5297297297297298, 39.76196420912173, 0.18502119081455226
refl:, 1.2074072828613354, 0.5226666666666667, 39.12922312098109, 0.1871223337651491
refl:, 1.2074072828613354, 0.5157894736842105, 38.51854903626071, 0.18904214560856375
refl:, 1.2074072828613354, 0.509090909090909, 37.928677376424204, 0.190789263829168
refl:, 1.2074072828613354, 0.5025641025641026, 37.35844816099257, 0.192373003931677
refl:, 1.2074072828613354, 0.4962025316455696, 36.806794523773725, 0.19380521484150348
refl:, 1.2074072828613354, 0.49, 36.272732803336254, 0.19509368153538312
refl:, 1.2074072828613354, 0.4839506172839506, 35.755353950220794, 0.19624859910356954
refl:, 1.2074072828613354, 0.47804878048780486, 35.253816041990945, 0.19727952544337504
refl:, 1.2074072828613354, 0.47228915662650606, 34.76733773550362, 0.19818701041137332
refl:, 1.2074072828613354, 0.4666666666666667, 34.295192516152845, 0.19897646054575882
refl:, 1.2074072828613354, 0.4611764705882353, 33.83670362811667, 0.19965889435875267
refl:, 1.2074072828613354, 0.4558139534883721, 33.391239589169516, 0.20023796589376308
refl:, 1.2074072828613354, 0.4505747126436782, 32.95821020943854, 0.20072106469864034
refl:, 1.2074072828613354, 0.4454545454545454, 32.537063046367535, 0.2011073384968122
refl:, 1.2074072828613354, 0.44044943820224725, 32.12728023870626, 0.2013932866552523
refl:, 1.2074072828613354, 0.43555555555555553, 31.728375671037586, 0.20159478768644698
refl:, 1.2074072828613354, 0.4307692307692308, 31.33989242755132, 0.20171671281423958
refl:, 1.2074072828613354, 0.4260869565217391, 30.96140049976025, 0.20176796051792995
refl:, 1.2074072828613354, 0.421505376344086, 30.592494717856905, 0.20173546160628877
refl:, 1.2074072828613354, 0.41702127659574467, 30.23279287960783, 0.2015919333831556
refl:, 1.2074072828613354, 0.4126315789473684, 29.88193405422147, 0.20140043913698885
refl:, 1.2074072828613354, 0.4083333333333333, 29.539577041619506, 0.20116952947891878
refl:, 1.2074072828613354, 0.4041237113402062, 29.2053989700849, 0.2008832766544487
refl:, 1.2074072828613354, 0.4, 28.879094017427605, 0.200502480384686
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.1392134275610875 / 0.1392134275610875 = 1.0
field decay(t = 100.01): 1.3144200424788948e-11 / 0.1392134275610875 = 9.441761944278838e-11
run 0 finished at t = 100.01 (10001 timesteps)
-----------
Initializing structure...
Meep: using complex fields.
field decay(t = 50.01): 0.13921342793500532 / 0.13921342793500532 = 1.0
field decay(t = 100.01): 9.871965144118156e-11 / 0.13921342793500532 = 7.091244925544889e-10
run 0 finished at t = 100.01 (10001 timesteps)
refl:, 1.23100969126526, 0.8, 79.99999999999994, 0.028667884842985224
refl:, 1.23100969126526, 0.784, 74.82079670263325, 0.0021347453734021486
refl:, 1.23100969126526, 0.7686274509803922, 71.11813585396203, 0.0041164349187047915
refl:, 1.23100969126526, 0.7538461538461539, 68.12392493650243, 0.017504855468406976
refl:, 1.23100969126526, 0.739622641509434, 65.57213417097068, 0.03282494395652445
refl:, 1.23100969126526, 0.7259259259259259, 63.331955818023665, 0.04792423816150501
refl:, 1.23100969126526, 0.7127272727272727, 61.32721657609733, 0.06193140028783135
refl:, 1.23100969126526, 0.7, 59.508764013939604, 0.07458578833476998
refl:, 1.23100969126526, 0.6877192982456141, 57.84259890532987, 0.08617897772571746
refl:, 1.23100969126526, 0.6758620689655173, 56.30398980710504, 0.0966374263474978
refl:, 1.23100969126526, 0.664406779661017, 54.87424701948212, 0.10612404031633497
refl:, 1.23100969126526, 0.6533333333333333, 53.53881754577338, 0.1147736073308378
refl:, 1.23100969126526, 0.6426229508196721, 52.2860932123638, 0.12261878935484406
refl:, 1.23100969126526, 0.632258064516129, 51.10662968356778, 0.1297955968814688
refl:, 1.23100969126526, 0.6222222222222222, 49.992614952243095, 0.1363423963503902
refl:, 1.23100969126526, 0.6124999999999999, 48.937495982303304, 0.14232561295216176
refl:, 1.23100969126526, 0.6030769230769231, 47.93570930843333, 0.1477973393651342
refl:, 1.23100969126526, 0.593939393939394, 46.98248211140486, 0.1527979150330299
refl:, 1.23100969126526, 0.5850746268656717, 46.07368235996428, 0.15737616810594365
refl:, 1.23100969126526, 0.5764705882352942, 45.205703916101704, 0.16157175886307418
refl:, 1.23100969126526, 0.5681159420289855, 44.37537706820699, 0.16542702710683707
refl:, 1.23100969126526, 0.56, 43.57989789517631, 0.16896483740231436
refl:, 1.23100969126526, 0.552112676056338, 42.81677180338089, 0.17222042834178233
refl:, 1.23100969126526, 0.5444444444444444, 42.083767886770914, 0.17520079945615433
refl:, 1.23100969126526, 0.536986301369863, 41.378881661291835, 0.1779352569426369
refl:, 1.23100969126526, 0.5297297297297298, 40.70030435652892, 0.18044129732794095
refl:, 1.23100969126526, 0.5226666666666667, 40.046397397863686, 0.182734107792854
refl:, 1.23100969126526, 0.5157894736842105, 39.41567103836258, 0.18483861166276766
refl:, 1.23100969126526, 0.509090909090909, 38.806766338771496, 0.18675601688212784
refl:, 1.23100969126526, 0.5025641025641026, 38.218439871701385, 0.1884996414001091
refl:, 1.23100969126526, 0.4962025316455696, 37.64955065969264, 0.19008080645459632
refl:, 1.23100969126526, 0.49, 37.099048958377544, 0.1915079168724047
refl:, 1.23100969126526, 0.4839506172839506, 36.56596657389412, 0.19279885897563126
refl:, 1.23100969126526, 0.47804878048780486, 36.049408464086156, 0.19395530230339666
refl:, 1.23100969126526, 0.47228915662650606, 35.54854542021205, 0.19498186511670068
refl:, 1.23100969126526, 0.4666666666666667, 35.062607663064846, 0.19588322966927102
refl:, 1.23100969126526, 0.4611764705882353, 34.59087921692149, 0.19666472127374038
refl:, 1.23100969126526, 0.4558139534883721, 34.13269294834013, 0.19734835806450404
refl:, 1.23100969126526, 0.4505747126436782, 33.6874261758209, 0.197930496667696
refl:, 1.23100969126526, 0.4454545454545454, 33.25449677173539, 0.19840503576945626
refl:, 1.23100969126526, 0.44044943820224725, 32.83335969046958, 0.19877805012361263
refl:, 1.23100969126526, 0.43555555555555553, 32.42350386700291, 0.19904737688450852
refl:, 1.23100969126526, 0.4307692307692308, 32.0244494386133, 0.19925404758581897
refl:, 1.23100969126526, 0.4260869565217391, 31.63574524940934, 0.1993954429097289
refl:, 1.23100969126526, 0.421505376344086, 31.2569666032255, 0.19942422718572894
refl:, 1.23100969126526, 0.41702127659574467, 30.887713235292672, 0.1993498431053417
refl:, 1.23100969126526, 0.4126315789473684, 30.527607477190394, 0.19917542826890433
refl:, 1.23100969126526, 0.4083333333333333, 30.176292593038497, 0.19903411145233238
refl:, 1.23100969126526, 0.4041237113402062, 29.83343126780691, 0.19890224385206828
refl:, 1.23100969126526, 0.4, 29.498704231103652, 0.1984945592268268
In [9]:
plt.figure(dpi=200)
plt.pcolormesh(kxs, wvls, Rmeep, cmap='hot', shading='gouraud', vmin=0, vmax=Rmeep.max())
plt.axis([kxs[0,0], kxs[0,-1], wvl_min, wvl_max])
plt.yticks([t for t in np.linspace(0.4,0.8,5)])
plt.xlabel(r"Bloch-periodic wavevector (k$_x$/2π)")
plt.ylabel("wavelength (μm)")
plt.title("reflectance (meep)")
cbar = plt.colorbar()
cbar.set_ticks([t for t in np.linspace(0,0.4,5)])
cbar.set_ticklabels(["{:.1f}".format(t) for t in np.linspace(0,0.4,5)])
In [10]:
plt.figure(dpi=200)
plt.pcolormesh(thetas, wvls, Rmeep, cmap='hot', shading='gouraud', vmin=0, vmax=Rmeep.max())
plt.axis([thetas.min(), thetas.max(), wvl_min, wvl_max])
plt.xticks([t for t in range(0,100,20)])
plt.yticks([t for t in np.linspace(0.4,0.8,5)])
plt.xlabel("angle of incident planewave (degrees)")
plt.ylabel("wavelength (μm)")
plt.title("reflectance (meep)")
cbar = plt.colorbar()
cbar.set_ticks([t for t in np.linspace(0,0.4,5)])
cbar.set_ticklabels(["{:.1f}".format(t) for t in np.linspace(0,0.4,5)])
In [11]:
n1=1
n2=3.5

# compute angle of refracted planewave in medium n2
# for incident planewave in medium n1 at angle theta_in
theta_out = lambda theta_in: math.asin(n1*math.sin(theta_in)/n2)

# compute Fresnel reflectance for P-polarization in medium n2
# for incident planewave in medium n1 at angle theta_in
Rfresnel = lambda theta_in: math.fabs((n1*math.cos(theta_out(theta_in))-n2*math.cos(theta_in))/(n1*math.cos(theta_out(theta_in))+n2*math.cos(theta_in)))**2

Ranalytic = np.empty((nfreq, theta_in.size))
for m in range(wvl.size):
    for n in range(theta_in.size):
        Ranalytic[m,n] = Rfresnel(math.radians(thetas[m,n]))

plt.figure(dpi=200)
plt.pcolormesh(thetas, wvls, Ranalytic, cmap='hot', shading='gouraud', vmin=0, vmax=Ranalytic.max())
plt.axis([thetas.min(), thetas.max(), wvl_min, wvl_max])
plt.xticks([t for t in range(0,100,20)])
plt.yticks([t for t in np.linspace(0.4,0.8,5)])
plt.xlabel("angle of incident planewave (degrees)")
plt.ylabel("wavelength (μm)")
plt.title("reflectance (analytic)")
cbar = plt.colorbar()
cbar.set_ticks([t for t in np.linspace(0,0.4,5)])
cbar.set_ticklabels(["{:.1f}".format(t) for t in np.linspace(0,0.4,5)])