Merge pull request #422 from danieljprice/moddump_radtoLTE

src/utils/moddump_rad_to_LTE.f90

@@ -0,0 +1,57 @@
+!--------------------------------------------------------------------------!
+! The Phantom Smoothed Particle Hydrodynamics code, by Daniel Price et al. !
+! Copyright (c) 2007-2023 The Authors (see AUTHORS)                        !
+! See LICENCE file for usage and distribution conditions                   !
+! http://phantomsph.bitbucket.io/                                          !
+!--------------------------------------------------------------------------!
+module moddump
+ !
+ ! Convert a radiative dump into a non-radiative dump, assuming LTE (i.e., ieos=12)
+ ! INSTRUCTIONS: 1) Set gmw (mean molecular weight) in the code below to the
+ !                  same value assumed in the radiative dump
+ !               1) Make this moddump with a radiative setup (e.g. SETUP=radstar),
+ !                  otherwise the rad array is not read.
+ !               2) Run this moddump
+ !               3) Write/modify the infile to remove radiation-related options,
+ !                  make sure ieos=12 and mu is correct
+ !               4) Recompile Phantom with the desired non-radiative setup (e.g.
+ !                  SETUP=star).
+ !
+ ! :References: None
+ !
+ ! :Owner: Mike Lau
+ !
+ ! :Runtime parameters: None
+ !
+ ! :Dependencies: 
+ !
+ implicit none
+
+ contains
+
+subroutine modify_dump(npart,npartoftype,massoftype,xyzh,vxyzu)
+ use dim, only:do_radiation
+ use io,  only:fatal
+ use eos, only:ieos,gmw
+ use part,only:rad,iradxi
+ integer, intent(inout) :: npart
+ integer, intent(inout) :: npartoftype(:)
+ real,    intent(inout) :: massoftype(:)
+ real,    intent(inout) :: xyzh(:,:),vxyzu(:,:)
+ integer                :: i
+
+ if (.not. do_radiation) call fatal("moddump_rad_to_LTE","Not compiled with radiation")
+ do i=1,npart
+    if (isnan(rad(iradxi,i))) call fatal("moddump_rad_to_LTE","rad array contains NaNs")
+    vxyzu(4,i) = vxyzu(4,i) + rad(iradxi,i)
+ enddo
+
+ ieos = 12
+ gmw = 0.6   ! CHANGE MU HERE for writing into infile
+ print*,'mu has been changed to',gmw  ! mu should not change from what was assumed with radiation
+
+end subroutine modify_dump
+
+end module moddump
+
+