diff --git a/ConSol/ConSol/Source/AxisymmEuler_std_noSwirl.cpp b/ConSol/ConSol/Source/AxisymmEuler_std_noSwirl.cpp
index bd25f48..44b30ae 100644
--- a/ConSol/ConSol/Source/AxisymmEuler_std_noSwirl.cpp
+++ b/ConSol/ConSol/Source/AxisymmEuler_std_noSwirl.cpp
@@ -1,35 +1,35 @@
 //
 // Created by Fabian Moenkeberg on 2020-02-13.
 //
 
 #include "AxisymmEuler_std_noSwirl.hpp"
 
 std::vector<double> AxisymmEuler_std_noSwirl::source(std::vector<std::vector<double>> &U, MeshBase *mesh, double t, double dt, int iCell){
-    std::vector<double> ret(U.size(),0);
+    std::vector<double> ret(U[iCell].size(),0);
 
 //    long double rrho, rp, rm1, rm2, rE;
 //    rrho = U[iCell][0];
 //    rm1 = U[iCell][1];
 //    rm2 = U[iCell][2];
 //    rE = U[iCell][3];
 //    rp = (gamma - 1)*(rE - 0.5*(rm1*rm1 + rm2*rm2)/rrho);
 //
 //    std::vector<double> r = mesh->MidPoint(iCell);
 //    ret[2]  = rp/r[1];
     long double rho, p, m1, m2, E;
     rho = U[iCell][0];
     m1 = U[iCell][1];
     m2 = U[iCell][2];
     E = U[iCell][3];
     p = (gamma - 1)*(E - 0.5*(m1*m1 + m2*m2)/rho);
 
-    if (iCell == 2225){
-        std::cout << iCell;
-    }
+//    if (iCell == 2225){
+//        std::cout << iCell;
+//    }
     std::vector<double> r = mesh->MidPoint(iCell);
-    ret[0]  = -m1/r[1];
+    ret[0]  = -m2/r[1];
     ret[1]  = -m1*m2/rho/r[1];
-    ret[2]  = -m1*m1/rho/r[1];
-    ret[3]  = -m1/rho*(E+p)/r[1];
+    ret[2]  = -m2*m2/rho/r[1];
+    ret[3]  = -m2/rho*(E+p)/r[1];
     return ret;
 }
\ No newline at end of file