Skip to content
Snippets Groups Projects
Commit bef5fd2b authored by tomrink's avatar tomrink
Browse files

snapshot...

parent a160c239
No related branches found
No related tags found
No related merge requests found
...@@ -4,7 +4,7 @@ import pickle ...@@ -4,7 +4,7 @@ import pickle
import matplotlib.pyplot as plt import matplotlib.pyplot as plt
import os import os
from util.util import get_time_tuple_utc, GenericException, add_time_range_to_filename, is_night, is_day, \ from util.util import get_time_tuple_utc, GenericException, add_time_range_to_filename, is_night, is_day, \
get_grid_values_all, check_oblique, make_times, find_bin_index, get_timestamp get_grid_values_all, check_oblique, make_times, find_bin_index, get_timestamp, homedir
from aeolus.datasource import CLAVRx, CLAVRx_VIIRS, GOESL1B from aeolus.datasource import CLAVRx, CLAVRx_VIIRS, GOESL1B
import h5py import h5py
import re import re
...@@ -24,7 +24,7 @@ dir_fmt = '%Y_%m_%d_%j' ...@@ -24,7 +24,7 @@ dir_fmt = '%Y_%m_%d_%j'
ds_dct = {} ds_dct = {}
goes_ds_dct = {} goes_ds_dct = {}
pirep_file = '/home/rink/data/pireps/pireps_20180101_20200331.csv' pirep_file = homedir+'data/pirep/pireps_20180101_20200331.csv'
# --- CLAVRx Radiometric parameters and metadata ------------------------------------------------ # --- CLAVRx Radiometric parameters and metadata ------------------------------------------------
l1b_ds_list = ['temp_10_4um_nom', 'temp_11_0um_nom', 'temp_12_0um_nom', 'temp_13_3um_nom', 'temp_3_75um_nom', l1b_ds_list = ['temp_10_4um_nom', 'temp_11_0um_nom', 'temp_12_0um_nom', 'temp_13_3um_nom', 'temp_3_75um_nom',
...@@ -57,9 +57,9 @@ ds_range = ['actual_range' for i in range(23)] + [None for i in range(3)] ...@@ -57,9 +57,9 @@ ds_range = ['actual_range' for i in range(23)] + [None for i in range(3)]
# An example file for accessing and copying metadata # An example file for accessing and copying metadata
a_clvr_file = '/home/rink/data/clavrx/clavrx_OR_ABI-L1b-RadC-M3C01_G16_s20190020002186.level2.nc' a_clvr_file = homedir+'data/clavrx/clavrx_OR_ABI-L1b-RadC-M3C01_G16_s20190020002186.level2.nc'
# VIIRS # VIIRS
#a_clvr_file = '/home/rink/data/clavrx/clavrx_snpp_viirs.A2019071.0000.001.2019071061610.uwssec_B00038187.level2.h5' #a_clvr_file = homedir+'data/clavrx/clavrx_snpp_viirs.A2019071.0000.001.2019071061610.uwssec_B00038187.level2.h5'
# Location of files for tile/FOV extraction # Location of files for tile/FOV extraction
icing_files = [f for f in glob.glob('/data/Personal/rink/icing_ml/icing_2*_DAY.h5')] icing_files = [f for f in glob.glob('/data/Personal/rink/icing_ml/icing_2*_DAY.h5')]
...@@ -81,7 +81,7 @@ train_params_night = ['cld_height_acha', 'cld_geo_thick', 'supercooled_cloud_fra ...@@ -81,7 +81,7 @@ train_params_night = ['cld_height_acha', 'cld_geo_thick', 'supercooled_cloud_fra
'cld_press_acha', 'cld_reff_acha', 'cld_opd_acha', 'cloud_phase', 'cloud_mask'] 'cld_press_acha', 'cld_reff_acha', 'cld_opd_acha', 'cloud_phase', 'cloud_mask']
def setup(pirep_file='/home/rink/data/pireps/pireps_20180101_20200331.csv'): def setup(pirep_file=pirep_file):
ice_dict, no_ice_dict, neg_ice_dict = pirep_icing(pirep_file) ice_dict, no_ice_dict, neg_ice_dict = pirep_icing(pirep_file)
return ice_dict, no_ice_dict, neg_ice_dict return ice_dict, no_ice_dict, neg_ice_dict
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment