flowchart TD
    A(["ts2clm()"]) --> B[/INPUT:<br>date, temperature/]
    B --> C{Checks?}
    C --> |Yes| D[/OUTPUT:<br>date, temperature/]
    C --> |No| E([End])
    D --> F(["make_whole_fast()"])
    
    G(["make_whole_fast()"]) --> H[/INPUT:<br>date, temperature/]
    H --> I[/OUTPUT:<br>doy, date, temperature/]
    I --> J(["clim_spread()"])
    
    K(["clim_spread()"]) --> L[/INPUT:<br>doy, date, temperature/]
    L ---> M[spread doy as rows<br>spread year as cols<br>grow doy by windowHalfWidth]
    M --> N[/OUTPUT:<br>temperature in<br>doy x year matrix/]
    N --> O(["clim_calc_cpp()"])
    
    P(["clim_calc_cpp()"]) --> Q[/INPUT:<br>temperature in<br>doy x year matrix/]
    Q --> 

ts2clm()

ts2clm() accepts a dataframe with date (x = t) and temperature (y = temp). Additional function arguments include:

1. climatologyPeriod Required. To this argument should be passed two values (see example below). The first value should be the chosen date for the start of the climatology period, and the second value the end date of said period. This chosen period (preferably 30 years in length) is then used to calculate the seasonal cycle and the extreme value threshold.
2. maxPadLength Specifies the maximum length of days over which to interpolate (pad) missing data (specified as NA) in the input temperature time series; i.e., any consecutive blocks of NAs with length greater than maxPadLength will be left as NA. The default is FALSE. Set as an integer to interpolate. Setting maxPadLength to TRUE will return an error.
3. windowHalfWidth Width of sliding window about day-of-year (to one side of the center day-of-year) used for the pooling of values and calculation of climatology and threshold percentile. Default is 5 days, which gives a window width of 11 days centred on the 6th day of the series of 11 days.
4. pctile Threshold percentile (%) for detection of events (MHWs). Default is 90th percentile. Should the intent be to use these threshold data for MCSs, set pctile = 10. Or some other low value.
5. smoothPercentile Boolean switch selecting whether to smooth the climatology and threshold percentile time series with a moving average of smoothPercentileWidth. Default is TRUE.
6. clmOnly Choose to calculate and return only the climatologies. The default is FALSE.
7. var This argument has been introduced to allow the user to choose if the variance of the seasonal signal per doy should be calculated. The default of FALSE will prevent the calculation, potentially increasing speed of calculations on gridded data and reducing the size of the output. The variance was initially introduced as part of the standard output from Hobday et al. (2016), but few researchers use it and so it is generally regarded now as unnecessary.
8. roundClm This argument allows the user to choose how many decimal places the seas and thresh outputs will be rounded to. Default is 4. To prevent rounding set roundClm = FALSE. This argument may only be given numeric values or FALSE.