[mesa-users] Fwd: Implementing new mass loss receipt

[Bill Paxton]

Hi,

Can someone else help Ehsan with this -- I'm buried at the moment and can't get to it.

Thanks,
Bill




Begin forwarded message:

> From: "Ehsan Moravveji" <moravveji at iasbs.ac.ir>
> Date: September 17, 2012 9:54:02 AM PDT
> To: "Bill Paxton" <paxton at kitp.ucsb.edu>
> Subject: Implementing new mass loss receipt
> 
> Hello Bill,
> There is a very recent mass loss receipt for Wolf-Rayet stars of WC-type
> by Sander et al. (2012).
> I liked to implement it from my run_star_extras.
> To this end, I wrote the following subroutine in the my personalized
> user_wolf_rayet.f file:
> 
>  ! Semi-empirical mass loss rate "only" for WC stars based on
>  ! Eq. 6 in Sander A., Hamman, W.R. and Todt H. 2012, A&A, 540, 144
>  ! Note to set optical depth tau = 20. They give:
>  ! log(Mdot [Msun/yr]) = (0.83 +/- 0.11) * log(L/Lsun) - 9.10 +/- 0.57
>  subroutine eval_sander_wind(s,w)
>  type (star_info), pointer :: s
>  real(dp), intent(out) :: w
>  real(dp) :: log10w
>  integer :: i20
> 
> !   include 'formats'
>  i20 = minloc(abs(s% tau(:) - 20d0), dim=1)
> 
>  log10w = 0.83d0 * dlog10(s% L(i20) / lsol) - 9.10d0
>  w = 10d0**log10w
> 
>  write(*,*) ' eval_sander_wind: log_Mdot = ', log10w
> 
>  end subroutine eval_sander_wind
> 
> 
> and then in my extras_finish_step set:
>       s% other_wind => eval_sander_wind
>       write(*,*), 'WC Mdot rate = ', dlog10(s% star_mdot)
> which gives me NaN ...
> 
> To solve it "crudely", I instead do the following:
>       call eval_sander_wind(s, mdot)
>       s% star_mdot = mdot
> which works fine.
> 
> Do you admit the latter method?
> 
> Merci :-D
> Ehsan.
> 
> 
> 
> -- 
> This message has been scanned for viruses and
> dangerous content by MailScanner, and is
> believed to be clean.
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20120917/117cdf7c/attachment.html>