[mesa-users] Hook in HR when hitting max_age

Bill Paxton paxton at kitp.ucsb.edu
Fri Jan 31 17:53:52 EST 2014


Hi Josiah,

So smaller max timesteps make the problem worse!!!!

smaller dt => worse approx for eps_grav?

groan.......

-B





On Jan 31, 2014, at 2:45 PM, Josiah Schwab wrote:

> 
> paxton at kitp.ucsb.edu writes:
> 
>> For now, interpolating to the max age is probably the best work around
>> -- with checks to make sure the results are stable when you go to a
>> smaller max timestep (convergence, convergence, convergence).  ;p
> 
> I did a little test over lunch.
> 
> First I ran Radek's inlist as is, which takes 1714 steps and then
> terminates at the max_age limit.  This has a max_timestep of 1e7 yr.
> 
> Then, I restarted from step 1700 and 
>  + evolved past the max_age with a max_timestep of 1e7 yr
>  + evolved to the max_age with a max_timestep of 3e6 yr
>  + evolved to the max_age with a max_timestep of 1e6 yr
>  + evolved to the max_age with a max_timestep of 1e5 yr
> 
> I've attached an HR diagram showing the results.
> 
> Which is just to say that this is all about having a smaller timestep,
> rather than reaching the last timestep.
> 
> Josiah
> 
> <histories.pdf>





More information about the Mesa-users mailing list