geodynamics / mineos

Mineos computes synthetic seismograms in a spherically symmetric non-rotating Earth by summing normal modes.
http://geodynamics.org/cig/software/mineos/
GNU General Public License v2.0
24 stars 19 forks source link

minos_bran produces NaNs #3

Open Chunli-Dai opened 10 years ago

Chunli-Dai commented 10 years ago

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.