The purpose of PPart_LW is to perform radiative transfer computations using k-distribution data sets (that can be produced by Kdistribution, or any other code). This code is based on analytic radiative transfer solutions that are valid only in the case of a emitting and absorbing, but non-scattering, atmosphere. This code should therefore be used, in the case of the terrestrial atmosphere, for clear-sky LW computations only.
Download and extract the archive (
tar -zxvf ppart_lw.tgz); then
move into the
PPart_LW directory, and use the
all command to compile. You should use the
command in the case you want to recompile from scratch only (when include files
are modified). If only source files are modified, the
command will recompile modified files.
PPart_LW has been parallelized, which means you will have to launch it using
mpirun -np <NPROCS> ./PPart.exe
<NPROCS> the number of processes requested for the
computation; MacOS users might encounter a error where the OS considers the
user requested a number of processes higher than recommended. You can bypass it
mpirun -oversubscribe -np <NPROCS> ./PPart.exe
Prerequistes: only a fortran compiler is needed (no external libraries). the gfortran compiler has been used for development, but other common compilers should work too (ifort, pgfortran, etc.)
The code can run either using a pre-computed (and provided) k-distribution
data set among 5 standard atmospheric profiles, or from a user-defined spectral
data file. This user-defined spectral data file should be name
ecrad_opt_prop.txt, and located into the
directory. The format of this file can be found in the
/Doc/gas_opt_prop.pdf documentation file, and a example
/Doc/input.for source file provides a fortran subroutine for
reading this file. A example
ecrad_opt_prop.txt file is provided
/data directory (for the Mid-Latitude Summer profile).
Simulation input is provided through the
options.in files (located in the main
options.infile provides the possibility to use either a (provided) standard atmospheric data set or a user-defined spectral data set (in this case, the
/data/ecrad_opt_prop.txtfile has to be found). It provides the possibility to use a specular or diffuse reflective ground, and the possibility to perform the spectral integration over a limited number of spectral intervals. Then the user should specify whether a angular integration should be performed or a single-direction radiative transfer computation is required. Finally, the user can disable the computation of a NER matrix for every spectral interval.
data.infile lets the user specify the spectral integration domain (if a limited spectral domain should be used), and the direction to use for a single-direction radiative transfer computation. Finally, the user should specify the index of the standard atmospheric profile to use, when the corresponding option has been selected.
Results are located into the
/results directory; several
gnuplot scripts are provided in order to visualise various results. When a NER
matrix has been computed for every spectral interval, a script named
anim.bash can be found in the
directory; this script will produce a animated gif of these NER matrices (the
spectral interval evolving with time). This script requires the GraphicsMagick
PPart_LW is free software released under the GPLv2+ license: GNU GPL version 2 or later. You can freely study, modify or extend it. You are also welcome to redistribute it under certain conditions; refer to the license for details.