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]:



    


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



    











In [3]:



    


temps = npAllSprls['bin_mid']



    











In [4]:



    


temp_lte0 = npAllSprls['bin_mid'][:20]
temp_gt0 = npAllSprls['bin_mid'][20:]
temp_gte5 = npAllSprls['bin_mid'][24:]

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





In [ ]:



    


ntMed_TZ_lte0 = npAllSprls['binNt_median_TZ'][:20]
finiteIx = np.isfinite(np.log10(ntMed_TZ_lte0))
pf_ntMed_TZ_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(ntMed_TZ_lte0[finiteIx]),  1))
ntMed_TZ_lte0_dT = 1 - (10**pf_ntMed_TZ_lte0(-6)/10**pf_ntMed_TZ_lte0(-5)) # 1 deg C change percentage
        
print(ntMed_TZ_lte0_dT)

# fig, ax = plt.subplots(figsize=(16,20))
# ax.plot(ntMed_TZ_lte0,temp_lte0,'b-')
# ax.plot(ntMed_TZ_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
# ax.plot(10**(pf_ntMed_TZ_lte0(temp_lte0)),temp_lte0,'r-')

# ax.set_ylim(-18.5,5)
# ax.set_title('TZ\t\t$dN_t/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
# ax.set_xscale('log',nonposx='mask')
# ax.invert_yaxis()
# ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
# ax.tick_params(axis='both', which='major', labelsize=22)
# ax.grid(which='both')



    











In [ ]:



    


ntMed_SR_lte0 = npAllSprls['binNt_median_SR'][:20]
finiteIx = np.isfinite(np.log10(ntMed_SR_lte0))
pf_ntMed_SR_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(ntMed_SR_lte0[finiteIx]),  1))
ntMed_SR_lte0_dT = 1 - (10**pf_ntMed_SR_lte0(-6)/10**pf_ntMed_SR_lte0(-5)) # 1 deg C change percentage
        
print(ntMed_SR_lte0_dT)

# fig, ax = plt.subplots(figsize=(16,20))
# ax.plot(ntMed_SR_lte0,temp_lte0,'b-')
# ax.plot(ntMed_SR_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
# ax.plot(10**(pf_ntMed_SR_lte0(temp_lte0)),temp_lte0,'r-')

# ax.set_ylim(-18.5,5)
# ax.set_title('ESR\t\t$dN_t/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
# ax.set_xscale('log',nonposx='mask')
# ax.invert_yaxis()
# ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
# ax.tick_params(axis='both', which='major', labelsize=22)
# ax.grid(which='both')



    











In [ ]:



    


ntMed_AR_lte0 = npAllSprls['binNt_median_AR'][:20]
finiteIx = np.isfinite(np.log10(ntMed_AR_lte0))
pf_ntMed_AR_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(ntMed_AR_lte0[finiteIx]),  1))
ntMed_AR_lte0_dT = 1 - (10**pf_ntMed_AR_lte0(-6)/10**pf_ntMed_AR_lte0(-5)) # 1 deg C change percentage
        
print(ntMed_AR_lte0_dT)

# fig, ax = plt.subplots(figsize=(16,20))
# ax.plot(ntMed_AR_lte0,temp_lte0,'b-')
# ax.plot(ntMed_AR_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
# ax.plot(10**(pf_ntMed_AR_lte0(temp_lte0)),temp_lte0,'r-')

# ax.set_ylim(-18.5,5)
# ax.set_title('AR\t\t$dN_t/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
# ax.set_xscale('log',nonposx='mask')
# ax.invert_yaxis()
# ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
# ax.tick_params(axis='both', which='major', labelsize=22)
# ax.grid(which='both')

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





In [ ]:



    


ntMed_TZ_gt0 = npAllSprls['binNt_median_TZ'][20:]
finiteIx = np.isfinite(np.log10(ntMed_TZ_gt0))
pf_ntMed_TZ_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],np.log10(ntMed_TZ_gt0[finiteIx]),  1))
ntMed_TZ_gt0_dT = 1 - (10**pf_ntMed_TZ_gt0(2)/10**pf_ntMed_TZ_gt0(3)) # 1 deg C change percentage
        
print(ntMed_TZ_gt0_dT)

# fig, ax = plt.subplots(figsize=(16,20))
# ax.plot(ntMed_TZ_gt0,temp_gt0,'b-')
# ax.plot(ntMed_TZ_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
# ax.plot(10**(pf_ntMed_TZ_gt0(temp_gt0)),temp_gt0,'r-')

# ax.set_ylim(-5,20)
# ax.set_title('TZ\t\t$dN_t/dT$ for $T > 0^{\circ}C$',fontsize=24)
# ax.set_xscale('log',nonposx='mask')
# ax.invert_yaxis()
# ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
# ax.tick_params(axis='both', which='major', labelsize=22)
# ax.grid(which='both')



    











In [ ]:



    


ntMed_SR_gt0 = npAllSprls['binNt_median_SR'][20:]
finiteIx = np.isfinite(np.log10(ntMed_SR_gt0))
pf_ntMed_SR_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],np.log10(ntMed_SR_gt0[finiteIx]),  1))
ntMed_SR_gt0_dT = 1 - (10**pf_ntMed_SR_gt0(2)/10**pf_ntMed_SR_gt0(3)) # 1 deg C change percentage
        
print(ntMed_SR_gt0_dT)

# fig, ax = plt.subplots(figsize=(16,20))
# ax.plot(ntMed_SR_gt0,temp_gt0,'b-')
# ax.plot(ntMed_SR_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
# ax.plot(10**(pf_ntMed_SR_gt0(temp_gt0)),temp_gt0,'r-')

# ax.set_ylim(-5,20)
# ax.set_title('ESR\t\t$dN_t/dT$ for $T > 0^{\circ}C$',fontsize=24)
# ax.set_xscale('log',nonposx='mask')
# ax.invert_yaxis()
# ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
# ax.tick_params(axis='both', which='major', labelsize=22)
# ax.grid(which='both')



    











In [ ]:



    


ntMed_AR_gt0 = npAllSprls['binNt_median_AR'][20:]
finiteIx = np.isfinite(np.log10(ntMed_AR_gt0))
pf_ntMed_AR_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],np.log10(ntMed_AR_gt0[finiteIx]),  1))
ntMed_AR_gt0_dT = 1 - (10**pf_ntMed_AR_gt0(2)/10**pf_ntMed_AR_gt0(3)) # 1 deg C change percentage
        
print(ntMed_AR_gt0_dT)

# fig, ax = plt.subplots(figsize=(16,20))
# ax.plot(ntMed_AR_gt0,temp_gt0,'b-')
# ax.plot(ntMed_AR_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
# ax.plot(10**(pf_ntMed_AR_gt0(temp_gt0)),temp_gt0,'r-')

# ax.set_ylim(-5,20)
# ax.set_title('AR\t\t$dN_t/dT$ for $T > 0^{\circ}C$',fontsize=24)
# ax.set_xscale('log',nonposx='mask')
# ax.invert_yaxis()
# ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
# ax.tick_params(axis='both', which='major', labelsize=22)
# ax.grid(which='both')

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





In [ ]:



    


twcMed_TZ_lte0 = npAllSprls['binTWC_median_TZ'][:20]
finiteIx = np.isfinite(np.log10(twcMed_TZ_lte0))
pf_twcMed_TZ_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(twcMed_TZ_lte0[finiteIx]),  1))
twcMed_TZ_lte0_dT = 1 - (10**pf_twcMed_TZ_lte0(-6)/10**pf_twcMed_TZ_lte0(-5)) # 1 deg C change percentage
        
print(twcMed_TZ_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(twcMed_TZ_lte0,temp_lte0,'b-')
ax.plot(twcMed_TZ_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(10**(pf_twcMed_TZ_lte0(temp_lte0)),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('TZ\t\t$dTWC/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.set_xscale('log',nonposx='mask')
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


twcMed_SR_lte0 = npAllSprls['binTWC_median_SR'][:20]
finiteIx = np.isfinite(np.log10(twcMed_SR_lte0))
pf_twcMed_SR_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(twcMed_SR_lte0[finiteIx]),  1))
twcMed_SR_lte0_dT = 1 - (10**pf_twcMed_SR_lte0(-6)/10**pf_twcMed_SR_lte0(-5)) # 1 deg C change percentage
        
print(twcMed_SR_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(twcMed_SR_lte0,temp_lte0,'b-')
ax.plot(twcMed_SR_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(10**(pf_twcMed_SR_lte0(temp_lte0)),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('ESR\t\t$dTWC/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.set_xscale('log',nonposx='mask')
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


twcMed_AR_lte0 = npAllSprls['binTWC_median_AR'][:20]
finiteIx = np.isfinite(np.log10(twcMed_AR_lte0))
pf_twcMed_AR_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],np.log10(twcMed_AR_lte0[finiteIx]),  1))
twcMed_AR_lte0_dT = 1 - (10**pf_twcMed_AR_lte0(-6)/10**pf_twcMed_AR_lte0(-5)) # 1 deg C change percentage
        
print(twcMed_AR_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(twcMed_AR_lte0,temp_lte0,'b-')
ax.plot(twcMed_AR_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(10**(pf_twcMed_AR_lte0(temp_lte0)),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('AR\t\t$dTWC/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.set_xscale('log',nonposx='mask')
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')

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





In [ ]:



    


twcMed_TZ_gt0 = npAllSprls['binTWC_median_TZ'][20:]
finiteIx = np.isfinite(np.log10(twcMed_TZ_gt0))
pf_twcMed_TZ_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],np.log10(twcMed_TZ_gt0[finiteIx]),  1))
twcMed_TZ_gt0_dT = 1 - (10**pf_twcMed_TZ_gt0(2)/10**pf_twcMed_TZ_gt0(3)) # 1 deg C change percentage
        
print(twcMed_TZ_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(twcMed_TZ_gt0,temp_gt0,'b-')
ax.plot(twcMed_TZ_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(10**(pf_twcMed_TZ_gt0(temp_gt0)),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('TZ\t\t$dTWC/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.set_xscale('log',nonposx='mask')
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


twcMed_SR_gt0 = npAllSprls['binTWC_median_SR'][20:]
finiteIx = np.isfinite(np.log10(twcMed_SR_gt0))
pf_twcMed_SR_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],np.log10(twcMed_SR_gt0[finiteIx]),  1))
twcMed_SR_gt0_dT = 1 - (10**pf_twcMed_SR_gt0(2)/10**pf_twcMed_SR_gt0(3)) # 1 deg C change percentage
        
print(twcMed_SR_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(twcMed_SR_gt0,temp_gt0,'b-')
ax.plot(twcMed_SR_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(10**(pf_twcMed_SR_gt0(temp_gt0)),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('ESR\t\t$dTWC/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.set_xscale('log',nonposx='mask')
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


twcMed_AR_gt0 = npAllSprls['binTWC_median_AR'][20:]
finiteIx = np.isfinite(np.log10(twcMed_AR_gt0))
pf_twcMed_AR_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],np.log10(twcMed_AR_gt0[finiteIx]),  1))
twcMed_AR_gt0_dT = 1 - (10**pf_twcMed_AR_gt0(2)/10**pf_twcMed_AR_gt0(3)) # 1 deg C change percentage
        
print(twcMed_AR_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(twcMed_AR_gt0,temp_gt0,'b-')
ax.plot(twcMed_AR_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(10**(pf_twcMed_AR_gt0(temp_gt0)),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('AR\t\t$dTWC/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.set_xscale('log',nonposx='mask')
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')

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





In [ ]:



    


dmmMed_TZ_lte0 = npAllSprls['binDmm_median_TZ'][:20]

finiteIx = np.isfinite(dmmMed_TZ_lte0)
pf_dmmMed_TZ_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],dmmMed_TZ_lte0[finiteIx],  1))
dmmMed_TZ_lte0_dT = pf_dmmMed_TZ_lte0(-5) - pf_dmmMed_TZ_lte0(-6)
        
print(dmmMed_TZ_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(dmmMed_TZ_lte0,temp_lte0,'b-')
ax.plot(dmmMed_TZ_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_dmmMed_TZ_lte0(temp_lte0),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('TZ\t\t$dD_{mm}/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


dmmMed_SR_lte0 = npAllSprls['binDmm_median_SR'][:20]

finiteIx = np.isfinite(dmmMed_SR_lte0)
pf_dmmMed_SR_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],dmmMed_SR_lte0[finiteIx],  1))
dmmMed_SR_lte0_dT = pf_dmmMed_SR_lte0(-5) - pf_dmmMed_SR_lte0(-6)
        
print(dmmMed_SR_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(dmmMed_SR_lte0,temp_lte0,'b-')
ax.plot(dmmMed_SR_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_dmmMed_SR_lte0(temp_lte0),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('ESR\t\t$dD_{mm}/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


dmmMed_AR_lte0 = npAllSprls['binDmm_median_AR'][:20]

finiteIx = np.isfinite(dmmMed_AR_lte0)
pf_dmmMed_AR_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],dmmMed_AR_lte0[finiteIx],  1))
dmmMed_AR_lte0_dT = pf_dmmMed_AR_lte0(-5) - pf_dmmMed_AR_lte0(-6)
        
print(dmmMed_AR_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(dmmMed_AR_lte0,temp_lte0,'b-')
ax.plot(dmmMed_AR_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_dmmMed_AR_lte0(temp_lte0),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('AR\t\t$dD_{mm}/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')

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





In [ ]:



    


dmmMed_TZ_gt0 = npAllSprls['binDmm_median_TZ'][20:]

finiteIx = np.isfinite(dmmMed_TZ_gt0)
pf_dmmMed_TZ_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],dmmMed_TZ_gt0[finiteIx],  1))
dmmMed_TZ_gt0_dT = pf_dmmMed_TZ_gt0(6) - pf_dmmMed_TZ_gt0(5)
        
print(dmmMed_TZ_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(dmmMed_TZ_gt0,temp_gt0,'b-')
ax.plot(dmmMed_TZ_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_dmmMed_TZ_gt0(temp_gt0),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('TZ\t\t$dD_{mm}/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


dmmMed_SR_gt0 = npAllSprls['binDmm_median_SR'][20:]

finiteIx = np.isfinite(dmmMed_SR_gt0)
pf_dmmMed_SR_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],dmmMed_SR_gt0[finiteIx],  1))
dmmMed_SR_gt0_dT = pf_dmmMed_SR_gt0(6) - pf_dmmMed_SR_gt0(5)
        
print(dmmMed_SR_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(dmmMed_SR_gt0,temp_gt0,'b-')
ax.plot(dmmMed_SR_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_dmmMed_SR_gt0(temp_gt0),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('ESR\t\t$dD_{mm}/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


dmmMed_AR_gt0 = npAllSprls['binDmm_median_AR'][20:]

finiteIx = np.isfinite(dmmMed_AR_gt0)
pf_dmmMed_AR_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],dmmMed_AR_gt0[finiteIx],  1))
dmmMed_AR_gt0_dT = pf_dmmMed_AR_gt0(6) - pf_dmmMed_AR_gt0(5)
        
print(dmmMed_AR_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(dmmMed_AR_gt0,temp_gt0,'b-')
ax.plot(dmmMed_AR_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_dmmMed_AR_gt0(temp_gt0),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('AR\t\t$dD_{mm}/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')

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





In [ ]:



    


ARatMed_TZ_lte0 = npAllSprls['binARat_median_TZ'][:20]

finiteIx = np.isfinite(ARatMed_TZ_lte0)
pf_ARatMed_TZ_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],ARatMed_TZ_lte0[finiteIx],  1))
ARatMed_TZ_lte0_dT = pf_ARatMed_TZ_lte0(-5) - pf_ARatMed_TZ_lte0(-6)
        
print(ARatMed_TZ_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(ARatMed_TZ_lte0,temp_lte0,'b-')
ax.plot(ARatMed_TZ_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_ARatMed_TZ_lte0(temp_lte0),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('TZ\t\t$dAR/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


ARatMed_SR_lte0 = npAllSprls['binARat_median_SR'][:20]

finiteIx = np.isfinite(ARatMed_SR_lte0)
pf_ARatMed_SR_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],ARatMed_SR_lte0[finiteIx],  1))
ARatMed_SR_lte0_dT = pf_ARatMed_SR_lte0(-5) - pf_ARatMed_SR_lte0(-6)
        
print(ARatMed_SR_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(ARatMed_SR_lte0,temp_lte0,'b-')
ax.plot(ARatMed_SR_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_ARatMed_SR_lte0(temp_lte0),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('ESR\t\t$dAR/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


ARatMed_AR_lte0 = npAllSprls['binARat_median_AR'][:20]

finiteIx = np.isfinite(ARatMed_AR_lte0)
pf_ARatMed_AR_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],ARatMed_AR_lte0[finiteIx],  1))
ARatMed_AR_lte0_dT = pf_ARatMed_AR_lte0(-5) - pf_ARatMed_AR_lte0(-6)
        
print(ARatMed_AR_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(ARatMed_AR_lte0,temp_lte0,'b-')
ax.plot(ARatMed_AR_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_ARatMed_AR_lte0(temp_lte0),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('AR\t\t$dAR/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')

$dAR/dT$ for $T > 0^{\circ}C$





In [ ]:



    


ARatMed_TZ_gt0 = npAllSprls['binARat_median_TZ'][20:]

finiteIx = np.isfinite(ARatMed_TZ_gt0)
pf_ARatMed_TZ_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],ARatMed_TZ_gt0[finiteIx],  1))
ARatMed_TZ_gt0_dT = pf_ARatMed_TZ_gt0(6) - pf_ARatMed_TZ_gt0(5)
        
print(ARatMed_TZ_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(ARatMed_TZ_gt0,temp_gt0,'b-')
ax.plot(ARatMed_TZ_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_ARatMed_TZ_gt0(temp_gt0),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('TZ\t\t$dAR/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


ARatMed_SR_gt0 = npAllSprls['binARat_median_SR'][20:]

finiteIx = np.isfinite(ARatMed_SR_gt0)
pf_ARatMed_SR_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],ARatMed_SR_gt0[finiteIx],  1))
ARatMed_SR_gt0_dT = pf_ARatMed_SR_gt0(6) - pf_ARatMed_SR_gt0(5)
        
print(ARatMed_SR_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(ARatMed_SR_gt0,temp_gt0,'b-')
ax.plot(ARatMed_SR_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_ARatMed_SR_gt0(temp_gt0),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('ESR\t\t$dAR/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


ARatMed_AR_gt0 = npAllSprls['binARat_median_AR'][20:]

finiteIx = np.isfinite(ARatMed_AR_gt0)
pf_ARatMed_AR_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],ARatMed_AR_gt0[finiteIx],  1))
ARatMed_AR_gt0_dT = pf_ARatMed_AR_gt0(6) - pf_ARatMed_AR_gt0(5)
        
print(ARatMed_AR_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(ARatMed_AR_gt0,temp_gt0,'b-')
ax.plot(ARatMed_AR_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_ARatMed_AR_gt0(temp_gt0),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('AR\t\t$dAR/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')

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





In [ ]:



    


rhMed_TZ_lte0 = npAllSprls['binRH_median_TZ'][:20]

finiteIx = np.isfinite(rhMed_TZ_lte0)
pf_rhMed_TZ_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],rhMed_TZ_lte0[finiteIx],  1))
rhMed_TZ_lte0_dT = pf_rhMed_TZ_lte0(-5) - pf_rhMed_TZ_lte0(-6)
        
print(rhMed_TZ_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(rhMed_TZ_lte0,temp_lte0,'b-')
ax.plot(rhMed_TZ_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_rhMed_TZ_lte0(temp_lte0),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('TZ\t\t$dRH/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


rhMed_SR_lte0 = npAllSprls['binRH_median_SR'][:20]

finiteIx = np.isfinite(rhMed_SR_lte0)
pf_rhMed_SR_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],rhMed_SR_lte0[finiteIx],  1))
rhMed_SR_lte0_dT = pf_rhMed_SR_lte0(-5) - pf_rhMed_SR_lte0(-6)
        
print(rhMed_SR_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(rhMed_SR_lte0,temp_lte0,'b-')
ax.plot(rhMed_SR_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_rhMed_SR_lte0(temp_lte0),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('ESR\t\t$dRH/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


rhMed_AR_lte0 = npAllSprls['binRH_median_AR'][:20]

finiteIx = np.isfinite(rhMed_AR_lte0)
pf_rhMed_AR_lte0 = np.poly1d(np.polyfit(temp_lte0[finiteIx],rhMed_AR_lte0[finiteIx],  1))
rhMed_AR_lte0_dT = pf_rhMed_AR_lte0(-5) - pf_rhMed_AR_lte0(-6)
        
print(rhMed_AR_lte0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(rhMed_AR_lte0,temp_lte0,'b-')
ax.plot(rhMed_AR_lte0,temp_lte0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_rhMed_AR_lte0(temp_lte0),temp_lte0,'r-')

ax.set_ylim(-18.5,5)
ax.set_title('AR\t\t$dRH/dT$ for $T \leq 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')

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





In [ ]:



    


rhMed_TZ_gt0 = npAllSprls['binRH_median_TZ'][20:]

finiteIx = np.isfinite(rhMed_TZ_gt0)
pf_rhMed_TZ_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],rhMed_TZ_gt0[finiteIx],  1))
rhMed_TZ_gt0_dT = pf_rhMed_TZ_gt0(6) - pf_rhMed_TZ_gt0(5)
        
print(rhMed_TZ_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(rhMed_TZ_gt0,temp_gt0,'b-')
ax.plot(rhMed_TZ_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_rhMed_TZ_gt0(temp_gt0),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('TZ\t\t$dRH/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


rhMed_SR_gt0 = npAllSprls['binRH_median_SR'][20:]

finiteIx = np.isfinite(rhMed_SR_gt0)
pf_rhMed_SR_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],rhMed_SR_gt0[finiteIx],  1))
rhMed_SR_gt0_dT = pf_rhMed_SR_gt0(6) - pf_rhMed_SR_gt0(5)
        
print(rhMed_SR_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(rhMed_SR_gt0,temp_gt0,'b-')
ax.plot(rhMed_SR_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_rhMed_SR_gt0(temp_gt0),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('ESR\t\t$dRH/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')



    











In [ ]:



    


rhMed_AR_gt0 = npAllSprls['binRH_median_AR'][20:]

finiteIx = np.isfinite(rhMed_AR_gt0)
pf_rhMed_AR_gt0 = np.poly1d(np.polyfit(temp_gt0[finiteIx],rhMed_AR_gt0[finiteIx],  1))
rhMed_AR_gt0_dT = pf_rhMed_AR_gt0(6) - pf_rhMed_AR_gt0(5)
        
print(rhMed_AR_gt0_dT)

fig, ax = plt.subplots(figsize=(16,20))
ax.plot(rhMed_AR_gt0,temp_gt0,'b-')
ax.plot(rhMed_AR_gt0,temp_gt0,'bo',markeredgecolor='white',markeredgewidth=1,markersize=12)
ax.plot(pf_rhMed_AR_gt0(temp_gt0),temp_gt0,'r-')

ax.set_ylim(-5,20)
ax.set_title('AR\t\t$dRH/dT$ for $T > 0^{\circ}C$',fontsize=24)
ax.invert_yaxis()
ax.set_ylabel('Temperature ($^{\circ}$C)',fontsize=24)
ax.tick_params(axis='both', which='major', labelsize=22)
ax.grid(which='both')

...

...

...

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





In [ ]:



    


ntMed_TZ_lte0 = npAllSprls['binNt_median_TZ'][:20]
ntMedAvg_TZ_lte0 = np.nanmean(ntMed_TZ_lte0)
print('{}'.format(ntMedAvg_TZ_lte0))

ntMed_SR_lte0 = npAllSprls['binNt_median_SR'][:20]
ntMedAvg_SR_lte0 = np.nanmean(ntMed_SR_lte0)
print('{}'.format(ntMedAvg_SR_lte0))

ntMed_AR_lte0 = npAllSprls['binNt_median_AR'][:20]
ntMedAvg_AR_lte0 = np.nanmean(ntMed_AR_lte0)
print('{}'.format(ntMedAvg_AR_lte0))

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





In [ ]:



    


ntMed_TZ_gte5 = npAllSprls['binNt_median_TZ'][24:]
ntMedAvg_TZ_gte5 = np.nanmean(ntMed_TZ_gte5)
print('{}'.format(ntMedAvg_TZ_gte5))

ntMed_SR_gte5 = npAllSprls['binNt_median_SR'][24:]
ntMedAvg_SR_gte5 = np.nanmean(ntMed_SR_gte5)
print('{}'.format(ntMedAvg_SR_gte5))

ntMed_AR_gte5 = npAllSprls['binNt_median_AR'][24:]
ntMedAvg_AR_gte5 = np.nanmean(ntMed_AR_gte5)
print('{}'.format(ntMedAvg_AR_gte5))

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





In [ ]:



    


twcMed_TZ_lte0 = npAllSprls['binTWC_median_TZ'][:20]
twcMedAvg_TZ_lte0 = np.nanmean(twcMed_TZ_lte0)
print('{}'.format(twcMedAvg_TZ_lte0))

twcMed_SR_lte0 = npAllSprls['binTWC_median_SR'][:20]
twcMedAvg_SR_lte0 = np.nanmean(twcMed_SR_lte0)
print('{}'.format(twcMedAvg_SR_lte0))

twcMed_AR_lte0 = npAllSprls['binTWC_median_AR'][:20]
twcMedAvg_AR_lte0 = np.nanmean(twcMed_AR_lte0)
print('{}'.format(twcMedAvg_AR_lte0))

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





In [ ]:



    


twcMed_TZ_gte5 = npAllSprls['binTWC_median_TZ'][24:]
twcMedAvg_TZ_gte5 = np.nanmean(twcMed_TZ_gte5)
print('{}'.format(twcMedAvg_TZ_gte5))

twcMed_SR_gte5 = npAllSprls['binTWC_median_SR'][24:]
twcMedAvg_SR_gte5 = np.nanmean(twcMed_SR_gte5)
print('{}'.format(twcMedAvg_SR_gte5))

twcMed_AR_gte5 = npAllSprls['binTWC_median_AR'][24:]
twcMedAvg_AR_gte5 = np.nanmean(twcMed_AR_gte5)
print('{}'.format(twcMedAvg_AR_gte5))

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





In [ ]:



    


DmmMed_TZ_lte0 = npAllSprls['binDmm_median_TZ'][:20]
DmmMedAvg_TZ_lte0 = np.nanmean(DmmMed_TZ_lte0)
print('{}'.format(DmmMedAvg_TZ_lte0))

DmmMed_SR_lte0 = npAllSprls['binDmm_median_SR'][:20]
DmmMedAvg_SR_lte0 = np.nanmean(DmmMed_SR_lte0)
print('{}'.format(DmmMedAvg_SR_lte0))

DmmMed_AR_lte0 = npAllSprls['binDmm_median_AR'][:20]
DmmMedAvg_AR_lte0 = np.nanmean(DmmMed_AR_lte0)
print('{}'.format(DmmMedAvg_AR_lte0))

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





In [ ]:



    


DmmMed_TZ_gte5 = npAllSprls['binDmm_median_TZ'][24:]
DmmMedAvg_TZ_gte5 = np.nanmean(DmmMed_TZ_gte5)
print('{}'.format(DmmMedAvg_TZ_gte5))

DmmMed_SR_gte5 = npAllSprls['binDmm_median_SR'][24:]
DmmMedAvg_SR_gte5 = np.nanmean(DmmMed_SR_gte5)
print('{}'.format(DmmMedAvg_SR_gte5))

DmmMed_AR_gte5 = npAllSprls['binDmm_median_AR'][24:]
DmmMedAvg_AR_gte5 = np.nanmean(DmmMed_AR_gte5)
print('{}'.format(DmmMedAvg_AR_gte5))

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





In [ ]:



    


ARatMed_TZ_lte0 = npAllSprls['binARat_median_TZ'][:20]
ARatMedAvg_TZ_lte0 = np.nanmean(ARatMed_TZ_lte0)
print('{}'.format(ARatMedAvg_TZ_lte0))

ARatMed_SR_lte0 = npAllSprls['binARat_median_SR'][:20]
ARatMedAvg_SR_lte0 = np.nanmean(ARatMed_SR_lte0)
print('{}'.format(ARatMedAvg_SR_lte0))

ARatMed_AR_lte0 = npAllSprls['binARat_median_AR'][:20]
ARatMedAvg_AR_lte0 = np.nanmean(ARatMed_AR_lte0)
print('{}'.format(ARatMedAvg_AR_lte0))

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





In [ ]:



    


ARatMed_TZ_gte5 = npAllSprls['binARat_median_TZ'][24:]
ARatMedAvg_TZ_gte5 = np.nanmean(ARatMed_TZ_gte5)
print('{}'.format(ARatMedAvg_TZ_gte5))

ARatMed_SR_gte5 = npAllSprls['binARat_median_SR'][24:]
ARatMedAvg_SR_gte5 = np.nanmean(ARatMed_SR_gte5)
print('{}'.format(ARatMedAvg_SR_gte5))

ARatMed_AR_gte5 = npAllSprls['binARat_median_AR'][24:]
ARatMedAvg_AR_gte5 = np.nanmean(ARatMed_AR_gte5)
print('{}'.format(ARatMedAvg_AR_gte5))

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





In [ ]:



    


rhMed_TZ_lte0 = npAllSprls['binRH_median_TZ'][:20]
rhMedAvg_TZ_lte0 = np.nanmean(rhMed_TZ_lte0)
print('{}'.format(rhMedAvg_TZ_lte0))

rhMed_SR_lte0 = npAllSprls['binRH_median_SR'][:20]
rhMedAvg_SR_lte0 = np.nanmean(rhMed_SR_lte0)
print('{}'.format(rhMedAvg_SR_lte0))

rhMed_AR_lte0 = npAllSprls['binRH_median_AR'][:20]
rhMedAvg_AR_lte0 = np.nanmean(rhMed_AR_lte0)
print('{}'.format(rhMedAvg_AR_lte0))

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





In [ ]:



    


rhMed_TZ_gte5 = npAllSprls['binRH_median_TZ'][24:]
rhMedAvg_TZ_gte5 = np.nanmean(rhMed_TZ_gte5)
print('{}'.format(rhMedAvg_TZ_gte5))

rhMed_SR_gte5 = npAllSprls['binRH_median_SR'][24:]
rhMedAvg_SR_gte5 = np.nanmean(rhMed_SR_gte5)
print('{}'.format(rhMedAvg_SR_gte5))

rhMed_AR_gte5 = npAllSprls['binRH_median_AR'][24:]
rhMedAvg_AR_gte5 = np.nanmean(rhMed_AR_gte5)
print('{}'.format(rhMedAvg_AR_gte5))

IQR Crap





In [5]:



    


top = 1
bot = 4
topTemp = np.where(temps == top)[0][0]
botTemp = np.where(temps == bot)[0][0]+1



    











In [6]:



    


Nt_stdv_SR_lte0 = npAllSprls['binNt_stdv_SR'][:20]
Nt_stdv_SR_gt0lt5 = npAllSprls['binNt_stdv_SR'][20:24]
Nt_stdv_SR_gte5 = npAllSprls['binNt_stdv_SR'][24:]
Nt_stdv_SR_custom = npAllSprls['binNt_stdv_SR'][topTemp:botTemp]
print('StdDev T<=0:\n{}'.format(Nt_stdv_SR_lte0))
print('Min:\t{:.2e}'.format(np.nanmin(Nt_stdv_SR_lte0)))
print('Max:\t{:.2e}'.format(np.nanmax(Nt_stdv_SR_lte0)))
print('Mean:\t{:.2e}'.format(np.nanmean(Nt_stdv_SR_lte0)))
print('Median:\t{:.2e}\n'.format(np.nanmedian(Nt_stdv_SR_lte0)))

print('StdDev 0<T<5:\n{}'.format(Nt_stdv_SR_gt0lt5))
print('Min:\t{:.2e}'.format(np.nanmin(Nt_stdv_SR_gt0lt5)))
print('Max:\t{:.2e}'.format(np.nanmax(Nt_stdv_SR_gt0lt5)))
print('Mean:\t{:.2e}'.format(np.nanmean(Nt_stdv_SR_gt0lt5)))
print('Median:\t{:.2e}\n'.format(np.nanmedian(Nt_stdv_SR_gt0lt5)))

print('StdDev T>=5:\n{}'.format(Nt_stdv_SR_gte5))
print('Min:\t{:.2e}'.format(np.nanmin(Nt_stdv_SR_gte5)))
print('Max:\t{:.2e}'.format(np.nanmax(Nt_stdv_SR_gte5)))
print('Mean:\t{:.2e}'.format(np.nanmean(Nt_stdv_SR_gte5)))
print('Median:\t{:.2e}\n'.format(np.nanmedian(Nt_stdv_SR_gte5)))

print('StdDev {}<=T<={}:\n{}'.format(temps[topTemp],temps[botTemp],Nt_stdv_SR_custom))
print('Min:\t{:.2e}'.format(np.nanmin(Nt_stdv_SR_custom)))
print('Max:\t{:.2e}'.format(np.nanmax(Nt_stdv_SR_custom)))
print('Mean:\t{:.2e}'.format(np.nanmean(Nt_stdv_SR_custom)))
print('Median:\t{:.2e}\n'.format(np.nanmedian(Nt_stdv_SR_custom)))



    


















    






StdDev T<=0:
[       nan        nan 0.04255942 0.03093237 0.0982552  0.04033663
 0.05893312 0.06876636 0.10510475 0.09186531 0.06776093 0.05946396
 0.0505815  0.04950092 0.04046586 0.0281963  0.0198174  0.0147025
 0.01114647 0.00630322]
Min:	6.30e-03
Max:	1.05e-01
Mean:	4.91e-02
Median:	4.60e-02

StdDev 0<T<5:
[0.0066027  0.00366156 0.00239544 0.00149231]
Min:	1.49e-03
Max:	6.60e-03
Mean:	3.54e-03
Median:	3.03e-03

StdDev T>=5:
[8.22425880e-04 7.41081409e-04 7.60571750e-04 5.40307729e-04
 4.37779936e-04 5.95536724e-04 5.41883658e-04 4.13352406e-04
 4.92408110e-04 1.80143958e-04 1.97504133e-04 1.34181881e-04
 1.42614792e-04 7.81164925e-05 3.06177455e-05            nan]
Min:	3.06e-05
Max:	8.22e-04
Mean:	4.07e-04
Median:	4.38e-04

StdDev 1.0<=T<=5.0:
[0.0066027  0.00366156 0.00239544 0.00149231]
Min:	1.49e-03
Max:	6.60e-03
Mean:	3.54e-03
Median:	3.03e-03


















In [8]:



    


Nt_rStdv_SR_lte0 = npAllSprls['binNt_stdv_SR'][:20]/npAllSprls['binNt_mean_SR'][:20]
Nt_rStdv_SR_gt0lt5 = npAllSprls['binNt_stdv_SR'][20:24]/npAllSprls['binNt_mean_SR'][20:24]
Nt_rStdv_SR_gte5 = npAllSprls['binNt_stdv_SR'][24:]/npAllSprls['binNt_mean_SR'][24:]
Nt_rStdv_SR_custom = npAllSprls['binNt_stdv_SR'][topTemp:botTemp]/npAllSprls['binNt_mean_SR'][topTemp:botTemp]
print('relStdDev T<=0:\n{}'.format(Nt_rStdv_SR_lte0))
print('Min:\t{:.3f}'.format(np.nanmin(Nt_rStdv_SR_lte0)))
print('Max:\t{:.3f}'.format(np.nanmax(Nt_rStdv_SR_lte0)))
print('Mean:\t{:.3f}'.format(np.nanmean(Nt_rStdv_SR_lte0)))
print('Median:\t{:.3f}\n'.format(np.nanmedian(Nt_rStdv_SR_lte0)))

print('relStdDev 0<T<5:\n{}'.format(Nt_rStdv_SR_gt0lt5))
print('Min:\t{:.3f}'.format(np.nanmin(Nt_rStdv_SR_gt0lt5)))
print('Max:\t{:.3f}'.format(np.nanmax(Nt_rStdv_SR_gt0lt5)))
print('Mean:\t{:.3f}'.format(np.nanmean(Nt_rStdv_SR_gt0lt5)))
print('Median:\t{:.3f}\n'.format(np.nanmedian(Nt_rStdv_SR_gt0lt5)))

print('relStdDev T>=5:\n{}'.format(Nt_rStdv_SR_gte5))
print('Min:\t{:.3f}'.format(np.nanmin(Nt_rStdv_SR_gte5)))
print('Max:\t{:.3f}'.format(np.nanmax(Nt_rStdv_SR_gte5)))
print('Mean:\t{:.3f}'.format(np.nanmean(Nt_rStdv_SR_gte5)))
print('Median:\t{:.3f}\n'.format(np.nanmedian(Nt_rStdv_SR_gte5)))

print('relStdDev {}<=T<={}:\n{}'.format(temps[topTemp],temps[botTemp],Nt_rStdv_SR_custom))
print('Min:\t{:.3f}'.format(np.nanmin(Nt_rStdv_SR_custom)))
print('Max:\t{:.3f}'.format(np.nanmax(Nt_rStdv_SR_custom)))
print('Mean:\t{:.3f}'.format(np.nanmean(Nt_rStdv_SR_custom)))
print('Median:\t{:.3f}\n'.format(np.nanmedian(Nt_rStdv_SR_custom)))



    


















    






relStdDev T<=0:
[       nan        nan 0.17883354 0.40376687 0.92748839 0.4626925
 0.65272504 0.68421449 0.61537768 1.04441698 1.1390607  0.94053537
 0.82175297 1.01448959 1.03856811 1.02146908 1.21593557 1.1494163
 1.12746275 0.95731595]
Min:	0.179
Max:	1.216
Mean:	0.855
Median:	0.949

relStdDev 0<T<5:
[1.19671089 1.16705165 1.31075156 1.15462805]
Min:	1.155
Max:	1.311
Mean:	1.207
Median:	1.182

relStdDev T>=5:
[0.9342157  1.05874622 1.11223146 1.03345458 0.89860678 0.98674203
 1.28128558 1.06111426 1.06724571 0.78456412 0.8367635  0.78789877
 0.76762208 0.7501262  0.44298434        nan]
Min:	0.443
Max:	1.281
Mean:	0.920
Median:	0.934

relStdDev 1.0<=T<=5.0:
[1.19671089 1.16705165 1.31075156 1.15462805]
Min:	1.155
Max:	1.311
Mean:	1.207
Median:	1.182


















In [ ]:



    


Nt_iqr_SR_lte0 = npAllSprls['binNt_75pct_SR'][:20] - npAllSprls['binNt_25pct_SR'][:20]
Nt_iqr_SR_gt0lt5 = npAllSprls['binNt_75pct_SR'][20:24] - npAllSprls['binNt_25pct_SR'][20:24]
Nt_iqr_SR_gte5 = npAllSprls['binNt_75pct_SR'][24:] - npAllSprls['binNt_25pct_SR'][24:]
Nt_iqr_SR_custom = npAllSprls['binNt_75pct_SR'][topTemp:botTemp] - npAllSprls['binNt_25pct_SR'][topTemp:botTemp]
print('IQR T<=0:\n{}'.format(Nt_iqr_SR_lte0))
print('Min:\t{:.2e}'.format(np.nanmin(Nt_iqr_SR_lte0)))
print('Max:\t{:.2e}'.format(np.nanmax(Nt_iqr_SR_lte0)))
print('Mean:\t{:.2e}'.format(np.nanmean(Nt_iqr_SR_lte0)))
print('Median:\t{:.2e}\n'.format(np.nanmedian(Nt_iqr_SR_lte0)))

print('IQR 0<T<5:\n{}'.format(Nt_iqr_SR_gt0lt5))
print('Min:\t{:.2e}'.format(np.nanmin(Nt_iqr_SR_gt0lt5)))
print('Max:\t{:.2e}'.format(np.nanmax(Nt_iqr_SR_gt0lt5)))
print('Mean:\t{:.2e}'.format(np.nanmean(Nt_iqr_SR_gt0lt5)))
print('Median:\t{:.2e}\n'.format(np.nanmedian(Nt_iqr_SR_gt0lt5)))

print('IQR T>=5:\n{}'.format(Nt_iqr_SR_gte5))
print('Min:\t{:.2e}'.format(np.nanmin(Nt_iqr_SR_gte5)))
print('Max:\t{:.2e}'.format(np.nanmax(Nt_iqr_SR_gte5)))
print('Mean:\t{:.2e}'.format(np.nanmean(Nt_iqr_SR_gte5)))
print('Median:\t{:.2e}\n'.format(np.nanmedian(Nt_iqr_SR_gte5)))

print('IQR {}<=T<={}:\n{}'.format(temps[topTemp],temps[botTemp],Nt_iqr_SR_custom))
print('Min:\t{:.2e}'.format(np.nanmin(Nt_iqr_SR_custom)))
print('Max:\t{:.2e}'.format(np.nanmax(Nt_iqr_SR_custom)))
print('Mean:\t{:.2e}'.format(np.nanmean(Nt_iqr_SR_custom)))
print('Median:\t{:.2e}\n'.format(np.nanmedian(Nt_iqr_SR_custom)))

Content source: stechma2/PECAN-MP-Analysis

Similar notebooks:

notebook.community | gallery | about