call daily_mean_solar (lat, time, cosz, fracday, rrsun)
call daily_mean_solar (lat, time_since_ae, cosz, fracday, rrsun)
call daily_mean_solar (lat, time, cosz, solar)
call daily_mean_solar (lat, time_since_ae, cosz, solar)

The first option (used in conjunction with time_manager_mod) generates the real variable time_since_ae from the time_type input time, and then calls daily_mean_solar with this real input (option 2). The third and fourth options correspond to the first and second and are used with then spectral 2-layer model, where only cosz and solar are desired as output. These routines generate dummy arguments and then call option 2, where the calculation is done.

The time of year is set by

real, intent(in) :: time_since_ae

With time_type input, the time of year is extracted from

type(time_type), intent(in) :: time

Separate routines exist within this interface for scalar, 1D or 2D input and output fields:

real, intent(in), dimension(:,:) :: lat
real, intent(in), dimension(:)   :: lat
real, intent(in)                 :: lat
real, intent(out), dimension(:,:) :: cosz, fracday
real, intent(out), dimension(:)   :: cosz, fracday
real, intent(out)                 :: cosz, fracday

Parameters

[in]	<lat>	Latitudes of model grid points [radians]
[in]	<time_since_ae>	Time of year; autumnal equinox = 0.0, one year = 2 * pi [radians]
[in]	<time>	Time at which astronomical values are desired (time_type variable) [seconds, days]
[out]	<cosz>	Cosine of solar zenith angle, set to zero when entire period is in darkness [dimensionless]
[out]	<fracday>	Daylight fraction of time interval [dimensionless]
[out]	<rrsun>	Earth-Sun distance (r) relative to semi-major axis of orbital ellipse (a):(a/r)**2 [dimensionless]
[out]	<solar>	shortwave flux factor: cosine of zenith angle * daylight fraction / (earth-sun distance squared) [dimensionless]