|
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) |
Return to
addradial.py
CVS log
Up to [Development] / JSOC / proj / globalhs / sosh