bugfix in one IC movie calculation

magic-sph · Sep 9, 2024 · b640b20 · b640b20
1 parent 9828a47
commit b640b20
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/store_movie_IC.f90 b/src/store_movie_IC.f90
@@ -349,8 +349,8 @@ subroutine store_movie_frame_IC(bICB,b_ic,db_ic,ddb_ic,aj_ic,dj_ic)
                            frames(n_o)=        help*O_sin_theta(nThetaR) * &
                            &    ( cBrB(nThetaR,nPhi0)*BtB(nThetaR,nPhi0) - &
                            &      cBtB(nThetaR,nPhi0)*BrB(nThetaR,nPhi0) )
-                           frames(n_o)=            help*O_sin_theta(nThetaR) * &
-                           &    ( cBrB(nThetaR,nPhi180)*BtB(nThetaR,nPhi180) - &
+                           frames(n_o+n_field_size)=help*O_sin_theta(nThetaR) * &
+                           &    ( cBrB(nThetaR,nPhi180)*BtB(nThetaR,nPhi180) -  &
                            &      cBtB(nThetaR,nPhi180)*BrB(nThetaR,nPhi180) )
                         end do
                      else