[mesa-users] New wind scheme, Grafener et al.

[niloofar afsari]

Hi,
I have written a subroutine that implements a mass loss recipe from "Grafener, G. & Hamann, W.-R. 2008, A&A 482, 945" .I'm copying it here, in case anyone needs this mass loss prescription. 

 
         subroutine eval_grafener_wind(w)            real(dp), intent(inout) :: w            real(dp) :: w1, logMdot, gamma_edd, xsurf, beta, gammazero      xsurf = surface_h1     gamma_edd= exp10_cr(-4.813d0)*(1+xsurf)*(L1/Lsun)*(Msun/M1)            beta= 1.727d0 + 0.250d0*log10_cr(Z/Zsolar)            gammazero= 0.326d0 - 0.301d0*log10_cr(Z/Zsolar)- 0.045d0*log10_cr(Z/Zsolar)**2             logMdot = &           + 10.046d0 &           + beta*log10_cr(gamma_edd-gammazero) &           - 3.5d0*log10_cr(T1) &                          + 0.42d0*log10_cr(L1/Lsun) &                          - 0.45d0*xsurf   w = exp10_cr(logMdot)                      if (dbg) write(*,*) 'grafener wind', w                     end subroutine eval_grafener_wind
Nilou
 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20141201/e589c514/attachment.html>