(chore) Added density in boundary layer output file

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()