![[ICO]](/icons/blank.gif){kind=icon}
![[PARENTDIR]](/icons/back.gif){kind=icon}
![[ ]](/icons/compressed.gif){kind=icon}
![[TXT]](/icons/text.gif){kind=icon}
Routine for calculating the mass-loss rate of M-type AGB stars at solar metallicityDescription:Here we provide a simple routine in FORTRAN for calculating the mass-loss rate as a function of stellar parameters (mass, luminosity, and effective temperature), based on data from the extensive grid of DARWIN models for M-type AGB stars (Bladh et al. 2019, A&A). The resulting mass-loss rates are calculated by interpolation between the current stellar mass (M), the effective temperature (T), and the stellar luminosity (log(L)). The remaining parameters are selected from the best fit to observed wind properties of nearby AGB stars (log(nd/nH)=-14.5 and Δup = 2.0). This mass-loss routine can be modified to include other mass-loss formulas or a lower limit for the mass-loss rate.Validity:This subroutine is valid for stellar masses up to ~3.25 solar masses and effective temperatures above 2200K. The upper limit of luminosity depends on the current steller mass and varies between log(L/L☉)=4.15 and log(L/L☉)=4.60.M < 0.75M☉ is set to M = 0.75M☉ M > 3.00M☉ is set to M = 3.00M☉ T < 2200K is set to T = 2400K In this mass-loss routine the mass-loss rate outside the luminosity range of the grid is assumed to scale linearly with luminosity. This may not be appropriate for luminosities that are much higher than the upper luminosity of the model grid. The mass-loss routine therefore includes two flags (MF and LF) that signal if the current mass or luminosity is outside the range of the model grid. If MF/LF=0 then the mass and luminosity are within the grid range. If MF/LF=-1 then the mass and luminosity values are outside the range of the grid and the resulting mass-loss rate might not be valid. Mass-loss routine command lines:The model data is read in with the line: CALL GRIDDATA("mgrid_solar.dat").The mass-loss rate for a specific combination of stellar parameters (M, L, T) is given by the line: MLR = MASSLOSS(MI, M, log(L), T, MF, LF), where MF and LF are warning flags. MI, M, log(L), and T are real*4 variables and MF and LF are integer variables. The initial mass MI is not used in the current interpolation but is included for implementating alternative mass-loss descriptions. File Summary: mlr_mstar_solar.zip---------------------------------------------------------------------------------------------------- FileName Explanations ---------------------------------------------------------------------------------------------------- mlr_mstar_solar.f Mass-loss routine in FORTRAN77 mlrmod_mstar.inc Mass-loss routine module mgrid_solar.dat The input data, based on the extensive model grid in Bladh et al. 2019 ---------------------------------------------------------------------------------------------------- Author's address: sara.bladh@physics.uu.se ==================================================================================================== (End) Sara Bladh [Uppsala university, Sweden] 08-May-2019 |
| Name | Last modified | Size | Description |
|---|---|---|---|---|
| Parent Directory | - | ||
| mlr_mstar_solar.zip | 2019-05-23 17:38 | 171K | |
| HEADER.html | 2019-05-27 13:54 | 3.5K | |