[mesa-users] Taking smaller time steps

[Liron Mcley]

Thank you for the help!
I added a few lines in extras_check_model (run_star_extras) and it worked
perfectly.

I have another small question - I'm using a few of the s% x_ctrl variables
and I would like to add them to the profile, for example, I want to see s%
x_ctrl(1) for each time step. So I tried adding this:

call do_val(i, 'L_conv_max_per_mass', s% x_ctrl(1))

at the end of the do_val() lines in do_profile_info(...) which is in
star/private/profile.f . But it's not doing anything and I can't see why. I
also tried putting the line inbetween the do_val() lines and nothing.

Thank you,
Liron



On 1 June 2014 19:21, Bill Paxton <paxton at kitp.ucsb.edu> wrote:

> Hi,
>
> Take a look at the built-in timestep controls by looking in
> star/defaults/controls.defaults -- search for the string "max_timestep ="
> to find the section.
>
> For terminal output concerning the selection of timesteps, turn on some or
> all of these:
>
>       report_why_dt_limits = .true.
>          ! if true, produce terminal output about choice of timestep
>       report_all_dt_limits = .false.
>          ! if true, produce terminal output about all influences for
> choice of timestep
>       report_hydro_dt_info = .false.
>          ! if true, produce terminal output about choice of timestep based
> on varcontrol_target
>       report_dX_nuc_drop_dt_limits = .false.
>          ! if true, report timestep limits from drop in abundance from
> nuclear reactions
>
> For the nitty-gritty details, star/private/timestep.f is the place to
> look.   But before you start rewriting it, let us know!  ;D
>
> Rather than change code in mesa, you can use your run_star_extras to write
> your own code to set the next timestep in extras_check_model.
>
> Cheers,
> Bill
>
>
>
>
>
> On Jun 1, 2014, at 5:01 AM, Liron Mcley wrote:
>
> > Hi everyone,
> >
> > I went through the mailing list and found that theoretically it is
> possible to make mesa take smaller time steps by modifying a few lines in
> /star/private/timestep.f . But this isn't working for me at all.
> > I put,
> >    write (*,*) 'time_step (s)=', s% dt_next / secyer
> > at the beginning of timestep_controller just to see if I get this line
> when mesa runs, but nothing. Can anyone please help me ?
> >
> > Thank you,
> > Liron Mcley
> >
> ------------------------------------------------------------------------------
> > Time is money. Stop wasting it! Get your web API in 5 minutes.
> > www.restlet.com/download
> >
> http://p.sf.net/sfu/restlet_______________________________________________
> > mesa-users mailing list
> > mesa-users at lists.sourceforge.net
> > https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140629/3dcd1068/attachment.html>