snapshot...

modules/util/infer_cloud_products.py

@@ -59,9 +59,12 @@ def infer_cloud_products(clvrx_path, output_dir, full_disk=True, satellite='GOES
         cld_opd_2 = np.where(cld_frac == 2, cld_opd_2 / 0.50, cld_opd_2)
         cld_opd_2 = np.where(cld_frac == 3, cld_opd_2 / 0.75, cld_opd_2)
 
-        cld_opd_1[np.invert(np.logical_and(cld_opd_1 < 0.0, cld_opd_1 > 160.0))] = -1.0
+        cld_opd_1[np.logical_and(cld_opd_1 < 0.0, cld_opd_1 > 160.0)] = -1.0
         cld_opd_1[np.isnan(opd)] = -1.0
 
+        cld_opd_2[np.logical_and(cld_opd_2 < 0.0, cld_opd_2 > 160.0)] = -1.0
+        cld_opd_2[np.isnan(opd)] = -1.0
+
         write_cld_prods_file_nc4(clvrx_str_time, out_file, cld_frac, cld_opd_1, cld_opd_2, x_rad, y_rad, None, None,
                                 satellite=satellite, domain=domain, has_time=True)
         t1 = time.time()