[mesa-users] other_wind

[Bill Paxton]

Hi Chris,

Good question -- let's dig into the code to find the answer.

grep other_wind in star/private.
look in winds.f90 where it is called:
                      call s% other_wind(s% id, L1, M1, R1, T1, wind, ierr)
find how M1 is set
M1 = M_phot
find declaration of M_phot as arg to set_mdot
         real(dp), intent(in) :: L_phot, M_phot, T_phot ! photosphere values (cgs)
to double check, find where set_mdot is called to see what it is given as arg M_phot.
grep lead to the call in evolve.f90
         call set_mdot(s, s% L_phot*Lsun, s% mstar, s% Teff, ierr)

So it turns out that M1 is actually s% mstar.  
That's the mass coordinate at the outer edge of the outer cell.
But in spite of the name M_phot in winds.f90, it is not the photosphere mass unless the photosphere is at the surface of the model.

Cheers,
Bill




On Jun 15, 2015, at 11:59 AM, Chris Peacock wrote:

> HI all, 
> 
> I am using a custom subroutine of the form " *_other_wind", and I noticed that the default variables included in the subroutine are Msurf, Lsurf, Tsurf, Rsurf,  w and ierr.
>  Before I proceed any further with my model, I would like to double check that the variable 'Msurf' refers to the outermost mass shell (in my case, the outer atmosphere) of the star.  
> 
> Thanks
> 
> Chris
> ------------------------------------------------------------------------------
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users