             }
 #endif
             complex<float> this_f1c = -qOverm * intVecArray(thisT, this_e1_dot_gradvf0);
-
+            //if (iX == write_iX && iP == write_iP) {
+            //        for(int i=0; i<nSteps;i++){
+            //            std::cout<<"this_e1_dot_gradvf0[i]: "<<this_e1_dot_gradvf0[i]<<std::endl;
+            //        }
+            //        std::cout<<"this_f1c: "<<this_f1c<<std::endl;
+            //}
 #if LOWMEM_ORBIT_WRITE >= 1
             if (iX == write_iX && iP == write_iP) {
 

            float v0z_i = ThisParticleList[iP].v_c3;

            float h = dv * Ze;

#pragma omp critical // "atomic" does not work for complex numbers
            {

                B1_XYZ = hanningWeight[i] * exp(-_i * wrf * thisT[i]) * getE1orB1_XYZ(thisParticle_XYZ, &r[0], &b1_CYL[0], r.size(), nPhi);
                thisB1c_XYZ[i] = B1_XYZ * (1 - thisParticle_XYZ.status);

                //if (iX == write_iX && iP == write_iP) {
                //    std::cout<<"E1_XYZ: "<<E1_XYZ<<std::endl;
                //    std::cout<<"status: "<<thisParticle_XYZ.status<<std::endl;
                //    std::cout<<"hanningWeight: "<<hanningWeight[i]<<std::endl;
                //}

#if DEBUG_MOVE >= 2
                std::cout << "thisE1c[i].c1: " << thisE1c_XYZ[i].c1 << std::endl;
                std::cout << "thisE1c[i].c2: " << thisE1c_XYZ[i].c2 << std::endl;
