[mesa-users] Using other_wind.f

[Jieun Choi]

Hi everyone,

I'm experimenting with a different wind scheme which ties the regular
Blöcker wind to an exponential factor involving the envelope mass. However,
it seems like MESA only allows you to access Lsurf, Msurf, Rsurf, and Tsurf
instead of the full star_info structure (from which I'd get the envelope
mass). Is there some way to get around this?

other_wind is called in winds.f using

else if (scheme == 'other') then
>     if (dbg) write(*,*) 'call other_wind'
>     call s% other_wind(s% id, L1, M1, R1, T1, wind, ierr)
>     if (ierr /= 0) return


and I copied the following from other_wind.f to use as a template in
run_star_extras.f

      subroutine null_other_wind(id, Lsurf, Msurf, Rsurf, Tsurf, w, ierr)
>          use star_def
>          integer, intent(in) :: id
>          real(dp), intent(in) :: Lsurf, Msurf, Rsurf, Tsurf ! surface
> values (cgs)
>          ! NOTE: surface is outermost cell. not necessarily at photosphere.
>          ! NOTE: don't assume that vars are set at this point.
>          ! so if you want values other than those given as args,
>          ! you should use values from s% xh(:,:) and s% xa(:,:) only.
>          ! rather than things like s% Teff or s% lnT(:) which have not
> been set yet.
>          real(dp), intent(out) :: w ! wind in units of Msun/year (value is
> >= 0)
>          integer, intent(out) :: ierr
>          w = 0
>          ierr = 0
>       end subroutine null_other_wind


Thanks,
Jieun
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20130105/7c0130aa/attachment.html>