made the wedge simulation work, maybe

Author: FilipDominec

example_metamaterial_wedge/batch.sh

@@ -0,0 +1,13 @@
+#!/bin/bash
+
+## example:
+## ./batch.sh radius=10u resolution=4u simtime=30p
+
+../scatter.py model=RodArray $@
+../effparam.py
+last_simulation_name=`cat last_simulation_name.dat`
+rm -f effparam.dat; ln -s effparam/${last_simulation_name}_effparam.dat effparam.dat
+
+../wedge_refraction.py $@
+last_simulation_name=`cat last_simulation_name.dat`
+../scripts_postpro/plot_FA.py $last_simulation_name/wedge_at_x0.000e+00_z*.h5

scripts_postpro/plot_FA.py

@@ -35,7 +35,7 @@
 h5file = h5py.File(filename, "r")
 print "Found datasets:", h5file.keys()
 time1 = time.time()
-data  = np.array(h5file['ex.r'])
+data  = np.array(h5file['ex.r']) * (1+0j)
 data += np.array(h5file['ex.i']) * 1j
 print "Loaded dataset with shape:", data.shape, 'in %04d s.' % (time.time()-time1)
 try:
@@ -54,6 +54,8 @@
 freq    = np.fft.fftshift(freq) #+ freq[len(freq)/2]
 Efy     = np.fft.fftshift(Efy)
 
+print("debug", len(t), (t[1]-t[0]), freq, Efy) 
+
 ## Fourier transform along the spatial axis
 kT      = np.fft.fftfreq(len(y), d=(y[1]-y[0]))        # calculate the frequency axis with proper spacing
 Ef      = np.fft.fft(Efy, axis=0) / len(y) * 2*np.pi     # calculate the FFT values
@@ -73,7 +75,7 @@ def fftshift2(arr): return np.vstack([arr[len(arr)/2:,:], arr[:len(arr)/2,:]])
     contours = plt.contourf(t, y, Et, cmap=matplotlib.cm.RdBu, extend='both') # levels=np.arange(0.,1,.01), 
 else:
     toplot = np.log10(np.abs(Ef))
-    contours = plt.contourf(freq, kT, toplot,  cmap=matplotlib.cm.gist_earth, levels=np.linspace(np.min(toplot)*.8+np.max(toplot)*.2,np.max(toplot),200) ,extend='both') #  
+    contours = plt.contourf(freq, kT, toplot,  cmap=matplotlib.cm.gist_earth, levels=np.linspace(np.min(toplot)*.5+np.max(toplot)*.5,np.max(toplot),200) ,extend='both') #  
     #contours = plt.contourf(freq, kT, np.abs(Ef), cmap=matplotlib.cm.gist_earth, extend='both') # levels=np.arange(0.,1,.01), 
     plt.plot([0, maxfreq], [0, 0], c='w',lw=.5)
     plt.plot([0, maxfreq], [0, maxfreq/c], c='w',lw=.5)

scripts_postpro/plot_TY.py

@@ -33,7 +33,7 @@
 h5file = h5py.File(filename, "r")
 print "Found datasets:", h5file.keys()
 time1 = time.time()
-data  = np.array(h5file['ex.r'])
+data  = np.array(h5file['ex.r']) * (1+0j)
 data += np.array(h5file['ex.i']) * 1j
 print "Loaded dataset with shape:", data.shape, 'in %04d s.' % (time.time()-time1)
 try:

wedge_refraction.py

@@ -34,7 +34,7 @@
 #import meep
 
 class Wedge_model(meep_utils.AbstractMeepModel): #{{{
-    """  Array of circular metallic wires along electric field
+    """  Array of circular wires parallel to the electric field
     FD 2013-07-13
     """
     def cell_centers(self):
@@ -96,7 +96,6 @@ def where_wire(self, r):
 
 
 # Model selection
-import model_simple_structures
 model_param = meep_utils.process_param(sys.argv[1:])
 model = Wedge_model(**model_param)
 
@@ -161,10 +160,10 @@ def complex_vec(self, vec):   ## Note: the 'vec' coordinates are _relative_ to t
 if "snapshote" in model.comment:
     slices += [meep_utils.Slice(model=model, field=f, components=(meep.Ex, meep.Ey, meep.Ez), at_t=np.inf, name='SnapshotE')]
 
-slices += [meep_utils.Slice(field=f, components=(meep.Ex), min_timestep=.1e-12, name='wedge',
-                volume=meep.volume( 
-                        meep.vec(0, -model.size_y/2+model.pml_thickness, model.size_z/2-model.pml_thickness), 
-                        meep.vec(0,  model.size_y/2-model.pml_thickness, model.size_z/2-model.pml_thickness)),
+                #volume=meep.volume( 
+                        #meep.vec(0, -model.size_y/2+model.pml_thickness, model.size_z/2-model.pml_thickness), 
+                        #meep.vec(0,  model.size_y/2-model.pml_thickness, model.size_z/2-model.pml_thickness)),
+slices += [meep_utils.Slice(field=f, components=(meep.Ex), min_timestep=.1e-12, name='wedge', at_x=0, at_z=model.size_z/2-model.pml_thickness,
                 model=model, outputdir=model.simulation_name, pad=model.pml_thickness, outputhdf=True, outputvtk=True, outputgif=True)]
 
         #pad = model.pml_thickness