meaning of failure to produce any output
Posted: Mon May 18, 2020 2:55 am
Hi all,
I have a general question. I'm trying to compute just the fundamental pressure mode for a large number of MESA models, so cycling over a GYRE template inlist restricted to a certain range (in muhz, which is not the best unit for this kind of star but is the one I'm most familiar with--I've made sure the max and min values are right).
GYRE only produces output for about a third of my models, despite very similar conditions, and I can't figure out why. I've tried changing parameters in MESA that should increase the resolution of the model, lowering the timestep in MESA to make sure I'm getting the stopping point I want, etc.
Is there something obviously wrong with this inlist? Attached
In cases where no output is generated, the terminal reports
gyre [5.1]
----------
OpenMP Threads : 16
Input filename : gyre_template.out
Model Init
----------
Reading from FGONG file
File name 12.0m.A1.90.Z0.02.eta0.8.endage15.6000Myr_forMC_seismicTol.gyre
File version 1300
Read 671 points
No need to add central point
Mode Search
-----------
Mode parameters
l : 0
m : 0
Building frequency scan
added scan interval : 0.1009E-03 -> 0.1009E-02 (10000 points, INVERSE)
Building x grid
Found inner turning points, x range ****** -> 1.0000
Adding 3 global point(s) in iteration 1
Adding 3 global point(s) in iteration 2
Adding 2 global point(s) in iteration 3
Adding 1 global point(s) in iteration 4
Adding 0 global point(s) in iteration 5
Final grid has 1 segment(s) and 680 point(s):
Segment 1 : x range 0.0000 -> 1.0000 (1 -> 680)
Starting search (adiabatic)
Root bracketing
Time elapsed : 1.932 s
Root Solving
l m n_pg n_p n_g Re(omega) Im(omega) chi n_iter
Time elapsed : 0.000 s
Basically I don't know what [nothing] means, or what's going wrong (or why a not-insignificant portion of these models work just fine). Any insight is appreciated
I have a general question. I'm trying to compute just the fundamental pressure mode for a large number of MESA models, so cycling over a GYRE template inlist restricted to a certain range (in muhz, which is not the best unit for this kind of star but is the one I'm most familiar with--I've made sure the max and min values are right).
GYRE only produces output for about a third of my models, despite very similar conditions, and I can't figure out why. I've tried changing parameters in MESA that should increase the resolution of the model, lowering the timestep in MESA to make sure I'm getting the stopping point I want, etc.
Is there something obviously wrong with this inlist? Attached
In cases where no output is generated, the terminal reports
gyre [5.1]
----------
OpenMP Threads : 16
Input filename : gyre_template.out
Model Init
----------
Reading from FGONG file
File name 12.0m.A1.90.Z0.02.eta0.8.endage15.6000Myr_forMC_seismicTol.gyre
File version 1300
Read 671 points
No need to add central point
Mode Search
-----------
Mode parameters
l : 0
m : 0
Building frequency scan
added scan interval : 0.1009E-03 -> 0.1009E-02 (10000 points, INVERSE)
Building x grid
Found inner turning points, x range ****** -> 1.0000
Adding 3 global point(s) in iteration 1
Adding 3 global point(s) in iteration 2
Adding 2 global point(s) in iteration 3
Adding 1 global point(s) in iteration 4
Adding 0 global point(s) in iteration 5
Final grid has 1 segment(s) and 680 point(s):
Segment 1 : x range 0.0000 -> 1.0000 (1 -> 680)
Starting search (adiabatic)
Root bracketing
Time elapsed : 1.932 s
Root Solving
l m n_pg n_p n_g Re(omega) Im(omega) chi n_iter
Time elapsed : 0.000 s
Basically I don't know what [nothing] means, or what's going wrong (or why a not-insignificant portion of these models work just fine). Any insight is appreciated