diff --git a/modules/util/viirs_surfrad.py b/modules/util/viirs_surfrad.py
index 3376c845e58476bfae47e5d66c5d6ec8a872acb1..a9e71a0337059d963046ed14ec4a6ab22025e5c9 100644
--- a/modules/util/viirs_surfrad.py
+++ b/modules/util/viirs_surfrad.py
@@ -243,72 +243,6 @@ def run(h5f, param_s, train_tiles, valid_tiles, lbl_param_s, lbl_train_tiles, lb
     return cnt_total, cnt_kept
 
 
-def scan(directory):
-
-    data_src = CLAVRx_VIIRS(directory)
-    files = data_src.flist
-
-    for idx, file in enumerate(files):
-        h5f = h5py.File(file, 'r')
-        ts = data_src.ftimes[idx][0]
-        try:
-            solzen = get_grid_values_all(h5f, 'solar_zenith_angle')
-        except Exception as e:
-            # print(e)
-            h5f.close()
-            continue
-
-        # if is_day(solzen) and moon_phase(ts):
-        if is_night(solzen) and moon_phase(ts):
-            print(file)
-        h5f.close()
-
-
-def scan_for_location(txt_file, lon_range=[111.0, 130.0], lat_range=[14.0, 32.0]):
-    with open(txt_file) as file:
-        for idx, fpath in enumerate(file):
-            fpath = fpath.strip()
-            h5f = h5py.File(fpath, 'r')
-            try:
-                lon_s = get_grid_values_all(h5f, 'longitude', stride=4)
-                lat_s = get_grid_values_all(h5f, 'latitude', stride=4)
-                c_lon, c_lat = lon_s[406, 400], lat_s[406, 400]
-                if (lon_range[0] < c_lon < lon_range[1]) and (lat_range[0] < c_lat < lat_range[1]):
-                    print(fpath)
-
-            except Exception as e:
-                # print(e)
-                h5f.close()
-                continue
-
-
-def test_nlcomp(file):
-    h5f = h5py.File(file, 'r')
-
-    cld_phs = get_grid_values_all(h5f, 'cloud_phase', scale_factor_name=None, range_name=None)
-    keep_0 = np.invert(np.isnan(cld_phs))
-
-    reff = get_grid_values_all(h5f, 'cld_reff_nlcomp')
-    keep_1 = np.invert(np.isnan(reff))
-
-    opd = get_grid_values_all(h5f, 'cld_opd_nlcomp')
-    keep_2 = np.invert(np.isnan(opd))
-
-    cld_dz = get_grid_values_all(h5f, 'cld_geo_thick')
-    keep_3 = np.logical_and(np.invert(np.isnan(cld_dz)), cld_dz > 5.0)
-
-    keep = keep_0 & keep_1 & keep_2 & keep_3
-
-    cld_phs = cld_phs[keep]
-    reff = reff[keep]
-    opd = opd[keep]
-    cld_dz = cld_dz[keep]
-
-    lwc_c, iwc_c = compute_lwc_iwc(cld_phs, reff, opd, cld_dz)
-
-    return lwc_c, iwc_c
-
-
 # def run_mean_std(directory):
 #
 #     data_dct = {name: [] for name in mod_res_params}