[mesa-users] MLT problem

Kenny Van kvan at ualberta.ca
Mon Jan 11 14:30:32 EST 2016


Hi Bill,

I'll try that before sending over the files I used. Unfortunately I have a
slightly modified version of MESA 7763 that adds additional acceleration
terms to the atmospheric calculations.

Thanks

On Mon, Jan 11, 2016 at 11:51 AM, Bill Paxton <paxton at kitp.ucsb.edu> wrote:

> Hi Kenny,
>
> I'd suggest doing a small patch to the mlt.f90 code.  Where you find the
> test
>             if (is_bad_num(d_Bcubed_dvb(mlt_dgrada))) then
> instead of printing out the error messages you get now, change it to do
>                ierr = -1
>                return
> Then the calling routine in mesa/star will know it has to do a retry.
>
> This is a leftover test that I should have removed before the release.
> Sorry 'bout that.
>
> -Bill
>
>
> On Jan 11, 2016, at 10:24 AM, Kenny Van wrote:
>
> > Hi everyone,
> >
> > So I've run into some issues with the code breaking. From the output of
> the error it seems like the issue occurs in the mlt.f90 file in the private
> directory inside of MLT. I have been trying to diagnose the problem myself
> and I've run into a bit of a wall. Looking at the error output I think the
> output is either in how the calculation is done for dA_dvb, which produces
> NaN. This error then propagates to d_Bcubed_dvb which is calculated using
> dA_dvb. If I understand the theory correctly, the A value, a0, diff_grads
> are all non negative values. The other place where this may be going wrong
> is with the d_diff_grads_dvb value. While I have some intuition on the
> diff_grads value, I'm not sure if d_diff_grads_dvb can be negative.
> Hopefully someone has some insight on what I can do to try and fix this
> error.
> >
> >
> >                                d_Bcubed_dvb(mlt_dgrada)
>      NaN
> >                            d_diff_grads_dvb(mlt_dgrada)
>  -1.0000000000000000D+00
> >                                      dA_dvb(mlt_dgrada)
>      NaN
> >
> >
> > Thanks
> >
> ------------------------------------------------------------------------------
> > Site24x7 APM Insight: Get Deep Visibility into Application Performance
> > APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month
> > Monitor end-to-end web transactions and take corrective actions now
> > Troubleshoot faster and improve end-user experience. Signup Now!
> >
> http://pubads.g.doubleclick.net/gampad/clk?id=267308311&iu=/4140_______________________________________________
> > 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/20160111/47ec9c70/attachment.html>


More information about the Mesa-users mailing list