[mesa-users] New public release 9793
ebauer at physics.ucsb.edu
Wed May 31 14:51:27 EDT 2017
It looks like my mail client may have done something weird with my attempt to attach “changes.txt”, so I’m just copying all that inline here to make sure everyone can see it.
- use_se_output: .false.
- set_Z_all_HELM: .false.
+ relax_initial_angular_momentum: .false.
+ max_steps_to_relax_angular_momentum: 1000
+ timescale_for_relax_angular_momentum: 1d-10
+ max_dt_for_relax_angular_momentum: 1d-9
+ num_timescales_for_relax_angular_momentum: 1000
+ relax_angular_momentum_filename: ''
+ remove_initial_center_by_entropy: 0
+ remove_center_by_entropy: 0
+ fallback_check_total_energy: .true.
+ remove_initial_surface_by_v_surf_km_s: 0
+ remove_surface_by_v_surf_km_s: 0
+ set_HELM_OPAL_Zs: .false.
+ Z_all_OPAL: 0.04d0
+ set_HELM_ion_neutral_blends: .false.
+ max_logRho_neutral_HELM: 3d0
+ logT_neutral_HELM: 5.0d0
+ logT_ion_HELM: 6.5d0
+ eosDT_use_linear_interp_for_X: .false.
* Z_all_HELM: 1d99 => 0.04d0
* min_x_for_keep: 1d-23 => 1d-5
* min_x_for_n: 1d-18 => 1d-4
* min_x_for_add: 1d-17 => 1d-4
* kappa_file_prefix: 'gn93' => 'gs98'
* kappa_CO_prefix: 'gn93_co' => 'gs98_co'
- use_eosDE_get: .true.
- helmeos_max_logRho_neutral: -1d99
- helmeos_logT_neutral: -1d99
- helmeos_logT_ion: -1d99
+ center_R_lower_limit: -1
+ mach1_mass_upper_limit: -1
+ use_fixed_L_for_BB_outer_BC: .false.
+ fixed_L_for_BB_outer_BC: 0
+ max_logT_for_net: 10.2d0
+ cgs_thermal_diffusion_eta_full_on: 0d0
+ cgs_thermal_diffusion_eta_full_off: 2d0
+ do_Ne22_sedimentation_heating: .false.
+ diffusion_use_full_net: .false.
+ eps_Ne22_sedimentation_factor: 1
+ tol_bad_max_correction: 0d0
+ dvdt_visc_factor: 1d0
+ use_artificial_viscosity: .false.
+ artificial_viscosity_Q_shift: -1
+ post_shock_viscosity_decay_factor: -1
+ pre_shock_viscosity_decay_factor: -1
+ art_visc_full_on_logRho_ge_this: -99
+ art_visc_full_off_logRho_le_this: -99
+ force_timestep_min: 0
+ force_timestep_min_years: 0
+ force_timestep_min_factor: 2d0
+ limit_for_max_E_residual: 1d99
+ hard_limit_for_max_E_residual: 1d99
* diffusion_use_cgs_solver: .false. => .true.
* diffusion_v_max: 1d-5 => 1d-3
* diffusion_gamma_full_off: 50 => 175
* diffusion_gamma_full_on: 40 => 150
* kap_phot_factor_for_kap_floor: 1d0 => -1d0
* Gamma_lnS_eps_grav_full_on: 160d0 => 150d0
* Gamma_lnS_eps_grav_full_off: 130d0 => 120d0
* shock_spread_linear: -1d0 => 0
* shock_spread_quadratic: -1d0 => 1d-3
- trace_binary_rlo: .true.
+ show_binary_log_description_at_start: .true.
+ relax_primary_to_th_eq: .false.
+ log_Lnuc_div_L_for_relax_primary_to_th_eq: 0.005
+ min_age_for_relax_primary_to_th_eq: 1d2
+ max_steps_for_relax_primary_to_th_eq: 1000
+ no_history_during_relax_primary_to_th_eq: .true.
+ reset_age_for_relax_primary_to_th_eq: .true.
+ tsync_for_relax_primary_to_th_eq: 1
+ history_interval: 5
+ terminal_interval: 1
+ write_header_frequency: 10
+ extra_binary_terminal_output_file: ''
+ report_rlo_solver_progress: .false.
* limit_retention_by_mdot_edd: .true. => .false.
* max_tries_to_achieve: 0 => 20
> On May 31, 2017, at 11:40 AM, Evan Bauer <ebauer at physics.ucsb.edu> wrote:
> Hello MESA Users,
> We are happy to announce a new MESA public release version: 9793. Release notes below.
> Evan (on behalf of Pablo, Rich, Josiah, Rob, Aaron, Frank, and Bill)
> MESA 9793 - Release Notes
> Binary module terminal output has been significantly changed to provide more useful information at runtime.
> Element Diffusion
> Option added to simplify the user interface for diffusion classes (diffusion_use_full_net = .true.). This assigns each isotope in the current net to its own diffusion class, circumventing the need to manually specify each isotope that serves as a class representative.
> Default solver for the Burgers equations has been changed to properly account for electron degeneracy in white dwarfs. The old solver based on the work of Thoul et al. (1994) for the Sun is still available, but is no longer the default.
> The new diffusion solver (diffusion_use_cgs_solver = .true.) now accounts for thermal diffusion terms in the Burgers equations. It includes options to turn off these terms smoothly as a function of electron degeneracy, since they are invalid and usually negligible for WD interiors. With these thermal diffusion terms included for non-degenerate material, this solver is appropriate for use on the main sequence as well as in WDs. Thus, in most circumstances, the user no longer needs to worry about switching between solvers for different types of diffusion problems.
> The new diffusion solver also supports radiative levitation now, where before this required using the solver based on Thoul et al. (1994) (as revised by Haili Hu to include extra driving terms from radiative levitation). The user interface remains unchanged.
> Note that radiative levitation is limited to isotopes that are supported by OP data for calculating radiative accelerations, so one must be careful to avoid using a class structure involving isotopes that aren’t supported. In particular, radiative levitation with the new diffusion_use_full_net option will cause problems for any net including isotopes not supported by OP.
> “diffusion_gamma_full_on/off” has been changed to 150-175 as default, meaning diffusion occurs for all material that has not crystallized.
> An option that accounts for Ne22 sedimentation heating in old, cooling WDs is now available.
> Minor bug fixed in calculating the ionization states for diffusion classes. This only affected diffusion, and only mattered when background material for diffusion was not hydrogen.
> Fixed a bug that could cause timesteps to crash unnecessarily due to stale diff_iters value. (Thanks to Bill Wolf for bug report!)
> GYRE has been upgraded to 5.0rc6. There are significant changes to GYRE’s input files in this updated version; users should consult example input files, and/or the online documentation, for details. https://bitbucket.org/rhdtownsend/gyre/wiki/Home <https://bitbucket.org/rhdtownsend/gyre/wiki/Home>
> Pgstar plots now have the option to call a decorator function, set in a users run_star_extras file in a similar way to the other_* functions, this allows additional customization options for all pgplots. See $MESA_DIR/star/other/pgstar_decorator.f90 for an example.
> Network plots now have the option to show a colorbar for the abundances.
> The default version of JINA reaclib provided with MESA has been updated to reaclib v2.2. To restore the previous default, use
> jina_reaclib_filename = ‘jina_reaclib_results_20130213default2'
> in your star_job inlist. When switching between different reaclib versions, users must either: run $MESA_DIR/empty_caches or set rates_cache_dir in their star_job inlist to a local folder. Otherwise mesa will not use any rates from the new reaclib version that have been cached from a previous version. Running $MESA_DIR/clean will also clear the caches.
> It is now much simpler for a user to specify a custom table providing a weak reaction rate. See the test suite case "custom_rates" for examples of providing user-specified reaction rate tables.
> Bug fix for redo_conv_for_dr_lt_mixing_length so that it now applies to central convection zones.
> The default opacity files all use gs98 relative metal fractions. (Previously, this was a mix of gn93 and gs98). As a reminder, these choices are controlled by the star_job options
> The star_job options associated with the blend (in Z) between OPAL and HELM and the blend (in T) between ionized and neutral HELM have been refactored. See the documentation for the controls:
> While "stella" is now part of the mesa release, it is not yet ready for general use. Please wait for the next release to begin using it.
> A file containing a summary of the changes to inlist defaults is attached. Some internal and experimental controls are not included.
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the Mesa-users