Add files via upload

Author: cococastano

etch_model_version3_cellular_automata.py

@@ -60,15 +60,16 @@
 #recipe_steps = {'step01':{'bosch':15,'iso':100,'cycles':7},
 ##                'step02':{'bosch':300,'iso':None,'cycles':300},
 #                'step03':{'bosch':None,'iso':100,'cycles':1}}
-recipe_steps = {'step01':{'bosch':20,'iso':125,'cycles':7},
+recipe_steps = {'step01':{'bosch':15,'iso':100,'cycles':4},
 #                'step02':{'bosch':300,'iso':None,'cycles':300},
-                'step03':{'bosch':None,'iso':100,'cycles':1}}
+                'step03':{'bosch':None,'iso':800,'cycles':1}}
 
 
 # load mask
 im_dir = 'C:/Users/nicas/Documents/E241-MicroNanoFab/masks/'
-im_file = 'mask5_R5_C3_v0.png'
-pixel_um_conv = 49.291/80.4384
+im_file = 'python_model_fil_sq.png'
+pixel_um_conv = 276/100  # px/um
+gap = 249/pixel_um_conv
 # 151.37/100  # for R2_C2
 # for R5_C3: 49.291/80.4384  # px/um
 
@@ -85,7 +86,7 @@
 out_file = im_dir + im_file[:dummy_i] + '_out' + im_file[dummy_i:]
 cv2.imwrite(out_file, conts_im)
 
-cell_size = 5 # microns
+cell_size = 8 # microns
 wafer_thickness = 500  # microns
 
 t_start = 0 # seconds
@@ -102,13 +103,19 @@
 def vert_rate(z):
     a = 0.141
     b = 0.0007
-    return a*np.exp(-b*z)
+#    return a*np.exp(-b*z)
+    return (0.14-0.02) + z*(0.03/500)
+#    d = a*np.exp(-b*z)
+#    return d + d/2* np.cos(z * 2*np.pi/10)
 def horiz_rate(z):
-    return 0.8#0.8*vert_rate(z)
+    r = (0.65-0.1)+z*(0.1/500)
+    return r*vert_rate(z)
+#    return 0.6 + 0.2* np.cos(z * 2*np.pi/10)
 
 #vert_rate = 8.5/60  # um/s
 def bosch_vert_step(z):
-    return 0.84 - 0.1/500*z
+    return (0.84-0.1) + 0.1/500*z
+#    return 0.84 + 0.4* np.cos(z * 2*np.pi/10)
 
 #bosch_vert_step = 0.84  # um/step
 
@@ -268,8 +275,10 @@ def bosch_vert_step(z):
 ax = fig.add_subplot(111, projection='3d')
 #ax = plt.subplot(111,projection='polar')
 plot_flag = True
-
+#
+#dummy_steps = loop_steps[loop_steps.index(step)-1:]
 for step_i, step in enumerate(loop_steps,start=1):
+#for step_i, step in enumerate(dummy_steps,start=loop_steps.index(step)):
 
     #plot some results
     if step_i % (int(len(loop_steps)/10)) == 0:
@@ -616,7 +625,7 @@ def bosch_vert_step(z):
 
 ######
 
-pts = etch_grid[loop_steps[-1]][0]
+pts = etch_grid['step03_iso01_isotime520']#loop_steps[-1]][0]
 states = etch_grid[loop_steps[-1]][1]
 iden = etch_grid[loop_steps[-1]][2]
 
@@ -653,13 +662,15 @@ def bosch_vert_step(z):
 plotter.add_scalar_bar(title='z_height',height=0.08,width=0.4,
                        position_x=0.01,position_y=0.1)
 
+pts = etch_grid[loop_steps[-1]][0]
 
-   
-cells = np.asarray(exposed_cells.points)
-#obj = make_grid([cells],cell_size)
+pts = etch_grid['step03_iso01_isotime520'][0]#loop_steps[-1]][0]   
+#cells = np.asarray(exposed_cells.points)
+obj = make_grid([pts],cell_size)
+
+#obj = pv.read('C:/Users/nicas/Documents/E241-MicroNanoFab/codes/etch_model_version5_1/exposed_obj.vtk')
+#obj = make_grid([np.array(obj.points)],cell_size)
 
-obj = pv.read('C:/Users/nicas/Documents/E241-MicroNanoFab/codes/etch_model_version5_1/neigh_obj.vtk')
-obj = make_grid([np.array(obj.points)],cell_size)
 slices = obj.slice(normal=[1,1,0])
 plotter = pv.BackgroundPlotter(window_size=[1024, 768])
 plotter.add_mesh(obj, show_edges=False,

etch_sim_utilities.py

@@ -242,7 +242,9 @@ def make_grid(containers,cell_size):
         point_count = 0
         pts = np.array([])
         if type(container) == set: container = list(container)
+        n_pts = len(container)
         for i_cell, cell_center in enumerate(container):
+            if i_cell%1000 == 0: print('building cell %i of %i' %(i_cell,n_pts))
             cells.append(8)
             for p in range(8):
                 cells.append(point_count)
@@ -400,7 +402,7 @@ def compute_normals(points, use_nn=False, flat=False,
 #    ax = fig.add_subplot(111, projection='3d')
     if type(points) == dict:
         pts = [list(cell) for cell in points]
-        cloud = make_cloud([exposed_cells])[0]
+        cloud = make_cloud([points])[0]
         cloud = np.vstack((cloud,pts))
         knns = NearestNeighbors(n_neighbors=use_nn).fit(cloud)
         dists, indices = knns.kneighbors(cloud, return_distance=True)
@@ -424,7 +426,7 @@ def compute_normals(points, use_nn=False, flat=False,
             unit_norms[pt] = unit_norm
 
     elif points.shape[0] == 3:
-        cloud = make_cloud([exposed_cells])[0]
+        cloud = make_cloud([points])[0]
         knns = NearestNeighbors(n_neighbors=use_nn).fit(cloud)
         dists, indices = knns.kneighbors(cloud, return_distance=True)
         pt_i = np.where((cloud[:,0] == points[0]) & (cloud[:,1] == points[1]) 
@@ -577,12 +579,19 @@ def compute_normals_for_cells(exposed_cells,removed_cells,
                                 normal_vect[1] += round(y_v - y,3)
                                 normal_vect[2] += round(z_v - z,3)
             if np.linalg.norm(normal_vect) == 0: 
+                # if this happens it is a point in space with a spehere of 
+                # neighbors making the normal_vect magnitude 0 so jus assign a 
+                # dummy normal
+                normal_vects[cell] = np.array([1,0,0])
                 alt_method_flag = True
-                alt_method_points[cell] = []
             else:
                 normal_vects[cell] = (normal_vect / np.linalg.norm(normal_vect))
     if alt_method_flag == True:
-        alt_method_points = compute_normals(alt_method_points, use_nn=8)
+        print('COMPUTING NORMAL W/ PT CLOUD')
+        try:
+            alt_method_points = compute_normals(alt_method_points, use_nn=5)
+        except:
+            print('\tEXCEPTION')
 
     print('\t\t\tassigning normals to %i cells' % n_normals)
     for cell in list(exposed_cells.keys()):