u_wd Function

public pure function u_wd(rho)

Degenerate specific internal energy as a function of density

See Benz W., Bowers R.L., Cameron A.G.W., Press W.H., 1990, APJ, 348, 647., eqs.(2.5)

FT 19.12.2022

Arguments

Type IntentOptional Attributes Name
double precision, intent(in) :: rho

Return Value doubleprecision


Calls

proc~~u_wd~~CallsGraph proc~u_wd u_wd proc~g_wd g_wd proc~u_wd->proc~g_wd proc~f_wd f_wd proc~g_wd->proc~f_wd

Called by

proc~~u_wd~~CalledByGraph proc~u_wd u_wd proc~test_wd_eos_cgs test_wd_eos_cgs proc~test_wd_eos_cgs->proc~u_wd

Contents

Variables

Source Code

u_wd

Variables

Type Visibility Attributes Name Initial
double precision, private, parameter :: rho_min = 1.D-35
double precision, private :: x

Source Code

  DOUBLE PRECISION PURE FUNCTION u_wd(rho)

    !********************************************
    !
    !# Degenerate specific internal energy as a
    !  function of density
    !
    !  See [Benz W., Bowers R.L., Cameron A.G.W., Press W.H., 1990, APJ, 348, 647.](https://ui.adsabs.harvard.edu/abs/1990ApJ...348..647B/abstract){:target="_blank"}, eqs.(2.5)
    !
    !  FT 19.12.2022
    !
    !********************************************

    IMPLICIT NONE

    DOUBLE PRECISION, INTENT(IN):: rho

    DOUBLE PRECISION, PARAMETER:: rho_min= 1.D-35

    DOUBLE PRECISION:: x

    IF(rho == 0)THEN
      u_wd= 0.D0
      RETURN
    ENDIF

    x= (rho/b_wd)**third

    u_wd= a_wd*g_wd(x)/rho

  END FUNCTION u_wd