ramses garvity

src/camr.cc

@@ -412,8 +412,8 @@ template <class T> int CAmr::loadData(uns::CParticles<T> * particles,
                   if (nvarg==ndim+1) { // there is pot/ax/ay/az
                     offset=1;          //  shift to ax
                   }
-                  for (int k=0; k<nvarg; k++) {
-                    particles->acc.push_back(varg[(offset+k)*ngrida*twotondim+ind*ngrida+i]);
+                  for (int k=offset; k<nvarg; k++) {
+                    particles->acc.push_back(varg[k*ngrida*twotondim+ind*ngrida+i]);
                   }
                   if (ndim<3) { // ndim = 2 , set az=0.0
                     particles->acc.push_back(0.0);