diff --git a/modules/icing/util.py b/modules/icing/util.py
index 2797be110bebb8b290f9e3935b90dfdc289b41dd..4bdcfc8c59168320b7191a66b60ce36bbe6328ee 100644
--- a/modules/icing/util.py
+++ b/modules/icing/util.py
@@ -12,6 +12,20 @@ from aeolus.datasource import CLAVRx, CLAVRx_VIIRS, GOESL1B, CLAVRx_H08
 import h5py
 import datetime
 
+# 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 get_training_parameters(day_night='DAY', l1b_andor_l2='both'):
     if day_night == 'DAY':