Merge branch 'master' of github.com:danieljprice/phantom into star

danieljprice · Jan 30, 2024 · 41c8c06 · 41c8c06
2 parents 2042002 + 3063500
commit 41c8c06
Show file tree

Hide file tree

Showing 454 changed files with 7,340 additions and 1,937 deletions.
diff --git a/.github/workflows/mpi.yml b/.github/workflows/mpi.yml
@@ -19,7 +19,7 @@ jobs:
           - yes
         ntasks:
           - 1
-          - 4
+          - 2
         input: # [SETUP, phantom_tests]
           - ['test', '']
           - ['testkd', '']
@@ -43,7 +43,7 @@ jobs:
 
     env:
       OMP_STACKSIZE: 512M
-      OMP_NUM_THREADS: 4
+      OMP_NUM_THREADS: 2
 
     steps:
 

diff --git a/.mailmap b/.mailmap
@@ -21,12 +21,13 @@ Rebecca Nealon <[email protected]> Nealon <[email protected]
 Alex Pettitt <[email protected]> Alex <[email protected]>
 Alex Pettitt <[email protected]> <[email protected]>
 Alex Pettitt <[email protected]> <[email protected]>
-Terrence Tricco <[email protected]>
-<[email protected]> <[email protected]>
-<[email protected]> <tricco@solo.(none)>
-<[email protected]> <tricco@rin.(none)>
-<[email protected]> <[email protected]>
-<[email protected]> <[email protected]>
+<[email protected]> <[email protected]>
+<[email protected]> <[email protected]>
+<[email protected]> <tricco@solo.(none)>
+<[email protected]> <tricco@rin.(none)>
+<[email protected]> <[email protected]>
+<[email protected]> <[email protected]>
+Terrence Tricco <[email protected]> <[email protected]> 
 James Wurster <[email protected]> James Wurster <[email protected]>
 James Wurster <[email protected]> jameswurster <[email protected]>
 James Wurster <[email protected]> jameswurster <[email protected]>
@@ -46,12 +47,12 @@ Stéven Toupin <[email protected]> stoupin <devnull@localhost>
 Guillaume Laibe <[email protected]> glaibe <[email protected]>
 Guillaume Laibe <[email protected]> glaibe <[email protected]>
 Alice Cerioli <[email protected]> ALICE CERIOLI <[email protected]>
+Alice Cerioli <[email protected]> <[email protected]>
 Thomas Reichardt <[email protected]>
 Thomas Reichardt <[email protected]> Thomas  Reichardt <[email protected]>
 Thomas Reichardt <[email protected]> Mr Thomas Reichardt <[email protected]>
 Roberto Iaconi <[email protected]> Roberto Iaconi <rxi552@raijin6.(none)>
 Roberto Iaconi <[email protected]> Roberto Iaconi <[email protected]>
-Alice Cerioli <[email protected]> <[email protected]>
 Daniel Mentiplay <[email protected]> <[email protected]>
 Daniel Mentiplay <[email protected]> <[email protected]>
 Daniel Mentiplay <[email protected]> <[email protected]>
@@ -72,10 +73,15 @@ Enrico Ragusa <[email protected]> <[email protected]>
 Enrico Ragusa <[email protected]> <[email protected]>
 Kieran Hirsh <[email protected]> <[email protected]>
 Giulia Ballabio <[email protected]> Giulia Ballabio <[email protected]>
+Mats Esseldeurs <[email protected]>
+Mats Esseldeurs <[email protected]> <[email protected]>
+Mats Esseldeurs <[email protected]> <[email protected]>
 Lionel Siess <[email protected]>
 Lionel Siess <[email protected]> <[email protected]>
 Lionel Siess <[email protected]> <[email protected]>
 Lionel Siess <[email protected]> <[email protected]>
+Mats Esseldeurs <[email protected]>
+Mats Esseldeurs <[email protected]> <[email protected]>
 David Liptai <[email protected]> <[email protected]>
 David Liptai <[email protected]> <[email protected]>
 David Liptai <[email protected]> <[email protected]>
@@ -85,8 +91,13 @@ Fangyi (Fitz) Hu <[email protected]> Fitz-Hu <[email protected]
 Fangyi (Fitz) Hu <[email protected]> root <root@Heartie的小电脑.localdomain>
 Fangyi (Fitz) Hu <[email protected]> root <root@Heartie的小电脑.localdomain>
 Fangyi (Fitz) Hu <[email protected]> fitzHu <[email protected]>
+Fangyi (Fitz) Hu <[email protected]> Fitz Hu <[email protected]>
+Fangyi (Fitz) Hu <[email protected]> fhu <[email protected]>
 Megha Sharma <[email protected]> Megha Sharma <[email protected]>
 Megha Sharma <[email protected]> megha sharma <[email protected]>
+Megha Sharma <[email protected]> Megha Sharma <[email protected]>
+Megha Sharma <[email protected]> Megha Sharma <[email protected]>
+Megha Sharma <[email protected]> Megha Sharma <[email protected]> 
 Mike Lau <[email protected]> Mike Lau <[email protected]>
 Elisabeth Borchert <[email protected]> emborchert <[email protected]>
 Ward Homan <[email protected]> ward <[email protected]>
@@ -103,3 +114,9 @@ Sahl Rowther <[email protected]> <[email protected]>
 Sahl Rowther <[email protected]> sahl95 <[email protected]>
 Caitlyn Hardiman <[email protected]> caitlynhardiman <[email protected]>
 Amena Faruqi <[email protected]> <[email protected]>
+Amena Faruqi <[email protected]> Amena Faruqi <[email protected]>
+Alison Young <[email protected]> Alison Young <[email protected]>
+Simone Ceppi <[email protected]> Simone Ceppi <[email protected]>
+Mats Esseldeurs <[email protected]> mats esseldeurs <[email protected]>
+Mats Esseldeurs <[email protected]> MatsEsseldeurs <[email protected]>
+Nicolás Cuello <[email protected]> Nicolas Cuello <[email protected]>
diff --git a/AUTHORS b/AUTHORS
@@ -11,33 +11,30 @@ Conrad Chan <[email protected]>
 James Wurster <[email protected]>
 David Liptai <[email protected]>
 Lionel Siess <[email protected]>
+Fangyi (Fitz) Hu <[email protected]>
 Daniel Mentiplay <[email protected]>
+Megha Sharma <[email protected]>
 Arnaud Vericel <[email protected]>
 Mark Hutchison <[email protected]>
-Fitz Hu <[email protected]>
-Megha Sharma <[email protected]>
 Rebecca Nealon <[email protected]>
 Elisabeth Borchert <[email protected]>
 Ward Homan <[email protected]>
 Christophe Pinte <[email protected]>
-Fangyi (Fitz) Hu <[email protected]>
-Megha Sharma <[email protected]>
-Terrence Tricco <[email protected]>
+Terrence Tricco <[email protected]>
+Simone Ceppi <[email protected]>
+Mats Esseldeurs <[email protected]>
 Mats Esseldeurs <[email protected]>
 Stephane Michoulier <[email protected]>
-Simone Ceppi <[email protected]>
-MatsEsseldeurs <[email protected]>
+Spencer Magnall <[email protected]>
 Caitlyn Hardiman <[email protected]>
 Enrico Ragusa <[email protected]>
-fhu <[email protected]>
 Sergei Biriukov <[email protected]>
 Cristiano Longarini <[email protected]>
 Giovanni Dipierro <[email protected]>
 Roberto Iaconi <[email protected]>
 Hauke Worpel <[email protected]>
+Amena Faruqi <[email protected]>
 Alison Young <[email protected]>
-Simone Ceppi <[email protected]>
-Amena Faruqi <[email protected]>
 Stephen Neilson <[email protected]>
 Martina Toscani <[email protected]>
 Benedetta Veronesi <[email protected]>
@@ -46,31 +43,26 @@ Simon Glover <[email protected]>
 Thomas Reichardt <[email protected]>
 Jean-François Gonzalez <[email protected]>
 Christopher Russell <[email protected]>
+Phantom benchmark bot <[email protected]>
 Alessia Franchini <[email protected]>
-Alex Pettitt <[email protected]>
 Jolien Malfait <[email protected]>
-Phantom benchmark bot <[email protected]>
-Kieran Hirsh <[email protected]>
+Alex Pettitt <[email protected]>
 Nicole Rodrigues <[email protected]>
-Amena Faruqi <[email protected]>
+Kieran Hirsh <[email protected]>
 David Trevascus <[email protected]>
+Nicolás Cuello <[email protected]>
+Farzana Meru <[email protected]>
 Chris Nixon <[email protected]>
-Megha Sharma <[email protected]>
-Nicolas Cuello <[email protected]>
+Miguel Gonzalez-Bolivar <[email protected]>
+Mike Lau <[email protected]>
 Benoit Commercon <[email protected]>
-Farzana Meru <[email protected]>
+Orsola De Marco <[email protected]>
 Giulia Ballabio <[email protected]>
-Joe Fisher <[email protected]>
 Maxime Lombart <[email protected]>
-Megha Sharma <[email protected]>
-Orsola De Marco <[email protected]>
-Terrence Tricco <[email protected]>
+Joe Fisher <[email protected]>
 Zachary Pellow <[email protected]>
 s-neilson <[email protected]>
-Alison Young <[email protected]>
 Cox, Samuel <[email protected]>
-Jorge Cuadra <[email protected]>
-Nicolás Cuello <[email protected]>
 Steven Rieder <[email protected]>
 Stéven Toupin <[email protected]>
-mats esseldeurs <[email protected]>
+Jorge Cuadra <[email protected]>
diff --git a/Makefile b/Makefile
@@ -4,9 +4,8 @@
 #
 #   See build/Makefile for the main Makefile
 #
-#   (c) 2007-2023 The Authors (see AUTHORS)
+#   (c) 2007-2024 The Authors (see AUTHORS)
 #
-#  $Id: Makefile,v 98b9fad01f38 2013/03/25 23:02:49 daniel $
 #----------------------------------------------------------------
 
 .PHONY: phantom

diff --git a/build/Makefile b/build/Makefile
@@ -1,6 +1,6 @@
 #--------------------------------------------------------------------------!
 # The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
-# Copyright (c) 2007-2023 The Authors (see AUTHORS)                        !
+# Copyright (c) 2007-2024 The Authors (see AUTHORS)                        !
 # See LICENCE file for usage and distribution conditions                   !
 # http://users.monash.edu.au/~dprice/phantom                               !
 #--------------------------------------------------------------------------!
@@ -14,14 +14,12 @@
 #   the SETUP variable
 #
 #  OWNER: Daniel Price
-#
-#  $Id: 2788b71b1c08e560e77dce9849c5cb24a668f4b9 $
 #+
 #--------------------------------------------------------------------------
 
 .KEEP_STATE:
 
-PHANTOM_VERSION_MAJOR=2023
+PHANTOM_VERSION_MAJOR=2024
 PHANTOM_VERSION_MINOR=0
 PHANTOM_VERSION_MICRO=0
 VERSION=$(PHANTOM_VERSION_MAJOR).$(PHANTOM_VERSION_MINOR).$(PHANTOM_VERSION_MICRO)
@@ -244,10 +242,6 @@ ifeq ($(NONIDEALMHD), yes)
     FPPFLAGS += -DNONIDEALMHD
 endif
 
-ifeq ($(H2CHEM), yes)
-    FPPFLAGS += -DH2CHEM
-endif
-
 ifeq ($(DISC_VISCOSITY), yes)
     FPPFLAGS += -DDISC_VISCOSITY
 endif
@@ -468,7 +462,8 @@ OBJDIR=obj
 # external forces
 #
 ifeq (X$(SRCPOTS), X)
-SRCPOTS= extern_corotate.f90 \
+SRCPOTS= extern_gr.F90 \
+         extern_corotate.f90 \
          extern_binary.f90 \
          extern_spiral.f90 \
          extern_lensethirring.f90 \
@@ -478,6 +473,7 @@ SRCPOTS= extern_corotate.f90 \
          extern_densprofile.f90 \
          extern_staticsine.f90 \
          extern_gwinspiral.f90 \
+         extern_geopot.f90 \
          externalforces.f90
 endif
 ifeq (X$(SRCPOT), X)
@@ -487,14 +483,15 @@ endif
 # metrics for GR
 #
 ifeq ($(GR),yes)
-    SRCPOT=extern_gr.F90 $(SRCPOTS:externalforces.f90=externalforces_gr.f90)
+    SRCPOT=$(SRCPOTS:externalforces.f90=externalforces_gr.f90)
 endif
+
 ifdef METRIC
     SRCMETRIC= metric_${METRIC}.f90
 else
     SRCMETRIC= metric_minkowski.f90
 endif
-SRCGR=inverse4x4.f90 $(SRCMETRIC) metric_tools.f90 utils_gr.f90
+SRCGR=inverse4x4.f90 einsteintk_utils.f90 $(SRCMETRIC) metric_tools.f90 utils_gr.f90 interpolate3D.f90 tmunu2grid.f90 
 #
 # chemistry and cooling
 #
@@ -547,6 +544,10 @@ SOURCES= physcon.f90 ${CONFIG} ${SRCKERNEL} io.F90 units.f90 \
          mf_write.f90 evolve.F90 utils_orbits.f90 utils_linalg.f90 \
          checksetup.F90 initial.F90
 
+# Needed as einsteintk_wrapper depends on initial
+ifeq ($(GR),yes)
+	SOURCES+=einsteintk_wrapper.f90
+endif
 OBJECTS1 = $(SOURCES:.f90=.o)
 OBJECTS = $(OBJECTS1:.F90=.o)
 
@@ -566,16 +567,14 @@ ifeq ($(UNAME), Darwin)
 endif
 
 
+ifeq ($(KROME), krome)
 	@echo ""
 	@echo "=============== CHEMISTRY ==============="
 	@echo ""
-ifeq ($(KROME), krome)
 	@echo "krome coupling status = enabled"
-else
-	@echo "krome coupling status = disabled"
-endif
 	@echo ""
 	@echo "========================================="
+endif
 	@sh ../scripts/phantom_version_gen.sh "$(IDFLAGS)"
 	@echo ""
 	@echo "The Phantom is here (in $(BINDIR)/phantom)"
@@ -618,14 +617,14 @@ SRCDUMP= physcon.f90 ${CONFIG} ${SRCKERNEL} utils_omp.F90 io.F90 units.f90 \
          utils_dumpfiles.f90 utils_vectors.f90 utils_mathfunc.f90 \
          utils_datafiles.f90 utils_filenames.f90 utils_system.f90 utils_tables.f90 datafiles.f90 gitinfo.f90 \
          centreofmass.f90 \
-         timestep.f90 ${SRCEOS} cullendehnen.f90 dust_formation.F90 \
+         timestep.f90 ${SRCEOS} cullendehnen.f90 dust_formation.f90 \
          ${SRCGR} ${SRCPOT} \
          memory.F90 \
          utils_sphNG.f90 \
          setup_params.f90 ${SRCFASTMATH} checkoptions.F90 \
          viscosity.f90 damping.f90 options.f90 checkconserved.f90 prompting.f90 ${SRCDUST} \
          ${SRCREADWRITE_DUMPS} \
-         utils_sort.f90 sort_particles.F90
+         utils_sort.f90 sort_particles.f90
 OBJDUMP1= $(SRCDUMP:.f90=.o)
 OBJDUMP= $(OBJDUMP1:.F90=.o)
 
@@ -1323,7 +1322,7 @@ getdims:
 	@echo $(MAXP)
 
 get_setup_opts:
-	@echo "${GR:yes=GR} ${METRIC} ${MHD:yes=MHD} ${NONIDEALMHD:yes=non-ideal} ${DUST:yes=dust} ${GRAVITY:yes=self-gravity} ${RADIATION:yes=radiation}  ${H2CHEM:yes=H2_Chemistry} ${DISC_VISCOSITY:yes=disc_viscosity} ${ISOTHERMAL:yes=isothermal} ${PERIODIC:yes=periodic}" | xargs | sed -e 's/ /, /g' -e 's/_/ /g'
+	@echo "${GR:yes=GR} ${METRIC} ${MHD:yes=MHD} ${NONIDEALMHD:yes=non-ideal} ${DUST:yes=dust} ${GRAVITY:yes=self-gravity} ${RADIATION:yes=radiation} ${DISC_VISCOSITY:yes=disc_viscosity} ${ISOTHERMAL:yes=isothermal} ${PERIODIC:yes=periodic}" | xargs | sed -e 's/ /, /g' -e 's/_/ /g'
 
 get_setup_file:
 	@echo "$(SETUPFILE)"

diff --git a/build/Makefile_defaults_gfortran b/build/Makefile_defaults_gfortran
@@ -12,7 +12,7 @@
 #  endif
 #
 FC= gfortran
-FFLAGS+= -O3 -Wall -Wno-unused-dummy-argument -frecord-marker=4 -gdwarf-2 \
+FFLAGS+= -O3 -Wall -Wno-unused-dummy-argument -frecord-marker=4 -g \
         -finline-functions-called-once -finline-limit=1500 -funroll-loops -ftree-vectorize \
         -std=f2008 -fall-intrinsics
 DBLFLAG= -fdefault-real-8 -fdefault-double-8

diff --git a/build/Makefile_defaults_ifort b/build/Makefile_defaults_ifort
@@ -15,8 +15,8 @@ KNOWN_SYSTEM=yes
 
 # for ifort version 18+ -openmp flag is obsolete
 IFORT_VERSION_MAJOR=${shell ifort -v 2>&1 | head -1 | cut -d' ' -f 3 | cut -d'.' -f 1}
-ifeq ($(shell [ $(IFORT_VERSION_MAJOR) -gt 17 ] && echo true),true)
-    OMPFLAGS= -qopenmp
+ifeq ($(shell [ $(IFORT_VERSION_MAJOR) -lt 17 ] && echo true),true)
+    OMPFLAGS= -openmp
 else
-    OMPFLAGS = -openmp
+    OMPFLAGS = -qopenmp
 endif
diff --git a/build/Makefile_setups b/build/Makefile_setups
@@ -775,7 +775,7 @@ ifeq ($(SETUP), dustystar)
     FPPFLAGS= -DDUST_NUCLEATION -DSTAR
     SETUPFILE= setup_star.f90
     MODFILE= utils_binary.f90 set_binary.f90 moddump_binary.f90
-    ANALYSIS= ${SRCNIMHD} utils_summary.o utils_omp.o ptmass.o energies.o analysis_common_envelope.f90 dust_formation.F90
+    ANALYSIS= ${SRCNIMHD} utils_summary.o utils_omp.o ptmass.o energies.o analysis_common_envelope.f90 dust_formation.f90
     KNOWN_SETUP=yes
     MAXP=10000000
     GRAVITY=yes
@@ -1036,6 +1036,25 @@ ifeq ($(SETUP), testgr)
     SETUPFILE= setup_grdisc.f90
 endif
 
+ifeq ($(SETUP), flrw)
+#   constant density FLRW cosmology with perturbations  
+ GR=yes
+ KNOWN_SETUP=yes
+ IND_TIMESTEPS=no
+ METRIC=et
+ SETUPFILE= setup_flrw.f90
+ PERIODIC=yes
+endif
+
+ifeq ($(SETUP), flrwpspec)
+#  FLRW universe using a CMB powerspectrum and the Zeldovich approximation
+ GR=yes
+ KNOWN_SETUP=yes
+ IND_TIMESTEPS=no
+ METRIC=et
+ SETUPFILE= setup_flrwpspec.f90
+ PERIODIC=yes
+endif
 ifeq ($(SETUP), default)
 #   default setup, uniform box
     KNOWN_SETUP=yes

diff --git a/docs/conf.py b/docs/conf.py
@@ -20,13 +20,13 @@
 # -- Project information -----------------------------------------------------
 
 project = 'Phantom'
-copyright = '2023 The Authors'
+copyright = '2024 The Authors'
 author = 'Daniel Price'
 
 # The short X.Y version
-version = '2023.0'
+version = '2024.0'
 # The full version, including alpha/beta/rc tags
-release = '2023.0.0'
+release = '2024.0.0'
 
 
 # -- General configuration ---------------------------------------------------

diff --git a/docs/eos-list.rst b/docs/eos-list.rst
@@ -7,8 +7,8 @@
 |           |                                                                                  |
 |           | where :math:`c_s^2 \equiv K` is a constant stored in the dump file header        |
 |           |                                                                                  |
-+-----------+----------------------------------------------------------------------------------+
-| 2         | **Adiabatic equation of state (code default)**                                   |
+|           |                                                                                  |
+|           | Adiabatic equation of state (code default)                                       |
 |           |                                                                                  |
 |           | :math:`P = (\gamma - 1) \rho u`                                                  |
 |           |                                                                                  |