I think there's a typo in the get_gage_q function in the Validation module (validation_confluence.py, line 217). If there is more than one index (which I'm guessing is rare if it happens), and there is a model_q value, you use the length of the time series to compare to this model_q, rather than the mean_q values (you compare it to glt rather than gmq). The line is:
index=np.array(index[np.argmin(np.abs(np.array(glt)-model_q))])
But I think it should be:
index=np.array(index[np.argmin(np.abs(np.array(gmq)-model_q))])
I think there's a typo in the get_gage_q function in the Validation module (validation_confluence.py, line 217). If there is more than one index (which I'm guessing is rare if it happens), and there is a model_q value, you use the length of the time series to compare to this model_q, rather than the mean_q values (you compare it to glt rather than gmq). The line is:
index=np.array(index[np.argmin(np.abs(np.array(glt)-model_q))])
But I think it should be:
index=np.array(index[np.argmin(np.abs(np.array(gmq)-model_q))])