danieljprice · alisonkyoung1 · Jan 18, 2023 · Jan 25, 2023 · Jan 27, 2023 · Jan 18, 2023
diff --git a/.gitignore b/.gitignore
@@ -36,4 +36,4 @@ build/test000
 *.cmod
 *.ilm
 *.stb
-
+eos_lom.dat
diff --git a/AUTHORS b/AUTHORS
@@ -16,6 +16,7 @@ Yann Bernard <[email protected]>
 Daniel Mentiplay <[email protected]>
 Megha Sharma <[email protected]>
 Arnaud Vericel <[email protected]>
+Alison Young <[email protected]>
 Mark Hutchison <[email protected]>
 Mats Esseldeurs <[email protected]>
 Rebecca Nealon <[email protected]>
@@ -34,7 +35,6 @@ Giovanni Dipierro <[email protected]>
 Roberto Iaconi <[email protected]>
 Amena Faruqi <[email protected]>
 Hauke Worpel <[email protected]>
-Alison Young <[email protected]>
 Stephen Nielson <[email protected]>
 Martina Toscani <[email protected]>
 Benedetta Veronesi <[email protected]>

diff --git a/build/Makefile b/build/Makefile
@@ -506,14 +506,15 @@ SRCCHEM= fs_data.f90 mol_data.f90 utils_spline.f90 \
          cooling_koyamainutsuka.f90 \
          cooling_ism.f90 \
          cooling_molecular.f90 \
+	 cooling_radapprox.f90\
          cooling_functions.f90 \
          cooling_solver.f90 \
          h2chem.f90 cooling.f90
 #
 # equations of state
 #
 SRCMESA= eos_mesa_microphysics.f90 eos_mesa.f90
-SRCEOS = eos_barotropic.f90 eos_stratified.f90 eos_piecewise.f90 ${SRCMESA} eos_shen.f90 eos_helmholtz.f90 eos_idealplusrad.f90 ionization.F90 eos_gasradrec.f90 eos_HIIR.f90 eos.f90
+SRCEOS = eos_barotropic.f90 eos_stratified.f90 eos_piecewise.f90 ${SRCMESA} eos_shen.f90 eos_helmholtz.f90 eos_idealplusrad.f90 ionization.F90 eos_gasradrec.f90 eos_HIIR.f90 eos_stamatellos.f90 eos.f90
 
 ifeq ($(HDF5), yes)
     SRCREADWRITE_DUMPS= utils_hdf5.f90 utils_dumpfiles_hdf5.f90 readwrite_dumps_common.f90 readwrite_dumps_fortran.F90 readwrite_dumps_hdf5.F90 readwrite_dumps.F90

diff --git a/data/eos/lombardi/README b/data/eos/lombardi/README
@@ -0,0 +1,10 @@
+The data tables for the equation of state and opacities for the modified Lombardi radiative cooling approximation are too large to be stored in the Phantom git repository.
+They will be downloaded automatically when you run the code
+
+or can be retrieved manually using wget from Zenodo, e.g.:
+
+wget https://zenodo.org/records/13842491/files/eos_lom.dat
+
+The files are:
+
+eos_lom.dat
diff --git a/docs/examples/selfgravity_gravitationalinstability.rst b/docs/examples/selfgravity_gravitationalinstability.rst
@@ -17,3 +17,10 @@ In the *setup* file there are the disc parameters, and in the *input* file it is
 - ieos = 2 , to choose an adiabatic equation of state
 - icooling = 3, to choose a β cooling prescription with constant β (icooling = 7 prescribes a varying β cooling with the radius)
 - beta_cool = #, to choose the value of β cooling
+
+To use the radiative cooling approximation of Young et al. (2024) use:
+
+- ieos = 23
+- icooling = 9
+
+See :doc:`Radiation hydrodynamics in phantom </physics/radiation>`.
diff --git a/docs/physics/eos-list.rst b/docs/physics/eos-list.rst
@@ -124,3 +124,9 @@
    |           | from HORMONE, Hirai+2020, as used in Lau+2022b                                   |
    |           |                                                                                  |
    +-----------+----------------------------------------------------------------------------------+
+   | 23        | **Tabulated EoS of Stamatellos et al. 2007 (includes opacities)**                |
+   |           |                                                                                  |
+   |           | Tabulated equation of state with opacities from Lombardi et al. 2015. For use    | 
+   |           |  with icooling = 9, the radiative cooling approximation (Young et al. 2024).     |													  |
+   |		   |																				  |
+   +-----------+----------------------------------------------------------------------------------+
diff --git a/docs/physics/radiation.rst b/docs/physics/radiation.rst
@@ -103,6 +103,28 @@ the star will just steadily cool. For red giants with a sink particle core a sim
 is to supply a constant luminosity input from the core (:doc:`sink heating <sinks>`). This is experimental.
 Another option would be to include a nuclear burning network (please somebody contribute this).
 
+
+Polytropic radiative cooling approximation
+-------------------------------------------
+This is an alternative to computationally expensive radiative transfer in regimes where radiative cooling is 
+important. This method estimates the optical depth for each particle and its equilibrium temperature. From 
+these the new temperature and internal energy is updated at each timestep. The method implemented here is 
+the "modified Lombardi" method of `Young et al. (2024) <https://ui.adsabs.harvard.edu/abs/2024MNRAS.531.1746Y>`__,
+ which was based on Stamatellos et al. (2007) and Lombardi et al. (2015). This method is designed for 
+ self-gravitating discs around a central star. Stellar heating is included from the most massive sink particle.
+
+Use icooling = 9 and ieos = 23 to use the tabulated equation of state which has the opacity tables required
+for the cooling calculation. The additional parameters are::
+
+	EOS_file =  myeos.dat   ! File containing tabulated EOS values
+	Lstar =       0.440    ! Luminosity of host star for calculating Tmin (Lsun)
+	Tfloor =       5.000    ! temperature floor (K); on if > 0
+
+
+N.B. This version does not currently include radiation from more than one sink particle. This version does not 
+couple with flux-limited diffusion at the moment.
+
+
 Irradiation from stars with phantom + MCFOST
 ---------------------------------------------
 In regimes where the radiation diffusion time is relatively short, anything not inside stars

diff --git a/src/main/H2regions.f90 b/src/main/H2regions.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/apr.f90 b/src/main/apr.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/apr_region.f90 b/src/main/apr_region.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/bondiexact.f90 b/src/main/bondiexact.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/bondiexact_gr.f90 b/src/main/bondiexact_gr.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/boundary.f90 b/src/main/boundary.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/boundary_dynamic.f90 b/src/main/boundary_dynamic.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/centreofmass.f90 b/src/main/centreofmass.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/checkconserved.f90 b/src/main/checkconserved.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/checkoptions.F90 b/src/main/checkoptions.F90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/checksetup.f90 b/src/main/checksetup.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/config.F90 b/src/main/config.F90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/cons2prim.f90 b/src/main/cons2prim.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!

diff --git a/src/main/cons2primsolver.f90 b/src/main/cons2primsolver.f90
@@ -1,6 +1,6 @@
 !--------------------------------------------------------------------------!
 ! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-! Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
+! Copyright (c) 2007-2025 The Authors (see AUTHORS)                        !
 ! See LICENCE file for usage and distribution conditions                   !
 ! http://phantomsph.github.io/                                             !
 !--------------------------------------------------------------------------!
-Original file line number
+Diff line change
@@ Expand Up / @@ -36,4 +36,4 @@ build/test000 @@
     *.cmod
     *.ilm
     *.stb
+    eos_lom.dat