         float tmpJzIm = imag(j1zc[iX]);
         nc_j1zc_re.putVar(&tmpJzRe);
         nc_j1zc_im.putVar(&tmpJzIm);
+
+        std::cout<<"j1xc[iX]: "<<j1xc[iX]<<std::endl;
     }
 
     // ProfilerStop();

                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 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;

            float v0z_i = ThisParticleList[iP].v_c3;

            float h = dv * Ze;


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