notebook.community

Edit and run 





In [1]:



    


import numpy as np
from matplotlib import pyplot as plt
import os
import warnings
warnings.filterwarnings("ignore",category=FutureWarning)
import xarray as xr
from matplotlib.ticker import MultipleLocator, FormatStrFormatter

warnings.filterwarnings("ignore",message='invalid value encountered in less_equal')

%matplotlib inline



    











In [2]:



    


npTypeZn = np.load('/Users/danstechman/GoogleDrive/PECAN-Data/mp-data/typeZone_stats.npz')



    











In [3]:



    


srCs1 = npTypeZn['srC_1JulS1'][()]
pfs1 = npTypeZn['pf_6JulS1'][()]
pfs2 = npTypeZn['pf_6JulS2'][()]

temps = npTypeZn['bin_mid']
temp_lte0 = temps[temps <= 0]
temp_gt0 = temps[temps > 0]
temp_gte5 = temps[temps >= 5]

lte0 = np.where(temps <= 0)[0]
gt0 = np.where(temps > 0)[0]
gte5 = np.where(temps >= 5)[0]
gt0lt5 = np.where((temps > 0) & (temps < 5))[0]

...

...

...

$dRH/dT$ for $T \leq 0^{\circ}C$





In [4]:



    


rhMed_SRt_lte0 = npTypeZn['binRH_median_SRt'][lte0]
finiteIx = np.isfinite(rhMed_SRt_lte0)
pf_rhMed_SRt_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],rhMed_SRt_lte0[finiteIx],  1))
rhMed_SRt_lte0_dT = pf_rhMed_SRt_lte0(-5) - pf_rhMed_SRt_lte0(-6) 
print(rhMed_SRt_lte0_dT)

rhMed_ARt_lte0 = npTypeZn['binRH_median_ARt'][lte0]
finiteIx = np.isfinite(rhMed_ARt_lte0)
pf_rhMed_ARt_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],rhMed_ARt_lte0[finiteIx],  1))
rhMed_ARt_lte0_dT = pf_rhMed_ARt_lte0(-5) - pf_rhMed_ARt_lte0(-6)     
print(rhMed_ARt_lte0_dT)

rhMed_SRp_lte0 = npTypeZn['binRH_median_SRp'][lte0]
finiteIx = np.isfinite(rhMed_SRp_lte0)
pf_rhMed_SRp_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],rhMed_SRp_lte0[finiteIx],  1))
rhMed_SRp_lte0_dT = pf_rhMed_SRp_lte0(-5) - pf_rhMed_SRp_lte0(-6) 
print(rhMed_SRp_lte0_dT)

rhMed_ARp_lte0 = npTypeZn['binRH_median_ARp'][lte0]
finiteIx = np.isfinite(rhMed_ARp_lte0)
pf_rhMed_ARp_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],rhMed_ARp_lte0[finiteIx],  1))
rhMed_ARp_lte0_dT = pf_rhMed_ARp_lte0(-5) - pf_rhMed_ARp_lte0(-6)     
print(rhMed_ARp_lte0_dT)


rh_srCs1_lte0 = srCs1['rh'][srCs1['tempC'] <= 0]
finiteIx_srCs1 = np.isfinite(rh_srCs1_lte0)
pf_rh_srCs1_lte0 = np.poly1d(np.polyfit(srCs1['tempC'][srCs1['tempC'] <= 0][finiteIx_srCs1],rh_srCs1_lte0[finiteIx_srCs1],  1))
rh_srCs1_lte0_dT = pf_rh_srCs1_lte0(-5) - pf_rh_srCs1_lte0(-6)
print('\n{}'.format(rh_srCs1_lte0_dT))

rh_pfs1_lte0 = pfs1['rh'][pfs1['tempC'] <= 0]
finiteIx_pfs1 = np.isfinite(rh_pfs1_lte0)
pf_rh_pfs1_lte0 = np.poly1d(np.polyfit(pfs1['tempC'][pfs1['tempC'] <= 0][finiteIx_pfs1],rh_pfs1_lte0[finiteIx_pfs1],  1))
rh_pfs1_lte0_dT = pf_rh_pfs1_lte0(-5) - pf_rh_pfs1_lte0(-6)
print('{}'.format(rh_pfs1_lte0_dT))

rh_pfs2_lte0 = pfs2['rh'][pfs2['tempC'] <= 0]
finiteIx_pfs2 = np.isfinite(rh_pfs2_lte0)
pf_rh_pfs2_lte0 = np.poly1d(np.polyfit(pfs2['tempC'][pfs2['tempC'] <= 0][finiteIx_pfs2],rh_pfs2_lte0[finiteIx_pfs2],  1))
rh_pfs2_lte0_dT = pf_rh_pfs2_lte0(-5) - pf_rh_pfs2_lte0(-6)
print('{}'.format(rh_pfs2_lte0_dT))



    


















    






-1.18167057924704
-2.5616459818238013
-0.9247557392883152
-0.15614178260902634

-1.0968525526386372
-1.5387395664021568
-1.7387309162027833

$dRH/dT$ for $T \geq 5^{\circ}C$





In [5]:



    


rhMed_SRt_gte5 = npTypeZn['binRH_median_SRt'][gte5]
finiteIx = np.isfinite(rhMed_SRt_gte5)
pf_rhMed_SRt_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],rhMed_SRt_gte5[finiteIx],  1))
rhMed_SRt_gte5_dT = pf_rhMed_SRt_gte5(6) - pf_rhMed_SRt_gte5(5) 
print(rhMed_SRt_gte5_dT)

rhMed_ARt_gte5 = npTypeZn['binRH_median_ARt'][gte5]
finiteIx = np.isfinite(rhMed_ARt_gte5)
pf_rhMed_ARt_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],rhMed_ARt_gte5[finiteIx],  1))
rhMed_ARt_gte5_dT = pf_rhMed_ARt_gte5(6) - pf_rhMed_ARt_gte5(5)     
print(rhMed_ARt_gte5_dT)

rhMed_SRp_gte5 = npTypeZn['binRH_median_SRp'][gte5]
finiteIx = np.isfinite(rhMed_SRp_gte5)
pf_rhMed_SRp_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],rhMed_SRp_gte5[finiteIx],  1))
rhMed_SRp_gte5_dT = pf_rhMed_SRp_gte5(6) - pf_rhMed_SRp_gte5(5) 
print(rhMed_SRp_gte5_dT)

rhMed_ARp_gte5 = npTypeZn['binRH_median_ARp'][gte5]
finiteIx = np.isfinite(rhMed_ARp_gte5)
pf_rhMed_ARp_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],rhMed_ARp_gte5[finiteIx],  1))
rhMed_ARp_gte5_dT = pf_rhMed_ARp_gte5(6) - pf_rhMed_ARp_gte5(5)     
print(rhMed_ARp_gte5_dT)


rh_srCs1_gte5 = srCs1['rh'][srCs1['tempC'] >= 5]
finiteIx_srCs1 = np.isfinite(rh_srCs1_gte5)
pf_rh_srCs1_gte5 = np.poly1d(np.polyfit(srCs1['tempC'][srCs1['tempC'] >= 5][finiteIx_srCs1],rh_srCs1_gte5[finiteIx_srCs1],  1))
rh_srCs1_gte5_dT = pf_rh_srCs1_gte5(6) - pf_rh_srCs1_gte5(5)
print('\n{}'.format(rh_srCs1_gte5_dT))

rh_pfs2_gte5 = pfs2['rh'][pfs2['tempC'] >= 5]
finiteIx_pfs2 = np.isfinite(rh_pfs2_gte5)
pf_rh_pfs2_gte5 = np.poly1d(np.polyfit(pfs2['tempC'][pfs2['tempC'] >= 5][finiteIx_pfs2],rh_pfs2_gte5[finiteIx_pfs2],  1))
rh_pfs2_gte5_dT = pf_rh_pfs2_gte5(6) - pf_rh_pfs2_gte5(5)
print('{}'.format(rh_pfs2_gte5_dT))



    


















    






-3.018034724246462
-0.8623992895150892
-2.5503460219079557
-2.970400886023114

-1.404971655512469
-0.2641469462220698

$d𝒜/dT$ for $T \leq 0^{\circ}C$





In [6]:



    


ARatMed_SRt_lte0 = npTypeZn['binARat_median_SRt'][lte0]
finiteIx = np.isfinite(ARatMed_SRt_lte0)
pf_ARatMed_SRt_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],ARatMed_SRt_lte0[finiteIx],  1))
ARatMed_SRt_lte0_dT = pf_ARatMed_SRt_lte0(-5) - pf_ARatMed_SRt_lte0(-6) 
print(ARatMed_SRt_lte0_dT)

ARatMed_ARt_lte0 = npTypeZn['binARat_median_ARt'][lte0]
finiteIx = np.isfinite(ARatMed_ARt_lte0)
pf_ARatMed_ARt_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],ARatMed_ARt_lte0[finiteIx],  1))
ARatMed_ARt_lte0_dT = pf_ARatMed_ARt_lte0(-5) - pf_ARatMed_ARt_lte0(-6)     
print(ARatMed_ARt_lte0_dT)

ARatMed_SRp_lte0 = npTypeZn['binARat_median_SRp'][lte0]
finiteIx = np.isfinite(ARatMed_SRp_lte0)
pf_ARatMed_SRp_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],ARatMed_SRp_lte0[finiteIx],  1))
ARatMed_SRp_lte0_dT = pf_ARatMed_SRp_lte0(-5) - pf_ARatMed_SRp_lte0(-6) 
print(ARatMed_SRp_lte0_dT)

ARatMed_ARp_lte0 = npTypeZn['binARat_median_ARp'][lte0]
finiteIx = np.isfinite(ARatMed_ARp_lte0)
pf_ARatMed_ARp_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],ARatMed_ARp_lte0[finiteIx],  1))
ARatMed_ARp_lte0_dT = pf_ARatMed_ARp_lte0(-5) - pf_ARatMed_ARp_lte0(-6)     
print(ARatMed_ARp_lte0_dT)


ARat_srCs1_lte0 = srCs1['ar'][srCs1['tempC'] <= 0]
finiteIx_srCs1 = np.isfinite(ARat_srCs1_lte0)
pf_ARat_srCs1_lte0 = np.poly1d(np.polyfit(srCs1['tempC'][srCs1['tempC'] <= 0][finiteIx_srCs1],ARat_srCs1_lte0[finiteIx_srCs1],  1))
ARat_srCs1_lte0_dT = pf_ARat_srCs1_lte0(-5) - pf_ARat_srCs1_lte0(-6)
print('\n{}'.format(ARat_srCs1_lte0_dT))

ARat_pfs1_lte0 = pfs1['ar'][pfs1['tempC'] <= 0]
finiteIx_pfs1 = np.isfinite(ARat_pfs1_lte0)
pf_ARat_pfs1_lte0 = np.poly1d(np.polyfit(pfs1['tempC'][pfs1['tempC'] <= 0][finiteIx_pfs1],ARat_pfs1_lte0[finiteIx_pfs1],  1))
ARat_pfs1_lte0_dT = pf_ARat_pfs1_lte0(-5) - pf_ARat_pfs1_lte0(-6)
print('{}'.format(ARat_pfs1_lte0_dT))

ARat_pfs2_lte0 = pfs2['ar'][pfs2['tempC'] <= 0]
finiteIx_pfs2 = np.isfinite(ARat_pfs2_lte0)
pf_ARat_pfs2_lte0 = np.poly1d(np.polyfit(pfs2['tempC'][pfs2['tempC'] <= 0][finiteIx_pfs2],ARat_pfs2_lte0[finiteIx_pfs2],  1))
ARat_pfs2_lte0_dT = pf_ARat_pfs2_lte0(-5) - pf_ARat_pfs2_lte0(-6)
print('{}'.format(ARat_pfs2_lte0_dT))



    


















    






-0.001587942923338137
0.004274667218122197
-0.003437811263710122
0.0037373435510775765

-0.0032427929360022145
-0.01720816122992541
0.0008080408010409657

$d𝒜/dT$ for $T \geq 5^{\circ}C$





In [7]:



    


ARatMed_SRt_gte5 = npTypeZn['binARat_median_SRt'][gte5]
finiteIx = np.isfinite(ARatMed_SRt_gte5)
pf_ARatMed_SRt_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],ARatMed_SRt_gte5[finiteIx],  1))
ARatMed_SRt_gte5_dT = pf_ARatMed_SRt_gte5(-5) - pf_ARatMed_SRt_gte5(-6) 
print(ARatMed_SRt_gte5_dT)

ARatMed_ARt_gte5 = npTypeZn['binARat_median_ARt'][gte5]
finiteIx = np.isfinite(ARatMed_ARt_gte5)
pf_ARatMed_ARt_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],ARatMed_ARt_gte5[finiteIx],  1))
ARatMed_ARt_gte5_dT = pf_ARatMed_ARt_gte5(-5) - pf_ARatMed_ARt_gte5(-6)     
print(ARatMed_ARt_gte5_dT)

ARatMed_SRp_gte5 = npTypeZn['binARat_median_SRp'][gte5]
finiteIx = np.isfinite(ARatMed_SRp_gte5)
pf_ARatMed_SRp_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],ARatMed_SRp_gte5[finiteIx],  1))
ARatMed_SRp_gte5_dT = pf_ARatMed_SRp_gte5(-5) - pf_ARatMed_SRp_gte5(-6) 
print(ARatMed_SRp_gte5_dT)

ARatMed_ARp_gte5 = npTypeZn['binARat_median_ARp'][gte5]
finiteIx = np.isfinite(ARatMed_ARp_gte5)
pf_ARatMed_ARp_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],ARatMed_ARp_gte5[finiteIx],  1))
ARatMed_ARp_gte5_dT = pf_ARatMed_ARp_gte5(-5) - pf_ARatMed_ARp_gte5(-6)     
print(ARatMed_ARp_gte5_dT)


ARat_srCs1_gte5 = srCs1['ar'][srCs1['tempC'] >= 5]
finiteIx_srCs1 = np.isfinite(ARat_srCs1_gte5)
pf_ARat_srCs1_gte5 = np.poly1d(np.polyfit(srCs1['tempC'][srCs1['tempC'] >= 5][finiteIx_srCs1],ARat_srCs1_gte5[finiteIx_srCs1],  1))
ARat_srCs1_gte5_dT = pf_ARat_srCs1_gte5(-5) - pf_ARat_srCs1_gte5(-6)
print('\n{}'.format(ARat_srCs1_gte5_dT))

ARat_pfs2_gte5 = pfs2['ar'][pfs2['tempC'] >= 5]
finiteIx_pfs2 = np.isfinite(ARat_pfs2_gte5)
pf_ARat_pfs2_gte5 = np.poly1d(np.polyfit(pfs2['tempC'][pfs2['tempC'] >= 5][finiteIx_pfs2],ARat_pfs2_gte5[finiteIx_pfs2],  1))
ARat_pfs2_gte5_dT = pf_ARat_pfs2_gte5(-5) - pf_ARat_pfs2_gte5(-6)
print('{}'.format(ARat_pfs2_gte5_dT))



    


















    






-0.00048045685671660454
-0.005911044800520404
0.0007528008198094405
-0.004135782685875555

-0.0010492083743433023
0.006686653702242973

$dD_{mm}/dT$ for $T \leq 0^{\circ}C$





In [8]:



    


DmmMed_SRt_lte0 = npTypeZn['binDmm_median_SRt'][lte0]
finiteIx = np.isfinite(DmmMed_SRt_lte0)
pf_DmmMed_SRt_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],DmmMed_SRt_lte0[finiteIx],  1))
DmmMed_SRt_lte0_dT = pf_DmmMed_SRt_lte0(-5) - pf_DmmMed_SRt_lte0(-6) 
print(DmmMed_SRt_lte0_dT)

DmmMed_ARt_lte0 = npTypeZn['binDmm_median_ARt'][lte0]
finiteIx = np.isfinite(DmmMed_ARt_lte0)
pf_DmmMed_ARt_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],DmmMed_ARt_lte0[finiteIx],  1))
DmmMed_ARt_lte0_dT = pf_DmmMed_ARt_lte0(-5) - pf_DmmMed_ARt_lte0(-6)     
print(DmmMed_ARt_lte0_dT)

DmmMed_SRp_lte0 = npTypeZn['binDmm_median_SRp'][lte0]
finiteIx = np.isfinite(DmmMed_SRp_lte0)
pf_DmmMed_SRp_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],DmmMed_SRp_lte0[finiteIx],  1))
DmmMed_SRp_lte0_dT = pf_DmmMed_SRp_lte0(-5) - pf_DmmMed_SRp_lte0(-6) 
print(DmmMed_SRp_lte0_dT)

DmmMed_ARp_lte0 = npTypeZn['binDmm_median_ARp'][lte0]
finiteIx = np.isfinite(DmmMed_ARp_lte0)
pf_DmmMed_ARp_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],DmmMed_ARp_lte0[finiteIx],  1))
DmmMed_ARp_lte0_dT = pf_DmmMed_ARp_lte0(-5) - pf_DmmMed_ARp_lte0(-6)     
print(DmmMed_ARp_lte0_dT)


Dmm_srCs1_lte0 = srCs1['Dmm'][srCs1['tempC'] <= 0]
finiteIx_srCs1 = np.isfinite(Dmm_srCs1_lte0)
pf_Dmm_srCs1_lte0 = np.poly1d(np.polyfit(srCs1['tempC'][srCs1['tempC'] <= 0][finiteIx_srCs1],Dmm_srCs1_lte0[finiteIx_srCs1],  1))
Dmm_srCs1_lte0_dT = pf_Dmm_srCs1_lte0(-5) - pf_Dmm_srCs1_lte0(-6)
print('\n{}'.format(Dmm_srCs1_lte0_dT))

Dmm_pfs1_lte0 = pfs1['Dmm'][pfs1['tempC'] <= 0]
finiteIx_pfs1 = np.isfinite(Dmm_pfs1_lte0)
pf_Dmm_pfs1_lte0 = np.poly1d(np.polyfit(pfs1['tempC'][pfs1['tempC'] <= 0][finiteIx_pfs1],Dmm_pfs1_lte0[finiteIx_pfs1],  1))
Dmm_pfs1_lte0_dT = pf_Dmm_pfs1_lte0(-5) - pf_Dmm_pfs1_lte0(-6)
print('{}'.format(Dmm_pfs1_lte0_dT))

Dmm_pfs2_lte0 = pfs2['Dmm'][pfs2['tempC'] <= 0]
finiteIx_pfs2 = np.isfinite(Dmm_pfs2_lte0)
pf_Dmm_pfs2_lte0 = np.poly1d(np.polyfit(pfs2['tempC'][pfs2['tempC'] <= 0][finiteIx_pfs2],Dmm_pfs2_lte0[finiteIx_pfs2],  1))
Dmm_pfs2_lte0_dT = pf_Dmm_pfs2_lte0(-5) - pf_Dmm_pfs2_lte0(-6)
print('{}'.format(Dmm_pfs2_lte0_dT))



    


















    






0.02645123839009278
0.004151348039215685
0.033522727272727426
-0.03279532967032961

0.028867281296104164
0.05367702263745433
-0.06788525356751629

$dD_{mm}/dT$ for $T \geq 5^{\circ}C$





In [9]:



    


DmmMed_SRt_gte5 = npTypeZn['binDmm_median_SRt'][gte5]
finiteIx = np.isfinite(DmmMed_SRt_gte5)
pf_DmmMed_SRt_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],DmmMed_SRt_gte5[finiteIx],  1))
DmmMed_SRt_gte5_dT = pf_DmmMed_SRt_gte5(-5) - pf_DmmMed_SRt_gte5(-6) 
print(DmmMed_SRt_gte5_dT)

DmmMed_ARt_gte5 = npTypeZn['binDmm_median_ARt'][gte5]
finiteIx = np.isfinite(DmmMed_ARt_gte5)
pf_DmmMed_ARt_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],DmmMed_ARt_gte5[finiteIx],  1))
DmmMed_ARt_gte5_dT = pf_DmmMed_ARt_gte5(-5) - pf_DmmMed_ARt_gte5(-6)     
print(DmmMed_ARt_gte5_dT)

DmmMed_SRp_gte5 = npTypeZn['binDmm_median_SRp'][gte5]
finiteIx = np.isfinite(DmmMed_SRp_gte5)
pf_DmmMed_SRp_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],DmmMed_SRp_gte5[finiteIx],  1))
DmmMed_SRp_gte5_dT = pf_DmmMed_SRp_gte5(-5) - pf_DmmMed_SRp_gte5(-6) 
print(DmmMed_SRp_gte5_dT)

DmmMed_ARp_gte5 = npTypeZn['binDmm_median_ARp'][gte5]
finiteIx = np.isfinite(DmmMed_ARp_gte5)
pf_DmmMed_ARp_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],DmmMed_ARp_gte5[finiteIx],  1))
DmmMed_ARp_gte5_dT = pf_DmmMed_ARp_gte5(-5) - pf_DmmMed_ARp_gte5(-6)     
print(DmmMed_ARp_gte5_dT)


Dmm_srCs1_gte5 = srCs1['Dmm'][srCs1['tempC'] >= 5]
finiteIx_srCs1 = np.isfinite(Dmm_srCs1_gte5)
pf_Dmm_srCs1_gte5 = np.poly1d(np.polyfit(srCs1['tempC'][srCs1['tempC'] >= 5][finiteIx_srCs1],Dmm_srCs1_gte5[finiteIx_srCs1],  1))
Dmm_srCs1_gte5_dT = pf_Dmm_srCs1_gte5(-5) - pf_Dmm_srCs1_gte5(-6)
print('\n{}'.format(Dmm_srCs1_gte5_dT))

Dmm_pfs2_gte5 = pfs2['Dmm'][pfs2['tempC'] >= 5]
finiteIx_pfs2 = np.isfinite(Dmm_pfs2_gte5)
pf_Dmm_pfs2_gte5 = np.poly1d(np.polyfit(pfs2['tempC'][pfs2['tempC'] >= 5][finiteIx_pfs2],Dmm_pfs2_gte5[finiteIx_pfs2],  1))
Dmm_pfs2_gte5_dT = pf_Dmm_pfs2_gte5(-5) - pf_Dmm_pfs2_gte5(-6)
print('{}'.format(Dmm_pfs2_gte5_dT))



    


















    






0.005892857142857144
0.054469273743016744
0.002611607142857242
0.027397297891036898

-0.017299089699716497
0.020716517832148507

$dN_t/dT$ for $T \leq 0^{\circ}C$





In [10]:



    


NtMed_SRt_lte0 = npTypeZn['binNt_median_SRt'][lte0]
finiteIx = np.isfinite(np.log10(NtMed_SRt_lte0))
pf_NtMed_SRt_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(NtMed_SRt_lte0[finiteIx]),  1))
NtMed_SRt_lte0_dT = 1 - (10**pf_NtMed_SRt_lte0(-6)/10**pf_NtMed_SRt_lte0(-5)) # 1 deg C change percentage 
print(NtMed_SRt_lte0_dT)

NtMed_ARt_lte0 = npTypeZn['binNt_median_ARt'][lte0]
finiteIx = np.isfinite(np.log10(NtMed_ARt_lte0))
pf_NtMed_ARt_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(NtMed_ARt_lte0[finiteIx]),  1))
NtMed_ARt_lte0_dT = 1 - (10**pf_NtMed_ARt_lte0(-6)/10**pf_NtMed_ARt_lte0(-5)) # 1 deg C change percentage 
print(NtMed_ARt_lte0_dT)

NtMed_SRp_lte0 = npTypeZn['binNt_median_SRp'][lte0]
finiteIx = np.isfinite(np.log10(NtMed_SRp_lte0))
pf_NtMed_SRp_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(NtMed_SRp_lte0[finiteIx]),  1))
NtMed_SRp_lte0_dT = 1 - (10**pf_NtMed_SRp_lte0(-6)/10**pf_NtMed_SRp_lte0(-5)) # 1 deg C change percentage 
print(NtMed_SRp_lte0_dT)

NtMed_ARp_lte0 = npTypeZn['binNt_median_ARp'][lte0]
finiteIx = np.isfinite(np.log10(NtMed_ARp_lte0))
pf_NtMed_ARp_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(NtMed_ARp_lte0[finiteIx]),  1))
NtMed_ARp_lte0_dT = 1 - (10**pf_NtMed_ARp_lte0(-6)/10**pf_NtMed_ARp_lte0(-5)) # 1 deg C change percentage 
print(NtMed_ARp_lte0_dT)



Nt_srCs1_lte0 = srCs1['Nt'][srCs1['tempC'] <= 0]
finiteIx_srCs1 = np.isfinite(np.log10(Nt_srCs1_lte0))
pf_Nt_srCs1_lte0 = np.poly1d(np.polyfit(srCs1['tempC'][srCs1['tempC'] <= 0][finiteIx_srCs1],np.log10(Nt_srCs1_lte0[finiteIx_srCs1]),  1))
Nt_srCs1_lte0_dT = 1 - (10**pf_Nt_srCs1_lte0(-6)/10**pf_Nt_srCs1_lte0(-5))
print('\n{}'.format(Nt_srCs1_lte0_dT))

Nt_pfs1_lte0 = pfs1['Nt'][pfs1['tempC'] <= 0]
finiteIx_pfs1 = np.isfinite(np.log10(Nt_pfs1_lte0))
pf_Nt_pfs1_lte0 = np.poly1d(np.polyfit(pfs1['tempC'][pfs1['tempC'] <= 0][finiteIx_pfs1],np.log10(Nt_pfs1_lte0[finiteIx_pfs1]),  1))
Nt_pfs1_lte0_dT = 1 - (10**pf_Nt_pfs1_lte0(-6)/10**pf_Nt_pfs1_lte0(-5))
print('{}'.format(Nt_pfs1_lte0_dT))

Nt_pfs2_lte0 = pfs2['Nt'][pfs2['tempC'] <= 0]
finiteIx_pfs2 = np.isfinite(np.log10(Nt_pfs2_lte0))
pf_Nt_pfs2_lte0 = np.poly1d(np.polyfit(pfs2['tempC'][pfs2['tempC'] <= 0][finiteIx_pfs2],np.log10(Nt_pfs2_lte0[finiteIx_pfs2]),  1))
Nt_pfs2_lte0_dT = 1 - (10**pf_Nt_pfs2_lte0(-6)/10**pf_Nt_pfs2_lte0(-5))
print('{}'.format(Nt_pfs2_lte0_dT))



    


















    






-0.1954092655529147
-0.34557856874573867
-0.19645235544460693
-0.485966405122366

-0.1886079458737584
0.3505182543832218
0.17926551740473973

$dN_t/dT$ for $T \geq 5^{\circ}C$





In [11]:



    


NtMed_SRt_gte5 = npTypeZn['binNt_median_SRt'][gte5]
finiteIx = np.isfinite(np.log10(NtMed_SRt_gte5))
pf_NtMed_SRt_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(NtMed_SRt_gte5[finiteIx]),  1))
NtMed_SRt_gte5_dT = 1 - (10**pf_NtMed_SRt_gte5(2)/10**pf_NtMed_SRt_gte5(3)) # 1 deg C change percentage 
print(NtMed_SRt_gte5_dT)

NtMed_ARt_gte5 = npTypeZn['binNt_median_ARt'][gte5]
finiteIx = np.isfinite(np.log10(NtMed_ARt_gte5))
pf_NtMed_ARt_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(NtMed_ARt_gte5[finiteIx]),  1))
NtMed_ARt_gte5_dT = 1 - (10**pf_NtMed_ARt_gte5(2)/10**pf_NtMed_ARt_gte5(3)) # 1 deg C change percentage 
print(NtMed_ARt_gte5_dT)

NtMed_SRp_gte5 = npTypeZn['binNt_median_SRp'][gte5]
finiteIx = np.isfinite(np.log10(NtMed_SRp_gte5))
pf_NtMed_SRp_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(NtMed_SRp_gte5[finiteIx]),  1))
NtMed_SRp_gte5_dT = 1 - (10**pf_NtMed_SRp_gte5(2)/10**pf_NtMed_SRp_gte5(3)) # 1 deg C change percentage 
print(NtMed_SRp_gte5_dT)

NtMed_ARp_gte5 = npTypeZn['binNt_median_ARp'][gte5]
finiteIx = np.isfinite(np.log10(NtMed_ARp_gte5))
pf_NtMed_ARp_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(NtMed_ARp_gte5[finiteIx]),  1))
NtMed_ARp_gte5_dT = 1 - (10**pf_NtMed_ARp_gte5(2)/10**pf_NtMed_ARp_gte5(3)) # 1 deg C change percentage 
print(NtMed_ARp_gte5_dT)



Nt_srCs1_gte5 = srCs1['Nt'][srCs1['tempC'] >= 5]
finiteIx_srCs1 = np.isfinite(np.log10(Nt_srCs1_gte5))
pf_Nt_srCs1_gte5 = np.poly1d(np.polyfit(srCs1['tempC'][srCs1['tempC'] >= 5][finiteIx_srCs1],np.log10(Nt_srCs1_gte5[finiteIx_srCs1]),  1))
Nt_srCs1_gte5_dT = 1 - (10**pf_Nt_srCs1_gte5(2)/10**pf_Nt_srCs1_gte5(3))
print('\n{}'.format(Nt_srCs1_gte5_dT))

Nt_pfs2_gte5 = pfs2['Nt'][pfs2['tempC'] >= 5]
finiteIx_pfs2 = np.isfinite(np.log10(Nt_pfs2_gte5))
pf_Nt_pfs2_gte5 = np.poly1d(np.polyfit(pfs2['tempC'][pfs2['tempC'] >= 5][finiteIx_pfs2],np.log10(Nt_pfs2_gte5[finiteIx_pfs2]),  1))
Nt_pfs2_gte5_dT = 1 - (10**pf_Nt_pfs2_gte5(2)/10**pf_Nt_pfs2_gte5(3))
print('{}'.format(Nt_pfs2_gte5_dT))



    


















    






-0.1521055773193618
0.1992797946465693
-0.0235674028416768
-0.025856268971608642

-0.062407688092258606
-0.06551822285913622

$dTWC/dT$ for $T \leq 0^{\circ}C$





In [12]:



    


twcMed_SRt_lte0 = npTypeZn['binTWC_median_SRt'][lte0]
finiteIx = np.isfinite(np.log10(twcMed_SRt_lte0))
pf_twcMed_SRt_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(twcMed_SRt_lte0[finiteIx]),  1))
twcMed_SRt_lte0_dT = 1 - (10**pf_twcMed_SRt_lte0(-6)/10**pf_twcMed_SRt_lte0(-5)) # 1 deg C change percentage 
print(twcMed_SRt_lte0_dT)

twcMed_ARt_lte0 = npTypeZn['binTWC_median_ARt'][lte0]
finiteIx = np.isfinite(np.log10(twcMed_ARt_lte0))
pf_twcMed_ARt_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(twcMed_ARt_lte0[finiteIx]),  1))
twcMed_ARt_lte0_dT = 1 - (10**pf_twcMed_ARt_lte0(-6)/10**pf_twcMed_ARt_lte0(-5)) # 1 deg C change percentage 
print(twcMed_ARt_lte0_dT)

twcMed_SRp_lte0 = npTypeZn['binTWC_median_SRp'][lte0]
finiteIx = np.isfinite(np.log10(twcMed_SRp_lte0))
pf_twcMed_SRp_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(twcMed_SRp_lte0[finiteIx]),  1))
twcMed_SRp_lte0_dT = 1 - (10**pf_twcMed_SRp_lte0(-6)/10**pf_twcMed_SRp_lte0(-5)) # 1 deg C change percentage 
print(twcMed_SRp_lte0_dT)

twcMed_ARp_lte0 = npTypeZn['binTWC_median_ARp'][lte0]
finiteIx = np.isfinite(np.log10(twcMed_ARp_lte0))
pf_twcMed_ARp_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(twcMed_ARp_lte0[finiteIx]),  1))
twcMed_ARp_lte0_dT = 1 - (10**pf_twcMed_ARp_lte0(-6)/10**pf_twcMed_ARp_lte0(-5)) # 1 deg C change percentage 
print(twcMed_ARp_lte0_dT)



TWC_srCs1_lte0 = srCs1['twc'][srCs1['tempC'] <= 0]
finiteIx_srCs1 = np.isfinite(np.log10(TWC_srCs1_lte0))
pf_TWC_srCs1_lte0 = np.poly1d(np.polyfit(srCs1['tempC'][srCs1['tempC'] <= 0][finiteIx_srCs1],np.log10(TWC_srCs1_lte0[finiteIx_srCs1]),  1))
TWC_srCs1_lte0_dT = 1 - (10**pf_TWC_srCs1_lte0(-6)/10**pf_TWC_srCs1_lte0(-5))
print('\n{}'.format(TWC_srCs1_lte0_dT))

TWC_pfs1_lte0 = pfs1['twc'][pfs1['tempC'] <= 0]
finiteIx_pfs1 = np.isfinite(np.log10(TWC_pfs1_lte0))
pf_TWC_pfs1_lte0 = np.poly1d(np.polyfit(pfs1['tempC'][pfs1['tempC'] <= 0][finiteIx_pfs1],np.log10(TWC_pfs1_lte0[finiteIx_pfs1]),  1))
TWC_pfs1_lte0_dT = 1 - (10**pf_TWC_pfs1_lte0(-6)/10**pf_TWC_pfs1_lte0(-5))
print('{}'.format(TWC_pfs1_lte0_dT))

TWC_pfs2_lte0 = pfs2['twc'][pfs2['tempC'] <= 0]
finiteIx_pfs2 = np.isfinite(np.log10(TWC_pfs2_lte0))
pf_TWC_pfs2_lte0 = np.poly1d(np.polyfit(pfs2['tempC'][pfs2['tempC'] <= 0][finiteIx_pfs2],np.log10(TWC_pfs2_lte0[finiteIx_pfs2]),  1))
TWC_pfs2_lte0_dT = 1 - (10**pf_TWC_pfs2_lte0(-6)/10**pf_TWC_pfs2_lte0(-5))
print('{}'.format(TWC_pfs2_lte0_dT))



    


















    






-0.12906899245236514
-0.30914564298487535
-0.10913909680971368
-0.6064912918860832

-0.10012271479656398
0.3868639290341105
0.06244524378263461

$dTWC/dT$ for $T \geq 5^{\circ}C$





In [13]:



    


twcMed_SRt_gte5 = npTypeZn['binTWC_median_SRt'][gte5]
finiteIx = np.isfinite(np.log10(twcMed_SRt_gte5))
pf_twcMed_SRt_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(twcMed_SRt_gte5[finiteIx]),  1))
twcMed_SRt_gte5_dT = 1 - (10**pf_twcMed_SRt_gte5(2)/10**pf_twcMed_SRt_gte5(3)) # 1 deg C change percentage 
print(twcMed_SRt_gte5_dT)

twcMed_ARt_gte5 = npTypeZn['binTWC_median_ARt'][gte5]
finiteIx = np.isfinite(np.log10(twcMed_ARt_gte5))
pf_twcMed_ARt_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(twcMed_ARt_gte5[finiteIx]),  1))
twcMed_ARt_gte5_dT = 1 - (10**pf_twcMed_ARt_gte5(2)/10**pf_twcMed_ARt_gte5(3)) # 1 deg C change percentage 
print(twcMed_ARt_gte5_dT)

twcMed_SRp_gte5 = npTypeZn['binTWC_median_SRp'][gte5]
finiteIx = np.isfinite(np.log10(twcMed_SRp_gte5))
pf_twcMed_SRp_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(twcMed_SRp_gte5[finiteIx]),  1))
twcMed_SRp_gte5_dT = 1 - (10**pf_twcMed_SRp_gte5(2)/10**pf_twcMed_SRp_gte5(3)) # 1 deg C change percentage 
print(twcMed_SRp_gte5_dT)

twcMed_ARp_gte5 = npTypeZn['binTWC_median_ARp'][gte5]
finiteIx = np.isfinite(np.log10(twcMed_ARp_gte5))
pf_twcMed_ARp_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(twcMed_ARp_gte5[finiteIx]),  1))
twcMed_ARp_gte5_dT = 1 - (10**pf_twcMed_ARp_gte5(2)/10**pf_twcMed_ARp_gte5(3)) # 1 deg C change percentage 
print(twcMed_ARp_gte5_dT)



TWC_srCs1_gte5 = srCs1['twc'][srCs1['tempC'] >= 5]
finiteIx_srCs1 = np.isfinite(np.log10(TWC_srCs1_gte5))
pf_TWC_srCs1_gte5 = np.poly1d(np.polyfit(srCs1['tempC'][srCs1['tempC'] >= 5][finiteIx_srCs1],np.log10(TWC_srCs1_gte5[finiteIx_srCs1]),  1))
TWC_srCs1_gte5_dT = 1 - (10**pf_TWC_srCs1_gte5(2)/10**pf_TWC_srCs1_gte5(3))
print('\n{}'.format(TWC_srCs1_gte5_dT))

TWC_pfs2_gte5 = pfs2['twc'][pfs2['tempC'] >= 5]
finiteIx_pfs2 = np.isfinite(np.log10(TWC_pfs2_gte5))
pf_TWC_pfs2_gte5 = np.poly1d(np.polyfit(pfs2['tempC'][pfs2['tempC'] >= 5][finiteIx_pfs2],np.log10(TWC_pfs2_gte5[finiteIx_pfs2]),  1))
TWC_pfs2_gte5_dT = 1 - (10**pf_TWC_pfs2_gte5(2)/10**pf_TWC_pfs2_gte5(3))
print('{}'.format(TWC_pfs2_gte5_dT))



    


















    






-0.10504387438026885
0.36344029900128816
0.02615622898291725
0.06961735706367889

-0.06008546762710942
-0.021612761024605298

...

...

...

Average of $RH$ Median Values for $T \leq 0^{\circ}C$





In [14]:



    


rhMed_SRt_lte0 = npTypeZn['binRH_median_SRt'][lte0]
rhMedAvg_SRt_lte0 = np.nanmean(rhMed_SRt_lte0)
print('{}'.format(rhMedAvg_SRt_lte0))

rhMed_ARt_lte0 = npTypeZn['binRH_median_ARt'][lte0]
rhMedAvg_ARt_lte0 = np.nanmean(rhMed_ARt_lte0)
print('{}'.format(rhMedAvg_ARt_lte0))

rhMed_SRp_lte0 = npTypeZn['binRH_median_SRp'][lte0]
rhMedAvg_SRp_lte0 = np.nanmean(rhMed_SRp_lte0)
print('{}'.format(rhMedAvg_SRp_lte0))

rhMed_ARp_lte0 = npTypeZn['binRH_median_ARp'][lte0]
rhMedAvg_ARp_lte0 = np.nanmean(rhMed_ARp_lte0)
print('{}'.format(rhMedAvg_ARp_lte0))

print('\n{}'.format(np.nanmean(srCs1['rh'][srCs1['tempC'] <= 0])))
print('{}'.format(np.nanmean(pfs1['rh'][pfs1['tempC'] <= 0])))
print('{}'.format(np.nanmean(pfs2['rh'][pfs2['tempC'] <= 0])))



    


















    






108.7671504493976
93.79847944936004
93.19036683552969
81.00428921659646

105.74826728160639
104.25040758867199
100.81165250084683

Average of $RH$ Median Values for $T \geq 5^{\circ}C$





In [15]:



    


rhMed_SRt_gte5 = npTypeZn['binRH_median_SRt'][gte5]
rhMedAvg_SRt_gte5 = np.nanmean(rhMed_SRt_gte5)
print('{}'.format(rhMedAvg_SRt_gte5))

rhMed_ARt_gte5 = npTypeZn['binRH_median_ARt'][gte5]
rhMedAvg_ARt_gte5 = np.nanmean(rhMed_ARt_gte5)
print('{}'.format(rhMedAvg_ARt_gte5))

rhMed_SRp_gte5 = npTypeZn['binRH_median_SRp'][gte5]
rhMedAvg_SRp_gte5 = np.nanmean(rhMed_SRp_gte5)
print('{}'.format(rhMedAvg_SRp_gte5))

rhMed_ARp_gte5 = npTypeZn['binRH_median_ARp'][gte5]
rhMedAvg_ARp_gte5 = np.nanmean(rhMed_ARp_gte5)
print('{}'.format(rhMedAvg_ARp_gte5))


print('\n{}'.format(np.nanmean(srCs1['rh'][srCs1['tempC'] >= 5])))
print('{}'.format(np.nanmean(pfs2['rh'][pfs2['tempC'] >= 5])))



    


















    






67.40310741377964
49.61811222785604
55.53799957744915
70.62118295146568

86.32477512414326
73.91955286562525

Average of $𝒜$ Median Values for $T \leq 0^{\circ}C$





In [16]:



    


ARatMed_SRt_lte0 = npTypeZn['binARat_median_SRt'][lte0]
ARatMedAvg_SRt_lte0 = np.nanmean(ARatMed_SRt_lte0)
print('{}'.format(ARatMedAvg_SRt_lte0))

ARatMed_ARt_lte0 = npTypeZn['binARat_median_ARt'][lte0]
ARatMedAvg_ARt_lte0 = np.nanmean(ARatMed_ARt_lte0)
print('{}'.format(ARatMedAvg_ARt_lte0))

ARatMed_SRp_lte0 = npTypeZn['binARat_median_SRp'][lte0]
ARatMedAvg_SRp_lte0 = np.nanmean(ARatMed_SRp_lte0)
print('{}'.format(ARatMedAvg_SRp_lte0))

ARatMed_ARp_lte0 = npTypeZn['binARat_median_ARp'][lte0]
ARatMedAvg_ARp_lte0 = np.nanmean(ARatMed_ARp_lte0)
print('{}'.format(ARatMedAvg_ARp_lte0))


print('\n{}'.format(np.nanmean(srCs1['ar'][srCs1['tempC'] <= 0])))
print('{}'.format(np.nanmean(pfs1['ar'][pfs1['tempC'] <= 0])))
print('{}'.format(np.nanmean(pfs2['ar'][pfs2['tempC'] <= 0])))



    


















    






0.45835055164812705
0.4852672636480371
0.428212426327026
0.5197895773833812

0.44799971054405946
0.5004027153536509
0.4968649817362637

Average of $𝒜$ Median Values for $T \geq 5^{\circ}C$





In [17]:



    


ARatMed_SRt_gte5 = npTypeZn['binARat_median_SRt'][gte5]
ARatMedAvg_SRt_gte5 = np.nanmean(ARatMed_SRt_gte5)
print('{}'.format(ARatMedAvg_SRt_gte5))

ARatMed_ARt_gte5 = npTypeZn['binARat_median_ARt'][gte5]
ARatMedAvg_ARt_gte5 = np.nanmean(ARatMed_ARt_gte5)
print('{}'.format(ARatMedAvg_ARt_gte5))

ARatMed_SRp_gte5 = npTypeZn['binARat_median_SRp'][gte5]
ARatMedAvg_SRp_gte5 = np.nanmean(ARatMed_SRp_gte5)
print('{}'.format(ARatMedAvg_SRp_gte5))

ARatMed_ARp_gte5 = npTypeZn['binARat_median_ARp'][gte5]
ARatMedAvg_ARp_gte5 = np.nanmean(ARatMed_ARp_gte5)
print('{}'.format(ARatMedAvg_ARp_gte5))


print('\n{}'.format(np.nanmean(srCs1['ar'][srCs1['tempC'] >= 5])))
print('{}'.format(np.nanmean(pfs2['ar'][pfs2['tempC'] >= 5])))



    


















    






0.7174562125214204
0.7586141443451795
0.7047048649663207
0.7474183274670622

0.7006623227845545
0.7013811249070705

Average of $D_{mm}$ Median Values for $T \leq 0^{\circ}C$





In [18]:



    


DmmMed_SRt_lte0 = npTypeZn['binDmm_median_SRt'][lte0]
DmmMedAvg_SRt_lte0 = np.nanmean(DmmMed_SRt_lte0)
print('{}'.format(DmmMedAvg_SRt_lte0))

DmmMed_ARt_lte0 = npTypeZn['binDmm_median_ARt'][lte0]
DmmMedAvg_ARt_lte0 = np.nanmean(DmmMed_ARt_lte0)
print('{}'.format(DmmMedAvg_ARt_lte0))

DmmMed_SRp_lte0 = npTypeZn['binDmm_median_SRp'][lte0]
DmmMedAvg_SRp_lte0 = np.nanmean(DmmMed_SRp_lte0)
print('{}'.format(DmmMedAvg_SRp_lte0))

DmmMed_ARp_lte0 = npTypeZn['binDmm_median_ARp'][lte0]
DmmMedAvg_ARp_lte0 = np.nanmean(DmmMed_ARp_lte0)
print('{}'.format(DmmMedAvg_ARp_lte0))


print('\n{}'.format(np.nanmean(srCs1['Dmm'][srCs1['tempC'] <= 0])))
print('{}'.format(np.nanmean(pfs1['Dmm'][pfs1['tempC'] <= 0])))
print('{}'.format(np.nanmean(pfs2['Dmm'][pfs2['tempC'] <= 0])))



    


















    






1.132291666666667
1.3161764705882353
1.343181818181818
1.272596153846154

1.1755813953488374
0.8638059701492538
0.9490131578947371

Average of $D_{mm}$ Median Values for $T \geq 5^{\circ}C$





In [19]:



    


DmmMed_SRt_gte5 = npTypeZn['binDmm_median_SRt'][gte5]
DmmMedAvg_SRt_gte5 = np.nanmean(DmmMed_SRt_gte5)
print('{}'.format(DmmMedAvg_SRt_gte5))

DmmMed_ARt_gte5 = npTypeZn['binDmm_median_ARt'][gte5]
DmmMedAvg_ARt_gte5 = np.nanmean(DmmMed_ARt_gte5)
print('{}'.format(DmmMedAvg_ARt_gte5))

DmmMed_SRp_gte5 = npTypeZn['binDmm_median_SRp'][gte5]
DmmMedAvg_SRp_gte5 = np.nanmean(DmmMed_SRp_gte5)
print('{}'.format(DmmMedAvg_SRp_gte5))

DmmMed_ARp_gte5 = npTypeZn['binDmm_median_ARp'][gte5]
DmmMedAvg_ARp_gte5 = np.nanmean(DmmMed_ARp_gte5)
print('{}'.format(DmmMedAvg_ARp_gte5))


print('\n{}'.format(np.nanmean(srCs1['Dmm'][srCs1['tempC'] >= 5])))
print('{}'.format(np.nanmean(pfs2['Dmm'][pfs2['tempC'] >= 5])))



    


















    






1.149166666666667
0.7697916666666667
1.1895833333333334
0.6065104166666666

1.3301724137931035
0.6817708333333335

Average of $N_t$ Median Values for $T \leq 0^{\circ}C$





In [20]:



    


ntMed_SRt_lte0 = npTypeZn['binNt_median_SRt'][lte0]
ntMedAvg_SRt_lte0 = np.nanmean(ntMed_SRt_lte0)
print('{}'.format(ntMedAvg_SRt_lte0))

ntMed_ARt_lte0 = npTypeZn['binNt_median_ARt'][lte0]
ntMedAvg_ARt_lte0 = np.nanmean(ntMed_ARt_lte0)
print('{}'.format(ntMedAvg_ARt_lte0))

ntMed_SRp_lte0 = npTypeZn['binNt_median_SRp'][lte0]
ntMedAvg_SRp_lte0 = np.nanmean(ntMed_SRp_lte0)
print('{}'.format(ntMedAvg_SRp_lte0))

ntMed_ARp_lte0 = npTypeZn['binNt_median_ARp'][lte0]
ntMedAvg_ARp_lte0 = np.nanmean(ntMed_ARp_lte0)
print('{}'.format(ntMedAvg_ARp_lte0))


print('\n{}'.format(np.nanmean(srCs1['Nt'][srCs1['tempC'] <= 0])))
print('{}'.format(np.nanmean(pfs1['Nt'][pfs1['tempC'] <= 0])))
print('{}'.format(np.nanmean(pfs2['Nt'][pfs2['tempC'] <= 0])))



    


















    






0.06258528720591569
0.023096177853716697
0.0058735591003069906
0.004183855613670667

0.028064582846677667
0.011991338084514589
0.009249095830266816

Average of $N_t$ Median Values for $T \geq 5^{\circ}C$





In [21]:



    


ntMed_SRt_gte5 = npTypeZn['binNt_median_SRt'][gte5]
ntMedAvg_SRt_gte5 = np.nanmean(ntMed_SRt_gte5)
print('{}'.format(ntMedAvg_SRt_gte5))

ntMed_ARt_gte5 = npTypeZn['binNt_median_ARt'][gte5]
ntMedAvg_ARt_gte5 = np.nanmean(ntMed_ARt_gte5)
print('{}'.format(ntMedAvg_ARt_gte5))

ntMed_SRp_gte5 = npTypeZn['binNt_median_SRp'][gte5]
ntMedAvg_SRp_gte5 = np.nanmean(ntMed_SRp_gte5)
print('{}'.format(ntMedAvg_SRp_gte5))

ntMed_ARp_gte5 = npTypeZn['binNt_median_ARp'][gte5]
ntMedAvg_ARp_gte5 = np.nanmean(ntMed_ARp_gte5)
print('{}'.format(ntMedAvg_ARp_gte5))


print('\n{}'.format(np.nanmean(srCs1['Nt'][srCs1['tempC'] >= 5])))
print('{}'.format(np.nanmean(pfs2['Nt'][pfs2['tempC'] >= 5])))



    


















    






0.0003004871395136338
0.00029442873122654655
8.454664557204488e-05
4.806236524058642e-05

0.00042241950085798936
0.00044106008822810175

Average of $TWC$ Median Values for $T \leq 0^{\circ}C$





In [22]:



    


twcMed_SRt_lte0 = npTypeZn['binTWC_median_SRt'][lte0]
twcMedAvg_SRt_lte0 = np.nanmean(twcMed_SRt_lte0)
print('{}'.format(twcMedAvg_SRt_lte0))

twcMed_ARt_lte0 = npTypeZn['binTWC_median_ARt'][lte0]
twcMedAvg_ARt_lte0 = np.nanmean(twcMed_ARt_lte0)
print('{}'.format(twcMedAvg_ARt_lte0))

twcMed_SRp_lte0 = npTypeZn['binTWC_median_SRp'][lte0]
twcMedAvg_SRp_lte0 = np.nanmean(twcMed_SRp_lte0)
print('{}'.format(twcMedAvg_SRp_lte0))

twcMed_ARp_lte0 = npTypeZn['binTWC_median_ARp'][lte0]
twcMedAvg_ARp_lte0 = np.nanmean(twcMed_ARp_lte0)
print('{}'.format(twcMedAvg_ARp_lte0))


print('\n{}'.format(np.nanmean(srCs1['twc'][srCs1['tempC'] <= 0])))
print('{}'.format(np.nanmean(pfs1['twc'][pfs1['tempC'] <= 0])))
print('{}'.format(np.nanmean(pfs2['twc'][pfs2['tempC'] <= 0])))



    


















    






0.6016717978447375
0.39216965399803094
0.10035667504800738
0.15351949434372475

0.4268076666204515
0.09170073436034976
0.0944918775240399

Average of $TWC$ Median Values for $T \geq 5^{\circ}C$





In [23]:



    


twcMed_SRt_gte5 = npTypeZn['binTWC_median_SRt'][gte5]
twcMedAvg_SRt_gte5 = np.nanmean(twcMed_SRt_gte5)
print('{}'.format(twcMedAvg_SRt_gte5))

twcMed_ARt_gte5 = npTypeZn['binTWC_median_ARt'][gte5]
twcMedAvg_ARt_gte5 = np.nanmean(twcMed_ARt_gte5)
print('{}'.format(twcMedAvg_ARt_gte5))

twcMed_SRp_gte5 = npTypeZn['binTWC_median_SRp'][gte5]
twcMedAvg_SRp_gte5 = np.nanmean(twcMed_SRp_gte5)
print('{}'.format(twcMedAvg_SRp_gte5))

twcMed_ARp_gte5 = npTypeZn['binTWC_median_ARp'][gte5]
twcMedAvg_ARp_gte5 = np.nanmean(twcMed_ARp_gte5)
print('{}'.format(twcMedAvg_ARp_gte5))


print('\n{}'.format(np.nanmean(srCs1['twc'][srCs1['tempC'] >= 5])))
print('{}'.format(np.nanmean(pfs2['twc'][pfs2['tempC'] >= 5])))



    


















    






0.049387467473419475
0.03243342144856784
0.017052943160593987
0.003157143602797839

0.09687523702310702
0.018950228302530937

...

...

...

Mean $ℛ_{RH}$ for $T \leq 0^{\circ}C$





In [24]:



    


RH_rStdv_SRt_lte0 = npTypeZn['binRH_stdv_SRt'][lte0]/npTypeZn['binRH_mean_SRt'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_SRt_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_SRt_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_SRt_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_SRt_lte0)*100))

RH_rStdv_ARt_lte0 = npTypeZn['binRH_stdv_ARt'][lte0]/npTypeZn['binRH_mean_ARt'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_ARt_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_ARt_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_ARt_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_ARt_lte0)*100))

RH_rStdv_SRp_lte0 = npTypeZn['binRH_stdv_SRp'][lte0]/npTypeZn['binRH_mean_SRp'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_SRp_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_SRp_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_SRp_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_SRp_lte0)*100))

RH_rStdv_ARp_lte0 = npTypeZn['binRH_stdv_ARp'][lte0]/npTypeZn['binRH_mean_ARp'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_ARp_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_ARp_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_ARp_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_ARp_lte0)*100))



    


















    






**Mean:	3.7

**Mean:	9.5

**Mean:	5.5

**Mean:	12.8

Mean $ℛ_{RH}$ for $T \geq 5^{\circ}C$





In [25]:



    


RH_rStdv_SRt_gte5 = npTypeZn['binRH_stdv_SRt'][gte5]/npTypeZn['binRH_mean_SRt'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_SRt_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_SRt_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_SRt_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_SRt_gte5)*100))

RH_rStdv_ARt_gte5 = npTypeZn['binRH_stdv_ARt'][gte5]/npTypeZn['binRH_mean_ARt'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_ARt_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_ARt_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_ARt_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_ARt_gte5)*100))

RH_rStdv_SRp_gte5 = npTypeZn['binRH_stdv_SRp'][gte5]/npTypeZn['binRH_mean_SRp'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_SRp_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_SRp_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_SRp_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_SRp_gte5)*100))

RH_rStdv_ARp_gte5 = npTypeZn['binRH_stdv_ARp'][gte5]/npTypeZn['binRH_mean_ARp'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_ARp_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_ARp_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_ARp_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_ARp_gte5)*100))



    


















    






**Mean:	18.8

**Mean:	5.4

**Mean:	6.0

**Mean:	8.9

Mean $ℛ_{𝒜}$ for $T \leq 0^{\circ}C$





In [26]:



    


ARat_rStdv_SRt_lte0 = npTypeZn['binARat_stdv_SRt'][lte0]/npTypeZn['binARat_mean_SRt'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_SRt_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_SRt_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_SRt_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_SRt_lte0)*100))

ARat_rStdv_ARt_lte0 = npTypeZn['binARat_stdv_ARt'][lte0]/npTypeZn['binARat_mean_ARt'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_ARt_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_ARt_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_ARt_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_ARt_lte0)*100))

ARat_rStdv_SRp_lte0 = npTypeZn['binARat_stdv_SRp'][lte0]/npTypeZn['binARat_mean_SRp'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_SRp_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_SRp_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_SRp_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_SRp_lte0)*100))

ARat_rStdv_ARp_lte0 = npTypeZn['binARat_stdv_ARp'][lte0]/npTypeZn['binARat_mean_ARp'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_ARp_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_ARp_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_ARp_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_ARp_lte0)*100))



    


















    






**Mean:	3.7

**Mean:	3.8

**Mean:	3.8

**Mean:	7.0

Mean $ℛ_{𝒜}$ for $T \geq 5^{\circ}C$





In [27]:



    


ARat_rStdv_SRt_gte5 = npTypeZn['binARat_stdv_SRt'][gte5]/npTypeZn['binARat_mean_SRt'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_SRt_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_SRt_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_SRt_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_SRt_gte5)*100))

ARat_rStdv_ARt_gte5 = npTypeZn['binARat_stdv_ARt'][gte5]/npTypeZn['binARat_mean_ARt'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_ARt_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_ARt_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_ARt_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_ARt_gte5)*100))

ARat_rStdv_SRp_gte5 = npTypeZn['binARat_stdv_SRp'][gte5]/npTypeZn['binARat_mean_SRp'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_SRp_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_SRp_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_SRp_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_SRp_gte5)*100))

ARat_rStdv_ARp_gte5 = npTypeZn['binARat_stdv_ARp'][gte5]/npTypeZn['binARat_mean_ARp'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_ARp_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_ARp_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_ARp_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_ARp_gte5)*100))



    


















    






**Mean:	9.9

**Mean:	1.7

**Mean:	9.5

**Mean:	4.9

Mean $ℛ_{D_{mm}}$ for $T \leq 0^{\circ}C$





In [28]:



    


Dmm_rStdv_SRt_lte0 = npTypeZn['binDmm_stdv_SRt'][lte0]/npTypeZn['binDmm_mean_SRt'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_SRt_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_SRt_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_SRt_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_SRt_lte0)*100))

Dmm_rStdv_ARt_lte0 = npTypeZn['binDmm_stdv_ARt'][lte0]/npTypeZn['binDmm_mean_ARt'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_ARt_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_ARt_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_ARt_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_ARt_lte0)*100))

Dmm_rStdv_SRp_lte0 = npTypeZn['binDmm_stdv_SRp'][lte0]/npTypeZn['binDmm_mean_SRp'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_SRp_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_SRp_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_SRp_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_SRp_lte0)*100))

Dmm_rStdv_ARp_lte0 = npTypeZn['binDmm_stdv_ARp'][lte0]/npTypeZn['binDmm_mean_ARp'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_ARp_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_ARp_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_ARp_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_ARp_lte0)*100))



    


















    






**Mean:	16.1

**Mean:	10.8

**Mean:	16.0

**Mean:	20.1

Mean $ℛ_{D_{mm}}$ for $T \geq 5^{\circ}C$





In [29]:



    


Dmm_rStdv_SRt_gte5 = npTypeZn['binDmm_stdv_SRt'][gte5]/npTypeZn['binDmm_mean_SRt'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_SRt_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_SRt_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_SRt_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_SRt_gte5)*100))

Dmm_rStdv_ARt_gte5 = npTypeZn['binDmm_stdv_ARt'][gte5]/npTypeZn['binDmm_mean_ARt'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_ARt_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_ARt_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_ARt_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_ARt_gte5)*100))

Dmm_rStdv_SRp_gte5 = npTypeZn['binDmm_stdv_SRp'][gte5]/npTypeZn['binDmm_mean_SRp'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_SRp_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_SRp_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_SRp_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_SRp_gte5)*100))

Dmm_rStdv_ARp_gte5 = npTypeZn['binDmm_stdv_ARp'][gte5]/npTypeZn['binDmm_mean_ARp'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_ARp_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_ARp_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_ARp_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_ARp_gte5)*100))



    


















    






**Mean:	27.3

**Mean:	13.2

**Mean:	24.3

**Mean:	17.3

Mean $ℛ_{N_t}$ for $T \leq 0^{\circ}C$





In [30]:



    


Nt_rStdv_SRt_lte0 = npTypeZn['binNt_stdv_SRt'][lte0]/npTypeZn['binNt_mean_SRt'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_SRt_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_SRt_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_SRt_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_SRt_lte0)*100))

Nt_rStdv_ARt_lte0 = npTypeZn['binNt_stdv_ARt'][lte0]/npTypeZn['binNt_mean_ARt'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_ARt_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_ARt_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_ARt_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_ARt_lte0)*100))

Nt_rStdv_SRp_lte0 = npTypeZn['binNt_stdv_SRp'][lte0]/npTypeZn['binNt_mean_SRp'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_SRp_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_SRp_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_SRp_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_SRp_lte0)*100))

Nt_rStdv_ARp_lte0 = npTypeZn['binNt_stdv_ARp'][lte0]/npTypeZn['binNt_mean_ARp'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_ARp_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_ARp_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_ARp_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_ARp_lte0)*100))



    


















    






**Mean:	79.3

**Mean:	74.3

**Mean:	48.8

**Mean:	76.6

Mean $ℛ_{N_t}$ for $T \geq 5^{\circ}C$





In [31]:



    


Nt_rStdv_SRt_gte5 = npTypeZn['binNt_stdv_SRt'][gte5]/npTypeZn['binNt_mean_SRt'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_SRt_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_SRt_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_SRt_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_SRt_gte5)*100))

Nt_rStdv_ARt_gte5 = npTypeZn['binNt_stdv_ARt'][gte5]/npTypeZn['binNt_mean_ARt'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_ARt_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_ARt_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_ARt_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_ARt_gte5)*100))

Nt_rStdv_SRp_gte5 = npTypeZn['binNt_stdv_SRp'][gte5]/npTypeZn['binNt_mean_SRp'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_SRp_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_SRp_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_SRp_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_SRp_gte5)*100))

Nt_rStdv_ARp_gte5 = npTypeZn['binNt_stdv_ARp'][gte5]/npTypeZn['binNt_mean_ARp'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_ARp_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_ARp_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_ARp_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_ARp_gte5)*100))



    


















    






**Mean:	86.1

**Mean:	43.0

**Mean:	45.1

**Mean:	19.1

Mean $ℛ_{TWC}$ for $T \leq 0^{\circ}C$





In [32]:



    


TWC_rStdv_SRt_lte0 = npTypeZn['binTWC_stdv_SRt'][lte0]/npTypeZn['binTWC_mean_SRt'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_SRt_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_SRt_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_SRt_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_SRt_lte0)*100))

TWC_rStdv_ARt_lte0 = npTypeZn['binTWC_stdv_ARt'][lte0]/npTypeZn['binTWC_mean_ARt'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_ARt_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_ARt_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_ARt_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_ARt_lte0)*100))

TWC_rStdv_SRp_lte0 = npTypeZn['binTWC_stdv_SRp'][lte0]/npTypeZn['binTWC_mean_SRp'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_SRp_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_SRp_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_SRp_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_SRp_lte0)*100))

TWC_rStdv_ARp_lte0 = npTypeZn['binTWC_stdv_ARp'][lte0]/npTypeZn['binTWC_mean_ARp'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_ARp_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_ARp_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_ARp_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_ARp_lte0)*100))



    


















    






**Mean:	49.2

**Mean:	53.4

**Mean:	42.0

**Mean:	91.6

Mean $ℛ_{TWC}$ for $T \geq 5^{\circ}C$





In [33]:



    


TWC_rStdv_SRt_gte5 = npTypeZn['binTWC_stdv_SRt'][gte5]/npTypeZn['binTWC_mean_SRt'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_SRt_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_SRt_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_SRt_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_SRt_gte5)*100))

TWC_rStdv_ARt_gte5 = npTypeZn['binTWC_stdv_ARt'][gte5]/npTypeZn['binTWC_mean_ARt'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_ARt_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_ARt_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_ARt_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_ARt_gte5)*100))

TWC_rStdv_SRp_gte5 = npTypeZn['binTWC_stdv_SRp'][gte5]/npTypeZn['binTWC_mean_SRp'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_SRp_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_SRp_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_SRp_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_SRp_gte5)*100))

TWC_rStdv_ARp_gte5 = npTypeZn['binTWC_stdv_ARp'][gte5]/npTypeZn['binTWC_mean_ARp'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_ARp_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_ARp_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_ARp_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_ARp_gte5)*100))



    


















    






**Mean:	78.7

**Mean:	44.0

**Mean:	79.6

**Mean:	33.8

Content source: stechma2/PECAN-MP-Analysis

Similar notebooks:

notebook.community | gallery | about