link_to_sample = 200
sample_das = [201, 24048, 14834, 14268, 12097, 8035, 5022] # sample drainage areas
## GENERATE AND SET PRECIPITATION TIME SERIES
total_t = 1000.* 365.25 * 24 # 1,000 years of rainfall generated.
thresh = 0.5 # Intensity threshold (lower limit, mm/hr)
PD = PrecipitationDistribution(mean_storm_duration = 10.75,
mean_interstorm_duration=433.58,
mean_storm_depth = 9.62,
total_t=total_t)
storm_arr = np.array(PD.get_storm_time_series())
intensity_threshold, = np.where(storm_arr[:, 2] > thresh)
interstorm_durs = (storm_arr[intensity_threshold][:,0][np.arange(1,
len(storm_arr[intensity_threshold]) - 1)] -
storm_arr[intensity_threshold][:, 1][np.arange(0,
len(storm_arr[intensity_threshold]) - 2)])
durations = (storm_arr[intensity_threshold][:,1] -
storm_arr[intensity_threshold][:,0])
durations_s = [x * 3600. for x in durations]
intensities = storm_arr[intensity_threshold][:,2]
## SET INITIAL FIELDS, ARRAYS AND ITERATORS
lowRvar_PrecipDist = PrecipitationDistribution(mean_storm_duration = 11.75,
mean_interstorm_duration = 146.25,
mean_storm_depth = 4.775,
total_t=total_t)
highRvar_PrecipDist = PrecipitationDistribution(mean_storm_duration = 10.75,
mean_interstorm_duration = 433.58,
mean_storm_depth = 9.62,
total_t=total_t)
thresh = 0.5
# Get actual time series for the high Rvar case
highRvar_storm_arr = np.array(highRvar_PrecipDist.get_storm_time_series())
highRvar_intensity_threshold, = np.where(highRvar_storm_arr[:, 2] > thresh)
highRvar_durations = highRvar_storm_arr[highRvar_intensity_threshold][:,1]-highRvar_storm_arr[highRvar_intensity_threshold][:,0]
highRvar_durations_s = [x *3600. for x in highRvar_durations]
highRvar_intensities = highRvar_storm_arr[highRvar_intensity_threshold][:,2]
highRvar_depth = highRvar_durations * highRvar_intensities
highRvar_durations_raw = highRvar_storm_arr[:,1]-highRvar_storm_arr[:,0]
highRvar_intensities_raw = highRvar_storm_arr[:,2]
highRvar_depth_raw = highRvar_durations_raw * highRvar_intensities_raw
# Get actual time series for the low Rvar case
lowRvar_storm_arr = np.array(lowRvar_PrecipDist.get_storm_time_series())
lowRvar_intensity_threshold, = np.where(lowRvar_storm_arr[:, 2] > thresh)
lowRvar_durations = lowRvar_storm_arr[lowRvar_intensity_threshold][:,1]-lowRvar_storm_arr[lowRvar_intensity_threshold][:,0]
lowRvar_durations_s = [x *3600. for x in lowRvar_durations]
lowRvar_intensities = lowRvar_storm_arr[lowRvar_intensity_threshold][:,2]