snapshot...

modules/icing/util.py

@@ -93,20 +93,24 @@ def run_icing_predict(clvrx_dir='/Users/tomrink/data/clavrx/RadC/', output_dir=h
                 nav = get_navigation(satellite, domain)
                 lons_2d, lats_2d, x_rad, y_rad = get_lon_lat_2d_mesh(nav, ll, cc, offset=int(8 / res_fac))
 
+            bt_fld_name = 'temp_10_4um_nom'
+
             ancil_data_dct, _, _ = make_for_full_domain_predict(h5f, name_list=
-                                ['solar_zenith_angle', 'sensor_zenith_angle', 'cld_height_acha', 'cld_geo_thick', 'temp_10_4um_nom'],
+                                ['solar_zenith_angle', 'sensor_zenith_angle', 'cld_height_acha', 'cld_geo_thick', bt_fld_name],
                                 satellite=satellite, domain=domain, res_fac=res_fac)
         elif satellite == 'VIIRS':
             data_dct, ll, cc, lats_2d, lons_2d = make_for_full_domain_predict_viirs_clavrx(h5f, name_list=train_params, res_fac=res_fac)
+
+            bt_fld_name = 'temp_11_0um_nom'
             ancil_data_dct, _, _, _, _ = make_for_full_domain_predict_viirs_clavrx(h5f, name_list=
-                                ['solar_zenith_angle', 'sensor_zenith_angle', 'cld_height_acha', 'cld_geo_thick', 'temp_11_0um_nom'],
+                                ['solar_zenith_angle', 'sensor_zenith_angle', 'cld_height_acha', 'cld_geo_thick', bt_fld_name],
                                 res_fac=res_fac)
 
         num_elems, num_lines = len(cc), len(ll)
 
         satzen = ancil_data_dct['sensor_zenith_angle']
         solzen = ancil_data_dct['solar_zenith_angle']
-        bt_10_4 = ancil_data_dct['temp_10_4um_nom']
+        bt_10_4 = ancil_data_dct[bt_fld_name]
         day_idxs = []
         nght_idxs = []
         avg_bt = []