[mesa-users] MLT problem

[Bill Paxton]

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