reimers-rgb

A made-up case for checking suspected bugs by running YREC: testsuite/Test_m1p0feh+0p0GS98_base_HeIgnite 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-rgb: massloss.f bug 1, mdot.f bug 2, htimer.f bug 1.
From YREC 7be4133, testsuite/Test_m1p0feh+0p0GS98_base_HeIgnite: a 1 Msun star from the ZAMS up the red giant branch to helium ignition.
Changed (each with its source):
LMDOT = .TRUE.: the case's own comment: "Reimers' parameterization for mass loss in evolved giants ... Reimers are currently enabled"
LREIMER = .TRUE.: the Reimers switch (nmlpar_mod.f), default coefficient CREIM
and the output names. Nothing else.

reimers-rgb.nml1 (0 setting lines changed from testsuite/Test_m1p0feh+0p0GS98_base_HeIgnite.nml1; download the whole file)

Unchanged.

reimers-rgb.nml2 (2 setting lines changed from testsuite/Test_m1p0feh+0p0GS98_base_HeIgnite.nml2; download the whole file)

@@ YREC's line 241, this file's line 247 @@
  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: the case's own comment: "Reimers' parameterization for mass loss in evolved giants ... Reimers are currently enabled"
  LSOLWIND = .FALSE. ! Remove mass loss scaled to the solar wind. Not currently implemented.
 ! LTSLIR = .FALSE.
 ! TSLITOL = 0.01
@@ YREC's line 252, this file's line 258 @@
  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