diff --git a/modules/util/viirs_l1b_l2.py b/modules/util/viirs_l1b_l2.py
index d7e072a1f19768e40035f95b83806bd123f6aa50..c85ad7776a31f6474117afc38cb006b0eb5cd139 100644
--- a/modules/util/viirs_l1b_l2.py
+++ b/modules/util/viirs_l1b_l2.py
@@ -1,6 +1,6 @@
 import numpy as np
 import h5py
-from util.util import get_grid_values, get_grid_values_all, is_night
+from util.util import get_grid_values, get_grid_values_all, is_night, compute_lwc_iwc
 import glob
 import os
 from pathlib import Path
@@ -240,6 +240,18 @@ def scan(directory):
         print(num_files, num_keep)
 
 
+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)
+    reff = get_grid_values_all(h5f, 'cld_reff_nlcomp')
+    opd = get_grid_values_all(h5f, 'cld_opd_nlcomp')
+    cld_dz = get_grid_values_all(h5f, 'cld_geo_thick')
+
+    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}