snapshot...

modules/util/viirs_surfrad.py

@@ -81,7 +81,7 @@ def process_cld_opd(grd_k):
     keep_cld = np.where(keep, (0.1 < grd_k), False)
     frac_cld = np.sum(keep_cld)/num_keep
     # if not (0.50 < frac_cld < 0.85):
-    if not (0.60 < frac_cld < 1.0):
+    if not (0.40 < frac_cld < 1.0):
         return None
     grd_k = np.where(np.invert(keep), 0, grd_k)  # Convert NaN to 0
     return grd_k
@@ -151,6 +151,7 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
 
                     h, b = np.histogram(valid_i.flatten(), bins=20, range=[0.0, 160.0])
                     hist_accum_valid_i += h
+
                     h, b = np.histogram(valid_m.flatten(), bins=20, range=[0.0, 160.0])
                     hist_accum_valid_m += h
 
@@ -164,6 +165,7 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
 
                     h, b = np.histogram(train_i.flatten(), bins=20, range=[0.0, 160.0])
                     hist_accum_train_i += h
+
                     h, b = np.histogram(train_m.flatten(), bins=20, range=[0.0, 160.0])
                     hist_accum_train_m += h
 
@@ -192,6 +194,7 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
 
         h, b = np.histogram(valid_i.flatten(), bins=20, range=[0.0, 160.0])
         hist_accum_valid_i += h
+
         h, b = np.histogram(valid_m.flatten(), bins=20, range=[0.0, 160.0])
         hist_accum_valid_m += h
 
@@ -205,6 +208,7 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
 
         h, b = np.histogram(train_i.flatten(), bins=20, range=[0.0, 160.0])
         hist_accum_train_i += h
+
         h, b = np.histogram(train_m.flatten(), bins=20, range=[0.0, 160.0])
         hist_accum_train_m += h
 
@@ -219,8 +223,11 @@ def run_all(directory, out_directory, day_night='ANY', pattern='clavrx_*.nc', st
     print('** total_num_train_samples, total_num_valid_samples: ', total_num_train_samples, total_num_valid_samples)
     print('--------------------------------------------------')
     print(hist_accum_train_i)
+    print('----------')
     print(hist_accum_train_m)
+    print('----------')
     print(hist_accum_valid_i)
+    print('----------')
     print(hist_accum_valid_m)