diff --git a/data/Fritz_anav.rda b/data/Fritz_anav.rda
new file mode 100644
index 0000000..b48786a
Binary files /dev/null and b/data/Fritz_anav.rda differ
diff --git a/data/SKIRTOR_anav.rda b/data/SKIRTOR_anav.rda
new file mode 100644
index 0000000..57cc656
Binary files /dev/null and b/data/SKIRTOR_anav.rda differ
diff --git a/man/Fritz_avan.Rd b/man/Fritz_avan.Rd
new file mode 100644
index 0000000..ab900a4
--- /dev/null
+++ b/man/Fritz_avan.Rd
@@ -0,0 +1,82 @@
+\name{Fritz_anav}
+\alias{Fritz_anav}
+\docType{data}
+\title{
+Fritz 2006 AGN Template Data
+}
+\description{
+The full extended Fritz 2006 AGN model as per Feltre 2012, for the averaged viewing angle. This is then interpolated by \code{AGNinterp} to create arbitrary AGN at non gridded points.
+}
+\usage{data("Fritz_anav")}
+\format{
+  The format is:\cr
+List of 8\cr
+ $ Wave  : num [1:178] 10 11.2 12.6 14.1 15.8 ...\cr
+ $ Aspec   : num [1:3, 1:4, 1:5, 1:8, 1:5, 1:178] 0.0672 0.0698 0.0727 0.064 0.064 ...\cr
+  ..- attr(*, "dimnames")=List of 6\cr
+  .. ..$ : chr [1:3] "60" "100" "140"\cr
+  .. ..$ : chr [1:4] "0" "2" "4" "6"\cr
+  .. ..$ : chr [1:5] "-1" "-0.75" "-0.5" "-0.25" ...\cr
+  .. ..$ : chr [1:8] "0.1" "0.3" "0.6" "1" ...\cr
+  .. ..$ : chr [1:5] "10" "30" "60" "100" ...\cr
+  .. ..$ : NULL\cr
+ $ params: chr [1:5] "ct" "al" "be" "ta" ...\cr
+ $ ct    : num [1:3] 60 100 140\cr
+ $ al    : num [1:4] 0 2 4 6\cr
+ $ be    : num [1:5] -1 -0.75 -0.5 -0.25 0\cr
+ $ ta    : num [1:8] 0.1 0.3 0.6 1 2 3 6 10\cr
+ $ rm    : num [1:5] 10 30 60 100 150\cr
+}
+\details{
+The below is taken directly from the included documentation:
+
+A composition of power-laws, with different spectral indexes as a function of the wavelength, was used to model the primary source (i.e. accretion disk SED) SED taken from the file pwlbb05.dat.
+
+Standard, Galactic composition of dust, with graphite and astronomical silicates. The grain size distribution is a power law with index -3.5. Scattering and absorption coefficients taken from Laor and Draine (1993). Number of different species: 2 (Graphite and Silicate only).
+
+The inner radius is set by the sublimation temperature of graphite grains: T = 1500 K. The dust sublimation radius is Rmin = 1.3*L46^0.5*T^(-2.8) pc [Barvainis (1987)] where L46=L/(1e46 erg/s), and T is the sublimation temperature. The dust density throughout the torus is expressed by the following: rho(r,theta)=A*r^(beta)exp(-gamma*cos(theta)), where A is a normalization costant, beta and gamma are, in the name of the model, be and al, respectively.
+
+Each AGN model has been normalised so that the total luminosity fo the source AGN integrates to 1.
+
+\describe{
+\item{Wave}{
+  The wavelength in Angstroms [178 elements].
+}
+  \item{Aspec}{
+  6 dimensional array containing (in order) dimension of \option{ct} (opening angle of torus in degrees); \option{al} (gamma dust parameter); \option{be} (beta dust parameter); \option{ta} (optical depth tau); \option{rm} (outer to inner torus radius ratio), and the final column is the normalised luminosity per unit wavelength [178 elements].
+}
+\item{params}{
+  Character vector; order list of parameters in \option{Aspec} model array.
+}
+\item{ct}{
+  Numeric vector; available grid values of \option{ct}. Opening angle of torus in degrees.
+}
+\item{al}{
+  Numeric vector; available grid values of \option{al}. Gamma dust parameter.
+}
+\item{be}{
+  Numeric vector; available grid values of \option{be}. Beta dust parameter.
+}
+\item{ta}{
+  Numeric vector; available grid values of \option{ta}. Optical depth tau.
+}
+\item{rm}{
+  Numeric vector; available grid values of \option{rm}. Outer to inner torus radius ratio.
+}
+}
+}
+\source{
+Downloaded the file ct20/ct40/ct60 from http://www.crya.unam.mx/gente/j.fritz/JFhp/AGN_models.html.
+}
+\references{
+Feltre et al, 2012, MNRAS, 426, 120 \cr
+Fritz et al, 2006, MNRAS, 366, 767
+}
+\examples{
+data(Fritz_anav)
+str(Fritz_anav)
+
+plot(Fritz_anav$Wave, Fritz_anav$Aspec["100","4","-0.5","2","60",],
+type='l', log='xy', xlab='Wave / Ang', ylab='Lum / L/A')
+}
+\keyword{datasets}
diff --git a/man/SKIRTOR_avan.Rd b/man/SKIRTOR_avan.Rd
new file mode 100644
index 0000000..870673a
--- /dev/null
+++ b/man/SKIRTOR_avan.Rd
@@ -0,0 +1,31 @@
+\docType{data}
+\name{SKIRTOR_anav}
+\alias{SKIRTOR_anav}
+\title{SKIRTOR angle-averaged AGN Template}
+\source{
+Stalevski et al, 2012, MNRAS, 420, 2756 \cr
+Stalevski et al, 2016, MNRAS, 458, 2288
+}
+\usage{data("SKIRTOR_anav")}
+
+\description{
+  A data frame representing the SKIRTOR AGN template used for SED fitting, viewing angle averaged.
+}
+
+\format{
+A data frame containing the SKIRTOR model. The format is:
+
+\describe{
+  \item{Wave}{Numeric vector; wavelengths in microns.}
+  \item{Aspec}{Numeric array [1:5, 1:4, 1:4, 1:8, 1:3, 1:1, 1:132]; the SED grid. Has named dimensions corresponding to model parameters.}
+  \item{params}{Character vector of parameter names: \code{"ta"}, \code{"p"}, \code{"q"}, \code{"ct"}, \code{"rm"}, \code{"Mcl"}.}
+  \item{ta}{Numeric scalar; Average edge-on optical depth at 9.7 micron.}
+  \item{p}{Numeric scalar; Index for the radial density gradient in the SKIRTOR AGN model.}
+  \item{q}{Numeric scalar; Index for the angular density gradient in the SKIRTOR AGN model.}
+  \item{ct}{Numeric scalar; Opening angle of torus in degrees.}
+  \item{rm}{Numeric scalar; Ratio of outer radius to inner radius in the SKIRTOR AGN model.}
+  \item{Mcl}{Numeric scalar; clumpiness parameter.}
+}
+}
+
+\keyword{datasets}