diff --git a/modules/contrail/util.py b/modules/contrail/util.py
index f31389b0e416df2f29b20d63e83643562fd7f2fd..d2c6bb38223bf7247d93a305215ce708cc31fb0a 100644
--- a/modules/contrail/util.py
+++ b/modules/contrail/util.py
@@ -75,6 +75,7 @@ def extract(mask_image, image_ts, clavrx_path):
     for bin_num in np.unique(binned_indexes):
         bins_dict[bin_num] = np.where(binned_indexes == bin_num)[0]
 
+    # MetPy for shearing deformation, static stability, vertical shear?, others...
     # This does point by point computation of model parameters for each contrail pixel
     voxel_dict = {key: [] for key in bins_dict.keys()}
     for key in bins_dict.keys():