diff --git a/modules/util/util.py b/modules/util/util.py
index a3b94576be9dcdec198d1075d756f68ae36ab6a8..24e7c7fca5e094e1d7369298f23deaa87a1b0354 100644
--- a/modules/util/util.py
+++ b/modules/util/util.py
@@ -662,6 +662,21 @@ def make_for_full_domain_predict(h5f, name_list=None, satellite='GOES16', domain
 rho_water = 1.
 rho_ice = 0.917
 
+# real(kind=real4), parameter:: Rho_Water = 1.0    !g / m ^ 3
+# real(kind=real4), parameter:: Rho_Ice = 0.917    !g / m ^ 3
+#
+# !--- compute
+# cloud
+# water
+# path
+# if (Iphase == 0) then
+# Cwp_Dcomp(Elem_Idx, Line_Idx) = 0.55 * Tau * Reff * Rho_Water
+# Lwp_Dcomp(Elem_Idx, Line_Idx) = 0.55 * Tau * Reff * Rho_Water
+# else
+# Cwp_Dcomp(Elem_Idx, Line_Idx) = 0.667 * Tau * Reff * Rho_Ice
+# Iwp_Dcomp(Elem_Idx, Line_Idx) = 0.667 * Tau * Reff * Rho_Ice
+# endif
+
 
 def compute_lwc_iwc(iphase, reff, opd, geo_dz):
     xy_shape = iphase.shape