[mesa-users] MLT problem
kvan at ualberta.ca
Tue Jan 19 11:47:31 EST 2016
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.
> On Mon, Jan 11, 2016 at 11:51 AM, Bill Paxton <paxton at kitp.ucsb.edu>
>> Hi Kenny,
>> I'd suggest doing a small patch to the mlt.f90 code. Where you find the
>> 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
>> 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.
>> 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
>> > d_Bcubed_dvb(mlt_dgrada)
>> > d_diff_grads_dvb(mlt_dgrada)
>> > dA_dvb(mlt_dgrada)
>> > 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!
>> > 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...
More information about the Mesa-users