[mesa-users] Debugging WD cooling in MESA

[Héctor MR <hec...@pi...>]

  Dear all,

  I have been following Bill Wolf's excellent guide for debugging in MESA (
http://wmwolf.github.io/projects/mesa_debugging/) in what concerns to
cooling WDs because I am encountering convergence problems. When I try to
cool the star, in this case for a maximum age of 1 Gyr, my computations
stop around ~0.327 Gyr due to some problems in the envelope. The luminosity
suddenly wiggles and goes down, whilst gradT_sub_grada and mixing_type
lurch as well. I suspect it is related to

          okay_to_reduce_gradT_excess = .true.

          ! these are for calculation of efficiency boosted gradT
          gradT_excess_f1 = 1d-4
          gradT_excess_f2 = 1d-2
          gradT_excess_age_fraction = 0.9d0

          ! turn on full all the time
          gradT_excess_lambda1 = 1.0
          gradT_excess_beta1 = 0.5
          gradT_excess_lambda2 = 1.0
          gradT_excess_beta2 = 0.5
          gradT_excess_dlambda = 1
          gradT_excess_dbeta = 1


  But I have been unable to solve it (turning Ledoux off seems not to help
either). Does anyone know what is going on? I just do not want to feel
compelled to cut the models via set_initial_tau_factor. I am attaching the
inlist, the nuclear net, the initial model and some plots.

  Thanks in advance!


  --
  Héctor