PyREC scaffolding (from YREC development, 7be4133)/ Suspected bugs / custom-mixGlossary

custom-mix

A made-up case for checking suspected bugs by running YREC: testsuite/Test_solarAAG21_dif 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:

custom-mix: parmin.f bug 2.
From YREC 7be4133, testsuite/Test_solarAAG21_dif: a solar model with the AAG21 mixture and diffusion.
Changed (each with its source):
ISETMIX = 2: the case's own "Custom mixture controls" block, and docs/source/namelist_control.md lines 173-185
ISETISO = 2: as ISETMIX
the case's commented "Custom mixture controls" block, uncommented with its own values
and the output names. Nothing else.

custom-mix.nml1 (16 setting lines changed from testsuite/Test_solarAAG21_dif.nml1; download the whole file)

@@ YREC's line 8, this file's line 15 @@
  ! These options adjust their relative abundances. It is recommended to tie these to preset mixtures used in opacity tables.
  ! This should be done at the start of a model run, and only work if LFIRST=True.
  ! Option to adjust the CNO mixture to be consistent with the mixture used in the opacity tables
- ISETMIX=1 ! 0 = do not adjust 1 = use specified mixture (AMIX below) 2 = custom mixture (variables at the end of the namelist)
+ ISETMIX = 2 ! bug check: the case's own "Custom mixture controls" block, and docs/source/namelist_control.md lines 173-185
  ! AMIX = 'AAG21' ! current implemented solar mixtures are GS98,AAG21,M22P and M22M
  ! Option to adjust the D/He3/Li6/Li7/Be9 abundances and the CNO isotope ratios to be consistent with a specified solar reference
- ISETISO = 1! 0 = do not adjust 1 = use specified mixture (AISO below) 2 = custom mixture (variables at the end of the namelist)
+ ISETISO = 2 ! bug check: as ISETMIX
  ! AISO = 'L21' ! Lodders 2021 is currently the only supported mixture, for a different one use the custom option
 
  ! LCALS = .FALSE. ! Toggle to auto-calibrate solar Y and alpha to match the solar R and L at the age of the Sun.
@@ YREC's line 112, this file's line 119 @@
 ! Custom mixture controls - CNO enabled if ISETMIX = 2, Light elements if ISETISO=2
 
 ! Total CNO (isotope ratios controlled with ISETISO)
-! FRAC_C = 0.172148 ! Mass fraction of metals in the form of carbon; this is the GS98 value
-! FRAC_N = 0.050426 ! Mass fraction of metals in the form of nitrogen; this is the GS98 value
-! FRAC_O = 0.468195 ! Mass fraction of metals in the form of oxygen; this is the GS98 value
+ FRAC_C = 0.172148 ! Mass fraction of metals in the form of carbon; this is the GS98 value
+ FRAC_N = 0.050426 ! Mass fraction of metals in the form of nitrogen; this is the GS98 value
+ FRAC_O = 0.468195 ! Mass fraction of metals in the form of oxygen; this is the GS98 value
 
 ! Isotope ratios and light elements
-! R12_13 = 88.27 ! C12/C13 ratio, L21 mix (88.26 L25)
-! R14_15 = 411.9 ! N14/N15 ratio, L21 mix (412.3 L25) - placeholder,N15 not currently used
-! R16_17 = 471.4 ! O16/O17 ratio, L21 mix (2644 L25)- placeholder,O17 not currently used
-! R16_18 = 2693.0 ! O16/O18 ratio, L21 mix (468.2 L25)
+ R12_13 = 88.27 ! C12/C13 ratio, L21 mix (88.26 L25)
+ R14_15 = 411.9 ! N14/N15 ratio, L21 mix (412.3 L25) - placeholder,N15 not currently used
+ R16_17 = 471.4 ! O16/O17 ratio, L21 mix (2644 L25)- placeholder,O17 not currently used
+ R16_18 = 2693.0 ! O16/O18 ratio, L21 mix (468.2 L25)
 
 ! Light element abundances
-! XH2_INI = 2.781e-5 ! Initial D mass fraction, L21 mix (2.780e-5 L25)
-! XHE3_INI = 3.461e-5 ! Initial D mass fraction, L21 mix (3.443e-5 L25)
-! XLI6_INI = 7.187e-10 ! Initial D mass fraction, L21 mix (7.830e-10 L25)
-! XLI7_INI = 1.025e-8 ! Initial D mass fraction, L21 mix (1.112e-8 L25)
-! XBE9_INI = 1.595e-10 ! Initial D mass fraction, L21 mix (1.810e-10 L25)
-! XB10_INI = 1.002e-9 ! Initial D mass fraction, L21 mix (1.099e-9 L25) - placeholder, not currently used
-! XB11_INI = 4.405e-9 ! Initial D mass fraction, L21 mix (4.885e-9 L25) - placeholder, not currently used
+ XH2_INI = 2.781e-5 ! Initial D mass fraction, L21 mix (2.780e-5 L25)
+ XHE3_INI = 3.461e-5 ! Initial D mass fraction, L21 mix (3.443e-5 L25)
+ XLI6_INI = 7.187e-10 ! Initial D mass fraction, L21 mix (7.830e-10 L25)
+ XLI7_INI = 1.025e-8 ! Initial D mass fraction, L21 mix (1.112e-8 L25)
+ XBE9_INI = 1.595e-10 ! Initial D mass fraction, L21 mix (1.810e-10 L25)
+ XB10_INI = 1.002e-9 ! Initial D mass fraction, L21 mix (1.099e-9 L25) - placeholder, not currently used
+ XB11_INI = 4.405e-9 ! Initial D mass fraction, L21 mix (4.885e-9 L25) - placeholder, not currently used
 
  $END

custom-mix.nml2 (0 setting lines changed from testsuite/Test_solarAAG21_dif.nml2; download the whole file)

Unchanged.