Skip to content
Snippets Groups Projects
Commit 70a01b34 authored by Sebastiano Bernuzzi's avatar Sebastiano Bernuzzi
Browse files

Added the example output for multiple radii in prerelease_example1

parent 3823e0e8
No related branches found
No related tags found
No related merge requests found
......@@ -120,14 +120,14 @@ wm = mwaves(path = bamsim['sim-folder'], code = 'bam', filenames = dfiles,
ignore_negative_m=True)
# Write .txt files in CoRe format
# Here dump to file only largest radius
# NOTE: mwaves() currently assumes every extraction radius has the same
# modes. Something to be fixed.
for (l,m) in wm.modes:
psilm = wm.get(var='Psi4',l=l, m=m)
psilm.write_to_txt('Psi4', bamsim['pre-release-folder'])
hlm = wm.get(l=l, m=m)
hlm.write_to_txt('h', bamsim['pre-release-folder'])
# modes.
for r in wm.radii:
for (l,m) in wm.modes:
psilm = wm.get(var='Psi4',l=l, m=m, r=r)
psilm.write_to_txt('Psi4', bamsim['pre-release-folder'])
hlm = wm.get(l=l, m=m)
hlm.write_to_txt('h', bamsim['pre-release-folder'])
# Compute and write energetics
wm.energetics(bamsim['massA'], bamsim['massB'], bamsim['madm'], bamsim['jadm'],
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment