first try adding watches disclaimer to plot

hyswap/plots.py

@@ -235,6 +235,7 @@ def plot_raster_hydrograph(df_formatted, ax=None,
 def plot_duration_hydrograph(percentiles_by_day, df, data_col, doy_col,
                              pct_list=[0, 5, 10, 25, 75, 90, 95, 100],
                              data_label=None, ax=None,
+                             disclaimer=False,
                              title="Duration Hydrograph",
                              ylab="Discharge, in Cubic Feet per Second",
                              xlab="Month", colors=None, **kwargs):
@@ -259,6 +260,11 @@ def plot_duration_hydrograph(percentiles_by_day, df, data_col, doy_col,
     ax : matplotlib.axes.Axes, optional
         Axes to plot on. If not provided, a new figure and axes will be
         created.
+    disclaimer : bool, optional
+        If True, displays the disclaimer 'For some streams, flow
+        statistics may have been computed from mixed regulated
+        and unregulated flows; this can affect depictions of flow
+        conditions.' below the x-axis.
     title : str, optional
         Title for the plot. If not provided, the default title will be
         'Duration Hydrograph'.
@@ -313,6 +319,11 @@ def plot_duration_hydrograph(percentiles_by_day, df, data_col, doy_col,
         label = df[data_col].name
     else:
         label = data_label
+    # Add disclaimer if True
+    if disclaimer is True:
+        txt = 'For some streams, flow statistics may have been computed from mixed regulated \nand unregulated flows; this can affect depictions of flow conditions.' # noqa: E501
+    else:
+        txt = ''
     # get colors
     if colors is None:
         colors = ["#e37676", "#e8c285", "#dbf595", "#a1cc9f",
@@ -354,6 +365,8 @@ def plot_duration_hydrograph(percentiles_by_day, df, data_col, doy_col,
     ax.set_ylabel(ylab)
     ax.set_yscale("log")
     ax.set_title(title)
+    # disclaimer
+    ax.text(0.5, 0.05, txt, ha='center', color='red')
     # get y-axis ticks and convert to comma-separated strings
     yticks = ax.get_yticks()
     yticklabels = [f'{int(y):,}' for y in yticks]