! getopal95-2: getopal95.f bug 2 (empty OP-table rows above log T = 8.0).
! From YREC 7be4133, examples/run_from_zahb_to_tahb/ZAHB_m1p0feh+0p0_TAHB,
! a 1 Msun core helium-burning model from the ZAHB to central Y = 0.0001.
! Changed: FLIV95, from legacy/opal95/GS98.OP17 to op05/AAG2021zx0.0187.tab,
! the table input/opacity/atomic/README recommends for the AAG21 mixture
! (and says the hydrogen-burning tables are also required for helium
! burning); the output names. Nothing else: the composition stays the case's.
 $CONTROL  

 DESCRIP(1) = "model5.1 example. Core helium burning, no diffusion, no rotation"
 DESCRIP(2) = "Gray Atm OPAL06+SCV EOS GS98" ! Second line label for output (optional)
 
 ISETMIX=0 ! 0 = do not adjust 1 = use specified mixture (AMIX below) 2 = custom mixture (see parmin.f for variables)
 AMIX = 'GS98'
 ISETISO = 0! 0 = do not adjust 1 = use specified mixture (AISO below) 2 = custom mixture (see parmin.f for variables)

 LCALS = .FALSE. ! Toggle to auto-calibrate solar Y and alpha to match the solar R and L at the age of the Sun. Does not enforce correct Z/X for mixture.
 LCALSOLZX = .TRUE.
 CALSOLZX = 0.02292d0 
 TOLL = 3.0D-6 ! Tolerance target for L in the calibration routines.
 TOLR = 3.0D-6 ! Tolerance target for R in the calibration routines.
 CALSOLAGE = 4.568D9 ! Tolerance solar age in the calibration routines.
  
 NUMRUN = 3

 KINDRN(1) = 2 ! 1 = evolve 2= rescale (zero timestep,MS or CHeB); 3 = rescale and evolve (pre-MS, no zero timestep)
 LFIRST(1) = .TRUE. ! T = use stored starting model F = use result of prior run
 CMIXLA(1) = 1.70767823 ! 10/28/26: solar calibration (base); was 1.707873
 RSCLX(1) = 0.71961003 ! 10/28/26: solar calibration (base); was 0.7196200
 RSCLZ(1) = 0.01649314 ! 10/28/26: solar calibration (base); was 0.016492
 NMODLS(1) = 10

 KINDRN(2) = 2 ! rescale
 LFIRST(2) = .FALSE.
 CMIXLA(2) = 1.70767823 ! 10/28/26: solar calibration (base); was 1.707873
 RSCLX(2) = 0.71961003 ! 10/28/26: solar calibration (base); was 0.7196200
 RSCLZ(2) = 0.01649314 ! 10/28/26: solar calibration (base); was 0.016492
 RSCLM(2) = -1 ! If positive, rescale to this mass (Msun)
 RSCLCM(2) = 0.46157 ! If positive, rescale to this core mass (CHeB only)
 NMODLS(2) = 10

 KINDRN(3) = 1
 LFIRST(3) = .FALSE.
 NMODLS(3) = 1000000
 XENV0A(3) = 0.71961003  ! 10/28/26: solar calibration (base); was 0.7196200 -- Envelope abundance label
 ZENV0A(3) = 0.01649314 ! 10/28/26: solar calibration (base); was 0.016492
 CMIXLA(3) = 1.70767823 ! 10/28/26: solar calibration (base); was 1.707873
 END_YCEN(3) = 0.0001 ! Run stops if NMODLS(I) is reached or this central helium is reached (TAHB)

 LSTORE = .TRUE. ! If T, store model structures at specified points to the .store file. 
 LSTATM = .FALSE. ! If T, ppend an atmosphere to the stored model file.
 LSTENV = .FALSE. ! If T, append an envelope to the stored model file. If LSTENV and LSTATM = F, no information is output beyond the fitting point.
 LSTMOD = .TRUE. ! Set the same as LSTORE. May be redundant.
 LSTPHYS = .TRUE. ! If T, output microphysics (e.g. epsilon, kappa) to store file.
 LSTROT = .FALSE. ! If T, append rotation data to store file.
 NPRTMOD = 1000 ! Output detailed model structures every NPRTMOD models.
 NPRTPT = 1 ! Output details for every NPRTPT points (if 1000 points and NPRTPT=10, 100 lines are output).
 LSTPCH = .TRUE. ! If T, output details of the final model to the .store file.

 LTRACK = .TRUE. ! If T, global properties are always output to the .track file. This should almost always be on.
 ITRVER = 0 ! Toggle for different track file formats. 0 (one line per model) strongly recommended.
 LRWSH = .FALSE. ! If T, the short file retains only the details from the last model. This minimizes output file size.
 LCORR = .TRUE. ! If true, output verbose details on model convergence to .short.  
 LPULSE = .FALSE. ! If T, generate pulsation output. LSTATM AND LSTENV are needed to get data for those regions (highly recommended if LPULSE is T)  
  
 LALEX06 = .TRUE. ! If T, use Ferguson+2005 molecular opacities (recommended)
 FALEX06 = "{YREC_INPUT}/opacity/molecular/alex06/alexmol06gs98.tab" ! Molecular opacity table path  
 LOPAL95 = .TRUE. ! If T, use atomic opacities in the OP/OPAL format (recommended).
 ZOPAL951 = 0.01649314  ! 10/28/26: solar calibration (base); was 0.016492 -- Generate a surface mixture table with this Z. Speeds up models without	 diffusion.
 FLIV95 = "{YREC_INPUT}/opacity/atomic/op05/AAG2021zx0.0187.tab" ! Atomic opacity table path
 LPUREZ = .TRUE. 
 FPUREZ = "{YREC_INPUT}/opacity/atomic/lanl/PURECO.DBGLAOL" ! Pure Y/Z  opacity table path

 FFIRST = "{YREC_START}/zahb/100mSOLpreZAHB.last"
 FLAST = "{YREC_OUTPUT}/getopal95-2.first" ! The last converged model is stored here.
 FMODPT = "{YREC_OUTPUT}/getopal95-2.full" ! No longer used.
 FSTOR = "{YREC_OUTPUT}/getopal95-2.store" ! Model snapshots during the run are stored here.
 FTRACK = "{YREC_OUTPUT}/getopal95-2.track" ! Global variables as a function of time are stored here.
 FSHORT = "{YREC_OUTPUT}/getopal95-2.short" ! Details about the numerics of the model run are stored here.
 FPMOD = "{YREC_OUTPUT}/getopal95-2.pmod" ! Pulsation output for the interior model.
 FPENV = "{YREC_OUTPUT}/getopal95-2.penv"! Pulsation output for the envelope model.
 FPATM = "{YREC_OUTPUT}/getopal95-2.atm" ! Pulsation output for the atmosphere model.
 FSNU = "{YREC_OUTPUT}/getopal95-2.snu" ! No longer used (formally, output if LSNU = T).
 FSCOMP = "{YREC_OUTPUT}/getopal95-2.excomp" ! No longer used (light elements vs. time, now included in track file)

 FOPALE06 = "{YREC_INPUT}/eos/opal2006/EOSOPAL06Z0.016492" ! OPAL 2006 EoS table path.
 FcondOpacP = "{YREC_INPUT}/opacity/conductive/potekhin/condall06.d" ! Conductive opacity table path.
 FATM = "{YREC_INPUT}/atmos/kurucz/atmk1990p00.tab" ! Model atmosphere look-up table for Kurucz and Kurucz/Castelli.
 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
 FSCVZ = "{YREC_INPUT}/eos/scv/z_tab_i.dat" ! SCZ EoS table for Z (generated from YREC EoS, needed for consistency)
 FFERMI = "{YREC_INPUT}/eos/yale/FERMI.TAB" ! Fermi integral tables for partial degeracy
   
 $END