Getting k-distributions for radiative transfer in Venus Atmosphere in longwave region (10 - 6000 cm-1). After running the executables one can obtain the next k-term in a given spectral band.
Note
- obtaining k-terms is an iterative process and requires running this programs several times.
- resulting parameterizations are strongly affected by input atmospheric profile and chosen spectroscopy (look in PTTables repository).
- Ubuntu 20.04 (or other Linux distro)
- gfortran compiler
- matplotlib (for plots)
- precalculated spectral PT-tables located in
/srv/PT_TABLES/
directory. PT-tables should be available for CO2, SO2 and H2O atmospheric components.
To build, type: ./make.sh
in the terminal. Make sure you are in the root directory of the project.
-
Fill the input file "band_settings.txt". Example (for the first channel):
__1 M 10.0D0 200.0D0 H2 H2O_gas_profile.dat
-
Fill the file "channel_settings.txt". Set values of thresholds of absorption coeffcient for lower and upper atmospheric levels. Example:
0.0 1 1e15 50
-
Run the executables
-
Run Python scripts
plot_cooling.py
,plot_flux_down.py
, orplot_flux_up.py
to plot cooling rates and upward/downward fluxes, to compare the accuracy of the obtained k-term (current channel number inband_settings.txt
file). -
If needed, repeat steps 1-4 with updated input values in
channels_settings.txt
file to get more accurate k-term.