From 70a01b3452f55430c40e84855473ff3ce2cebba1 Mon Sep 17 00:00:00 2001
From: sebastiano bernuzzi <sebastiano.bernuzzi@gmail.com>
Date: Sat, 23 Oct 2021 13:35:01 +0200
Subject: [PATCH] Added the example output for multiple radii in
prerelease_example1
---
msc/coredb_prerelease_example1.py | 14 +++++++-------
1 file changed, 7 insertions(+), 7 deletions(-)
diff --git a/msc/coredb_prerelease_example1.py b/msc/coredb_prerelease_example1.py
index bd05518..40304d2 100644
--- a/msc/coredb_prerelease_example1.py
+++ b/msc/coredb_prerelease_example1.py
@@ -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'],
--
GitLab