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

[Bill Paxton]

Hi,

Thanks!  I've added it to the code, so it will be in the next release.

Cheers,
Bill



On Dec 1, 2014, at 8:37 AM, niloofar afsari wrote:

> 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
> 
>  
> ------------------------------------------------------------------------------
> Download BIRT iHub F-Type - The Free Enterprise-Grade BIRT Server
> from Actuate! Instantly Supercharge Your Business Reports and Dashboards
> with Interactivity, Sharing, Native Excel Exports, App Integration & more
> Get technology previously reserved for billion-dollar corporations, FREE
> http://pubads.g.doubleclick.net/gampad/clk?id=157005751&iu=/4140/ostg.clktrk_______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20141201/9c54cc1a/attachment.html>