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]:
# npRIJall = np.load('/Users/danstechman/GoogleDrive/PECAN-Data/mp-data/esr_RIJcmp_stats.npz')
npRIJall = np.load('/Users/danstechman/GoogleDrive/PECAN-Data/mp-data/esr_RIJcmp_stats_excldU4S7.npz')

In [3]:
temps = npRIJall['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_ALLrij_lte0 = npRIJall['binRH_median_ALLrij'][lte0]
finiteIx = np.isfinite(np.log10(RHMed_ALLrij_lte0))
pf_RHMed_ALLrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],RHMed_ALLrij_lte0[finiteIx],  1))
RHMed_ALLrij_lte0_dT = pf_RHMed_ALLrij_lte0(-5) - pf_RHMed_ALLrij_lte0(-6)
print('{}'.format(RHMed_ALLrij_lte0_dT))

RHMed_ALLnrij_lte0 = npRIJall['binRH_median_ALLnrij'][lte0]
finiteIx = np.isfinite(np.log10(RHMed_ALLnrij_lte0))
pf_RHMed_ALLnrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],RHMed_ALLnrij_lte0[finiteIx],  1))
RHMed_ALLnrij_lte0_dT = pf_RHMed_ALLnrij_lte0(-5) - pf_RHMed_ALLnrij_lte0(-6)
print('{}\n\n'.format(RHMed_ALLnrij_lte0_dT))


RHMed_SRrij_lte0 = npRIJall['binRH_median_SRrij'][lte0]
finiteIx = np.isfinite(np.log10(RHMed_SRrij_lte0))
pf_RHMed_SRrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],RHMed_SRrij_lte0[finiteIx],  1))
RHMed_SRrij_lte0_dT = pf_RHMed_SRrij_lte0(-5) - pf_RHMed_SRrij_lte0(-6)
print('{}'.format(RHMed_SRrij_lte0_dT))

RHMed_SRnrij_lte0 = npRIJall['binRH_median_SRnrij'][lte0]
finiteIx = np.isfinite(np.log10(RHMed_SRnrij_lte0))
pf_RHMed_SRnrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],RHMed_SRnrij_lte0[finiteIx],  1))
RHMed_SRnrij_lte0_dT = pf_RHMed_SRnrij_lte0(-5) - pf_RHMed_SRnrij_lte0(-6)
print('{}'.format(RHMed_SRnrij_lte0_dT))


-1.4040355661747412
-1.3088817107362871


-1.523725319680409
-1.2417910983042049

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


In [5]:
RHMed_ALLrij_gte5 = npRIJall['binRH_median_ALLrij'][gte5]
finiteIx = np.isfinite(np.log10(RHMed_ALLrij_gte5))
pf_RHMed_ALLrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],RHMed_ALLrij_gte5[finiteIx],  1))
RHMed_ALLrij_gte5_dT = pf_RHMed_ALLrij_gte5(6) - pf_RHMed_ALLrij_gte5(5)
print('{}'.format(RHMed_ALLrij_gte5_dT))

RHMed_ALLnrij_gte5 = npRIJall['binRH_median_ALLnrij'][gte5]
finiteIx = np.isfinite(np.log10(RHMed_ALLnrij_gte5))
pf_RHMed_ALLnrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],RHMed_ALLnrij_gte5[finiteIx],  1))
RHMed_ALLnrij_gte5_dT = pf_RHMed_ALLnrij_gte5(6) - pf_RHMed_ALLnrij_gte5(5)
print('{}\n\n'.format(RHMed_ALLnrij_gte5_dT))


RHMed_SRrij_gte5 = npRIJall['binRH_median_SRrij'][gte5]
finiteIx = np.isfinite(np.log10(RHMed_SRrij_gte5))
pf_RHMed_SRrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],RHMed_SRrij_gte5[finiteIx],  1))
RHMed_SRrij_gte5_dT = pf_RHMed_SRrij_gte5(6) - pf_RHMed_SRrij_gte5(5)
print('{}'.format(RHMed_SRrij_gte5_dT))

RHMed_SRnrij_gte5 = npRIJall['binRH_median_SRnrij'][gte5]
finiteIx = np.isfinite(np.log10(RHMed_SRnrij_gte5))
pf_RHMed_SRnrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],RHMed_SRnrij_gte5[finiteIx],  1))
RHMed_SRnrij_gte5_dT = pf_RHMed_SRnrij_gte5(6) - pf_RHMed_SRnrij_gte5(5)
print('{}'.format(RHMed_SRnrij_gte5_dT))


-3.202165434613221
-3.0425909310103236


-2.62017262284229
-3.116374160743547

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


In [6]:
ARatMed_ALLrij_lte0 = npRIJall['binARat_median_ALLrij'][lte0]
finiteIx = np.isfinite(np.log10(ARatMed_ALLrij_lte0))
pf_ARatMed_ALLrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],ARatMed_ALLrij_lte0[finiteIx],  1))
ARatMed_ALLrij_lte0_dT = pf_ARatMed_ALLrij_lte0(-5) - pf_ARatMed_ALLrij_lte0(-6) 
print('{}'.format(ARatMed_ALLrij_lte0_dT))

ARatMed_ALLnrij_lte0 = npRIJall['binARat_median_ALLnrij'][lte0]
finiteIx = np.isfinite(np.log10(ARatMed_ALLnrij_lte0))
pf_ARatMed_ALLnrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],ARatMed_ALLnrij_lte0[finiteIx],  1))
ARatMed_ALLnrij_lte0_dT = pf_ARatMed_ALLnrij_lte0(-5) - pf_ARatMed_ALLnrij_lte0(-6)
print('{}\n\n'.format(ARatMed_ALLnrij_lte0_dT))


ARatMed_SRrij_lte0 = npRIJall['binARat_median_SRrij'][lte0]
finiteIx = np.isfinite(np.log10(ARatMed_SRrij_lte0))
pf_ARatMed_SRrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],ARatMed_SRrij_lte0[finiteIx],  1))
ARatMed_SRrij_lte0_dT = pf_ARatMed_SRrij_lte0(-5) - pf_ARatMed_SRrij_lte0(-6)
print('{}'.format(ARatMed_SRrij_lte0_dT))

ARatMed_SRnrij_lte0 = npRIJall['binARat_median_SRnrij'][lte0]
finiteIx = np.isfinite(np.log10(ARatMed_SRnrij_lte0))
pf_ARatMed_SRnrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],ARatMed_SRnrij_lte0[finiteIx],  1))
ARatMed_SRnrij_lte0_dT = pf_ARatMed_SRnrij_lte0(-5) - pf_ARatMed_SRnrij_lte0(-6)
print('{}'.format(ARatMed_SRnrij_lte0_dT))


0.00013050871029601785
-0.0019009404624312776


0.0014921020678541064
-0.0019777554671047137

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


In [7]:
ARatMed_ALLrij_gte5 = npRIJall['binARat_median_ALLrij'][gte5]
finiteIx = np.isfinite(np.log10(ARatMed_ALLrij_gte5))
pf_ARatMed_ALLrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],ARatMed_ALLrij_gte5[finiteIx],  1))
ARatMed_ALLrij_gte5_dT = pf_ARatMed_ALLrij_gte5(6) - pf_ARatMed_ALLrij_gte5(5)
print('{}'.format(ARatMed_ALLrij_gte5_dT))

ARatMed_ALLnrij_gte5 = npRIJall['binARat_median_ALLnrij'][gte5]
finiteIx = np.isfinite(np.log10(ARatMed_ALLnrij_gte5))
pf_ARatMed_ALLnrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],ARatMed_ALLnrij_gte5[finiteIx],  1))
ARatMed_ALLnrij_gte5_dT = pf_ARatMed_ALLnrij_gte5(6) - pf_ARatMed_ALLnrij_gte5(5)
print('{}\n\n'.format(ARatMed_ALLnrij_gte5_dT))


ARatMed_SRrij_gte5 = npRIJall['binARat_median_SRrij'][gte5]
finiteIx = np.isfinite(np.log10(ARatMed_SRrij_gte5))
pf_ARatMed_SRrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],ARatMed_SRrij_gte5[finiteIx],  1))
ARatMed_SRrij_gte5_dT = pf_ARatMed_SRrij_gte5(6) - pf_ARatMed_SRrij_gte5(5)
print('{}'.format(ARatMed_SRrij_gte5_dT))

ARatMed_SRnrij_gte5 = npRIJall['binARat_median_SRnrij'][gte5]
finiteIx = np.isfinite(np.log10(ARatMed_SRnrij_gte5))
pf_ARatMed_SRnrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],ARatMed_SRnrij_gte5[finiteIx],  1))
ARatMed_SRnrij_gte5_dT = pf_ARatMed_SRnrij_gte5(6) - pf_ARatMed_SRnrij_gte5(5)
print('{}'.format(ARatMed_SRnrij_gte5_dT))


0.007042545067728612
0.00017222091792234995


0.005418997106827805
0.00017004020394628672

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


In [8]:
DmmMed_ALLrij_lte0 = npRIJall['binDmm_median_ALLrij'][lte0]
finiteIx = np.isfinite(np.log10(DmmMed_ALLrij_lte0))
pf_DmmMed_ALLrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],DmmMed_ALLrij_lte0[finiteIx],  1))
DmmMed_ALLrij_lte0_dT = pf_DmmMed_ALLrij_lte0(-5) - pf_DmmMed_ALLrij_lte0(-6)
print('{}'.format(DmmMed_ALLrij_lte0_dT))

DmmMed_ALLnrij_lte0 = npRIJall['binDmm_median_ALLnrij'][lte0]
finiteIx = np.isfinite(np.log10(DmmMed_ALLnrij_lte0))
pf_DmmMed_ALLnrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],DmmMed_ALLnrij_lte0[finiteIx],  1))
DmmMed_ALLnrij_lte0_dT = pf_DmmMed_ALLnrij_lte0(-5) - pf_DmmMed_ALLnrij_lte0(-6)
print('{}\n\n'.format(DmmMed_ALLnrij_lte0_dT))


DmmMed_SRrij_lte0 = npRIJall['binDmm_median_SRrij'][lte0]
finiteIx = np.isfinite(np.log10(DmmMed_SRrij_lte0))
pf_DmmMed_SRrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],DmmMed_SRrij_lte0[finiteIx],  1))
DmmMed_SRrij_lte0_dT = pf_DmmMed_SRrij_lte0(-5) - pf_DmmMed_SRrij_lte0(-6)
print('{}'.format(DmmMed_SRrij_lte0_dT))

DmmMed_SRnrij_lte0 = npRIJall['binDmm_median_SRnrij'][lte0]
finiteIx = np.isfinite(np.log10(DmmMed_SRnrij_lte0))
pf_DmmMed_SRnrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],DmmMed_SRnrij_lte0[finiteIx],  1))
DmmMed_SRnrij_lte0_dT = pf_DmmMed_SRnrij_lte0(-5) - pf_DmmMed_SRnrij_lte0(-6)
print('{}'.format(DmmMed_SRnrij_lte0_dT))


0.007707688338493224
0.022123323013415686


0.02767857142857122
0.027096233230134104

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


In [9]:
DmmMed_ALLrij_gte5 = npRIJall['binDmm_median_ALLrij'][gte5]
finiteIx = np.isfinite(np.log10(DmmMed_ALLrij_gte5))
pf_DmmMed_ALLrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],DmmMed_ALLrij_gte5[finiteIx],  1))
DmmMed_ALLrij_gte5_dT = pf_DmmMed_ALLrij_gte5(6) - pf_DmmMed_ALLrij_gte5(5) 
print('{}'.format(DmmMed_ALLrij_gte5_dT))

DmmMed_ALLnrij_gte5 = npRIJall['binDmm_median_ALLnrij'][gte5]
finiteIx = np.isfinite(np.log10(DmmMed_ALLnrij_gte5))
pf_DmmMed_ALLnrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],DmmMed_ALLnrij_gte5[finiteIx],  1))
DmmMed_ALLnrij_gte5_dT = pf_DmmMed_ALLnrij_gte5(6) - pf_DmmMed_ALLnrij_gte5(5)
print('{}\n\n'.format(DmmMed_ALLnrij_gte5_dT))


DmmMed_SRrij_gte5 = npRIJall['binDmm_median_SRrij'][gte5]
finiteIx = np.isfinite(np.log10(DmmMed_SRrij_gte5))
pf_DmmMed_SRrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],DmmMed_SRrij_gte5[finiteIx],  1))
DmmMed_SRrij_gte5_dT = pf_DmmMed_SRrij_gte5(6) - pf_DmmMed_SRrij_gte5(5)
print('{}'.format(DmmMed_SRrij_gte5_dT))

DmmMed_SRnrij_gte5 = npRIJall['binDmm_median_SRnrij'][gte5]
finiteIx = np.isfinite(np.log10(DmmMed_SRnrij_gte5))
pf_DmmMed_SRnrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],DmmMed_SRnrij_gte5[finiteIx],  1))
DmmMed_SRnrij_gte5_dT = pf_DmmMed_SRnrij_gte5(6) - pf_DmmMed_SRnrij_gte5(5)
print('{}'.format(DmmMed_SRnrij_gte5_dT))


-0.007878787878787818
0.005022321428571397


0.02204545454545459
0.0019642857142856407

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


In [10]:
ntMed_ALLrij_lte0 = npRIJall['binNt_median_ALLrij'][lte0]
finiteIx = np.isfinite(np.log10(ntMed_ALLrij_lte0))
pf_ntMed_ALLrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(ntMed_ALLrij_lte0[finiteIx]),  1))
ntMed_ALLrij_lte0_dT = 1 - (10**pf_ntMed_ALLrij_lte0(-6)/10**pf_ntMed_ALLrij_lte0(-5)) # 1 deg C change percentage     
print('{}'.format(ntMed_ALLrij_lte0_dT))

ntMed_ALLnrij_lte0 = npRIJall['binNt_median_ALLnrij'][lte0]
finiteIx = np.isfinite(np.log10(ntMed_ALLnrij_lte0))
pf_ntMed_ALLnrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(ntMed_ALLnrij_lte0[finiteIx]),  1))
ntMed_ALLnrij_lte0_dT = 1 - (10**pf_ntMed_ALLnrij_lte0(-6)/10**pf_ntMed_ALLnrij_lte0(-5)) # 1 deg C change percentage      
print('{}\n\n'.format(ntMed_ALLnrij_lte0_dT))


ntMed_SRrij_lte0 = npRIJall['binNt_median_SRrij'][lte0]
finiteIx = np.isfinite(np.log10(ntMed_SRrij_lte0))
pf_ntMed_SRrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(ntMed_SRrij_lte0[finiteIx]),  1))
ntMed_SRrij_lte0_dT = 1 - (10**pf_ntMed_SRrij_lte0(-6)/10**pf_ntMed_SRrij_lte0(-5)) # 1 deg C change percentage    
print('{}'.format(ntMed_SRrij_lte0_dT))

ntMed_SRnrij_lte0 = npRIJall['binNt_median_SRnrij'][lte0]
finiteIx = np.isfinite(np.log10(ntMed_SRnrij_lte0))
pf_ntMed_SRnrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(ntMed_SRnrij_lte0[finiteIx]),  1))
ntMed_SRnrij_lte0_dT = 1 - (10**pf_ntMed_SRnrij_lte0(-6)/10**pf_ntMed_SRnrij_lte0(-5)) # 1 deg C change percentage    
print('{}'.format(ntMed_SRnrij_lte0_dT))


-0.12461149471419963
-0.23308133256876973


-0.22762985064809427
-0.22983764565753262

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


In [11]:
ntMed_ALLrij_gte5 = npRIJall['binNt_median_ALLrij'][gte5]
finiteIx = np.isfinite(np.log10(ntMed_ALLrij_gte5))
pf_ntMed_ALLrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(ntMed_ALLrij_gte5[finiteIx]),  1))
ntMed_ALLrij_gte5_dT = 1 - (10**pf_ntMed_ALLrij_gte5(-6)/10**pf_ntMed_ALLrij_gte5(-5)) # 1 deg C change percentage    
print('{}'.format(ntMed_ALLrij_gte5_dT))

ntMed_ALLnrij_gte5 = npRIJall['binNt_median_ALLnrij'][gte5]
finiteIx = np.isfinite(np.log10(ntMed_ALLnrij_gte5))
pf_ntMed_ALLnrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(ntMed_ALLnrij_gte5[finiteIx]),  1))
ntMed_ALLnrij_gte5_dT = 1 - (10**pf_ntMed_ALLnrij_gte5(-6)/10**pf_ntMed_ALLnrij_gte5(-5)) # 1 deg C change percentage 
print('{}\n\n'.format(ntMed_ALLnrij_gte5_dT))


ntMed_SRrij_gte5 = npRIJall['binNt_median_SRrij'][gte5]
finiteIx = np.isfinite(np.log10(ntMed_SRrij_gte5))
pf_ntMed_SRrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(ntMed_SRrij_gte5[finiteIx]),  1))
ntMed_SRrij_gte5_dT = 1 - (10**pf_ntMed_SRrij_gte5(-6)/10**pf_ntMed_SRrij_gte5(-5)) # 1 deg C change percentage
print('{}'.format(ntMed_SRrij_gte5_dT))


ntMed_SRnrij_gte5 = npRIJall['binNt_median_SRnrij'][gte5]
finiteIx = np.isfinite(np.log10(ntMed_SRnrij_gte5))
pf_ntMed_SRnrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(ntMed_SRnrij_gte5[finiteIx]),  1))
ntMed_SRnrij_gte5_dT = 1 - (10**pf_ntMed_SRnrij_gte5(-6)/10**pf_ntMed_SRnrij_gte5(-5)) # 1 deg C change percentage
print('{}'.format(ntMed_SRnrij_gte5_dT))


-0.07360393863120995
-0.11988040472671768


-0.09470030243816474
-0.11993443150071714

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


In [12]:
twcMed_ALLrij_lte0 = npRIJall['binTWC_median_ALLrij'][lte0]
finiteIx = np.isfinite(np.log10(twcMed_ALLrij_lte0))
pf_twcMed_ALLrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(twcMed_ALLrij_lte0[finiteIx]),  1))
twcMed_ALLrij_lte0_dT = 1 - (10**pf_twcMed_ALLrij_lte0(-6)/10**pf_twcMed_ALLrij_lte0(-5)) # 1 deg C change percentage
print('{}'.format(twcMed_ALLrij_lte0_dT))

twcMed_ALLnrij_lte0 = npRIJall['binTWC_median_ALLnrij'][lte0]
finiteIx = np.isfinite(np.log10(twcMed_ALLnrij_lte0))
pf_twcMed_ALLnrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(twcMed_ALLnrij_lte0[finiteIx]),  1))
twcMed_ALLnrij_lte0_dT = 1 - (10**pf_twcMed_ALLnrij_lte0(-6)/10**pf_twcMed_ALLnrij_lte0(-5)) # 1 deg C change percentage
print('{}\n\n'.format(twcMed_ALLnrij_lte0_dT))


twcMed_SRrij_lte0 = npRIJall['binTWC_median_SRrij'][lte0]
finiteIx = np.isfinite(np.log10(twcMed_SRrij_lte0))
pf_twcMed_SRrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(twcMed_SRrij_lte0[finiteIx]),  1))
twcMed_SRrij_lte0_dT = 1 - (10**pf_twcMed_SRrij_lte0(-6)/10**pf_twcMed_SRrij_lte0(-5)) # 1 deg C change percentage
print('{}'.format(twcMed_SRrij_lte0_dT))

twcMed_SRnrij_lte0 = npRIJall['binTWC_median_SRnrij'][lte0]
finiteIx = np.isfinite(np.log10(twcMed_SRnrij_lte0))
pf_twcMed_SRnrij_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(twcMed_SRnrij_lte0[finiteIx]),  1))
twcMed_SRnrij_lte0_dT = 1 - (10**pf_twcMed_SRnrij_lte0(-6)/10**pf_twcMed_SRnrij_lte0(-5)) # 1 deg C change percentage  
print('{}'.format(twcMed_SRnrij_lte0_dT))


-0.08673565091247282
-0.15621256309133225


-0.10263992831177249
-0.15813705798612143

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


In [13]:
TWCMed_ALLrij_gte5 = npRIJall['binTWC_median_ALLrij'][gte5]
finiteIx = np.isfinite(np.log10(TWCMed_ALLrij_gte5))
pf_TWCMed_ALLrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(TWCMed_ALLrij_gte5[finiteIx]),  1))
TWCMed_ALLrij_gte5_dT = 1 - (10**pf_TWCMed_ALLrij_gte5(-6)/10**pf_TWCMed_ALLrij_gte5(-5)) # 1 deg C change percentage
print('{}'.format(TWCMed_ALLrij_gte5_dT))

TWCMed_ALLnrij_gte5 = npRIJall['binTWC_median_ALLnrij'][gte5]
finiteIx = np.isfinite(np.log10(TWCMed_ALLnrij_gte5))
pf_TWCMed_ALLnrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(TWCMed_ALLnrij_gte5[finiteIx]),  1))
TWCMed_ALLnrij_gte5_dT = 1 - (10**pf_TWCMed_ALLnrij_gte5(-6)/10**pf_TWCMed_ALLnrij_gte5(-5)) # 1 deg C change percentage
print('{}\n\n'.format(TWCMed_ALLnrij_gte5_dT))


TWCMed_SRrij_gte5 = npRIJall['binTWC_median_SRrij'][gte5]
finiteIx = np.isfinite(np.log10(TWCMed_SRrij_gte5))
pf_TWCMed_SRrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(TWCMed_SRrij_gte5[finiteIx]),  1))
TWCMed_SRrij_gte5_dT = 1 - (10**pf_TWCMed_SRrij_gte5(-6)/10**pf_TWCMed_SRrij_gte5(-5)) # 1 deg C change percentage
print('{}'.format(TWCMed_SRrij_gte5_dT))

TWCMed_SRnrij_gte5 = npRIJall['binTWC_median_SRnrij'][gte5]
finiteIx = np.isfinite(np.log10(TWCMed_SRnrij_gte5))
pf_TWCMed_SRnrij_gte5 = np.poly1d(np.polyfit(temp_gte5[finiteIx],np.log10(TWCMed_SRnrij_gte5[finiteIx]),  1))
TWCMed_SRnrij_gte5_dT = 1 - (10**pf_TWCMed_SRnrij_gte5(-6)/10**pf_TWCMed_SRnrij_gte5(-5)) # 1 deg C change percentage
print('{}'.format(TWCMed_SRnrij_gte5_dT))


0.12088037344131808
-0.09249240649765778


0.12799962228442086
-0.09368190450609082

...

...

...

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


In [14]:
rhMed_ALLrij_lte0 = npRIJall['binRH_median_ALLrij'][lte0]
rhMedAvg_ALLrij_lte0 = np.nanmean(rhMed_ALLrij_lte0)
print('{}'.format(rhMedAvg_ALLrij_lte0))

rhMed_ALLnrij_lte0 = npRIJall['binRH_median_ALLnrij'][lte0]
rhMedAvg_ALLnrij_lte0 = np.nanmean(rhMed_ALLnrij_lte0)
print('{}'.format(rhMedAvg_ALLnrij_lte0))


rhMed_SRrij_lte0 = npRIJall['binRH_median_SRrij'][lte0]
rhMedAvg_SRrij_lte0 = np.nanmean(rhMed_SRrij_lte0)
print('{}'.format(rhMedAvg_SRrij_lte0))

rhMed_SRnrij_lte0 = npRIJall['binRH_median_SRnrij'][lte0]
rhMedAvg_SRnrij_lte0 = np.nanmean(rhMed_SRnrij_lte0)
print('{}'.format(rhMedAvg_SRnrij_lte0))


108.57388675064607
106.64797391025328
106.27009022648515
108.03117003650054

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


In [15]:
rhMed_ALLrij_gte5 = npRIJall['binRH_median_ALLrij'][gte5]
rhMedAvg_ALLrij_gte5 = np.nanmean(rhMed_ALLrij_gte5)
print('{}'.format(rhMedAvg_ALLrij_gte5))

rhMed_ALLnrij_gte5 = npRIJall['binRH_median_ALLnrij'][gte5]
rhMedAvg_ALLnrij_gte5 = np.nanmean(rhMed_ALLnrij_gte5)
print('{}'.format(rhMedAvg_ALLnrij_gte5))


rhMed_SRrij_gte5 = npRIJall['binRH_median_SRrij'][gte5]
rhMedAvg_SRrij_gte5 = np.nanmean(rhMed_SRrij_gte5)
print('{}'.format(rhMedAvg_SRrij_gte5))

rhMed_SRnrij_gte5 = npRIJall['binRH_median_SRnrij'][gte5]
rhMedAvg_SRnrij_gte5 = np.nanmean(rhMed_SRnrij_gte5)
print('{}'.format(rhMedAvg_SRnrij_gte5))


55.06341816264919
65.48604669267911
58.137902793819144
65.80391439435378

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


In [16]:
ARatMed_ALLrij_lte0 = npRIJall['binARat_median_ALLrij'][lte0]
ARatMedAvg_ALLrij_lte0 = np.nanmean(ARatMed_ALLrij_lte0)
print('{}'.format(ARatMedAvg_ALLrij_lte0))

ARatMed_ALLnrij_lte0 = npRIJall['binARat_median_ALLnrij'][lte0]
ARatMedAvg_ALLnrij_lte0 = np.nanmean(ARatMed_ALLnrij_lte0)
print('{}'.format(ARatMedAvg_ALLnrij_lte0))


ARatMed_SRrij_lte0 = npRIJall['binARat_median_SRrij'][lte0]
ARatMedAvg_SRrij_lte0 = np.nanmean(ARatMed_SRrij_lte0)
print('{}'.format(ARatMedAvg_SRrij_lte0))

ARatMed_SRnrij_lte0 = npRIJall['binARat_median_SRnrij'][lte0]
ARatMedAvg_SRnrij_lte0 = np.nanmean(ARatMed_SRnrij_lte0)
print('{}'.format(ARatMedAvg_SRnrij_lte0))


0.47271583818819046
0.4579316619840416
0.4682031571171745
0.4539074435369086

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


In [17]:
ARatMed_ALLrij_gte5 = npRIJall['binARat_median_ALLrij'][gte5]
ARatMedAvg_ALLrij_gte5 = np.nanmean(ARatMed_ALLrij_gte5)
print('{}'.format(ARatMedAvg_ALLrij_gte5))

ARatMed_ALLnrij_gte5 = npRIJall['binARat_median_ALLnrij'][gte5]
ARatMedAvg_ALLnrij_gte5 = np.nanmean(ARatMed_ALLnrij_gte5)
print('{}'.format(ARatMedAvg_ALLnrij_gte5))


ARatMed_SRrij_gte5 = npRIJall['binARat_median_SRrij'][gte5]
ARatMedAvg_SRrij_gte5 = np.nanmean(ARatMed_SRrij_gte5)
print('{}'.format(ARatMedAvg_SRrij_gte5))

ARatMed_SRnrij_gte5 = npRIJall['binARat_median_SRnrij'][gte5]
ARatMedAvg_SRnrij_gte5 = np.nanmean(ARatMed_SRnrij_gte5)
print('{}'.format(ARatMedAvg_SRnrij_gte5))


0.7003689676853355
0.7180264635387696
0.7202994101732528
0.7154642007173607

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


In [18]:
DmmMed_ALLrij_lte0 = npRIJall['binDmm_median_ALLrij'][lte0]
DmmMedAvg_ALLrij_lte0 = np.nanmean(DmmMed_ALLrij_lte0)
print('{}'.format(DmmMedAvg_ALLrij_lte0))

DmmMed_ALLnrij_lte0 = npRIJall['binDmm_median_ALLnrij'][lte0]
DmmMedAvg_ALLnrij_lte0 = np.nanmean(DmmMed_ALLnrij_lte0)
print('{}'.format(DmmMedAvg_ALLnrij_lte0))


DmmMed_SRrij_lte0 = npRIJall['binDmm_median_SRrij'][lte0]
DmmMedAvg_SRrij_lte0 = np.nanmean(DmmMed_SRrij_lte0)
print('{}'.format(DmmMedAvg_SRrij_lte0))

DmmMed_SRnrij_lte0 = npRIJall['binDmm_median_SRnrij'][lte0]
DmmMedAvg_SRnrij_lte0 = np.nanmean(DmmMed_SRnrij_lte0)
print('{}'.format(DmmMedAvg_SRnrij_lte0))


1.2656250000000002
1.2048611111111116
1.2129464285714289
1.146180555555556

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


In [19]:
DmmMed_ALLrij_gte5 = npRIJall['binDmm_median_ALLrij'][gte5]
DmmMedAvg_ALLrij_gte5 = np.nanmean(DmmMed_ALLrij_gte5)
print('{}'.format(DmmMedAvg_ALLrij_gte5))

DmmMed_ALLnrij_gte5 = npRIJall['binDmm_median_ALLnrij'][gte5]
DmmMedAvg_ALLnrij_gte5 = np.nanmean(DmmMed_ALLnrij_gte5)
print('{}'.format(DmmMedAvg_ALLnrij_gte5))


DmmMed_SRrij_gte5 = npRIJall['binDmm_median_SRrij'][gte5]
DmmMedAvg_SRrij_gte5 = np.nanmean(DmmMed_SRrij_gte5)
print('{}'.format(DmmMedAvg_SRrij_gte5))

DmmMed_SRnrij_gte5 = npRIJall['binDmm_median_SRnrij'][gte5]
DmmMedAvg_SRnrij_gte5 = np.nanmean(DmmMed_SRnrij_gte5)
print('{}'.format(DmmMedAvg_SRnrij_gte5))


0.9325000000000001
1.1645833333333335
0.9550000000000001
1.1858333333333335

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


In [20]:
ntMed_ALLrij_lte0 = npRIJall['binNt_median_ALLrij'][lte0]
ntMedAvg_ALLrij_lte0 = np.nanmean(ntMed_ALLrij_lte0)
print('{}'.format(ntMedAvg_ALLrij_lte0))

ntMed_ALLnrij_lte0 = npRIJall['binNt_median_ALLnrij'][lte0]
ntMedAvg_ALLnrij_lte0 = np.nanmean(ntMed_ALLnrij_lte0)
print('{}'.format(ntMedAvg_ALLnrij_lte0))


ntMed_SRrij_lte0 = npRIJall['binNt_median_SRrij'][lte0]
ntMedAvg_SRrij_lte0 = np.nanmean(ntMed_SRrij_lte0)
print('{}'.format(ntMedAvg_SRrij_lte0))

ntMed_SRnrij_lte0 = npRIJall['binNt_median_SRnrij'][lte0]
ntMedAvg_SRnrij_lte0 = np.nanmean(ntMed_SRnrij_lte0)
print('{}'.format(ntMedAvg_SRnrij_lte0))


0.04134617440842435
0.04181896482308808
0.058551538328716186
0.055421720122038005

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


In [21]:
ntMed_ALLrij_gte5 = npRIJall['binNt_median_ALLrij'][gte5]
ntMedAvg_ALLrij_gte5 = np.nanmean(ntMed_ALLrij_gte5)
print('{}'.format(ntMedAvg_ALLrij_gte5))

ntMed_ALLnrij_gte5 = npRIJall['binNt_median_ALLnrij'][gte5]
ntMedAvg_ALLnrij_gte5 = np.nanmean(ntMed_ALLnrij_gte5)
print('{}'.format(ntMedAvg_ALLnrij_gte5))


ntMed_SRrij_gte5 = npRIJall['binNt_median_SRrij'][gte5]
ntMedAvg_SRrij_gte5 = np.nanmean(ntMed_SRrij_gte5)
print('{}'.format(ntMedAvg_SRrij_gte5))

ntMed_SRnrij_gte5 = npRIJall['binNt_median_SRnrij'][gte5]
ntMedAvg_SRnrij_gte5 = np.nanmean(ntMed_SRnrij_gte5)
print('{}'.format(ntMedAvg_SRnrij_gte5))


0.0005943415996730159
0.00020179476694051184
0.0007706485595870881
0.00021790763608945301

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


In [22]:
twcMed_ALLrij_lte0 = npRIJall['binTWC_median_ALLrij'][lte0]
twcMedAvg_ALLrij_lte0 = np.nanmean(twcMed_ALLrij_lte0)
print('{}'.format(twcMedAvg_ALLrij_lte0))

twcMed_ALLnrij_lte0 = npRIJall['binTWC_median_ALLnrij'][lte0]
twcMedAvg_ALLnrij_lte0 = np.nanmean(twcMed_ALLnrij_lte0)
print('{}'.format(twcMedAvg_ALLnrij_lte0))


twcMed_SRrij_lte0 = npRIJall['binTWC_median_SRrij'][lte0]
twcMedAvg_SRrij_lte0 = np.nanmean(twcMed_SRrij_lte0)
print('{}'.format(twcMedAvg_SRrij_lte0))

twcMed_SRnrij_lte0 = npRIJall['binTWC_median_SRnrij'][lte0]
twcMedAvg_SRnrij_lte0 = np.nanmean(twcMed_SRnrij_lte0)
print('{}'.format(twcMedAvg_SRnrij_lte0))


0.6072602080434002
0.47996851906965904
0.5827800414795684
0.5490319345811208

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


In [23]:
twcMed_ALLrij_gte5 = npRIJall['binTWC_median_ALLrij'][gte5]
twcMedAvg_ALLrij_gte5 = np.nanmean(twcMed_ALLrij_gte5)
print('{}'.format(twcMedAvg_ALLrij_gte5))

twcMed_ALLnrij_gte5 = npRIJall['binTWC_median_ALLnrij'][gte5]
twcMedAvg_ALLnrij_gte5 = np.nanmean(twcMed_ALLnrij_gte5)
print('{}'.format(twcMedAvg_ALLnrij_gte5))


twcMed_SRrij_gte5 = npRIJall['binTWC_median_SRrij'][gte5]
twcMedAvg_SRrij_gte5 = np.nanmean(twcMed_SRrij_gte5)
print('{}'.format(twcMedAvg_SRrij_gte5))

twcMed_SRnrij_gte5 = npRIJall['binTWC_median_SRnrij'][gte5]
twcMedAvg_SRnrij_gte5 = np.nanmean(twcMed_SRnrij_gte5)
print('{}'.format(twcMedAvg_SRnrij_gte5))


0.04553579499109446
0.03778868209654382
0.05542281116684518
0.04168947453575408

...

...

...

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


In [ ]:
RH_rStdv_ALLrij_lte0 = npRIJall['binRH_stdv_ALLrij'][lte0]/npRIJall['binRH_mean_ALLrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_ALLrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_ALLrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_ALLrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_ALLrij_lte0)*100))

RH_rStdv_ALLnrij_lte0 = npRIJall['binRH_stdv_ALLnrij'][lte0]/npRIJall['binRH_mean_ALLnrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_ALLnrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_ALLnrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_ALLnrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_ALLnrij_lte0)*100))

RH_rStdv_SRrij_lte0 = npRIJall['binRH_stdv_SRrij'][lte0]/npRIJall['binRH_mean_SRrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_SRrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_SRrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_SRrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_SRrij_lte0)*100))

RH_rStdv_SRnrij_lte0 = npRIJall['binRH_stdv_SRnrij'][lte0]/npRIJall['binRH_mean_SRnrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_SRnrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_SRnrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_SRnrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_SRnrij_lte0)*100))

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


In [ ]:
RH_rStdv_ALLrij_gte5 = npRIJall['binRH_stdv_ALLrij'][gte5]/npRIJall['binRH_mean_ALLrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_ALLrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_ALLrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_ALLrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_ALLrij_gte5)*100))

RH_rStdv_ALLnrij_gte5 = npRIJall['binRH_stdv_ALLnrij'][gte5]/npRIJall['binRH_mean_ALLnrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_ALLnrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_ALLnrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_ALLnrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_ALLnrij_gte5)*100))

RH_rStdv_SRrij_gte5 = npRIJall['binRH_stdv_SRrij'][gte5]/npRIJall['binRH_mean_SRrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_SRrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_SRrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_SRrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_SRrij_gte5)*100))

RH_rStdv_SRnrij_gte5 = npRIJall['binRH_stdv_SRnrij'][gte5]/npRIJall['binRH_mean_SRnrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(RH_rStdv_SRnrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(RH_rStdv_SRnrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(RH_rStdv_SRnrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(RH_rStdv_SRnrij_gte5)*100))

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


In [ ]:
ARat_rStdv_ALLrij_lte0 = npRIJall['binARat_stdv_ALLrij'][lte0]/npRIJall['binARat_mean_ALLrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_ALLrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_ALLrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_ALLrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_ALLrij_lte0)*100))

ARat_rStdv_ALLnrij_lte0 = npRIJall['binARat_stdv_ALLnrij'][lte0]/npRIJall['binARat_mean_ALLnrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_ALLnrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_ALLnrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_ALLnrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_ALLnrij_lte0)*100))

ARat_rStdv_SRrij_lte0 = npRIJall['binARat_stdv_SRrij'][lte0]/npRIJall['binARat_mean_SRrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_SRrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_SRrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_SRrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_SRrij_lte0)*100))

ARat_rStdv_SRnrij_lte0 = npRIJall['binARat_stdv_SRnrij'][lte0]/npRIJall['binARat_mean_SRnrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_SRnrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_SRnrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_SRnrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_SRnrij_lte0)*100))

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


In [ ]:
ARat_rStdv_ALLrij_gte5 = npRIJall['binARat_stdv_ALLrij'][gte5]/npRIJall['binARat_mean_ALLrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_ALLrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_ALLrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_ALLrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_ALLrij_gte5)*100))

ARat_rStdv_ALLnrij_gte5 = npRIJall['binARat_stdv_ALLnrij'][gte5]/npRIJall['binARat_mean_ALLnrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_ALLnrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_ALLnrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_ALLnrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_ALLnrij_gte5)*100))

ARat_rStdv_SRrij_gte5 = npRIJall['binARat_stdv_SRrij'][gte5]/npRIJall['binARat_mean_SRrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_SRrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_SRrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_SRrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_SRrij_gte5)*100))

ARat_rStdv_SRnrij_gte5 = npRIJall['binARat_stdv_SRnrij'][gte5]/npRIJall['binARat_mean_SRnrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(ARat_rStdv_SRnrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(ARat_rStdv_SRnrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(ARat_rStdv_SRnrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(ARat_rStdv_SRnrij_gte5)*100))

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


In [ ]:
Dmm_rStdv_ALLrij_lte0 = npRIJall['binDmm_stdv_ALLrij'][lte0]/npRIJall['binDmm_mean_ALLrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_ALLrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_ALLrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_ALLrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_ALLrij_lte0)*100))

Dmm_rStdv_ALLnrij_lte0 = npRIJall['binDmm_stdv_ALLnrij'][lte0]/npRIJall['binDmm_mean_ALLnrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_ALLnrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_ALLnrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_ALLnrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_ALLnrij_lte0)*100))

Dmm_rStdv_SRrij_lte0 = npRIJall['binDmm_stdv_SRrij'][lte0]/npRIJall['binDmm_mean_SRrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_SRrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_SRrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_SRrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_SRrij_lte0)*100))

Dmm_rStdv_SRnrij_lte0 = npRIJall['binDmm_stdv_SRnrij'][lte0]/npRIJall['binDmm_mean_SRnrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_SRnrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_SRnrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_SRnrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_SRnrij_lte0)*100))

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


In [ ]:
Dmm_rStdv_ALLrij_gte5 = npRIJall['binDmm_stdv_ALLrij'][gte5]/npRIJall['binDmm_mean_ALLrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_ALLrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_ALLrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_ALLrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_ALLrij_gte5)*100))

Dmm_rStdv_ALLnrij_gte5 = npRIJall['binDmm_stdv_ALLnrij'][gte5]/npRIJall['binDmm_mean_ALLnrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_ALLnrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_ALLnrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_ALLnrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_ALLnrij_gte5)*100))

Dmm_rStdv_SRrij_gte5 = npRIJall['binDmm_stdv_SRrij'][gte5]/npRIJall['binDmm_mean_SRrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_SRrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_SRrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_SRrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_SRrij_gte5)*100))

Dmm_rStdv_SRnrij_gte5 = npRIJall['binDmm_stdv_SRnrij'][gte5]/npRIJall['binDmm_mean_SRnrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Dmm_rStdv_SRnrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Dmm_rStdv_SRnrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Dmm_rStdv_SRnrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Dmm_rStdv_SRnrij_gte5)*100))

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


In [ ]:
Nt_rStdv_ALLrij_lte0 = npRIJall['binNt_stdv_ALLrij'][lte0]/npRIJall['binNt_mean_ALLrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_ALLrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_ALLrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_ALLrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_ALLrij_lte0)*100))

Nt_rStdv_ALLnrij_lte0 = npRIJall['binNt_stdv_ALLnrij'][lte0]/npRIJall['binNt_mean_ALLnrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_ALLnrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_ALLnrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_ALLnrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_ALLnrij_lte0)*100))

Nt_rStdv_SRrij_lte0 = npRIJall['binNt_stdv_SRrij'][lte0]/npRIJall['binNt_mean_SRrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_SRrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_SRrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_SRrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_SRrij_lte0)*100))

Nt_rStdv_SRnrij_lte0 = npRIJall['binNt_stdv_SRnrij'][lte0]/npRIJall['binNt_mean_SRnrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_SRnrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_SRnrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_SRnrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_SRnrij_lte0)*100))

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


In [ ]:
Nt_rStdv_ALLrij_gte5 = npRIJall['binNt_stdv_ALLrij'][gte5]/npRIJall['binNt_mean_ALLrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_ALLrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_ALLrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_ALLrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_ALLrij_gte5)*100))

Nt_rStdv_ALLnrij_gte5 = npRIJall['binNt_stdv_ALLnrij'][gte5]/npRIJall['binNt_mean_ALLnrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_ALLnrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_ALLnrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_ALLnrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_ALLnrij_gte5)*100))

Nt_rStdv_SRrij_gte5 = npRIJall['binNt_stdv_SRrij'][gte5]/npRIJall['binNt_mean_SRrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_SRrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_SRrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_SRrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_SRrij_gte5)*100))

Nt_rStdv_SRnrij_gte5 = npRIJall['binNt_stdv_SRnrij'][gte5]/npRIJall['binNt_mean_SRnrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(Nt_rStdv_SRnrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(Nt_rStdv_SRnrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(Nt_rStdv_SRnrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(Nt_rStdv_SRnrij_gte5)*100))

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


In [ ]:
TWC_rStdv_ALLrij_lte0 = npRIJall['binTWC_stdv_ALLrij'][lte0]/npRIJall['binTWC_mean_ALLrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_ALLrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_ALLrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_ALLrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_ALLrij_lte0)*100))

TWC_rStdv_ALLnrij_lte0 = npRIJall['binTWC_stdv_ALLnrij'][lte0]/npRIJall['binTWC_mean_ALLnrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_ALLnrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_ALLnrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_ALLnrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_ALLnrij_lte0)*100))

TWC_rStdv_SRrij_lte0 = npRIJall['binTWC_stdv_SRrij'][lte0]/npRIJall['binTWC_mean_SRrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_SRrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_SRrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_SRrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_SRrij_lte0)*100))

TWC_rStdv_SRnrij_lte0 = npRIJall['binTWC_stdv_SRnrij'][lte0]/npRIJall['binTWC_mean_SRnrij'][lte0]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_SRnrij_lte0)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_SRnrij_lte0)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_SRnrij_lte0)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_SRnrij_lte0)*100))

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


In [ ]:
TWC_rStdv_ALLrij_gte5 = npRIJall['binTWC_stdv_ALLrij'][gte5]/npRIJall['binTWC_mean_ALLrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_ALLrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_ALLrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_ALLrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_ALLrij_gte5)*100))

TWC_rStdv_ALLnrij_gte5 = npRIJall['binTWC_stdv_ALLnrij'][gte5]/npRIJall['binTWC_mean_ALLnrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_ALLnrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_ALLnrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_ALLnrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_ALLnrij_gte5)*100))

TWC_rStdv_SRrij_gte5 = npRIJall['binTWC_stdv_SRrij'][gte5]/npRIJall['binTWC_mean_SRrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_SRrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_SRrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_SRrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_SRrij_gte5)*100))

TWC_rStdv_SRnrij_gte5 = npRIJall['binTWC_stdv_SRnrij'][gte5]/npRIJall['binTWC_mean_SRnrij'][gte5]
# print('Min:\t{:.1f}'.format(np.nanmin(TWC_rStdv_SRnrij_gte5)*100))
# print('Max:\t{:.1f}'.format(np.nanmax(TWC_rStdv_SRnrij_gte5)*100))
# print('Median:\t{:.1f}'.format(np.nanmedian(TWC_rStdv_SRnrij_gte5)*100))
print('**Mean:\t{:.1f}\n'.format(np.nanmean(TWC_rStdv_SRnrij_gte5)*100))