Fix dynamo benchmark data evaluation for Case 0

Author: hiromatsui

src/Fortran_libraries/MHD_src/sph_MHD/const_data_4_dynamobench.f90

@@ -85,16 +85,23 @@ subroutine const_dynamobench_data                                 &
 !
       if(bench%ipwr_ocore .gt. 0) then
         irank_copy = pwr%v_spectr(bench%ipwr_ocore)%irank_m
-        if(my_rank .eq. irank_copy)  then
-          call copy_kin_energy_4_dbench(bench%ipwr_ocore, pwr,          &
+        if(ipol%base%i_velo .gt. 0) then
+          if(my_rank .eq. irank_copy) then
+            call copy_kin_energy_4_dbench(bench%ipwr_ocore, pwr,        &
      &                                  bench%KE_bench)
-          call copy_mag_energy_4_dbench(bench%ipwr_ocore, pwr,          &
-     &                                  bench%ME_bench)
+          end if
+          call calypso_mpi_bcast_real(bench%KE_bench,                   &
+     &                                cast_long(ithree), irank_copy)
+        end if
+!
+        if(ipol%base%i_magne .gt. 0) then
+          if(my_rank .eq. irank_copy) then
+            call copy_mag_energy_4_dbench(bench%ipwr_ocore, pwr,        &
+     &                                    bench%ME_bench)
+          end if
+          call calypso_mpi_bcast_real(bench%ME_bench,                   &
+     &                                cast_long(ithree), irank_copy)
         end if
-        call calypso_mpi_bcast_real(bench%KE_bench, cast_long(ithree),  &
-     &                              irank_copy)
-        call calypso_mpi_bcast_real(bench%ME_bench, cast_long(ithree),  &
-     &                              irank_copy)
       end if
 !
       if(sph_MHD_bc%sph_bc_U%iflag_icb .eq. iflag_rotatable_ic) then

src/Fortran_libraries/MHD_src/sph_MHD/field_at_mid_equator.f90

@@ -221,7 +221,7 @@ subroutine cal_field_4_dynamobench                                &
       phi_prev(1:4) = phi_zero(1:4)
 !
       ave_drift = quad * (drift(1) + drift(2) + drift(3) + drift(4))
-      do nd = 1, 7
+      do nd = 1, d_circle%ntot_phys
         d_zero(0,nd) = quad * (d_zero(1,nd) + d_zero(2,nd)              &
      &                       + d_zero(3,nd) + d_zero(4,nd))
       end do