corerd.f

Go to the documentation of this file.

***
      REAL*8 FUNCTION corerd(KW,MC,M0,MFLASH)
*     
* A function to determine the radius of the core of a giant-like star.
* NOTE: this is out of date so rc should be obtained using HRDIAG!
* It is still OK to use but bear in mind that the core radius calculated
* for non-degenerate giant cores is only a rough estimate.
*
*     Author : C. A. Tout
*     Date :   26th February 1997
*     Updated 6/1/98 by J. Hurley
*
      INTEGER kw
      REAL*8 mc,mch,m0,mflash
      parameter(mch = 1.44d0)
*
* First do the black holes and neutron stars.
*
      IF(kw.EQ.14)THEN
         corerd = 4.24d-06*mc
      ELSEIF(kw.EQ.13)THEN
         corerd = 1.4d-05
*
* Main sequence stars.
*
      ELSEIF(kw.LE.1.OR.kw.EQ.7)THEN
         corerd = 0.d0
*
* Core-helium-burning stars, FAGB stars and non-degenerate giant cores.
*
      ELSEIF(kw.EQ.4.OR.kw.EQ.5.OR.(kw.LE.3.AND.m0.GT.mflash))THEN
         corerd = 0.2239d0*mc**0.62d0
*
* The degenerate giants and white dwarfs.
*
      ELSE
         corerd = 0.0115d0*sqrt(max(1.48204d-06,(mch/mc)**(2.d0/3.d0)
     &                                        - (mc/mch)**(2.d0/3.d0)))
* 
* Degenerate giants have hot subdwarf cores.
*
         IF(kw.LE.9) corerd = 5.d0*corerd
      ENDIF
*
      RETURN
      END
***