*******************************************************************************
BOXMOX PBL diurnal evolution case
Christoph Knote - 11/2015
*******************************************************************************

Simulates the diurnal evolution of chemistry within an idealized planetary
boundary layer (PBL).

Emissions ("urban" levels) are emitted with a diurnal cycle. Photolysis rates,
temperature and the height of the PBL have a diurnal cycle as well. Photolysis
rates representative for a clear sky day at 40 deg N. For relevant species, a
(time-independent) deposition sink is included. During the morning, the rising
boundary layer mixes PBL air with air above ("X-"trainment). At nightfall, the
collapsing boundary layer reduces the mixing volume, and the effect of emissions
is felt much stronger.

This "respirating" PBL pattern is repeated for several day, and a "steady-state"
(i.e., no accumulation of species over several days) is observed.

The "box" we imagine has a 1x1 m2 horizontal extent, and its vertical dimension
is from the ground up to the top of the boundary layer. This represents an
always and instantaneously well-mixed boundary layer, the resulting
concentrations are hence representative for mean boundary layer values.

*******************************************************************************

Set up to be used with the MOZART-4 mechanism. MOZART-T1 should work as well.

*******************************************************************************

Files used:
 * Deposition.csv
 * Emissions.csv
 * Environment.csv
 * InitialConditions.csv
 * PhotolysisRates.csv
 * Xtrainment.csv

And the BOXMOX.nml namelist.

plot.R is an R script (https://www.r-project.org) to make a nice plot of the
result ("pbl_diurnal_cycle.png"). Usage:

> Rscript plot.R

*******************************************************************************
