snapshot...

modules/icing/pirep_goes.py

@@ -40,15 +40,10 @@ ds_types = ['f4' for i in range(23)] + ['i4' for i in range(3)]
 a_clvr_file = '/home/rink/data/clavrx/clavrx_OR_ABI-L1b-RadC-M3C01_G16_s20190020002186.level2.nc'
 #a_clvr_file = '/Users/tomrink/data/clavrx/clavrx_OR_ABI-L1b-RadC-M3C01_G16_s20190020002186.level2.nc'
 
-data_dir = '/Users/tomrink/data/icing/'
-icing_files = ['icing_2018010600_2018033022_QC_DAY.h5', 'icing_2018040100_2018063021_QC_DAY.h5',
-               'icing_2018070101_2018093022_QC_DAY.h5', 'icing_2018100115_2018123022_QC_DAY.h5',
-               'icing_2019010100_2019033023_QC_DAY.h5']
+icing_files = [f for f in glob.glob('/home/rink/data/icing/icing*.h5')]
 icing_l1b_files = []
 
-no_icing_files = ['no_icing_2018010600_2018033022_DAY.h5', 'no_icing_2018040100_2018063022_DAY.h5',
-                  'no_icing_2018070100_2018093022_DAY.h5', 'no_icing_2018100115_2018123022_DAY.h5',
-                  'no_icing_2019010100_2019033023_DAY.h5']
+no_icing_files = [f for f in glob.glob('/home/rink/data/icing/no_icing*.h5')]
 no_icing_l1b_files = []
 
 train_params_day = ['cld_height_acha', 'cld_geo_thick', 'supercooled_cloud_fraction', 'cld_temp_acha', 'cld_press_acha',
@@ -805,7 +800,7 @@ def apply_qc_icing_pireps(icing_alt, cld_top_hgt, cld_phase, cld_opd, cld_mask,
     return mask, idxs
 
 
-def fov_extract(outfile='/Users/tomrink/fovs_out.h5', train_params=train_params_day):
+def fov_extract(outfile='/home/rink/fovs_out.h5', train_params=train_params_day):
     ice_times = []
     icing_int_s = []
     ice_lons = []
@@ -826,7 +821,7 @@ def fov_extract(outfile='/Users/tomrink/fovs_out.h5', train_params=train_params_
     num_ice = 0
     for fidx in range(len(icing_files)):
         fname = icing_files[fidx]
-        f = h5py.File(data_dir+fname, 'r')
+        f = h5py.File(fname, 'r')
         h5_s_icing.append(f)
 
         times = f['time'][:]
@@ -878,7 +873,7 @@ def fov_extract(outfile='/Users/tomrink/fovs_out.h5', train_params=train_params_
     num_no_ice = 0
     for fidx in range(len(no_icing_files)):
         fname = no_icing_files[fidx]
-        f = h5py.File(data_dir+fname, 'r')
+        f = h5py.File(fname, 'r')
         h5_s_no_icing.append(f)
 
         times = f['time']
@@ -989,7 +984,7 @@ def fov_extract(outfile='/Users/tomrink/fovs_out.h5', train_params=train_params_
     h5f_expl.close()
 
 
-def box_extract(outfile='/Users/tomrink/box_out.h5', train_params=train_params_day):
+def box_extract(outfile='/home/rink/box_out.h5', train_params=train_params_day):
     icing_int_s = []
     ice_time_s = []
     no_ice_time_s = []
@@ -1006,7 +1001,7 @@ def box_extract(outfile='/Users/tomrink/box_out.h5', train_params=train_params_d
 
     for fidx in range(len(icing_files)):
         fname = icing_files[fidx]
-        f = h5py.File(data_dir+fname, 'r')
+        f = h5py.File(fname, 'r')
         h5_s_icing.append(f)
 
         times = f['time'][:]
@@ -1039,7 +1034,7 @@ def box_extract(outfile='/Users/tomrink/box_out.h5', train_params=train_params_d
     num_no_ice = 0
     for fidx in range(len(no_icing_files)):
         fname = no_icing_files[fidx]
-        f = h5py.File(data_dir+fname, 'r')
+        f = h5py.File(fname, 'r')
         h5_s_no_icing.append(f)
 
         times = f['time']
@@ -1143,8 +1138,8 @@ def run_mean_std():
     mean_dct = {}
     std_dct = {}
 
-    ice_flist = [f for f in glob.glob('/home/rink/data/icing/', 'icing*.h5')]
-    no_ice_flist = [f for f in glob.glob('/home/rink/data/icing/', 'no_icing*.h5')]
+    ice_flist = [f for f in glob.glob('/home/rink/data/icing/icing*.h5')]
+    no_ice_flist = [f for f in glob.glob('/home/rink/data/icing/no_icing*.h5')]
 
     ice_h5f_lst = [h5py.File(f, 'r') for f in ice_flist]
     no_ice_h5f_lst = [h5py.File(f, 'r') for f in no_ice_flist]