Closes #7, updates .gitignore for hier data

Fixes VisibleDeprecationWarning by converting all instances of supp to an
int when loaded from the dict. Squelches matplotlib warning raised from
internal workings of pyplot.pause.

pyIB2d/.gitignore

@@ -1,5 +1,6 @@
 # Generated data
 viz_IB2d/
+hier_IB2d_data/
 
 # Byte-compiled / optimized / DLL files
 __pycache__/

pyIB2d/IBM_Blackbox/IBM_Driver.py

@@ -104,7 +104,7 @@ def main(struct_name, mu, rho, grid_Info, dt, T_FINAL, model_Info):
     Ly = grid_Info['Ly']   # Length of Eulerian grid in y-coordinate
     dx = grid_Info['dx']   # Spatial-size in x
     dy = grid_Info['dy']   # Spatial-size in y
-    supp = grid_Info['supp'] # Delta-function support
+    supp = int(grid_Info['supp']) # Delta-function support
 
     # PRINTING/PLOTTING INFO #
     pDump = grid_Info['pDump']              # Print (Plot) Dump interval

pyIB2d/IBM_Blackbox/Supp.py

@@ -80,7 +80,7 @@ def please_Move_Lagrangian_Point_Positions(u, v, xL_P, yL_P, xL_H, yL_H, x, y,\
     Ly =   grid_Info['Ly']
     dx =   grid_Info['dx']
     dy =   grid_Info['dy']
-    supp = grid_Info['supp']
+    supp = int(grid_Info['supp'])
     Nb =   grid_Info['Nb']
     ds =   grid_Info['ds']
 

pyIB2d/IBM_Blackbox/please_Find_Lagrangian_Forces_On_Eulerian_grid.py

@@ -88,7 +88,7 @@ def please_Find_Lagrangian_Forces_On_Eulerian_grid(dt, current_time, xLag, yLag,
     Ly =    grid_Info['Ly'] # Length of Eulerian grid in y-coordinate
     dx =    grid_Info['dx'] # Spatial-size in x
     dy =    grid_Info['dy'] # Spatial-size in y
-    supp =  grid_Info['supp'] # Delta-function support
+    supp =  int(grid_Info['supp']) # Delta-function support
     Nb =    grid_Info['Nb'] # # of Lagrangian pts. 
     ds =    grid_Info['ds'] # Lagrangian spacing
 
@@ -714,7 +714,7 @@ def give_Me_Delta_Function_Approximations_For_Force_Calc(x,y,grid_Info,xLag,yLag
     Ly =   grid_Info['Ly']
     dx =   grid_Info['dx']
     dy =   grid_Info['dy']
-    supp = grid_Info['supp']
+    supp = int(grid_Info['supp'])
     Nb =   grid_Info['Nb']
 
     # Find the indices of the points (xi, yj) where the 1D delta functions are 

pyIB2d/IBM_Blackbox/please_Plot_Results.py

@@ -30,6 +30,7 @@
 import matplotlib.pyplot as plt
 plt.rcParams['image.cmap'] = 'viridis'
 from math import sqrt
+import warnings
 
 ###########################################################################
 #
@@ -237,7 +238,9 @@ def please_Plot_Results(ds,X,Y,U,V,vort,uMag,p,chiX,chiY,\
     plt.box(on=True)
 
     plt.draw()
-    plt.pause(0.0001) #no idea why this is necessary, but it is
+    with warnings.catch_warnings():
+        warnings.simplefilter("ignore")
+        plt.pause(0.0001) #no idea why this is necessary, but it is
 
 
 

pyIB2d/IBM_Blackbox/please_Update_Fluid_Velocity.py

@@ -72,7 +72,7 @@ def please_Update_Fluid_Velocity(U, V, Fx, Fy, rho, mu, grid_Info, dt):
     Ly =   grid_Info['Ly']
     dx =   grid_Info['dx']
     dy =   grid_Info['dy']
-    supp = grid_Info['supp']
+    supp = int(grid_Info['supp'])
     Nb =   grid_Info['Nb']
     ds =   grid_Info['ds']