In [1]:
import numpy as np
import datetime as dt
In [2]:
flight = '20150709'
In [3]:
if flight == '20150617':
# Time/[temp at same time] when last evidence of ice was observed
mlBotTemp = np.array([np.nan, 2.3943, np.nan, 3.3489, 4.3656, 2.5243, 5.2405],dtype=float)
mlBotTime = np.array([np.nan, 17018, np.nan, 20530, 22632, 23810, 25750],dtype=float)
# S1 - No particles warmer than -5
# S2 - No particles between +2 and +8, so ML bottom actually warmer than +2
# S3 - No particles warmer than -4
# Time/[temp at same time] when first evidence of liquid water was observed
mlTopTemp = np.array([np.nan, 0.3578, np.nan, 0.9657, 0.4826, 0.1111, 0.3433],dtype=float)
mlTopTime = np.array([np.nan, 17074, np.nan, 20580, 22512, 23857, 25635],dtype=float)
# S1 - No particles warmer than -5
# S3 - No particles warmer than -4
elif flight == '20150620':
mlBotTemp = np.array([np.nan, 6.8400, np.nan, 5.9069, 5.7851, 4.6184, 4.7246],dtype=float)
mlBotTime = np.array([np.nan, 18647, np.nan, 21219, 22080, 24352, 28124],dtype=float)
# S1 - No particles warmer than -9
# S2 - Very little melting all the way through +7 (nothing recorded warmer)
# S3 - No particles warmer than -6
mlTopTemp = np.array([np.nan, 0.0282, np.nan, 0.1974, 0.1935, 2.2460, 0.0326],dtype=float)
mlTopTime = np.array([np.nan, 18582, np.nan, 21122, 22210, 24276, 28268],dtype=float)
# S1 - No particles warmer than -9
# S3 - No particles warmer than -6
elif flight == '20150701':
mlBotTemp = np.array([3.4335],dtype=float)
mlBotTime = np.array([22482],dtype=float)
mlTopTemp = np.array([0.7108],dtype=float)
mlTopTime = np.array([22375],dtype=float)
elif flight == '20150702':
mlBotTemp = np.array([np.nan, 1.5034, np.nan],dtype=float)
mlBotTime = np.array([np.nan, 15826, np.nan],dtype=float)
# S1 - No particles warmer than -13
# S3 - No particles between ~+1 and +9.59 - ML bottom somewhere in here
mlTopTemp = np.array([np.nan, 0.7152, np.nan],dtype=float)
mlTopTime = np.array([np.nan, 15852, np.nan],dtype=float)
# S1 - No particles warmer than -13
# S3 - Melting not observed at T<0.79; no particles warmer than this before +9.59
elif flight == '20150706':
mlBotTemp = np.array([np.nan, 1.9535, 3.7707, 4.9644, 1.9219, 3.5011, 3.0759, 3.6202],dtype=float)
mlBotTime = np.array([np.nan, 13099, 16228, 17252, 20851, 21872, 23420, 24383],dtype=float)
mlTopTemp = np.array([np.nan, 0.8741, 0.0271, 0.2957, 0.0111, 0.3134, 0.1221, 0.0554],dtype=float)
mlTopTime = np.array([np.nan, 13054, 16336, 17156, 20917, 21792, 23492, 24281],dtype=float)
elif flight == '20150709':
mlBotTemp = np.array([0.827, 2.76, 1.584, 2.2527, 2.8054, 4.7628, 3.823, 5.7988, 3.7277,
5.6487, 2.3316, 4.7067, 3.6857, 4.9586, 2.5909, 5.7844],dtype=float)
mlBotTime = np.array([9166, 10102, 10997, 11769, 13072, 14202, 14867, 15861, 16956, 17932,
19429, 20445, 21312, 22169, 23075, 23942],dtype=float)
mlTopTemp = np.array([0.0104, 0.0034, 0.0423, 0.1119, 0.4757, 3.4426, 1.3057, 2.5534, 0.4398,
2.0643, 0.0656, 1.5315, 0.6095, 2.740, 0.6454, 2.5692],dtype=float)
mlTopTime = np.array([9185, 10028, 11036, 11711, 13146, 14163, 14903, 15820, 17038,
17866, 19494, 20391, 21391, 22136, 23121, 23882],dtype=float)
In [4]:
mlBtime_hhmmss = np.floor(mlBotTime/3600)*10000 + np.floor(np.mod(mlBotTime,3600)/60)*100 + np.floor(np.mod(mlBotTime,60))
mlTtime_hhmmss = np.floor(mlTopTime/3600)*10000 + np.floor(np.mod(mlTopTime,3600)/60)*100 + np.floor(np.mod(mlTopTime,60))
In [5]:
print('Melting layer top times (hhmmss)')
for t in mlBtime_hhmmss:
if np.isnan(t):
print('\tN/A')
else:
print('\t{:.0f}'.format(t))
print('\nMelting layer bottom times (hhmmss)')
for t in mlTtime_hhmmss:
if np.isnan(t):
print('\tN/A')
else:
print('\t{:.0f}'.format(t))
print('\nMelting layer top temp (deg C)')
for t in mlTopTemp:
if np.isnan(t):
print('\tN/A')
else:
print('\t{:.2f}'.format(t))
print('\nMelting layer bottom temp (deg C)')
for t in mlBotTemp:
if np.isnan(t):
print('\tN/A')
else:
print('\t{:.2f}'.format(t))