|
version 1.3, 2019/02/24 18:34:10
|
version 1.4, 2019/03/04 17:33:22
|
|
|
|
| try: | try: |
| colorshift | colorshift |
| except: | except: |
| colorshift=0 |
colorshift=1 |
| | |
| try: | try: |
| animate | animate |
| Line 154 instr=sosh.catchc("Enter number of modes |
|
| Line 154 instr=sosh.catchc("Enter number of modes |
|
| if (instr == 'q'): | if (instr == 'q'): |
| sys.exit() | sys.exit() |
| nmodes=int(instr) | nmodes=int(instr) |
| for i in range(nmodes): |
for i in range(abs(nmodes)): |
| print("MODE #%i" % (i+1)) | print("MODE #%i" % (i+1)) |
| | |
| (l,m,n)=sosh.querylmn() |
(l,m,n)=sosh.querylmn(nmodes) |
| if l == -1: | if l == -1: |
| sys.exit() | sys.exit() |
| ntest = sosh.modeln[l==sosh.modell] | ntest = sosh.modeln[l==sosh.modell] |
| while (n not in ntest): | while (n not in ntest): |
| print("That n was not modelled for that l. Modelled values are in the range %i to %i. Try again." \ | print("That n was not modelled for that l. Modelled values are in the range %i to %i. Try again." \ |
| % (ntest.min(),ntest.max())) | % (ntest.min(),ntest.max())) |
| (l,m,n)=sosh.querylmn() |
(l,m,n)=sosh.querylmn(nmodes) |
| ntest = sosh.modeln[l==sosh.modell] | ntest = sosh.modeln[l==sosh.modell] |
| | |
| signedm=m | signedm=m |
| Line 218 for i in range(nmodes): |
|
| Line 218 for i in range(nmodes): |
|
| # freqlist.append(freq[l+signedm]*freqscale) | # freqlist.append(freq[l+signedm]*freqscale) |
| # amplist.append(modeparms[ind,3]) | # amplist.append(modeparms[ind,3]) |
| | |
| |
nmodes=abs(nmodes) |
| sumr=0.0 | sumr=0.0 |
| sumt=0.0 | sumt=0.0 |
| sump=0.0 | sump=0.0 |
|
|
|
| | |
| fig,im = sosh.drawfigure(var, fsize=figsize) | fig,im = sosh.drawfigure(var, fsize=figsize) |
| if (sosh.isave != 0): | if (sosh.isave != 0): |
| sosh.savefigure(ianimate=animate, label='_interior') |
sosh.savefigure(ianimate=animate, label='_interiorsum') |
| if (show != 0): | if (show != 0): |
| if (animate != 0): | if (animate != 0): |
| anim = animation.FuncAnimation(fig, sumanimate, frames=nframes, interval=interval) | anim = animation.FuncAnimation(fig, sumanimate, frames=nframes, interval=interval) |
|
|
|
| c=sosh.queryplotparms() | c=sosh.queryplotparms() |
| if (c == 'q'): | if (c == 'q'): |
| break | break |
| |
fs=sosh.catchc("Enter new freqscale: ", freqscale) |
| |
if (fs == 'q'): |
| |
break |
| |
else: |
| |
freqscale=float(fs) |
| |
nf=sosh.catchc("Enter new nframes: ", nframes) |
| |
if (nf == 'q'): |
| |
break |
| |
else: |
| |
nframes=int(nf) |
| |
sosh.nframes=nframes |
| | |
| | |
| | |
Return to
addradial.py
CVS log
Up to [Development] / JSOC / proj / globalhs / sosh