From 223ae9db16bb261780e50d44152ed9c3db15eb3d Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sun, 7 Mar 2021 11:23:01 -0600
Subject: [PATCH] snapshot...

---
 modules/icing/pirep_goes.py | 51 ++++++++++++++++++++++++++++++++++++-
 1 file changed, 50 insertions(+), 1 deletion(-)

diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py
index 772f9ac5..761f6655 100644
--- a/modules/icing/pirep_goes.py
+++ b/modules/icing/pirep_goes.py
@@ -76,8 +76,16 @@ def run():
     obs_times = []
     obs_lons = []
     obs_lats = []
+    temp_8_5um_nom = []
     temp_10_4um_nom = []
+    temp_11_0um_nom = []
+    temp_13_3um_nom = []
     cld_height_acha = []
+    sc_cld_frac = []
+    cld_tmp_acha = []
+    cld_opd_acha = []
+    cld_phase = []
+    solar_zenith_angle = []
 
     for idx, time in enumerate(time_keys):
         reports = ice_dict[time]
@@ -115,12 +123,53 @@ def run():
             gvals_a = get_grid_values(h5f, 'temp_10_4um_nom', ll[0], cc[0], 20)
             if gvals_a is None:
                 continue
+
             gvals_b = get_grid_values(h5f, 'cld_height_acha', ll[0], cc[0], 20)
             if gvals_b is None:
                 continue
 
+            gvals_c = get_grid_values(h5f, 'supercooled_cloud_fraction', ll[0], cc[0], 20)
+            if gvals_c is None:
+                continue
+
+            gvals_d = get_grid_values(h5f, 'cld_temp_acha', ll[0], cc[0], 20)
+            if gvals_c is None:
+                continue
+
+            gvals_e = get_grid_values(h5f, 'cld_opd_acha', ll[0], cc[0], 20)
+            if gvals_c is None:
+                continue
+
+            gvals_f = get_grid_values(h5f, 'cloud_phase', ll[0], cc[0], 20)
+            if gvals_f is None:
+                continue
+
+            gvals_g = get_grid_values(h5f, 'solar_zenith_angle', ll[0], cc[0], 20)
+            if gvals_g is None:
+                continue
+
+            gvals_h = get_grid_values(h5f, 'temp_8_5um_nom', ll[0], cc[0], 20)
+            if gvals_h is None:
+                continue
+
+            gvals_i = get_grid_values(h5f, 'temp_11_0um_nom', ll[0], cc[0], 20)
+            if gvals_i is None:
+                continue
+
+            gvals_j = get_grid_values(h5f, 'temp_13_3um_nom', ll[0], cc[0], 20)
+            if gvals_j is None:
+                continue
+
             obs_times.append(time)
             obs_lons.append(lon)
             obs_lats.append(lat)
             temp_10_4um_nom.append(gvals_a)
-            cld_height_acha.append(gvals_b)
\ No newline at end of file
+            cld_height_acha.append(gvals_b)
+            sc_cld_frac.append(gvals_c)
+            cld_tmp_acha.append(gvals_d)
+            cld_opd_acha.append(gvals_e)
+            cld_phase.append(gvals_f)
+            solar_zenith_angle.append(gvals_g)
+            temp_8_5um_nom.append(gvals_h)
+            temp_11_0um_nom.append(gvals_i)
+            temp_13_3um_nom.append(gvals_j)
\ No newline at end of file
-- 
GitLab