bubble-dynamics

commit a44793fa3b64c460c2cb675919424df4e888932b
parent d5cd12233e9f3fe0497e3122a592816312dc9e9a
Author: Ursula Anna Rasthofer <ursular@euler02.hpc-lca.ethz.ch>
Date:   Fri,  3 Feb 2017 19:12:13 +0100

generalized output for clouds such that all bubbles are written

Diffstat:
M
src/kernels/KMCluster_TY.h
 | 
10
+++++-----

1 file changed, 5 insertions(+), 5 deletions(-)
diff --git a/src/kernels/KMCluster_TY.h b/src/kernels/KMCluster_TY.h
@@ -157,18 +157,18 @@ public:
         // // const Real odat[6] = {U[0][0], U[1][0], U[4][0], U[2][0], U[5][0], U[3][0]};
         // const Real odat[7] = {U[0][0], U[7][0], U[1][0], pb[0], pb[1], pb[2], bd.pext(t)};
 
-        const size_t N = 102;
+        const size_t N = 2*_N+2;
         double out[N];
         double Vg = 0.0;
-        for (int i = 0; i < 50; ++i)
+        for (int i = 0; i < _N; ++i)
         {
             const double r = U[i][0];
             out[i] = r;
-            out[i+50] = bd.pBubble(U, i);
+            out[i+_N] = bd.pBubble(U, i);
             Vg += 4.0/3.0*r*r*r*M_PI;
         }
-        out[100] = bd.pext(t);
-        out[101] = Vg;
+        out[2*_N] = bd.pext(t);
+        out[2*_N+1] = Vg;
 
         m_dumper->write(step, t, dt, &out[0], N);