Commit 2867006b authored by Leonie Pick's avatar Leonie Pick

Checking values - no important changes to code.

parent 9ed855d7
This diff is collapsed.
......@@ -34,6 +34,8 @@ def Search_TargetEvents(HMC, HMC11y, HMC5d, dHMC, HTime, DTime, grid, Save):
scaling = HMC11y/min(HMC11y)
Index_thres1 = np.percentile(HMC,percentile)*scaling**(power)
#YearLow = DTime[np.argmin(scaling)]
#YearHigh = DTime[np.argmax(scaling)]
StormIndices = np.where(np.less(Index, Index_thres1))[0]
#print(np.percentile(HMC,percentile))
......
......@@ -86,8 +86,8 @@ def IndexDist(Time,YearsIndex,Storms,Kp_all,KpHours_all,HMC,Save):
HMC_step = 7.78
HMC_range = np.arange(10,480+HMC_step,HMC_step)
Index_range = Kp_range #HMC_range #Kp_range
step = Kp_step #HMC_step #Kp_step
Index_range = HMC_range #Kp_range
step = HMC_step #Kp_step
years = np.arange(1932,2016,1)
Percentages = np.zeros((len(years),len(Index_range)))
......@@ -123,11 +123,11 @@ def IndexDist(Time,YearsIndex,Storms,Kp_all,KpHours_all,HMC,Save):
HMCValues = -HMC[Storms[StormsBlock]]
KpValuesAll.extend(KpValues.tolist())
Hist[l,:], bin_edges = np.histogram(KpValues, bins = Index_range)
#KpValuesAll.extend(KpValues.tolist())
#Hist[l,:], bin_edges = np.histogram(KpValues, bins = Index_range)
#KpValuesAll.extend(HMCValues.tolist())
#Hist[l,:], bin_edges = np.histogram(HMCValues, bins = Index_range)
KpValuesAll.extend(HMCValues.tolist())
Hist[l,:], bin_edges = np.histogram(HMCValues, bins = Index_range)
Hist[l,:] *= 100/len(Storms1932)
......@@ -173,13 +173,13 @@ def IndexDist(Time,YearsIndex,Storms,Kp_all,KpHours_all,HMC,Save):
ax2.legend(loc=4,ncol=1,frameon=False,fontsize=10)
ax2.set_ylabel('Occurrence [%]')
ax2.set_xticklabels(np.arange(0,10,1))
ax2.set_xlabel('Kp')
ax2.set_ylim([0,10])
#ax2.set_xticklabels(np.arange(0,10,1))
#ax2.set_xlabel('Kp')
#ax2.set_ylim([0,10])
#ax2.set_xticklabels(np.around(Index_range[::3]/100,1))
#ax2.set_xlabel('-HMC/100 [nT]')
#ax2.set_ylim([0,20])
ax2.set_xticklabels(np.around(Index_range[::3]/100,1))
ax2.set_xlabel('-HMC/100 [nT]')
ax2.set_ylim([0,20])
#ax3.fill_between(np.arange(0,len(Kp_range)-1,1)+0.5,CumSum,CumSum2,color='maroon')
#ax3.set_ylim([0,105])
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment