reimers-3msun

A made-up case for checking suspected bugs by running YREC: testsuite/Test_m3p0feh+0p0GS98_std25_TAMS at 7be4133, with the changes below as diffs against that case (removed lines red, added green, the words that differ shaded darker); besides these, only the head and the output names differ. Its head says what it is:

reimers-3msun: massloss.f bug 1, mdot.f bug 2, dburn.f bug 1.
From YREC 7be4133, testsuite/Test_m3p0feh+0p0GS98_std25_TAMS: a 3 Msun main-sequence star.
Changed (each with its source):
LMDOT = .TRUE.: as in reimers-rgb: a 3 Msun track meant to reach the giant branch with Reimers mass loss starts on the main sequence
LREIMER = .TRUE.: the Reimers switch (nmlpar_mod.f), default coefficient CREIM
and the output names. Nothing else.

reimers-3msun.nml1 (0 setting lines changed from testsuite/Test_m3p0feh+0p0GS98_std25_TAMS.nml1; download the whole file)

Unchanged.

reimers-3msun.nml2 (2 setting lines changed from testsuite/Test_m3p0feh+0p0GS98_std25_TAMS.nml2; download the whole file)

@@ YREC's line 246, this file's line 252 @@
  LADOV = .FALSE. ! Experimental option to enforce adiabatic overshoot.
  LVFC = .FALSE. ! Implements FC as a variable, using the Zahn 1991 approach. Not recommended (disables FC setting.)
  LDIFAD = .FALSE. ! Treat meridional circulation in a diffusion+ advection approach (Zahn) vs purely diffusive (earlier)
- LMDOT = .FALSE. ! Allow mass loss (experimental code). Possible options are a constant mass loss rate (or gain); a scaled solar wind; or a Reimers' parameterization for mass loss in evolved giants. Constant Mdot (negative) and Reimers are currently enabled.
+ LMDOT = .TRUE. ! bug check: as in reimers-rgb: a 3 Msun track meant to reach the giant branch with Reimers mass loss starts on the main sequence
  LSOLWIND = .FALSE. ! Remove mass loss scaled to the solar wind. Not currently implemented.
 ! LTSLIR = .FALSE.
 ! TSLITOL = 0.01
@@ YREC's line 257, this file's line 263 @@
  MCORE = 40 ! Points to assign to new domain; not used now.
  ALFA = 1.5 ! Legacy exponent for angular momentum loss in Kawaler (1988) wind model.
 
+ LREIMER = .TRUE. ! bug check: the Reimers switch (nmlpar_mod.f), default coefficient CREIM
  $END