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

snapshot...

parent 6ff38478
No related branches found
No related tags found
No related merge requests found
...@@ -19,7 +19,8 @@ goes_ds_dct = {} ...@@ -19,7 +19,8 @@ goes_ds_dct = {}
ice_dict = None ice_dict = None
no_ice_dict = None no_ice_dict = None
time_keys = None time_keys = None
pirep_file = '/home/rink/data/pireps/pireps_2019010000_2019063023.csv' #pirep_file = '/home/rink/data/pireps/pireps_2019010000_2019063023.csv'
pirep_file = '/home/rink/data/pireps/pireps_20160101_20190630.csv'
ds_list = ['temp_8_5um_nom', 'temp_10_4um_nom', 'temp_11_0um_nom', 'temp_13_3um_nom', 'cld_height_acha', ds_list = ['temp_8_5um_nom', 'temp_10_4um_nom', 'temp_11_0um_nom', 'temp_13_3um_nom', 'cld_height_acha',
'supercooled_cloud_fraction', 'cld_temp_acha', 'cld_opd_acha', 'cloud_phase', 'solar_zenith_angle', 'supercooled_cloud_fraction', 'cld_temp_acha', 'cld_opd_acha', 'cloud_phase', 'solar_zenith_angle',
...@@ -147,25 +148,6 @@ def run(pirep_dct, outfile=None): ...@@ -147,25 +148,6 @@ def run(pirep_dct, outfile=None):
def analyze(ice_dct, no_ice_dct): def analyze(ice_dct, no_ice_dct):
ice_times = np.array(list(ice_dct.keys()))
no_ice_times = np.array(list(no_ice_dct.keys()))
dt_str = '201601010000'
dto = datetime.datetime.strptime(dt_str, '%Y%m%d%H%M').replace(tzinfo=timezone.utc)
t0 = dto.timestamp()
dt_str = '201906300000'
dto = datetime.datetime.strptime(dt_str, '%Y%m%d%H%M').replace(tzinfo=timezone.utc)
t1 = dto.timestamp()
tline = np.arange(t0, t1, 600)
hist_a, edges = np.histogram(no_ice_times, bins=tline)
print(len(np.nonzero(hist_a)[0]))
hist_b, edges = np.histogram(ice_times, bins=tline)
print(len(np.nonzero(hist_b)[0]))
print(np.sum(np.logical_and(hist_a > 0, hist_b > 0)))
last_file = None last_file = None
ice_files = [] ice_files = []
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment