"""
Capacitance sensor Calibration
"""
# 2972
y_cal = np.array([100, 500, 800, 1200, 1800, 2400, 3000, 3600, 3700])
x_cal = np.array([1901, 2176, 2393, 2668, 3095, 3496, 3914, 4330, 4403])
a_stage = cd.polyfit(x_cal, y_cal, 1)
coeff_cal = a_stage['polynomial']
slope = coeff_cal[0]
intercept = coeff_cal[1]
print coeff_cal
"""
read tank data
"""
block_1 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/2972/2972_015_001.CSV'
water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept)
block_2 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/2972/2972_015_002_22_08_2015.CSV'
water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept)
block_3 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/2972/2972_010_001.CSV'
water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept)
block_4 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/2972/2972_010_002.CSV'
water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept)
block_5 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/2972/2972_007_001_24_12_2015.CSV'
water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept)
# remove error value when logger is removed for downloading
# http://stackoverflow.com/a/16613835/2632856
water_level_4.drop(pd.Timestamp('2015-09-10 12:00:00'), inplace=True, axis=0)
# drop last and first value
for i in range(1, 6, 1):
eval(
__author__ = 'Tejas'
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import checkdam.checkdam as cd
from datetime import timedelta
"""
Read Check dam data
"""
block_1 = 'E:\Tejas\JK\\a.csv'
water_level_1 = cd.read_correct_ch_dam_data(block_1, 0.922, -1893)
block_2 = 'E:\Tejas\JK\\b.csv'
water_level_2 = cd.read_correct_ch_dam_data(block_2, 0.922, -1893)
block_3 = 'E:\Tejas\JK\\c.csv'
water_level_3 = cd.read_correct_ch_dam_data(block_2, 0.922, -1893)
for i in range(1, 3, 1):
eval("water_level_{0}.drop(water_level_{0}.tail(1).index, inplace=True, axis=0)".format(i))
eval("water_level_{0}.drop(water_level_{0}.head(1).index, inplace=True, axis=0)".format(i))
# for i in range(1, 11, 1):
# print "water_level_{0}".format(i)
# print eval("water_level_{0}.head()".format(i))
fig = plt.figure()
for i in range(1, 3, 1):
x = eval("water_level_{0}.index".format(i))
y = eval("water_level_{0}['stage(m)']".format(i))
plt.plot(x, y)
import matplotlib.colors as colors
from datetime import timedelta
import math
# latex parameters
rc('font', **{'family': 'sans-serif', 'sans-serif': ['Helvetica']})
rc('text', usetex=True)
plt.rc('text', usetex=True)
plt.rc('font', family='serif', size=18)
"""
Read Check dam data
"""
block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_001.CSV'
water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept)
block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_002.CSV'
water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept)
block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_003.CSV'
water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept)
block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_004.CSV'
water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept)
block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_005.CSV'
water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept)
block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_006.CSV'
water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept)
block_7 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_004_001.CSV'
water_level_7 = cd.read_correct_ch_dam_data(block_7, slope, intercept)
# water_level_7['stage(m)'] += 0.05
block_8 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_004_002_12_12_2014.CSV'
water_level_8 = cd.read_correct_ch_dam_data(block_8, slope, intercept)
# drop date time column
rain_df = rain_df.drop('Date_Time', 1)
"""
Check dam calibration
"""
y_cal = np.array([100, 400, 1000, 1600, 2250, 2750])
x_cal = np.array([1987, 2454, 3344, 4192, 5104, 5804])
a_stage = cd.polyfit(x_cal, y_cal, 1)
coeff_cal = a_stage['polynomial']
slope = coeff_cal[0]
intercept = coeff_cal[1]
"""
Read Check dam data
"""
block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_007_001.CSV'
water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept)
block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_007_002_25_8_14.CSV'
water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept)
block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_007_003.CSV'
water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept)
block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_007_004.CSV'
water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept)
block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_005_001.CSV'
water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept)
block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2526/2526_005_002_03_12_2014.CSV'
water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept)
for i in range(1, 7, 1):
eval("water_level_{0}.drop(water_level_{0}.tail(1).index, inplace=True, axis=0)".format(i))
eval("water_level_{0}.drop(water_level_{0}.head(1).index, inplace=True, axis=0)".format(i))
# raise SystemExit(0)
"""
Check dam calibration
"""
y_cal = np.array([100, 500, 1000, 2000, 2500, 3000])
x_cal = np.array([2018, 2761, 3492, 4924, 5609, 6320])
a_stage = cd.polyfit(x_cal, y_cal, 1)
coeff_cal = a_stage['polynomial']
slope = coeff_cal[0]
intercept = coeff_cal[1]
"""
Read Check dam data
"""
block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_001.CSV'
water_level_1 = cd.read_correct_ch_dam_data(block_1,
slope,
intercept,
stage_cutoff=stage_cutoff)
block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_002_25_8_14.CSV'
water_level_2 = cd.read_correct_ch_dam_data(block_2,
slope,
intercept,
stage_cutoff=stage_cutoff)
block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_003.CSV'
water_level_3 = cd.read_correct_ch_dam_data(block_3,
slope,
intercept,
stage_cutoff=stage_cutoff)
block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_004.CSV'
water_level_4 = cd.read_correct_ch_dam_data(block_4,
slope,
intercept,
__author__ = 'kiruba' import numpy as np import matplotlib.pyplot as plt import pandas as pd import itertools import checkdam.checkdam as cd # calibration y_cal = np.array([100, 1000, 2000, 3000, 4000, 5000]) x_cal = np.array([1875, 2516, 3212, 3901, 4605, 5280]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] interecept = coeff_cal[1] # depth correction based on capacitance block_1 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_012_002_21_08_2015.CSV' block_2 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_015_001.CSV' block_3 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_010_001_08_01_2015.CSV' block_1_df = cd.read_correct_ch_dam_data(block_1, calibration_slope=slope, calibration_intercept=interecept) block_2_df = cd.read_correct_ch_dam_data(block_2, calibration_slope=slope, calibration_intercept=interecept) block_3_df = cd.read_correct_ch_dam_data(block_3, calibration_slope=slope, calibration_intercept=interecept) fig = plt.figure() plt.plot(block_1_df.index, block_1_df['stage(m)'], 'g-') plt.plot(block_2_df.index, block_2_df['stage(m)'], 'g-') plt.plot(block_3_df.index, block_3_df['stage(m)'], 'g-') plt.show()
rain_df = rain_df.drop('Date_Time', 1)
"""
Check dam calibration
"""
y_cal = np.array([100, 400, 1000, 1500, 2000, 3000])
x_cal = np.array([1971, 2336, 3083, 3720, 4335, 5604])
a_stage = cd.polyfit(x_cal, y_cal, 1)
coeff_cal = a_stage['polynomial']
slope = coeff_cal[0]
intercept = coeff_cal[1]
"""
Read Check dam data
"""
block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_005_001.CSV'
water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept)
# block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_003_008.CSV'
# water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept)
# water_level_2 = water_level_2[ :"2014-08-23"]
# block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_001.CSV'
# water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept)
# block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_002_12_12_2014.CSV'
# water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept)
# block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_003_24_12_2014.CSV'
# water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept)
# block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_004_24_1_2014.CSV'
# water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept)
# fig = plt.figure()
# for i in range(1, 7, 1):
# x = eval("water_level_{0}.index".format(i))
# raise SystemExit(0) """ Check dam calibration """ y_cal = np.array([100, 500, 1000, 2000, 2500, 3000]) x_cal = np.array([2018, 2761, 3492, 4924, 5609, 6320]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] intercept = coeff_cal[1] """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept, stage_cutoff=stage_cutoff) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_002_25_8_14.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept, stage_cutoff=stage_cutoff) block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_003.CSV' water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept, stage_cutoff=stage_cutoff) block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_006_004.CSV' water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept, stage_cutoff=stage_cutoff) block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_002_001.CSV' water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept, stage_cutoff=stage_cutoff) block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_002_002_03_12_2014.CSV' water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept, stage_cutoff=stage_cutoff) block_7 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_002_003_11_12_2014.CSV' water_level_7 = cd.read_correct_ch_dam_data(block_7, slope, intercept, stage_cutoff=stage_cutoff) block_8 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_002_004_23_12_2014.CSV' water_level_8 = cd.read_correct_ch_dam_data(block_8, slope, intercept, stage_cutoff=stage_cutoff) block_9 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2510/2510_002_005_3_1_2015.CSV'
y_cal_1 = [100, 400, 1000, 1500, 2000, 3000] x_cal_1 = [1971, 2336, 3083, 3720, 4335, 5604] coeff = np.polyfit(x_cal_1, y_cal_1, 1) print coeff # print coeff[0]*3020 + coeff[1] # print coeff[0]*2810 + coeff[1] # raise SystemExit(0) fig = plt.figure() plt.plot(x_cal_1, y_cal_1, 'ro-') plt.show() raise SystemExit(0) """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_002.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept) block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_003.CSV' water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept) block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_004.CSV' water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept) block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_005.CSV' water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept) block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_008_006.CSV' water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept) block_7 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_004_001.CSV' water_level_7 = cd.read_correct_ch_dam_data(block_7, slope, intercept) # water_level_7['stage(m)'] += 0.05 block_8 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2525/2525_004_002_12_12_2014.CSV' water_level_8 = cd.read_correct_ch_dam_data(block_8, slope, intercept)
# drop date time column
rain_df = rain_df.drop('Date_Time', 1)
"""
Check dam calibration
"""
y_cal = np.array([100, 400, 1000, 1500, 2000, 3000])
x_cal = np.array([1971, 2336, 3083, 3720, 4335, 5604])
a_stage = cd.polyfit(x_cal, y_cal, 1)
coeff_cal = a_stage['polynomial']
slope = coeff_cal[0]
intercept = coeff_cal[1]
"""
Read Check dam data
"""
block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_005_001.CSV'
water_level_1 = cd.read_correct_ch_dam_data(block_1, slope, intercept)
# block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_003_008.CSV'
# water_level_2 = cd.read_correct_ch_dam_data(block_2, slope, intercept)
# water_level_2 = water_level_2[ :"2014-08-23"]
# block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_001.CSV'
# water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept)
# block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_002_12_12_2014.CSV'
# water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept)
# block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_003_24_12_2014.CSV'
# water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept)
# block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2665/2665_007_004_24_1_2014.CSV'
# water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept)
# fig = plt.figure()
# for i in range(1, 7, 1):
# x = eval("water_level_{0}.index".format(i))
x_cal_1 = np.array([1894, 2563, 3298, 4049, 4794, 5548]) a_stage_1 = cd.polyfit(x_cal_1, y_cal_1, 1) coeff_cal_1 = a_stage_1['polynomial'] slope_1 = coeff_cal_1[0] intercept_1 = coeff_cal_1[1] y_cal = np.array([100, 1000, 2000, 3000, 4000, 5000]) x_cal = np.array([1864, 2540, 3313, 4078, 4835, 5582]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] intercept = coeff_cal[1] """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/3075/3075_012_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope_1, intercept_1) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/3075/3075_012_00325_8_14.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope_1, intercept_1) # block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_001.CSV' # water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept) # block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_002_03_12_2014.CSV' # water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept) # block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_003_11_12_2014.CSV' # water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept) # block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_004_23_12_2014.CSV' # water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept) # block_7 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_005_03_01_2015.CSV' # water_level_7 = cd.read_correct_ch_dam_data(block_7, slope, intercept) # Take data upto july 17, reason : mouse bite water_level_2 = water_level_2[:'2014-07-17'] # print water_level_2.tail()
a_stage_1 = cd.polyfit(x_cal_1, y_cal_1, 1) coeff_cal_1 = a_stage_1['polynomial'] slope_1 = coeff_cal_1[0] intercept_1 = coeff_cal_1[1] y_cal = np.array([100, 1000, 2000, 3000, 4000, 5000]) x_cal = np.array([1864, 2540, 3313, 4078, 4835, 5582]) a_stage = cd.polyfit(x_cal, y_cal, 1) coeff_cal = a_stage['polynomial'] slope = coeff_cal[0] intercept = coeff_cal[1] """ Read Check dam data """ block_1 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/3075/3075_012_001.CSV' water_level_1 = cd.read_correct_ch_dam_data(block_1, slope_1, intercept_1) block_2 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/3075/3075_012_00325_8_14.CSV' water_level_2 = cd.read_correct_ch_dam_data(block_2, slope_1, intercept_1) # block_3 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_001.CSV' # water_level_3 = cd.read_correct_ch_dam_data(block_3, slope, intercept) # block_4 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_002_03_12_2014.CSV' # water_level_4 = cd.read_correct_ch_dam_data(block_4, slope, intercept) # block_5 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_003_11_12_2014.CSV' # water_level_5 = cd.read_correct_ch_dam_data(block_5, slope, intercept) # block_6 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_004_23_12_2014.CSV' # water_level_6 = cd.read_correct_ch_dam_data(block_6, slope, intercept) # block_7 = '/media/kiruba/New Volume/ACCUWA_Data/check_dam_water_level/2976/2976_019_005_03_01_2015.CSV' # water_level_7 = cd.read_correct_ch_dam_data(block_7, slope, intercept) # Take data upto july 17, reason : mouse bite water_level_2 = water_level_2[:'2014-07-17'] # print water_level_2.tail()
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
import itertools
import checkdam.checkdam as cd
# calibration
y_cal = np.array([100, 1000, 2000, 3000, 4000, 5000])
x_cal = np.array([1875, 2516, 3212, 3901, 4605, 5280])
a_stage = cd.polyfit(x_cal, y_cal, 1)
coeff_cal = a_stage['polynomial']
slope = coeff_cal[0]
interecept = coeff_cal[1]
# depth correction based on capacitance
block_1 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_012_002_21_08_2015.CSV'
block_2 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_015_001.CSV'
block_3 = '/media/kiruba/New Volume/ACCUWA_Data/lake_water_level/3055/3055_010_001_08_01_2015.CSV'
block_1_df = cd.read_correct_ch_dam_data(block_1,
calibration_slope=slope,
calibration_intercept=interecept)
block_2_df = cd.read_correct_ch_dam_data(block_2,
calibration_slope=slope,
calibration_intercept=interecept)
block_3_df = cd.read_correct_ch_dam_data(block_3,
calibration_slope=slope,
calibration_intercept=interecept)
fig = plt.figure()
plt.plot(block_1_df.index, block_1_df['stage(m)'], 'g-')
plt.plot(block_2_df.index, block_2_df['stage(m)'], 'g-')
plt.plot(block_3_df.index, block_3_df['stage(m)'], 'g-')
plt.show()
