diff --git a/modules/icing/pirep_goes.py b/modules/icing/pirep_goes.py index 9306f8fe7ccf9e3cfb2e312fa892a258f5634ff3..6c3709997d2ad3d83b3e32c5765f0b5421808cba 100644 --- a/modules/icing/pirep_goes.py +++ b/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]