version 1.2, 2019/02/22 22:26:26 version 1.3, 2019/02/24 18:34:10
 Line 181  for i in range(nmodes):
 Line 181  for i in range(nmodes):
(y,dy)=(plm[...,l-m],dplm[...,l-m])     (y,dy)=(plm[...,l-m],dplm[...,l-m])
print("Spherical harmonic retrieved.")     print("Spherical harmonic retrieved.")

if (i == 0):
if (m > 0):   if (m > 0):
phi[:,int(nx/2):nx]=(np.pi/4)/m     phi[:,int(nx/2):nx]=(np.pi/4)/m
phi[:,0:int(nx/2)]=(np.pi/4)/m + np.pi     phi[:,0:int(nx/2)]=(np.pi/4)/m + np.pi
 Line 199  for i in range(nmodes):
 Line 200  for i in range(nmodes):
fh=np.interp(r,sosh.rmesh,xih*np.sqrt(sosh.rho))   fh=np.interp(r,sosh.rmesh,xih*np.sqrt(sosh.rho))

idx = ((sosh.modeln == n) & (sosh.modell == l))   idx = ((sosh.modeln == n) & (sosh.modell == l))
freq = sosh.modelnu[idx]/1000.0    freq = sosh.modelnu[idx]*1000.0
#  ind = ((l==lmod) & (n==nmod)) #  ind = ((l==lmod) & (n==nmod))
#  freq = float(modeparms[ind,2])/1000.0 #  freq = float(modeparms[ind,2])/1000.0
rclist.append(fr)   rclist.append(fr)

Legend:
 Removed from v.1.2 changed lines Added in v.1.3