Retains the legacy calculation for the solar declination from the Ault et al. (2015) Matlab code. This returns the day of year based on the start of climatological spring, meant for use with the function daylength. That is, a DOY=1 is March 1. The function itself uses the calendar day of the year (from 1-366).
Usage
soldec(DOY)
Arguments
DOY
Scalar - day of year from 1 to 366.
Value
Scalar - Spring climatological calendar day of year for the location and calendar day of year matching the appropriate solar declination angle.
References
Ault, T.R., R. Zurita-Miller and M. Schwarz, 2015: A Matlab<U+00A9> toolbox for calculating spring indices from daily meteorological data. Computers and Geosciences, 83, DOI: 10.1016/j.cageo.2015.06.015