-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathpars_velocity.py
43 lines (32 loc) · 1.32 KB
/
pars_velocity.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
####################################################################
# Parameter file for computing and plotting velocity spectrograms. #
####################################################################
rotrates = ['0.0'] #Rotation rate cases (rad/s)
bouncetimes={} #Time at which bounce occurs (seconds)
bouncetimes['0.0'] = 0.300 #in seconds
bouncetimes['2.0'] = 0.316 #in seconds
#Spectrogram parameters
Window='bohman'
WindowWidth=35 #in ms
Scaling='spectrum'
Mode='complex' #'psd' for visualizing spectrograms, 'complex' for filtering and inverting spectrograms
#Which components of velocity to plot
plotv = ['velr','veltheta'] #,'veltheta'] #'velx','vely'
#Which rotrate to plot (only one at a time)
plotr = '0.0'
#Which points to plot (specify range of r (km) and theta (rad))
minr = 74.0
maxr = 86.0
mintheta = 0.9*pi/4
maxtheta = 1.1*3*pi/4#1.1*3*pi/4
#Plotting parameters
titles=[r'v_r',r'$v_\theta$'] #,r'$v_\theta$'] #r'$v_x$',r'$v_y$'
###################
# Filter parameters
###################
rotrates_to_filter = ['0.0']
############################################
#Spherical harmonic decomposition parameters
############################################
rotrate_to_decompose = '0.0' #one at a time for now
poloidal_numbers = range(0,11,1) #which angular momentum numbers l to compute spherical harmonics for