From 6b71dd47722765bca3946edd79d08ac23fb68ff9 Mon Sep 17 00:00:00 2001
From: Paul Dechamps <paul.dechamps@uliege.be>
Date: Mon, 1 Jul 2024 15:25:35 +0200
Subject: [PATCH] (chore) Added density in boundary layer output file
---
blast/utils.py | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/blast/utils.py b/blast/utils.py
index b32492a..93efffc 100644
--- a/blast/utils.py
+++ b/blast/utils.py
@@ -123,7 +123,7 @@ def getSolution(sections, write=False, toW=['x', 'xelm', 'theta', 'H', 'deltaSta
varKeys = ['theta', 'H', 'N', 'ue', 'Ct']
attrKeys = ['x', 'y', 'z', 'xoc', 'deltaStar', \
'cd', 'cf', 'Hstar', 'Hstar2', 'Hk', 'ctEq', 'us', 'delta',\
- 'vt', 'Ret']
+ 'vt', 'rhoe', 'Ret']
elemsKeys = ['xelm', 'yelm', 'zelm']
for iSec, sec in enumerate(sections):
@@ -194,8 +194,8 @@ def getSolution(sections, write=False, toW=['x', 'xelm', 'theta', 'H', 'deltaSta
f.write('# Stagnation point at index '+str(sec[0].x.size()-1)+'\n')
f.write('# Wake at index '+str(sec[0].x.size()+sec[1].x.size())+'\n')
- for key in toW:
- f.write(' {0:>21s}'.format(key))
+ for i, key in enumerate(toW):
+ f.write(' {0:>21s}'.format(str(i)+' '+key))
f.write('\n')
for i in range(len(sol[iSec]['x'])):
for key in toW:
@@ -229,4 +229,4 @@ def plot(cfg):
if 'legend' in cfg and cfg['legend'] is True: plt.legend(frameon=False)
if 'xlabel' in cfg: plt.xlabel(cfg['xlabel'])
if 'ylabel' in cfg: plt.ylabel(cfg['ylabel'])
- plt.show()
\ No newline at end of file
+ plt.show()
--
GitLab