diff --git a/fireInIce/1_process/src/cores_to_shp.py b/fireInIce/1_process/src/cores_to_shp.py
index 3110b32c3dd2e2f9998b4abc772463f78d36d2c7..61584e1697e96fd33fd4320194789e81562ef210 100644
--- a/fireInIce/1_process/src/cores_to_shp.py
+++ b/fireInIce/1_process/src/cores_to_shp.py
@@ -2,8 +2,8 @@ import pandas as pd
 import geopandas as gpd
 def main(core_data,core_shpfile):
     cd = pd.read_excel(core_data)
-    cd['Latitude'][2] = "58° 51’ 16.1” N" #datafix 2016, C3
-    cd['Longitude'][2] = "134° 10’ 31.1” W" #datafix 2016, C3
+    cd.loc[2, 'Latitude'] = "58° 51’ 16.1” N" #datafix 2016, C3
+    cd.loc[2, 'Longitude'] = "134° 10’ 31.1” W" #datafix 2016, C3
     
     core_name = []
     Lat = []
diff --git a/fireInIce/2_visualize/src/xs_process.py b/fireInIce/2_visualize/src/xs_process.py
index e3beb0a13f1a689d8fb39bdcafa480f2d8f1e962..c9dcb8990d39fdf44f97bdbb4fd530ddebf2f838 100644
--- a/fireInIce/2_visualize/src/xs_process.py
+++ b/fireInIce/2_visualize/src/xs_process.py
@@ -157,15 +157,6 @@ def main(
     buff_west, buff_east, buff_south, buff_north = buffer_shapefile(
         west, east, length, south, north, height, 25.0
     )
-    domain_geom = Polygon(
-        zip(
-            [buff_west, buff_west, buff_east, buff_east],
-            [buff_north, buff_south, buff_south, buff_north],
-        )
-    )
-    domain_polygon = gpd.GeoDataFrame(index=[0], crs=utm_crs, geometry=[domain_geom])
-    buff_length = buff_west - buff_east
-    buff_height = buff_north - buff_south
 
     # build interpolators to sampling data from dem and ice field
     extent = list(get_raster_extent(proj_ds))
diff --git a/fireInIce/Snakefile_glacial_scan b/fireInIce/Snakefile_glacial_scan
index 56ab1e7bcc26bc2207a31b4b480eafd04221f91a..099526520fef1518b86b1bd1b0b20cf24540f46c 100644
--- a/fireInIce/Snakefile_glacial_scan
+++ b/fireInIce/Snakefile_glacial_scan
@@ -1,5 +1,3 @@
-HUC8s = [19010301,19010304]
-
 sb_id = "5e472c3ee4b0ff554f6837bc"
 
 rule all: