sunbather.solveT

Functions

calc_HCratio(radheat, radcool, expcool, ...)	Calculates the ratio of total heating to total cooling.
calc_advection(r, rho, v, Te, mu)	Calculates advection heating/cooling (Linssen et al. 2024 Eq.
calc_cloc(radheat, radcool, expcool, ...)	Checks if there is a point in the atmosphere where we can use the construction algorithm.
calc_expansion(r, rho, v, Te, mu)	Calculates expansion cooling (Linssen et al. 2024 Eq.
check_converged(fc, HCratio, newTe, prevTe)	Checks whether the temperature profile is converged.
clean_converged_folder(folder)	Deletes all files in a folder that are not called "converged*".
constructTstruc(grid, newTe_relax, cloc, v, ...)	Proposes a new temperature profile based on a 'construction' algorithm, starting at the cloc and at higher altitudes.
get_new_Tstruc(old_Te, HCratio, fac)	Returns a new temperature profile based on a previous non-converged temperature profile and the associated heating/cooling imbalance.
make_converged_plot(altgrid, altmax, path, ...)	Makes a plot of the converged temperature profile, as well as the different heating/cooling rates.
make_rates_plot(altgrid, Te, newTe_relax, ...)	Makes a plot of the previous and newly proposed temperature profiles, as well as the different heating/cooling rates and their ratio based on the previous temperature profile.
relaxTstruc(grid, path, itno, Te, HCratio)	Proposes a new temperature profile using a 'relaxation' algorithm.
run_loop(path, itno, fc[, save_sp, maxit])	Solves for the nonisothermal temperature profile of a Parker wind profile through an iterative convergence scheme including Cloudy.
simtogrid(sim, grid)	Extracts various needed quantities from a Cloudy simulation and interpolates them onto the provided radius grid.