diff --git a/modules/contrail/util.py b/modules/contrail/util.py
index d73d6946c24033ad775f0ada16e816861fbe04cf..368b7f17f04cc8ccc27808db436a0066a7a801c3 100644
--- a/modules/contrail/util.py
+++ b/modules/contrail/util.py
@@ -103,10 +103,7 @@ def extract(mask_image, image_ts, clavrx_path):
     uwind_3d = uwind_3d.transpose('Pressure', 'Latitude', 'Longitude')
     vwind_3d = vwind_3d.transpose('Pressure', 'Latitude', 'Longitude')
     temp_3d = temp_3d.transpose('Pressure', 'Latitude', 'Longitude')
-
-    # uwind_2d = get_horz_layer(xr_dataset, 'u-wind', 500.0, lon_range=[lon_range[0], lon_range[1]], lat_range=[lat_range[0], lat_range[1]])
-    # vwind_2d = get_horz_layer(xr_dataset, 'v-wind', 500.0, lon_range=[lon_range[0], lon_range[1]], lat_range=[lat_range[0], lat_range[1]])
-    # shear_2d = shearing_deformation(uwind_2d, vwind_2d)
+    rh_3d = rh_3d.transpose('Pressure', 'Latitude', 'Longitude')
 
     horz_shear_3d = shearing_deformation(uwind_3d, vwind_3d)
     static_3d = static_stability(temp_3d.coords['Pressure'] * units.hPa, temp_3d)