small changes

README.md

@@ -1,5 +1,6 @@
-![output example](/doc/SiAgSi.png)
+![output example](/doc/OutputExample.png)
 **Output example:** Field distribution inside layered Si\Ag\Si sphere
+and Poynting vector distribution in Ag sphere with poweflow lines
 calculated with Scattnlay.
 
 How to use scattnlay

tests/python/field-Ag-flow.py

@@ -186,7 +186,7 @@ try:
     scale_ticks = np.linspace(min_tick, max_tick, 11)
 
     # Interpolation can be 'nearest', 'bilinear' or 'bicubic'
-    #ax.set_title('Eabs')
+    ax.set_title(r'$Re(E \times H)$')
     cax = ax.imshow(Eabs_data, interpolation = 'nearest', cmap = cm.jet,
                         origin = 'lower'
                         #, vmin = min_tick, vmax = max_tick
@@ -241,7 +241,7 @@ try:
     #     patch = patches.PathPatch(path, facecolor='none', lw=1, edgecolor='white')
     #     ax.add_patch(patch)
  
-    # plt.savefig("SiAgSi.png")
+    plt.savefig("Ag-flow.png")
     plt.draw()
 
     plt.show()

tests/python/field-SiAgSi-flow.py

@@ -118,7 +118,7 @@ print "m =", m
 
 npts = 241
 
-factor=2.5
+factor=2.2
 scan = np.linspace(-factor*x[0, 2], factor*x[0, 2], npts)
 
 coordX, coordZ = np.meshgrid(scan, scan)
@@ -209,7 +209,7 @@ try:
     from matplotlib.path import Path
     #import matplotlib.patches as patches
 
-    flow_total = 41
+    flow_total = 21
     for flow in range(0,flow_total):
         flow_x, flow_z = GetFlow(scale_x, scale_z, Ec, Hc,
                                  min(scale_x)+flow*(scale_x[-1]-scale_x[0])/(flow_total-1),