"""
Time-drift correlation
=======================

This example shows how to detect timestamp drift in irradiance
measurements using
:py:func:`solarpy.plotting.plot_time_drift_correlation`.
"""

# %%
import pvlib

import solarpy

# Read a year of 1-minute GHI measurements
data, meta = solarpy.iotools.read_t16(
    "https://raw.githubusercontent.com/AssessingSolar/solarpy/refs/heads/main/data/LYN_2023.csv",  # noqa: E501
    map_variables=True,
)

location = pvlib.location.Location(meta["latitude"], meta["longitude"])
clearsky = location.get_clearsky(data.index)
is_clearsky = pvlib.clearsky.detect_clearsky(
    data["ghi"], clearsky["ghi"], data.index
)

fig, ax = solarpy.plotting.plot_time_drift_correlation(
    times=data.index,
    ghi=data["ghi"],
    ghi_clear=clearsky["ghi"],
    is_clearsky=is_clearsky,
)