=======================
Version 1.0.0

2015/01
Delivering first working version to J.-L. Dufresne, LMD.


=======================
Version 1.1.0

2016/09
+ added: switches in "options.in" to control the separated computation of the following quantities:
H2O CIA, CO2 CIA, N2-N2 CIA, O2-O2 CIA, H2O AER continuum
+ added: on-screen messages to let the user know precisely what additionnal source of
opacity is computed or not (and why), and whether it is recorded or not.
+ added: switches in "options.in" to control the separate recording of the following quantities:
Rayleigh scattering, H2O CIA, CO2 CIA, N2-N2 CIA, O2-O2 CIA, H2O AER continuum, total optical thickness
+ added: scripts in the "results" folder in order to plot every opacity source, both as a contour plot
(as a function of wavenumber and altitude) and cumulated over the column (as a function of wavenumber).

=======================
Version 1.2.0

2016/10
+ corrected: adding the names of several record files in the "init" subroutine so that they are erased
at the beginning of a new computation.
+ corrected: the computation of the "omega" function that records the fraction of the atmosphere-space
exchange that is due to each gas layer (contribution from the ground was previously taken into account).
+ corrected: the computation of the emission altitude
+ added: the computation of the emission blackbody intensity, and the corresponding brightness temperature.
+ added: file results/ES_signal.txt that records the flux emitted by the ground, the ground-space net exchange,
the total optical depth between the ground and space, and the transmittivity between ground and space,
computed both analytically and from radiative transfer results.
+ modified: all visualization scripts now produce pdf files.

=======================
Version 1.3.0

2016/10
+ corrected: mean layer altitude and mean layer temperature have been set resp. to the altitude and the temperature
at the middle (geometric center) of the layer.
+ corrected: use of partial CO2 pressures instead of total atmospheric pressure when computing CO2 CIA
optical thickness.
+ corrected: a bug related to the detection of wavenumber intervals in the routine that computes
H2O AER continuum; was completely crashing the program.
+ added: a visualization script to plot the temperature profile
+ added: the computation of upward, downward and net fluxes (spectral signal and spectrally integrated)
at a given probe pressure level.
+ added: the possibility to specify a probe pressure level in "data.in"
+ added: a visualization script in order to plot the spectral signals of upward, downward and net fluxes
at the probe pressure level.

=======================
Version 1.4.0

2016/11
+ corrected: a bug introduced in the previous version. Temperatures at each atmospheric level were messed up.
+ corrected: detection of wavenumber intervals bug is still present. Changed previous modifications for more radical scheme.
+ corrected: a bug in "CIA_optical_thickness" subroutine, causing possible negative optical thicknesses
+ modified: the way optical thickness files are recorded, so that they are not overwritten in case of a multi-pass computation.

=======================
Version 1.5.0

2016/11
+ added: two parameters have been added in the "data.in" file, that let the user specify the ground and space temperatures.
These temperatures are now distinct from the temperature of the gas at ground and space levels, with the consequences
that a discontinuity of temperature between the surfaces and the atmosphere is possible, and that the user does no
longuer have to mess with the "composition.in" file in order to specify the temperature of space at the last level.
+ modified: the "data.in" input file: ground emissivity has been removed
+ added: a new input data file: data/hires_spectra/profile_parameters.in that contain the values of the ground temperature,
space temperature and ground emissivity. 

=======================
Version 1.6.0

2016/11
+ added: computation of average layer calorific capacities; the "data.in" has one more parameter.
+ added: the total radiative budget of each layer (both spectrally integrated and per wavelength) is
now also recorded in terms of cooling rates (K/day). A positive value stands for a heating of the layer.
+ added: a script to visualize the spectrally integrated total radiative budget as a function of altitude.
+ record files that have been added:
- results/cooling_rate_signal.txt: total cooling rate (K/day) as a function of altitude and wavenumber. Positive means heating.
- results/cooling_rate.txt: spectrally integrated total cooling rate (K/day) as a function of altitude.  Positive means heating.
- results/gas_ground_coolingrate_signal.txt: part of the total cooling rate (K/day) that is due to radiative exchanges
between the gas and ground, as  a function of altitude and wavenumber. Positive means heating. 
- results/gas_ground_coolingrate.txt: part of the spectrally integrated total cooling rate (K/day) that is due to radiative exchanges
between the gas and ground, as  a function of altitude. Positive means heating. 
- results/gas_space_coolingrate_signal.txt: part of the total cooling rate (K/day) that is due to radiative exchanges
between the gas and space, as  a function of altitude and wavenumber. Positive means heating. 
- results/gas_space_coolingrate.txt: part of the spectrally integrated total cooling rate (K/day) that is due to radiative exchanges
between the gas and space, as  a function of altitude. Positive means heating. 
- results/gas_gas_coolingrate_signal.txt: part of the total cooling rate (K/day) that is due to radiative exchanges
between the gas and the rest of the atmosphere, as  a function of altitude and wavenumber. Positive means heating. 
- results/gas_gas_coolingrate.txt: part of the spectrally integrated total cooling rate (K/day) that is due to radiative exchanges
between the gas and the rest of the atmosphere, as  a function of altitude. Positive means heating.

=======================
Version 2.0

2017/01
+ modified: all steps involving the computation of a additionnal optical thickness have been parallelized.
+ corrected: a bug that prevented the recording of total optical thicknesses when no additionnal opacity was computed.
+ added: program "flux_differences.exe" in the "results" folder. This program will compute the difference in flux signals
at the ground, at the probe level and at the top of the atmosphere, making subsequently possible to visualize these
differences. It will also compute the spectrally integrated difference (radiative forcing) at these 3 levels. A result folder
must exist within "results", containing all output files for two test cases, and then radiative forcings will be computed
as the difference of radiative transfer results recorded in these two folders.
The main source file "results/flux_differences.for" must be edited and the names of the two folders must be provided.
Then the program should be recompiled using the "f0" script command.
+ added: scripts "plot_dFground", "plot_dFprobe" and "plot_dFtoa" in the "results" folder; these scripts
make possible to visualize the difference in radiative net fluxes at ground level, at probe level and at the top of the atmosphere.
+ added: scripts in the "results" folder:
- results/copy_all_results.bash: use it in order to copy ALL result files in the folder of your choice. This script needs
one argument: the name of the folder you want to save output files into. If the folder already exists, the script will
prompt you but will do nothing (in order to prevent the overwriting of existing result files). You must therefore provide
the name of a non-existing folder.
-results/copy_light_results.bash: use this script in order to copy only the smallest result files in the folder of your
choice. The files that are copied have been chosen so that the total size of the folder does not exceed 300 Mb.
In particular, opacity files are not copied. If you want to keep a trace of all results and opacities, you should use the
"copy_all_results.bash" script, possibly in combination with the "compression.bash" script.
As for the previous script, this one takes only one argument: the name of the folder you want to copy files to. Use
a non-existing folder, as the script will create it for you.
-results/compression.bash: this script will compress (in .gz format) result files that exceed a few hundred megabytes.
This script will be transfered in the destination directory by both "copy_all_results.bash" and "copy_light_results.bash",
so you do not have to do it manually.
- results/plot_all: this script will produce all plots from available result files; it only produces the .pdf files, and
does not perform the visualization. This script is transfered in the destination directory by the "copy_all_results.bash" script. 
- results/plot_all_light: this script will produce a subset of plots from available result files; it only produces the .pdf files, and
does not perform the visualization. This script is transfered in the destination directory by the "copy_light_results.bash" script.
- results/plot_all_differences: this script produces the plots relative to radiative transfer results differences. It is transfered
into the destination directory by both "copy_all_results.bash" and "copy_light_results.bash" scripts.

=======================
Version 2.1

2017/02
+ corrected: a bug in the computation of the emission altitude; it was not devided by a factor 2.
+ modified: now the spectrally integrated NER matrix is plotted using a logarithmic range. The range is automatically set over 3
orders of magnitude (toping at the maximum possible NER). The problem is for negative values, that are not
supported by a logarithmic scale (!). So we need to use some tricks to overcome the possibilities of gnuplot. The new
"results/plot_si_NER" script will still produce the "si_NER.pdf" file, but in order to do so one new requirement
is that "imagemagick" and "python" are installed. You may have to edit the "results/crop_image" python script
in order to specify your version of python (in the first line).
