castelli-kttau5

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

castelli-kttau5: setups.f bug 1, kcsurfp.f bug 1, putstore.f bug 2.
From YREC 7be4133, testsuite/Test_solarGS98_kurucz: a solar model with Kurucz atmospheres.
Changed (each with its source):
KTTAU = 5: input/atmos/CastelliKurucz/README, line 15: "You NEED TO CHANGE KTTAU to 5" for these, "THE PREFERRED MODELS"
FATM = "{YREC_INPUT}/atmos/CastelliKurucz/atmk2004p00.tab": the solar-metallicity table of that folder
and the output names. Nothing else.

castelli-kttau5.nml1 (1 setting line changed from testsuite/Test_solarGS98_kurucz.nml1; download the whole file)

@@ YREC's line 84, this file's line 90 @@
 
 ! Equation of state tables
  FOPALE06 = "{YREC_INPUT}/eos/opal2006/EOSOPAL06Z0.016492" ! OPAL 2006 EoS table path.
- FATM = "{YREC_INPUT}/atmos/kurucz/atmk1990p00.tab" ! Model atmosphere look-up table for Kurucz and Kurucz/Castelli.
+ FATM = "{YREC_INPUT}/atmos/CastelliKurucz/atmk2004p00.tab" ! bug check: the solar-metallicity table of that folder
  FALLARD = "{YREC_INPUT}/atmos/allard/Nextgen2.all" ! Allard model atmosphere look-up table (solar Z only).
  FSCVH = "{YREC_INPUT}/eos/scv/h_tab_i.dat" ! SCZ EoS table for X
  FSCVHE = "{YREC_INPUT}/eos/scv/he_tab_i.dat" ! SCZ EoS table for Y

castelli-kttau5.nml2 (1 setting line changed from testsuite/Test_solarGS98_kurucz.nml2; download the whole file)

@@ YREC's line 104, this file's line 110 @@
  XSBE93 = 21000.0! Cross-section for p(9Be,alpha) (KeV barns)
 
 ! Atmospheric surface boundary condition options
- KTTAU = 3 ! Toggle switch. 0 = Eddington gray 1 = Krishnaswamy 2=legacy (don't use) 3=Kurucz 4=Allard. Note - for 3 and 4, switches to gray if outside table domain.
+ KTTAU = 5 ! bug check: input/atmos/CastelliKurucz/README, line 15: "You NEED TO CHANGE KTTAU to 5" for these, "THE PREFERRED MODELS"
 ! Atmosphere and envelope integrator options. Note - if a tabulated atmosphere is used, only an envelope integration is performed to the outer fitting point.
  TRIDT = .005 ! The code solves for a SBC between 3 models. This is the spacing between 2 of them in Log Teff. The third is a different L at the average Teff.
  TRIDL = .040 ! Spacing between the lower pair and upper point in Log L.