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