snapshot...

modules/icing/pirep_goes.py

@@ -642,11 +642,19 @@ def analyze2(filename, filename_l1b):
     plt.show()
 
 
+x_a = 10
+x_b = 30
+y_a = x_a
+y_b = x_b
+nx = ny = (x_b - x_a)
+nx_x_ny = nx * ny
+
+
 def run_daynight(filename, filename_l1b, day_night='ANY'):
     f = h5py.File(filename, 'r')
     f_l1b = h5py.File(filename_l1b, 'r')
 
-    solzen = f['solar_zenith_angle'][:, 10:30, 10:30]
+    solzen = f['solar_zenith_angle'][:, y_a:y_b, x_a:x_b]
     num_obs = solzen.shape[0]
 
     idxs = []
@@ -693,12 +701,6 @@ def run_daynight(filename, filename_l1b, day_night='ANY'):
 
 def run_qc(filename, filename_l1b, day_night='ANY'):
     pass_thresh_frac = 0.4
-    x_a = 10
-    x_b = 30
-    y_a = x_a
-    y_b = x_b
-    nx = ny = (x_b - x_a)
-    nx_x_ny = nx * ny
 
     f = h5py.File(filename, 'r')
     icing_alt = f['icing_altitude'][:]