Burn_nseAtDens.F90

Go to the documentation of this file.

!! NOTICE
!!  Copyright 2022 UChicago Argonne, LLC and contributors
!!
!!  Licensed under the Apache License, Version 2.0 (the "License");
!!  you may not use this file except in compliance with the License.
!!
!!  Unless required by applicable law or agreed to in writing, software
!!  distributed under the License is distributed on an "AS IS" BASIS,
!!  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
!!  See the License for the specific language governing permissions and
!!  limitations under the License.
!!
 
 
subroutine Burn_nseAtDens(qbar_nse,sumyi_nse,approxtemp,edot,Yedot, Ye, dens, emq)
  implicit none
  
  real, intent(IN) :: Ye, dens, emq
  real, intent(OUT) :: qbar_nse,sumyi_nse,approxtemp,edot,Yedot
  
! NOTE Ivo suggested values that are closer to sensible than just plain zero
!  This is representative of 100% Ni56 at a temperature of 4x10^9K
  qbar_nse = 8.64 ! MeV/nucleon
  sumyi_nse = 1./56.
  approxtemp = 4.0E9 !Kelvin
  edot = 0.0
  Yedot = 0.0
 
  return
end subroutine Burn_nseAtDens