Generated from /asc/asci2/site/flashcode/secure/release_4p6/source/physics/materialProperties/NSE/NSE_finalAtDens.F90 with ROBODoc v4.99.8 on Tue Oct 15 13:28:03 2019


[Functions] source/physics/materialProperties/NSE/NSE_finalAtDens





    call NSE_finalAtDens( real(OUT)   :: qbar_nse,
                          real(OUT)   :: sumyi_nse,
                          real(OUT)   :: approxtemp,
                          real(OUT)   :: edot,
                          real(OUT)   :: Yedot,
                          real(IN)    :: Ye,
                          real(IN)    :: dens,
                          real(IN)    :: emq)


 This function gives the NSE final state for a constant density burn
 via table lookup.  qbar is returned in MeV/nucleon.
 Only the compositional information (qbar, 1/Abar, edot, Yedot) should be
 used for hydrodynamics purposes.
 Accurate thermodynamic prperties of the NSE final state (rho, T) should be
 obtained by solving
        eint - qbar_nse = emq .
 This is important due to the limited accuracy of the table interpolation
 (i.e., the interpolated values cannot satisfy any constraints on their own.)


  outputs  (properties of the NSE state that we want)
    qbar_nse    --  Average binding energy per nucleon (averaged over nse)
    sumyi_nse   --  Sum over the abundances. (Y = X/A = X/(#of nucleons))
    approxtemp  --  approximate temperature of NSE final state (see note above)
    edot        --  Neutrino and anti-neutrino loss rates [ergs/g]
    Yedot       --  Change in electron fraction wrt time (neutronization rate)

  inputs  (these uniquely define the NSE final state)
    Ye          --  Electron fraction
    dens        --  Density in/out
    emq         --  "e-qbar" = Internal energy per gram minus nuclear binding energy
                    per nucleon.  The latter is converted to units of ergs/g by assuming
                    that each nucleon has mass 1 a.m.u
                    This is the total energy in initial and final state -- this quantity
                    is constant during an isochoric burn because there is no work done


   Takes place after burning is complete.  NSE expands after burning
   and the equilibrium state is a function of the density and the
   temperature.  Been burned, flame is done, but then the hot stuff
   expands.  This changes the composition and releases energy.


   NSE_interface : top-level interface file for description of subroutine function


   Dean Townsley, Alan Calder 2006-8
   original interpolation kernel by Alan Calder (2006)
   thanks to Flash Code Group (esp. Lynn Reid) for bugfixes

   Dean Townsley 2008  stub for when unit is not included