Runtime Error with GYRE v2.0

[ehsan]

Hi.
I just updated to the new version (2.0), and successfully completed "make test".
when I run GYRE over a given model (gs98-gs98.gyre in the attachments), it terminates with the following error:

Code: Select all

       1    -142       2     144    0.1619607532114381E-01    0.2093843124502976E-03    0.1704274604505912E-15    45
       1    -148       1     149    0.1579055484734884E-01    0.7165224201230905E-04    0.7503237960139201E-14    14
       1    -148       1     149    0.1579055484734884E-01    0.7165224201230628E-04    0.4965210610430295E-14    12
       1    -145       2     147    0.1587924138199240E-01    0.1854755422659538E-03    0.8947960780067702E-14    15
 ASSERT 'n_iter <= this%np%n_iter_max' failed at line 430 <gyre_nad_bvp:mode>:
 Too many iterations

The namelist that I use is

Code: Select all

&coeffs
        file = 'LOGS/gs98-gs98.GYRE'
        coeffs_type = 'EVOL'
        file_format = 'MESA'
        G = 6.6742800000e-08
/

&osc
        l = 1

/        
&num
        ivp_solver_type = 'MAGNUS_GL2'
        n_iter_max = 100
/

&scan
        grid_type = 'INVERSE'
        !freq_units = 'UHZ'
        freq_min = 0.01
        freq_max = 1.5
        n_freq = 2000
        
/

&shoot_grid
/

&recon_grid
/


&output
      summary_file = 'LOGS/gs98-gs98.h5'
      summary_item_list = 'l,n_p,n_g,omega,freq,E'
      mode_prefix = ''
      mode_item_list = 'M_star,R_star,L_star,m,rho,p,n,l,n_p,n_g,omega,freq,E,E_norm,W,x,V,As,U,c_1,Gamma_1,nabla_ad,delta,xi_r,xi_h,phip,dphip_dx,delS,delL,delp,delrho,delT,dE_dx,dW_dx,T'
      freq_units = 'UHZ'

/

I would be grateful to have this issue fixed, and kept informed.

Best regards.
Ehsan.

[rhtownsend]

Hi Ehsan --

It seems that your stellar model didn't get attached - please could you post it to this forum.

cheers,

Rich

[ehsan]

Hi Rich,
I re-attach it; let me know please if it does not arrive again.
I also noticed now that the forum does not upload .gyre ascii files. This is not bad to look at.
Additionally, I do not know how to attach multiple files to a message.

Best.
Ehsan.

[rhtownsend]

Hi Ehsan --

I've modified the forum settings to accept .mesa, .gyre, .fgong extensions, plus extensions for shell scripts and plain text files.

I'll have a look into the issues you're experiencing with this model.

cheers,

Rich

[ehsan]

This is an update on the same run.
With a factor of 5 times increase in n_iter_max to boost it to 500, I still get the same error.

Rich, can you please let me know what is being iterated? The root finding routine?
The fact that i get the error with an order of magnitude larger n_iter_max than the default points that there is something wrong with my settings, probably!

Ehsan.