Commit 38d04787 authored by Roland Haas's avatar Roland Haas
Browse files

POWER: reowork so that time is set as a real valued array

parent c5754a38
...@@ -450,13 +450,11 @@ def POWER(sim_path, radii, modes): ...@@ -450,13 +450,11 @@ def POWER(sim_path, radii, modes):
hTable = psi4ToStrain(mp_psi4_vars[i], f0) # table of strain hTable = psi4ToStrain(mp_psi4_vars[i], f0) # table of strain
time = hTable[:, 0] time = hTable[:, 0].real
h = hTable[:, 1] h = hTable[:, 1]
hplus = h.real hplus = h.real
hcross = h.imag hcross = h.imag
newhTable = np.column_stack((time, hplus, hcross)) newhTable = np.column_stack((time, hplus, hcross))
warnings.filterwarnings('ignore')
finalhTable = newhTable.astype(float)
strain.append(finalhTable) strain.append(finalhTable)
#------------------------------------------------------------------- #-------------------------------------------------------------------
...@@ -467,11 +465,9 @@ def POWER(sim_path, radii, modes): ...@@ -467,11 +465,9 @@ def POWER(sim_path, radii, modes):
# print(len(h_phase), "h_phase length") # print(len(h_phase), "h_phase length")
# print(len(time), "time length") # print(len(time), "time length")
angleTable = np.column_stack((time, h_phase)) ### start here angleTable = np.column_stack((time, h_phase)) ### start here
angleTable = angleTable.astype(float) ### b/c t is defined based on
phase.append(angleTable) ### time here phase.append(angleTable) ### time here
h_amp = np.absolute(h) h_amp = np.absolute(h)
ampTable = np.column_stack((time, h_amp)) ampTable = np.column_stack((time, h_amp))
ampTable = ampTable.astype(float)
amp.append(ampTable) amp.append(ampTable)
#---------------------------------------------------------------------- #----------------------------------------------------------------------
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment