Correct to eos=24 in docs and analysis_disc_stresses.f90

danieljprice · Jan 29, 2025 · c4b7f8e · c4b7f8e
1 parent 29e010f
commit c4b7f8e
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/docs/examples/selfgravity_gravitationalinstability.rst b/docs/examples/selfgravity_gravitationalinstability.rst
@@ -20,7 +20,7 @@ In the *setup* file there are the disc parameters, and in the *input* file it is
 
 To use the radiative cooling approximation of Young et al. (2024) use:
 
-- ieos = 23
+- ieos = 24
 - icooling = 9
 
-See :doc:`Radiation hydrodynamics in phantom </physics/radiation>`.
+See :doc:`Radiation hydrodynamics in phantom </physics/radiation>`.
diff --git a/src/utils/analysis_disc_stresses.f90 b/src/utils/analysis_disc_stresses.f90
@@ -74,7 +74,7 @@ subroutine do_analysis(dumpfile,numfile,xyzh,vxyzu,pmass,npart,time,iunit)
 
 ! Read analysis options
  call read_analysis_options
- if (ieos==23) call read_optab(eos_file,ierr)
+ if (ieos==24) call read_optab(eos_file,ierr)
 
  if (mhd) print*, 'This is an MHD dump: will calculate Maxwell Stress'
 
@@ -503,7 +503,7 @@ subroutine calc_stresses(npart,xyzh,vxyzu,pmass)
  call print_units
 
  sigma(:) = sigma(:)*umass/(udist*udist)
- if (ieos /= 23) then
+ if (ieos /= 24) then
     csbin(:) = csbin(:)*unit_velocity
  endif