correct E field advance in 1D

ALaDyn · Mar 16, 2015 · d4d23d6 · d4d23d6
1 parent f3b5a9f
commit d4d23d6
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/em_field.cpp b/em_field.cpp
@@ -887,6 +887,7 @@ void EM_FIELD::new_advance_E(CURRENT *current)
       EY   = &val[my_indice(edge,YGrid_factor, ZGrid_factor, 1, i,   j,   k,   N_grid[0], N_grid[1], N_grid[2], Ncomp)];
       EZ   = &val[my_indice(edge,YGrid_factor, ZGrid_factor, 2, i,   j,   k,   N_grid[0], N_grid[1], N_grid[2], Ncomp)];
 
+      *EX += dt*(                    - mygrid->den_factor*current->Jx(i, j, k) );
       *EY += dt*( - dxi*(BZ - BZ_XM) - mygrid->den_factor*current->Jy(i, j, k) );
       *EZ += dt*(   dxi*(BY - BY_XM) - mygrid->den_factor*current->Jz(i, j, k) );
 #else