giant3.f

Go to the documentation of this file.

      SUBROUTINE giant3(I,BODYI,W,Q,XN,QL)
*
*
*       Structure constants of giant star.
*       ----------------------------------
*
*       Theory of Rosemary Mardling, Ap. J. XX, YYY, 1995.
*       @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
*
      include 'common6.h'
      REAL*8  w(2),q(2),sw(2),mc
*     DATA  WW  /2.119,3.113,8.175,3.742,4.953,9.413/
*     DATA  QQ  /0.4909,0.4219,0.2372,0.4677,0.3560,0.1519/
      DATA  a0,a1,a2,a3  /0.944525,-0.392030,6.01655d-02,-3.34790d-03/
      DATA  b0,b1,b2  /-0.3789,1.481,-0.1018/
      DATA  c0,c1,c2,c3  /1.452104,3.923872,-11.88722,13.46106/
      DATA  e0,e1,e2,e3  /1.934977,2.214222,-4.855796,4.025394/
*
*
*       Set typical core mass of 0.3/0.5 Msun.
      mc = 0.3 + 0.1*float(kstar(i) - 3)
      IF(kstar(i).EQ.9) mc = 0.5*bodyi*zmbar
*
*       Form ratio of core mass and mass.
      sig = mc/(bodyi*zmbar)
*
*       Include safety check on mass ratio just in case.
      IF (sig.GT.0.9) sig = 0.9
*
*       Define mass, radius, envelope mass and luminosity in S.U.
      zm = bodyi*zmbar
      rsi = radius(i)*su
      zme = zm - mc
*       Obtain L^{1/3} from giant relation L = 1.98D+05*M_c^6.
      zl3 = 58.3*mc**2
*       Evaluate damping constant from Zahn's theory (R.M. 13/5/97).
*       FAC = (GM)^{1/2}*M_{env}^{1/3)/(R^{5/6}*L^{1/3}) = 8.48D+03 for S.U.
      ql = 8.48d+03*sqrt(zm)*(zme/zl3)**0.33/rsi**0.833
*
*       Set effective frequencies, overlap integrals and structure constants.
      DO 10 k = 1,2
*         K1 = 3*K - 2
      IF (k.EQ.1) THEN
              sw(k) = ((c3*sig + c2)*sig + c1)*sig + c0
          ELSE
              sw(k) = ((e3*sig + e2)*sig + e1)*sig + e0
      END IF
      w(k) = sw(k)**2
          q(k) = ((a3*sw(k) + a2)*sw(k) + a1)*sw(k) + a0
   10 CONTINUE
*
*       Evaluate new polytropic index.
      xn = (b2*sw(1) + b1)*sw(1) + b0
*     WRITE (24,20)  IC, IPAIR, KSTAR(J), MC/(BODYI*SMU), RSI, XN, QL
*  20 FORMAT (' GIANT:    IC KS K* MC/M R* n Q ',3I4,3F6.2,F7.1)
*     CALL FLUSH(24)
*
*       Include warning if n > 5.
*     IF (XN.GE.5.0) THEN
*         WRITE (6,30)  IC, IPAIR, KSTAR(J), CM(L,IC)/BODY(I),
*    &                  RADIUS(I)/RIN, XN, QL
*  30     FORMAT (' GIANT:    WARNING!    IC KS K* MC/M R/R0 n QL ',
*    &                                    3I4,F6.2,F6.1,F6.2,F7.1)
*     END IF

      RETURN
*
      END