SPARTACUS VEGETATION

Robin Hogan <r.j.hogan@ecmwf.int>
December 2017
http://www.met.reading.ac.uk/clouds/spartacus

This package provides a Matlab implementation of the "SPARTACUS"
algorithm for computing the interaction of shortwave (i.e. solar)
radiation with vegetation canopies, described in the following paper:

  Hogan, R. J., T. Quaife and R. Braghiere, 2017: Fast matrix
  treatment of 3D radiative transfer in vegetation canopies:
  SPARTACUS-Vegetation 1.1. Submitted to Geosci. Model Dev.

The software is released under the terms of the Apache License,
version 2.0 (see LICENSE file).  The Matlab files present are:

  do_all_plots.m - script to produce Figs. 2-5 of the paper (calls
    plot_sw_vegetation.m)

  plot_sw_vegetation.m - script to produce one figure (calls
    calc_sw_vegetation.m and plot_ref.m)

  plot_ref.m - function to plot reference Monte Carlo results for
    comparison with SPARTACUS

  calc_sw_vegetation.m - function to compute shortwave radiation
    properties of a vegetation canopy with the assumptions of the
    RAMI4PILPS intercomparison exercise (calls
    spartacus_sw_vegetation.m)

  spartacus_sw_vegetation.m - function that implements the SPARTACUS
    algorithm for an arbitrary number of layers and returns the fluxes
    at each layer interface (calls fast_expm.m)

  fast_expm.m - function that computes the matrix-exponential rapidly
    for 2x2 and 3x3 matrices with a specific form that arises in
    layers containing no vegetation

The RAMI4PILPS directory contains reference Monte Carlo results kindly
provided by Jean-Luc Widlowski, described in the following paper and
web site:

  Widlowski, J-L., B. Pinty, M. Clerici, Y. Dai, M. De Kauwe, K. de
  Ridder, A. Kallel, H. Kobayashi, T. Lavergne, W. Ni-Meister,
  A. Olchev, T. Quaife, S. Wang, W. Yang, Y. Yang and H. Yuan, 2011:
  RAMI4PILPS: An intercomparison of formulations for the partitioning
  of solar radiation in land surface models. J. Geophys. Res., 116,
  G02019, 25, DOI: 10.1029/2010JG001511.

http://rami-benchmark.jrc.ec.europa.eu/HTML/RAMI4PILPS/RAMI4PILPS.php
