snapshot...

modules/icing/pirep_goes.py

@@ -31,7 +31,8 @@ l1b_grd_dct = {name: [] for name in l1b_ds_list}
 
 ds_list = ['cld_height_acha', 'cld_geo_thick', 'cld_press_acha', 'sensor_zenith_angle', 'cloud_type', 'supercooled_prob_acha',
            'supercooled_cloud_fraction', 'cld_temp_acha', 'cld_opd_acha', 'cloud_phase', 'solar_zenith_angle',
-           'cloud_mask', 'cld_reff_acha']
+           'cloud_mask', 'cld_reff_acha', 'cld_reff_dcomp', 'cld_reff_dcomp_1', 'cld_reff_dcomp_2', 'cld_reff_dcomp_3',
+           'cld_opd_dcomp', 'cld_opd_dcomp_1', 'cld_opd_dcomp_2', 'cld_opd_dcomp_3', 'cld_cwp_dcomp']
 ds_grd_dct = {name: [] for name in ds_list}
 
 a_clvr_file = '/Users/tomrink/data/clavrx/clavrx_OR_ABI-L1b-RadC-M3C01_G16_s20190020002186.level2.nc'
@@ -112,12 +113,13 @@ def create_file(filename, data_dct, ds_list, lon_c, lat_c, time_s):
     # copy relevant attributes
     for ds_name in ds_list:
         hf5[ds_name].attrs.create('standard_name', data=h5f_expl[ds_name].attrs.get('standard_name'))
+        hf5[ds_name].attrs.create('long_name', data=h5f_expl[ds_name].attrs.get('long_name'))
         hf5[ds_name].attrs.create('units', data=h5f_expl[ds_name].attrs.get('units'))
 
     hf5.close()
 
 
-def run(pirep_dct, outfile=None):
+def run(pirep_dct, outfile=None, outfile_l1b=None):
     time_keys = list(pirep_dct.keys())
 
     nav = GEOSNavigation(sub_lon=-75.0, CFAC=5.6E-05, COFF=-0.101332, LFAC=-5.6E-05, LOFF=0.128212, num_elems=2500, num_lines=1500)
@@ -186,13 +188,15 @@ def run(pirep_dct, outfile=None):
     lat_c = np.array(lat_c)
     time_s = np.array(time_s)
 
-    create_file(outfile, data_dct, ds_list, lon_c, lat_c, time_s)
+    if outfile is not None:
+        create_file(outfile, data_dct, ds_list, lon_c, lat_c, time_s)
 
     data_dct = {}
     for ds_name in ds_list:
         data_dct[ds_name] = np.array(l1b_grd_dct[ds_name])
 
-    create_file(outfile, data_dct, l1b_ds_list, lon_c, lat_c, time_s)
+    if outfile_l1b is not None:
+        create_file(outfile_l1b, data_dct, l1b_ds_list, lon_c, lat_c, time_s)
 
 
 def analyze(ice_dct, no_ice_dct):
@@ -255,4 +259,8 @@ def analyze(ice_dct, no_ice_dct):
         filename = os.path.split(files[i])[1]
         so = re.search('_s\\d{11}', filename)
         dt_str = so.group()
-        print(dt_str[2:])
\ No newline at end of file
+        print(dt_str[2:])
+
+
+def apply_qc_icing_pireps():
+    pass
\ No newline at end of file