[mesa-users] MESA: doing something more than reporting...

Bill Paxton paxton at kitp.ucsb.edu
Tue Jan 18 21:45:04 EST 2011


Hi,

Turns out that some fool (me) left a debugging trap enabled.  ; - (
I'll put out a new release asap with that fixed.  
If you'd like to make a tiny patch, replace this file in star/private.
Then in star/test do
./mk
./export

Sorry for the inconvenience. 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: timestep.f
Type: application/octet-stream
Size: 74745 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20110118/dcda0a50/attachment.obj>
-------------- next part --------------


Cheers,
Bill


On Jan 18, 2011, at 5:34 PM, Evan O'Connor wrote:

> Hello all,
> 
> I just ran into something weird.  For the situation I am looking at, when I switch on 
> 
> report_all_dt_limits = .true.
> 
> the code crashes, if it is false then it is fine.  some of the error message is below.  
> 
>>>>>>>>>>>>>>>> 
>       1510   8.180391  1.244E+04   4.636411   4.755512  14.674456  11.547502   0.000000   0.012976   0.280000   0.477836  -3.242907    705      1
>   2.248222   3.082165   1.683802   3.145216   3.469588  -6.929051   3.126954   0.979459   0.000502   0.020000   0.502238   0.041347      2      1
> 1.1563E+07  19.104688   4.700862   4.088597   2.141492 -10.288952   0.000000   0.000728   0.004313   0.000396  1.993E-02  0.260E-04    varcontrol
> 
>                        dX n14  623    0.2006471728E-02    0.2000000000E-02    0.8955172923E+00    0.1445299224E-02
>                             s% D_mix(k)         623    0.0000000000000000E+00
>                                   X_old         623    1.4982242077306349E-02
>                              s% xa(j,k)         623    1.2975770349680954E-02
>                                      dX         623   -2.0064717276253950E-03
>                                   D_mix         618    0.0000000000000000E+00    4.7017707893630550E-03
>                                   D_mix         619    0.0000000000000000E+00    3.9011863067364982E-03
>                                   D_mix         620    0.0000000000000000E+00    3.1006018241099415E-03
>                                   D_mix         621    0.0000000000000000E+00    1.4994328588572581E-03
>                                   D_mix         622    0.0000000000000000E+00    7.4971642942862904E-04
>                                   D_mix         623    0.0000000000000000E+00    0.0000000000000000E+00
>                                   D_mix         624    1.2954886313713439E+13   -1.4994328588571471E-03
>                                   D_mix         625    1.6061308510682049E+13   -2.6642224508062681E-03
>                                   D_mix         626    1.8407963009267238E+13   -3.8290120427557223E-03
>                                   D_mix         627    2.0443522267985812E+13   -4.9938016347049544E-03
>                                   D_mix         628    2.2259202775794098E+13   -6.1585912266542975E-03
> 
> check_dX
>>>>>>>>>>>>>>>>>> 
> 
> For me this is occurring right after helium burning in a 15Msun star (when the Ne20 and N14 levels are changing).  For reference, the inlist_15_z2m2 is identical to mesa/start/test_suite/15_si_burn/inlist_15_z2m2 except this flag is changed and I do not load the test_suite model but rather start from pre-main sequence).  I've also seen this occur in the same place in stars up to 19Msun.
> 
> It seems like the extra reporting causes mesa to fail.  I'm not sure why, thoughts?
> 
> Evan
> 
> 
> 
> 
> ------------------------------------------------------------------------------
> Protect Your Site and Customers from Malware Attacks
> Learn about various malware tactics and how to avoid them. Understand 
> malware threats, the impact they can have on your business, and how you 
> can protect your company and customers by using code signing.
> http://p.sf.net/sfu/oracle-sfdevnl
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users



More information about the Mesa-users mailing list