Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

minos_bran produces NaNs #3

Open
Chunli-Dai opened this issue Jul 18, 2014 · 0 comments
Open

minos_bran produces NaNs #3

Chunli-Dai opened this issue Jul 18, 2014 · 0 comments

Comments

@Chunli-Dai
Copy link

When n=81, l=83, the results from minos_bran are NaNs.
This problem occurs when I use the following input:
input model file:
../models/prem_noocean.txt
output file:
prem_noocean_S
eigenfunction file (output):
eprem_noocean_S
enter eps and wgrav
1.E-12 200.
enter jcom (1=rad;2=tor;3=sph;4=ictor)
3
enter lmin,lmax,wmin,wmax,nmin,nmax
83 84 0. 200. 81 81

It seems that in line 2957 of minos_bran.f, the denominator of a2 is zero. The same case
for a1.
The NaN value of wtry is then given to we(2) (Line 209), which cause infinite loop for the next degree l=84.

I'd really appreciate it if someone can fix this problem. It'll help me a lot. Thanks.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant