Search found 84 matches

by warrick
Wed Jun 12, 2019 1:31 am
Forum: General Discussion
Topic: Using GYRE machinery to evaluate eigenvalues of Laplace's tidal equations
Replies: 4
Views: 4456

Re: Using GYRE machinery to evaluate eigenvalues of Laplace's tidal equations

Also, thanks for including the Python module for reading and interpolating in the eigenvalue data! Glad you've found it useful. I'm considering releasing the TAR stuff as a separate, stand-alone library, so that people can easily incorporate it in their science without needing the full GYRE install...
by warrick
Tue Jun 11, 2019 1:57 pm
Forum: General Discussion
Topic: Using GYRE machinery to evaluate eigenvalues of Laplace's tidal equations
Replies: 4
Views: 4456

Re: Using GYRE machinery to evaluate eigenvalues of Laplace's tidal equations

Aha, thanks. Is this documented somewhere that I should have found? I did a cursory glance through some GYRE-related papers but didn't find a description of details like this (presumably along with the TAR implementation). Also, thanks for including the Python module for reading and interpolating in...
by warrick
Tue Jun 11, 2019 4:06 am
Forum: General Discussion
Topic: Using GYRE machinery to evaluate eigenvalues of Laplace's tidal equations
Replies: 4
Views: 4456

Using GYRE machinery to evaluate eigenvalues of Laplace's tidal equations

Hi, I'm currently trying to reproduce some g-mode analysis that involves a formula with eigenvalues of Laplace's tidal equation. This is beyond my usual area of expertise but, given a vague notion that this is something GYRE needs for the Traditional Approximation of Rotation (TAR), I went digging a...
by warrick
Thu May 16, 2019 4:12 am
Forum: Bug reports
Topic: Different results from FGONG or AMDL version of model file
Replies: 3
Views: 7031

Re: Different results from FGONG or AMDL version of model file

Finally, to close this off for anywhere else who finds this thread, Rich already fixed this in the development version. I just cleaned up my input file to use Model S's value of G (6.67232e-8) and I get the same radial mode frequencies for the FGONG and the AMDL at the machine precision level.

W
by warrick
Thu May 16, 2019 3:30 am
Forum: Bug reports
Topic: Different results from FGONG or AMDL version of model file
Replies: 3
Views: 7031

Re: Different results from FGONG or AMDL version of model file

As usual, shortly after answering the question, I find the answer. This is almost certainly a bug in the calculation of V_2. e.g. in gyre_amdl_file.fpp, V_2 is set (for x > 0) by V_2 = V_g*Gamma_1 whereas I've inferred from the other (working) model readers that it should be V_2 = V_g*Gamma_1/x**2 T...
by warrick
Thu May 16, 2019 2:35 am
Forum: Bug reports
Topic: Different results from FGONG or AMDL version of model file
Replies: 3
Views: 7031

Different results from FGONG or AMDL version of model file

Hi, I seem to be getting different oscillation frequencies for a stellar model depending on whether I use an FGONG file or an AMDL file, converted using either ADIPLS' internal fgong converter (if you have MESA, it's at $MESA_DIR/adipls/adipack.c/adiajobs/fgong-amdl.d.f) or my own Python-based conve...
by warrick
Wed Mar 20, 2019 10:39 am
Forum: Bug reports
Topic: Fail to output the temperature when model_type is 'EVOL'
Replies: 4
Views: 5326

Re: Fail to output the temperature when model_type is 'EVOL'

Hi again Stefano, The following modifications will allow you to get T from FGONG files, at least in GYRE 5.2, which is what I'm using. The way I did this was in essence to compare $GYRE_DIR/src/model/gyre_fgong_file.fpp with $GYRE_DIR/src/model/gyre_mesa_file.fpp and add lines related to beta_rad to...
by warrick
Wed Mar 20, 2019 10:25 am
Forum: Bug reports
Topic: Fail to output the temperature when model_type is 'EVOL'
Replies: 4
Views: 5326

Re: Fail to output the temperature when model_type is 'EVOL'

Hi Stefano, The main issue looks like a genuine bug: I also get an error when I try to add T as an output option for an FGONG file. So the documentation is advertising something that doesn't appear to be available, unless we've both missed some extra parameter. From what I can see, the problem is th...
by warrick
Wed Mar 13, 2019 2:54 am
Forum: Bug reports
Topic: Problem reading FGONG file
Replies: 6
Views: 5664

Re: Problem reading FGONG file

Hi Rich, I'm responsible for producing this file, which I did with MESA r9575. I'm having a look through now and, though I can't find the modified inlist, I'm pretty sure I changed the MESA control `format_for_FGONG_data` to the old default ("1P5E16.9,x") to make the FGONG files play nicely with the...
by warrick
Wed Feb 06, 2019 6:27 am
Forum: Bug reports
Topic: Segfault when prune_modes= .true.
Replies: 7
Views: 4867

Re: Segfault when prune_modes= .true.

The problem with this is that not all the summary data is available at that point (since we haven't yet calculated all the modes). This is less of a problem for TXT output (because we write that file line-by-line), but for HDF5 output we would have to continually re-write the file as new modes are ...