diff --git a/modules/util/util.py b/modules/util/util.py
index 52c4c368e2b8210c15194e0a105b782c22f1b4bb..a3b94576be9dcdec198d1075d756f68ae36ab6a8 100644
--- a/modules/util/util.py
+++ b/modules/util/util.py
@@ -663,12 +663,13 @@ rho_water = 1.
 rho_ice = 0.917
 
 
-def compute_lwc_iwc(iphase, reff, opd):
+def compute_lwc_iwc(iphase, reff, opd, geo_dz):
     xy_shape = iphase.shape
 
     iphase = iphase.flatten()
     reff = reff.flatten()
     opd = opd.flatten()
+    geo_dz = geo_dz.flatten()
 
     lwp_dcomp = np.zeros(reff.shape[0])
     iwp_dcomp = np.zeros(reff.shape[0])
@@ -681,6 +682,9 @@ def compute_lwc_iwc(iphase, reff, opd):
     iwp_dcomp[ice] = 0.667 * opd[ice] * rho_ice * reff[ice]
     lwp_dcomp[no_ice] = 0.55 * opd[no_ice] * rho_water * reff[no_ice]
 
+    iwp_dcomp /= geo_dz
+    lwp_dcomp /= geo_dz
+
     lwp_dcomp = lwp_dcomp.reshape(xy_shape)
     iwp_dcomp = iwp_dcomp.reshape(xy_shape)