From 2e48e0042757b0d261eb140c004c5a47f1cbfd53 Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Fri, 23 Apr 2021 14:36:55 -0500
Subject: [PATCH] snapshot...
---
modules/icing/pirep_goes.py | 16 +++++++++-------
1 file changed, 9 insertions(+), 7 deletions(-)
diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 3c152108..28252cfd 100644
--- a/modules/icing/pirep_goes.py
+++ b/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'][:]
--
GitLab