rfast_surface

Modified version of the open source radiative transfer code created by Tyler Robinson and Arnaud Salvador: Robinson & Salvador 2023

Code has been modified to run forward models and retrievals using a user-defined albedo model (defined in user_models). Currently set-up for a five parameter surface parametrisation, further parameters will require small modifications in rfast_genspec_alb.py, rfast_retrieve_pp.py and rfast_analyze_pp.py. Under the folder agnostic_surface there are specific versions of genspec, user_models, retrieve, analyze, and spectra_gen to use for the agnostic surface configuration of the model. It is currently set for a linear parameterisation.

Description of changeable files
rfast_inputs.scr is the inputs file with user-chosen parameters
rfast_user_model.py script with models for the forward model, surface albedo parametrisation can be modified through surfalb function
rfast_genspec_alb.py is the forward modeling script
rfast_noise.py generates noisy spectra with user specified noise model
rfast_retrieve_pp.py main script for MCMC analysis - outputs a .h5 file with the chain data
rfast_analyze_pp.py first analysis script, generates corner plots, individual posteriors for land fractions, and fitted spectra
spectra_gen_pp.py second analysis script for plotting confidence intervals of fits

How to run

1. rfast_genspec_alb.py rfast_inputs.scr
2. rfast_noise.py rfast_inputs.scr
3. rfast_retrieve_pp.py rfast_inputs.scr
4. rfast_analyze_pp.py rfast_inputs.scr
5. spectra_gen_pp.py rfast_inputs.scr
   

This code relies on a co-located opacities file called hires_opacities. This and installation instructions can be found in the main rfast github.

Publications that use this code
https://doi.org/10.3847/1538-4357/adec69