Merge pull request #438 from danieljprice/misc-radiation

Few and minor cleanup to modules related to radiation
danieljprice · Jun 29, 2023 · b020e27 · b020e27
2 parents d3fbe44 + cb02c04
commit b020e27
Show file tree

Hide file tree

Showing 4 changed files with 3 additions and 16 deletions.
diff --git a/src/main/ptmass_heating.f90 b/src/main/ptmass_heating.f90
@@ -48,15 +48,4 @@ subroutine energ_sinkheat(nptmass,xyzmh_ptmass,xi,yi,zi,dudtheati)
 end subroutine energ_sinkheat
 
 
-!-----------------------------------------------------------------------
-!+
-!  Heating from nuclear burning
-!+
-!-----------------------------------------------------------------------
-! subroutine nuc_heating()
-
-
-! end subroutine nuc_heating
-
-
 end module ptmass_heating
diff --git a/src/main/radiation_implicit.f90 b/src/main/radiation_implicit.f90
@@ -461,8 +461,6 @@ subroutine fill_arrays(ncompact,ncompactlocal,npart,icompactmax,dt,xyzh,vxyzu,iv
        dvy = vxyzu(2,i) - vxyzu(2,j)
        dvz = vxyzu(3,i) - vxyzu(3,j)
 
-       rhomean = 0.5*(rhoi+rhoj)
-
        dvdotdr = dvx*dx + dvy*dy + dvz*dz
        dv = dvdotdr/dr
 
@@ -1005,7 +1003,7 @@ subroutine update_gas_radiation_energy(ivar,ijvar,vari,ncompact,npart,ncompactlo
     !
     EU0(1,i) = E1i
     EU0(2,i) = U1i
-    radprop(icv,i) = get_cv(rhoi,EU0(2,i),cv_type) ! should this be EU0(2,i)??
+    radprop(icv,i) = get_cv(rhoi,EU0(2,i),cv_type)
     radprop(ikappa,i) = get_kappa(iopacity_type,EU0(2,i),radprop(icv,i),rhoi)
 
     if (implicit_radiation_store_drad) then  ! use this for testing

diff --git a/src/setup/set_star.f90 b/src/setup/set_star.f90
@@ -697,7 +697,7 @@ subroutine read_options_star(star,need_iso,ieos,polyk,db,nerr,label)
     else
        star%isinkcore = .true.
        call read_inopt(star%input_profile,'input_profile'//trim(c),db,errcount=nerr)
-       call read_inopt(star%outputfilename,'outputfilename//trim(c)',db,errcount=nerr)
+       call read_inopt(star%outputfilename,'outputfilename'//trim(c),db,errcount=nerr)
        if (star%isoftcore==1) call read_inopt(star%isofteningopt,'isofteningopt'//trim(c),&
                                               db,errcount=nerr,min=0)
        if ((star%isofteningopt==1) .or. (star%isofteningopt==3)) then

diff --git a/src/utils/moddump_rad_to_LTE.f90 b/src/utils/moddump_rad_to_LTE.f90
@@ -6,7 +6,7 @@
 !--------------------------------------------------------------------------!
 module moddump
 !
-! moddump
+! Convert radiation dump to LTE dump (ieos=12)
 !
 ! :References: None
 !