[mesa-users] MLT problem

Kenny Van kvan at ualberta.ca
Tue Jan 19 11:47:31 EST 2016


Hi Bill,

I'm not sure I'm properly implementing the ierr = -1. If I understand
correctly, I would not only need to include ierr = -1 in the test of
(is_bad_num(d_Bcubed_dvb(mlt_dgrada))), but I would also need to define it
in the initial parameters of the subroutine correct?

On Mon, Jan 11, 2016 at 12:30 PM, Kenny Van <kvan at ualberta.ca> wrote:

> 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/20160119/41311537/attachment.html>


More information about the Mesa-users mailing list